Impact Factor 4.400

The 1st most cited open-access journal in Pharmacology & Pharmacy

Original Research ARTICLE

Front. Pharmacol., 06 July 2017 | https://doi.org/10.3389/fphar.2017.00446

Detection of Emodin Derived Glutathione Adduct in Normal Rats Administered with Large Dosage of Polygoni Multiflori Radix

  • State Key Laboratory of Natural Medicines, China Pharmaceutical University, Nanjing, China

Polygoni Multiflori Radix (PMR) has been commonly used as a tonic in China for centuries. PMR-associated hepatotoxicity has been drawing increasingly more attention in recent years in parallel with its wide utilization. Anthraquinones (AQs) are recognized as the main hepatotoxic components in PMR. However, the exact underlying mechanism of AQs poisoning is still not fully understood. Herein, we proposed a hypothesis that metabolic activation of AQs such as emodin was involved in PMR-induced liver injury, AQs followed to generate the electrophilic reactive metabolites and subsequently formed covalent adduct with cellular nucleophiles in the liver to exert hepatotoxicity. In the present study, the link of cytotoxicity of PMR in primary human hepatocytes and the depletion of glutathione (GSH) was investigated by MTT assay and UHPLC-QqQ-MS/MS analysis. The results showed that PMR depleted GSH and therefore induced cytotoxicity. Then, emodin-GSH adduct was identified in bile of liver injured rats after intragastric administration of PMR or emodin with the aid of UHPLC-QTOF-MS/MS method. Our findings not only provided confirmative evidence that the mechanism of hepatotoxicity induced by AQs in PMR involved key metabolic steps, but also revealed that emodin-GSH adduct had potential to be further developed as a sensitive and traceable biomarker for the assessment of PMR-induced liver injury.

Introduction

Polygoni Multiflori Radix (PMR), derived from the dried root of Polygonum multiflorum Thunb. (Family Polygonaceae), has been commonly used as a Chinese herbal medicine in clinical practice for health promotion and disease treatment for thousands of years, with the Chinese name “Heshouwu” (Lin et al., 2015a). It is traditionally valued for hair-blacking, liver and kidney-tonifying and anti-aging effects (Lin et al., 2015c). However, PMR-induced liver injury in clinic has been constantly reported in recent years (Yu et al., 2011). The earliest case of liver injury induced by PMR was recorded in 1996 in China (But et al., 1996). Since then, an increasing number of PMR-poisoning cases have been documented in many countries, including Australia, Italy, Holland, England, America, Korea, etc. (Park et al., 2001; Mazzanti et al., 2004; Panis et al., 2005; Cárdenas et al., 2006; Laird et al., 2008; Zhang et al., 2009; Jung et al., 2011). So far, a broad spectrum of chemicals such as stilbenes, anthraquinones (AQs), flavonoids, and phospholipids have been isolated from PMR (Lin et al., 2015b). Among these compounds, AQs including emodin, physcion, rhein, aloe-emodin, and chrysophanol, are commonly and conveniently believed to be primarily responsible for the PMR-associated hepatotoxicity, especially emodin, a predominant AQs occurring in PMR, is assumed to be one of the leading hepatotoxic component (Zhang et al., 2010; Li et al., 2012; Lin et al., 2015b; Lv et al., 2015; Ma et al., 2015). Nevertheless, the specific mechanism of AQs-induced liver injury is still unclear, which limits the establishment of effective approaches for diagnosis and treatment of liver injury induced by PMR.

It is well known that hepatotoxic effect of many xenobiotics requires metabolism (Tang et al., 1999; James et al., 2003; Peterson, 2013; Ruan et al., 2014). Hepatic metabolic activation of drugs to electrophilic reactive metabolites has been suggested to initiate the development of drug-induced liver toxicity (Guo, 2012). The reactive metabolites covalently bind to intracellular nucleophiles such as amino acids, glutathione (GSH), proteins, and nucleic acids, thereby triggering a series of pathologic alterations and causing liver injury (James et al., 2003; Fu et al., 2004; Lin et al., 2011). In general, conjugation with GSH is a major detoxification route for many xenobiotics (Chung et al., 1997; Chen et al., 2009), since this conjugation reaction protects cells against the harmful effects of reactive metabolites. On the other hand, the overproduction of electrophilic reactive metabolites will deplete GSH and eventually produce the toxicity (Kouzi et al., 1994; Shimizu et al., 2011). In this regard, the reactive metabolite-GSH conjugates have ever been developed as mechanism-based biomarkers to assess the toxicity of xenobiotics (Lin et al., 2011; Ruan et al., 2014).

Quinone compounds such as doxorubicin and menadione have been used as therapeutic drugs. However, clinical applications of these drugs are generally limited because of their hepatotoxicities. A prevailing mechanism of quinone-induced toxicity is that quinones as reactive electrophiles deplete GSH and covalently bind to proteins (Ishihara et al., 2011). AQs, structurally belonging to natural quinones, occur widely in many herbal drugs, including Rheum palmatum (Family Polygonaceae), Cassia occidentalis (Family Leguminosae), Morinda citrifolia (Family Rubiaceae), etc. Unfortunately, these AQ-containing drugs were documented to be toxic for human beings (Stadlbauer et al., 2008; Wang et al., 2009; Panigrahi et al., 2015). Despite many poisoning cases have been reported, there is still no clear consensus on the research of AQs-induced liver injury.

In the present study, for the purpose of verifying the depletion of GSH might be involved in PMR-induced liver injury, we tested the cytotoxic effect of PMR in relationship to its effect on the levels of intracellular GSH using MTT assay and ultra-high-performance liquid chromatography coupled to triple quadrupole mass spectrometry (UHPLC-QqQ-MS/MS) method. With repeated administration of 20 g/kg PMR for 3 weeks, significant elevations of serum total bilirubin (TBIL) and alkaline phosphatase (ALP) levels were observed, and massive necrosis was found in the liver with the histopathological analysis, indicating the liver injury induced by PMR. Furthermore, using ultra-high-performance liquid chromatography-quadrupole/time of flight mass spectrometry (UHPLC-QTOF-MS/MS) method, emodin-GSH conjugate was detected for the first time in the bile of liver injured rats after intragastric administration of PMR or emodin in corroboration of the mechanism proposed for AQs-induced liver injury initiated by metabolic activation.

Materials and Methods

Chemicals, Reagents and Animals

The decoction pieces of PMR were purchased from Bozhou traditional Chinese medicine market (Anhui Province, China). The samples were authenticated by Prof. Hui-Jun Li and deposited at State Key Laboratory of Natural Medicines (China Pharmaceutical University). Emodin was from Chengdu Must Bio-technology Co. (Chengdu, China). Dulbecco’s modified eagle medium (DMEM) was obtained from Gibco (Grand Island, NY, United States). Fetal bovine serum (FBS) was from Sijiqing (Hangzhou, China). GSH and GSSG assay kit and Bicinchoninic acid (BCA) assay kit were from Beyotime (Shanghai, China). Acetonitrile of HPLC grade was purchased from TEDIA Company (Fairfield, United States). Formic acid, acetic acid, and ammonium acetate (HPLC grade) were from ROE (Neward, New Castle, DE, United States). GSH and glutathione ethyl ester (internal standard, IS) were from Aladdin (Shanghai, China). Deionized water was prepared using a Milli-Q purification system produced by Millipore (Milford, MA, United States). All other reagents were analytical grade.

Male Sprague-Dawley (SD) rats, weighing 200-220 g, were provided by Sino-British Sippr/BK Lab Animal Ltd. (Shanghai, China). The animal studies were conducted in accordance with the Provision and General Recommendation of Chinese Experimental Animals Administration Legislation and the Instructive Notions with Respect to Caring for Laboratory Animals, and were approved by Animal Ethics Committee of China Pharmaceutical University and the Department of Science and Technology of Jiangsu Province [license number: SYXK (SU) 2016-0011]. The animals were housed under controlled conditions (temperature 22 ± 2°C, relative humidity 50 ± 10%) with a natural light-dark cycle for 1 week before the experiment was carried out.

Preparation of Herbal Samples

Polygoni Multiflori Radix decoction pieces were cracked into small chunks, then the samples were extracted thrice under reflux each with 10 times of 75% ethanol for 3 h. The three extracts were combined and concentrated under reduced pressure, lyophilized, producing the dried PMR extract with a yield of 20.78%. The contents of the main constituents in PMR extracts were quantified by UHPLC method (Ma et al., 2012). The analytical results showed that the concentrations of emodin, physcion and 2,3,5,4′-tetrahydroxystilbene-2-O-β-D-glucopyranoside in PMR extract were 0.16, 0.13, and 5.69%, respectively. The lyophilized powders were re-dissolved and dispersed in 0.5% carboxymethyl cellulose sodium salt (CMC-Na) aqueous solution for intragastric administration to rats.

Primary Human Hepatocytes Culture

Primary human hepatocytes were obtained from ScienCell (Carlsbad, CA, United States). The hepatocytes were recovered using the manufacturer’s protocol. Hepatocytes were grown in hepatocyte medium (ScienCell, Carlsbad, CA, United States) supplemented with 10% FBS, 1% hepatocyte growth supplement, 1% penicillin/streptomycin and maintained in a humidified atmosphere of 95% O2-5% CO2 at 37°C. The pH of the medium was maintained at 7.4. The hepatocytes were plated in 96-well round-bottom tissue culture plates to test cytotoxicity and plated in 6-well culture plates for GSH depletion assay. At the beginning of the experiment, when plated hepatocytes had reached about 85% confluence, the growth medium was removed from the wells, the hepatocytes were washed with PBS and fresh medium was used as base for treatment with the drugs.

Cytotoxicity Assay in Primary Human Hepatocytes

MTT dye uptake in primary human hepatocytes in the presence and absence of PMR extract were determined. In brief, primary human hepatocytes (5 × 103 cells/well) were seeded into 96-well plate and cultured for 24 h. The hepatocytes were then exposed to different concentrations of PMR extract. After 24 h of treatment, the media were removed and 100 μL MTT (0.5 mg/mL in medium) was added to each well. The plates were further incubated for another 4 h. Then, the media was removed from the wells and 100 μL of DMSO was added to each well. Finally the absorbance was taken at 490 nm after 10 min, on a micro plate reader. The percentage viability was calculated by comparing the absorbance of control and treated cells.

GSH Depletion Assay in Primary Human Hepatocytes

Cellular GSH was determined using a UHPLC-QqQ-MS/MS method. Primary human hepatocytes were exposed for 24 h to either DMSO (control) or to increasing concentrations of PMR extract in 6-well culture plates. At the end of the incubation, the cell samples were washed twice with PBS before repeated freezing-thawing for three times in acetic acid-ammonium acetate buffering solution (pH = 3.5) and then centrifuged at 13000 rpm for 15 min at 4°C. Methanol (600 μL) was added to the supernatant (200 μL) to precipitate the protein, the supernatant was dried under nitrogen, reconstituted using 100 μL of ammonium acetate solution, vortex-mixed for 1 min, and centrifuged at 13000 rpm for 10 min at 4°C. The supernatant was injected into the UHPLC-QqQ-MS/MS system for analysis. The protein concentration was determined with BCA assay kit, and the GSH levels were expressed as ng/mL per mg protein.

UHPLC-QqQ-MS/MS Analytical Condition and Method Validation

Samples were analyzed on an LC-MS/MS system consisting of an Agilent series 1290 UHPLC system and an Agilent 6460 triple quadrupole mass spectrometer (Agilent Technologies, Palo Alto, CA, United States). The separation was carried out on a Waters Xselect CSH C18 column (4.6 mm × 150 mm, 3.5 μm) at 30°C. Mobile phase adopted for this study was acetonitrile/water containing 0.1% (v/v) formic acid (10/90) and the flow rate was 0.6 mL/min. The mass spectrometer was operated in positive ion mode with electrospray ionization (ESI) source. The operating conditions were as follows: drying gas temperature, 325°C; drying gas flow, 10 L/min; sheath gas temperature, 300°C; sheath gas flow, 10 L/min; nebulizer pressure, 35 psi; capillary voltage, 3.5 kV. Quantitation was performed by multiple reaction monitoring (MRM) mode. Data collection and processing were conducted with MassHunter Workstation 05.00 (Agilent Technologies, United States).

Stock solution of GSH and IS were separately prepared in acetic acid-ammonium acetate buffering solution at the concentrations of 1 mg/mL. The GSH stock solution was diluted into a series of working solutions with the concentrations ranging from 5 to 20000 ng/mL, each containing 3.125 μg/mL of IS. Quality control (QC) samples with low, middle, and high concentrations were also prepared in the same way (10, 500, and 15000 ng/mL). All of the working solutions were stored at 4°C before use.

The UPLC-QqQ-MS/MS method was fully validated in terms of specificity, limit of quantification (LOQ), linearity, sensitivity, accuracy, precision, extraction recovery, matrix effect and stability, according to US FDA guidelines (United States Food and Drug Administration, 2001).

Animal Protocol

Male Sprague-Dawley rats were randomly divided into two groups (control group and treatment group in different periods, n = 5) by intragastric administration with vehicle (1 mL/100 g rat) and PMR extract (20 g/kg) for 3 weeks. Blood, bile and liver samples were harvested at 1st, 2nd, and 3rd week, respectively. After the latest administration of PMR, urine from 0 to 24 h was collected and stored at –80°C until analysis. Blood samples were drawn at 24 h via the post-ocular vein, the plasma was generated by centrifugation at 4000 rpm for 15 min at 4°C and stored at -80°C before use. The rats were anesthetized with 20% urethane (0.5 mL/100 g, intraperitoneal injection), then the common bile ducts were cannulated with PE-10 tubing. The bile from 0 to 12 h was collected and stored at -80°C before use. At the end of the experiment, rats were sacrificed by cervical dislocation, and the livers in each group (n = 5) were harvested for histopathological examination. In a parallel study, emodin was suspended in 0.5% CMC-Na aqueous solution. The rats were treated with a single dosage of either vehicle or emodin (200 mg/kg) as control for the identification of GSH adduct. Urine from 0 to 24 h and bile from 0 to 12 h were collected and stored at –80°C until use. Blood samples were drawn at 24 h via the post-ocular vein for analysis.

Biochemical Analysis and Histopathologic Examination

Blood samples collected at 24 h after the latest treatment of PMR were centrifuged at 4000 rpm for 15 min at 4°C. The alanine aminotransferase (ALT) activity, aspartate aminotransferase (AST) activity, total bile acid (TBA) concentration, TBIL concentration and ALP activity were determined by an automatic blood biochemical analyzer (Beckman Coulter LX20, United States). Liver function testing kits were supplied by Zhongda Hospital (Nanjing, China).

The rat liver tissues were fixed in 10% neutral buffered formaldehyde solution, paraffin-processed, and sectioned at 4 μm. The sections were stained with haematoxylin-eosin (HE) and examined for histopathological changes under the microscope (Olympus DX45, Japan).

Detection of Emodin-GSH Conjugate in Plasma, Urine and Bile

Methanol (threefolds) was added to plasma, urine and bile to precipitate the protein, then vortex-mixed for 1 min and centrifuged at 13000 rpm for 10 min. The supernatants of plasma and urine samples were dried under nitrogen and dissolved in 100 μL of methanol for UPLC-QTOF-MS/MS analysis. The supernatant (2 μL) of bile samples was directly injected into the UPLC-QTOF-MS/MS system for analysis.

Biological specimens were analyzed on an Agilent 6530 QTOF-MS/MS with ESI in negative-ion mode. The chromatographic separations were achieved on a Shimadzu VP-ODS column (250 mm × 4.6 mm, 5 μm). The mobile phase of water containing 0.1% (v/v) formic acid (solution A) and acetonitrile (solution B) was used with a gradient elution as follows: 0–55 min, 8–45% B; 55–57 min, 45–100% B; 57–65 min, 100% B. The UHPLC flow rate was 0.8 mL/min. The operating conditions for ESI mass spectrometry were as follows: drying gas temperature, 320°C; capillary voltage, 3500 V; drying gas flow rate, 10 L/min; nebulizing gas pressure, 35 psi; sheath gas temperature, 350°C; sheath gas flow, 11 L/min; fragmentor voltage, 125 V. The collision energies were 30 and 60 V. The molecular mass was accurately measured in negative-ion mode to analyze emodin-GSH conjugate. The characteristics of ion pairs for emodin-GSH conjugate were m/z 574 → 272.

Statistical Analysis

The statistical analysis was performed using a one-way ANOVA with the Dunnett’s post hoc test to determine the significance of the differences between the individual groups. A value of p < 0.05 was considered as statistically significant.

Results and Discussion

Effect of PMR on Hepatocytes Viability and Cellular GSH

MTT assay was used to estimate the cytotoxicity of PMR to primary human hepatocytes. As shown in Figure 1, cell viabilities were decreased in a dose-dependent manner after exposure to different doses of PMR (50, 150, 250, 350, and 400 μg/mL) for 24 h. The IC50 values of PMR were 331.30 μg/mL. These results revealed the in vitro hepatotoxicity of PMR.

FIGURE 1
www.frontiersin.org

FIGURE 1. Cell viabilities of primary human hepatocytes exposed to Polygoni Multiflori Radix (PMR) for 24 h. PMR inhibited the cell viabilities in a dose-dependent manner. Data are presented as means ± SD of three independent experiments. p < 0.5, ∗∗p < 0.01 vs. control.

As the cytotoxicity was observed with the treatment of PMR, we further investigated the relationship between GSH depletion and the cytotoxicity in primary human hepatocytes, with the aid of UHPLC-QqQ-MS/MS method. In order to achieve a good resolution, the chromatographic conditions were optimized. Finally, the MRM transitions were monitored at m/z 308.10 → 76.10 and m/z 336.10 → 76.10 for GSH and IS, respectively. The fragmentor voltage values set for GSH and IS were 110 and 100 V, and the collision energies were 25 and 30 V, respectively.

The representative UHPLC-QqQ-MS/MS MRM chromatograms of GSH and IS are shown in Figure 2. The calibration exhibited satisfactory linearity, ranging from 5 to 20000 ng/mL. The LOQ of GSH was 1.25 ng/mL. The intra- and inter-day precision and accuracy were almost less than 15%, indicating the accuracy and precision of the method were acceptable for the quantitative analysis. The extraction recoveries of GSH were 91.50–94.73%, respectively. The matrix effect was also found to be satisfactory.

FIGURE 2
www.frontiersin.org

FIGURE 2. Representative chromatograms of glutathione (GSH) and glutathione ethyl ester [internal standard (IS)] in primary human hepatocytes samples. (A) Blank hepatocytes sample; (B) blank hepatocytes sample spiked with GSH and IS; (C) hepatocytes sample treated with PMR. 1, GSH; 2, IS.

Hence, the intracellular GSH levels were quantified by the newly developed UHPLC-QqQ-MS/MS method and the quantitative results are illustrated in Figure 3. As shown in Figure 3, treatment with the concentration of 25 μg/mL PMR resulted in an increase in GSH level compared to control group, indicating a protective effect of PMR in low concentration, whereas dramatic GSH depletion was found to occur in the concentrations of 50, 100, 200, and 300 μg/mL. At the concentration of 50 μg/mL, PMR significantly decreased GSH level by 12.05% as compared to control group. All of these results suggested that GSH depletion might be a critical event in the mechanism of PMR-induced liver injury.

FIGURE 3
www.frontiersin.org

FIGURE 3. Depletion of GSH induced by PMR in primary human hepatocytes in vitro. PMR depleted cellular GSH in a dose-dependent manner when PMR concentration was greater than 50 μg/mL. Data are presented as means ± SD of three independent experiments. p < 0.05, ∗∗p < 0.01 vs. control.

Evaluation of Hepatotoxicity of PMR

Serum liver enzymes are the most commonly used biomarkers for diagnosis and assessment of hepatic diseases in clinic (Abboud and Kaplowitz, 2007). In the present study, five conventional biomarkers of liver injury, namely ALT activity, AST activity, TBA concentration, TBIL concentration and ALP activity, were examined to evaluate the liver injury induced by PMR. As can be seen in Figure 4, the levels of TBA did not show significant changes between the control group and the treatment groups. Intragastric gavage of PMR caused significant elevations of serum ALT and AST activities at the 1st week (serum ALT and AST: 103.00 and 504.67 U/L for treatment group; 52.40 and 96.25 U/L for control group). No obvious alterations of serum ALT and AST activities were observed at the 2nd and 3rd weeks (serum ALT and AST: 40.60 and 100.00 U/L at 2nd week; 45.25 and 92.25 U/L at 3rd week). The TBIL concentrations showed a time-dependent increase during the treatment of PMR, with significant alterations at 2nd and 3rd weeks compared to the control group (3.05 and 3.40 μmol/L at 2nd and 3rd weeks; 1.77 μmol/L for control group). Meanwhile, the administration of PMR for 3 weeks showed a significant elevation of ALP levels when compared to the control group (418.75 U/L at 3rd week; 205.00 U/L for control group). The histopathologic evaluations (Figure 5) demonstrated that a small amount of inflammatory cells were observed in livers of PMR-treated rats at 1st and 2nd weeks (Figures 5B,C), while liver cell necrosis could be found after 3 weeks of PMR treatment (Figure 5D), demonstrating the insensitivity of conventional biochemical parameters in plasma like ALT, AST, and TBA for the indicating function of PMR-induced liver injury. Remarkably, TBIL and ALP might be more sensitive biomarkers during the process of liver injury induced by PMR. These findings were consistent with previous report (Tu et al., 2015; Zhao et al., 2017). Additionally, the patients often take multi-herb preparation containing PMR (Lei et al., 2015), leading to the uncertainty whether the liver injury is truly induced by PMR. It is clinically necessary to identify the traceable biomarkers for the diagnosis of PMR-induced liver injury. Meanwhile, the GSH/GSSG ratio in PMR-administrated liver was analyzed according to the technical manual of the detection kit. As shown in Figure 6, the GSH/GSSG were significantly decreased during 3 weeks, indicating that PMR administration depleted GSH in livers of rats, and slight increases were observed at the 2nd and 3rd weeks when compared to the 1st week, this would be considered as a compensatory response of rats. The result revealed that GSH depletion might be involved in PMR-induced liver injury.

FIGURE 4
www.frontiersin.org

FIGURE 4. Effects of PMR on the levels of serum biochemical parameters. Total bilirubin (TBIL) concentrations were increased in a time-dependent manner, and a significant increase in plasma alkaline phosphatase (ALP) levels was observed after 3 weeks treatment of PMR, which might be more sensitive than alanine aminotransferase (ALT), aspartate aminotransferase (AST), and total bile acid (TBA). Data are expressed as mean ± SD (n = 5). p < 0.05, ∗∗p < 0.01 compared with control groups.

FIGURE 5
www.frontiersin.org

FIGURE 5. Histopathologic evaluation (H&E staining) of liver tissues obtained from rats treated with PMR (n = 5). Original magnification was 200× for each sample. Rats were treated with vehicle (A), PMR for 1 week (B), PMR for 2 weeks (C), PMR for 3 weeks (D). Bar = 200 μm. After 3 weeks treatment of PMR, necroses were observed in the livers of rats, indicating the liver injury induced by PMR. Inflammatory cell infiltration and necrotic region were indicated by solid and dotted arrows, respectively.

FIGURE 6
www.frontiersin.org

FIGURE 6. Hepatic GSH/GSSG ratio of PMR-treated rats. During the treatment of PMR, GSH/GSSG ratios were significantly decreased in livers of rats. p < 0.05 compared with the control group.

Identification of Emodin-GSH Adduct in Liver Injured Rats

Considering that quinones are easily subjected to metabolic activation (Wu et al., 2014), we speculate that AQs in PMR preferentially undergo phase I metabolism mediated by cytochromes P450s (CYPs) and form electrophilic epoxide intermediate metabolite, then covalently bind to nucleophilic GSH (He et al., 2015). Owing to the relative high abundance in PMR extract, emodin was exemplarily selected as model compound. Figure 7 is the proposed bioactivation route of emodin-GSH.

FIGURE 7
www.frontiersin.org

FIGURE 7. Proposed metabolic biotransformation of emodin and the generation of emodin-GSH adduct. CYPs mediated emodin bioactivation involved the formation of an electrophilic epoxide intermediate. Theoretically, the epoxidation reaction of anthracene ring occurred in multiple positions, resulting in the diversity of GSH conjugation of emodin.

A UHPLC-QTOF-MS/MS-based analytical approach was employed to identify emodin-GSH adduct in bio-samples (Figure 8). Because more abundant mass information were acquired in negative mode than in positive mode, therefore we chose negative ion mode for structural characterization and identification. Consequently, the target ion ([M-H]-) was predefined at m/z 574. One extensive peak with m/z 574.11 at retention time of 30.8 min was observed in the extract ion spectra of bile samples. No such peak was detected either in the control group or in the plasma and urine of rats after intragastric gavage of PMR or emodin (data not shown).

FIGURE 8
www.frontiersin.org

FIGURE 8. The representative UHPLC-QTOF-MS/MS chromatograms. (A,B) The typical total ion current (TIC), extract ion (EI) at m/z 574 in rat blank bile; (C,D) the typical TIC, EI at m/z 574 and MS/MS spectra in bile of rat treated with PMR; (E,F) the typical TIC, EI at m/z 574 and MS/MS spectra in bile of rat treated with emodin.

Taking the advantage of QTOF-MS/MS, the accurate mass values of this peak were within the range of m/z 574.1109–574.1172, the errors of accurate mass of the tested samples were from 0.18 to 6.08 ppm compared with the exact mass of the target adduct (calculated mass 574.1137 for C25H24O11N3S). Furthermore, the characteristic fragmentation ions at m/z 301.01, 272.08, 269.04, 254.07, and 241.04 were observed from the MS/MS spectra of the precursor ion ([M-H]-) (Figures 8D,F). The product ions at m/z 301.01 and 272.08 were two complementary ions, corresponding to the moieties of sulfur-substituted emodin and γ-glutamyl-dehydroalanyl-glycine (Dieckhaus et al., 2005), respectively. The ions at m/z 269.04 and 241.04 were generated by successive neutral losses of sulfur atom (32 Da) and carbon monoxide (28 Da) from the ion at m/z 301.01, the ion at m/z 254.07 was ascribed to the loss of water (18 Da) from at m/z 272.08. Thus, the detected peak with m/z 574.11 was supposed as the emodin-GSH adduct. The detailed fragmentation pathway is shown in Figure 9.

FIGURE 9
www.frontiersin.org

FIGURE 9. The proposed fracture way of emodin-GSH adduct.

Significance of Emodin-GSH Adduct Detected in Bile of PMR-Induced Liver Injured Rats

Drugs are generally metabolized by a variety of chemical processes involving oxidation, reduction, and hydrolysis (phase I reactions) or glucuronidation, sulfation, acetylation, and methylation (phase II reactions) (Wilkinson, 2005). AQs underwent both phase I and phase II metabolism. Several phase I metabolites had been identified (Song et al., 2009), the phase II conjugates, especially glucuronides were the predominant existing forms of AQs in vivo (Wu et al., 2014). The phase II conjugation reaction was recognized as a detoxification step through increasing water solubility of AQs and enhancing their removal. As we know, GSH is a major intracellular antioxidant which is critical for preserving normal cellular redox balance and protecting the liver against oxidative stress (Chen et al., 2009). However, nucleophilic addition to GSH plays a key role in the detoxification reaction due to the nucleophilicity of GSH. In the present work, the detection of emodin-GSH adduct in bile of PMR-induced liver injured rats suggested that a portion of emodin underwent a crucial metabolic pathway known as epoxidation via hepatic drug-metabolizing enzymes, such as CYPs toward formation of epoxide intermediate (Peterson, 2013; He et al., 2015). The resulting epoxide subsequently attacked nucleophiles such as GSH, proteins, DNA and membranes that eventually induced liver injury.

Very recently, the hydroquinone intermediates of emodin were characterized in liver microsomes and the electrophilicity of emodin was determined (Qin et al., 2016), which supported the mechanism proposed in the present research that AQs-induced liver injury initiated by metabolic activation. And certainly, the specific metabolism of emodin in vivo especially phase I metabolic pathway associated with the toxic effects warranted further confirmation.

According to the histopathologic evaluations, inflammatory cell infiltration were observed in livers of PMR-treated rats at 1st and 2nd weeks, liver cell necrosis was observed after 3 weeks administration of PMR (Figure 5D), the emodin-derived GSH adduct could be detected from the 1st week during the administration of PMR, strongly suggesting that it had potential as an early sensitive and traceable biomarker exposure to PMR. Further study is to synthesize sufficient amounts of the chemical entity of emodin-GSH adduct that enables a dynamic quantification analysis to reflect the status and severity of injury. Notably, the emodin-GSH adduct was not detected in the plasma and urine, a similar observation was made in another research (Yang et al., 2014), this could be explained by the fact that the conjugated emodin was mainly excreted via bile (Bachmann and Schlatter, 1981) or the undetectable concentration of emodin-GSH adduct in blood or urine. On the basis of the detection of emodin-GSH adduct, the derivatization and detection of emodin-GSH conjugates or emodin-protein adducts in blood or urine deserve further investigation for clinical detection.

Conclusion

To date, there is no clear consensus concerning to the hepatotoxic ingredients of PMR as well as the underlying mechanism. Some studies suggested that the hepatotoxicity of PMR was probably idiosyncratic in patients (Jung et al., 2011; Dong et al., 2013). A recent research demonstrated that combined treatment with lipopolysaccharide and PMR resulted in acute idiosyncratic liver injury in rats, and revealed that 2,3,5,4′-tetrahydroxy cis-stilbene-2-O-β-glucoside was closely associated with the idiosyncratic hepatotoxicity of PMR (Li et al., 2017). On the other hand, PMR would cause chronic hepatitis in general (Ma, 2013). The liver injury occurred from several weeks to several months after PMR administration. Our previous toxicokinetic studies about an excessive consumption of PMR disclosed that the 2,3,5,4′-tetrahydroxy trans-stilbene-2-O-β-glucoside could inhibit the glucuronidation of emodin in the phase II metabolism and bring about a chronically accumulative emodin exposure, ultimately reaching a threshold and resulting in hepatotoxicity (Ma et al., 2013, 2015). In this continued study of PMR-associated hepatotoxicity, we herein put forward a hypothesis that AQs followed metabolic activation to generate the electrophilic reactive metabolites and subsequently formed covalent adduct with cellular nucleophiles in the liver to exert hepatotoxicity. First of all, the relationship between hepatotoxicity and GSH depletion was established, which demonstrated the important role of GSH in regulating the PMR hepatotoxicity. Then, emodin-derived GSH adduct was identified in bile of PMR-treated rats, suggesting that the metabolic activation of AQs toward epoxide formation might be one of the prime triggers for liver injury. The detection of emodin-GSH conjugation also indicated that the GSH conjugate showed a potential to be developed as a sensitive and traceable biomarker for the diagnosis of PMR-induced liver injury. Taken together, this study will not only shed novel insight into revealing the underlying mechanism of PMR-induced liver injury, but also offer a significant clue for hepatotoxic research of other AQs-containing herbs.

Author Contributions

H-JL and PL conceived and designed the experiments; L-LJ and QY performed the experiments; L-LJ and D-SZ analyzed the data; Y-XF and QY contributed reagents/materials/analysis tools; H-JL and L-LJ wrote the paper.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (No. 81573562, 81322051), the Natural Science Foundation of Jiangsu Province (No. BK20151442) and the Project Funded by the Priority Academic Program Development (PAPD) of Jiangsu Higher Education Institutions.

References

Abboud, D. G., and Kaplowitz, N. (2007). Drug-induced liver injury. Drug Saf. 30, 277–294. doi: 10.2165/00002018-200730040-00001

CrossRef Full Text | Google Scholar

Bachmann, M., and Schlatter, C. (1981). Metabolism of [14C] emodin in the rat. Xenobiotica 3, 217–225. doi: 10.3109/00498258109045294

CrossRef Full Text | Google Scholar

But, P. P., Tomlinson, B., and Lee, K. L. (1996). Hepatitis related to the Chinese medicine Shou-Wu-Pian manufactured from Polygonum multiflorum. Vet. Hum. Toxicol. 38, 280–282.

PubMed Abstract | Google Scholar

Cárdenas, A., Restrepo, J. C., Sierra, F., and Correa, G. (2006). Acute hepatitis due to Shen-Min: a herbal product derived from Polygonum multiflorum. J. Clin. Gastroenterol. 40, 629–632. doi: 10.1097/00004836-200608000-00014

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Ji, L. L., Wang, H. T., and Wang, Z. T. (2009). Intracellular glutathione plays important roles in pyrrolizidine alkaloids-induced growth inhibition on hepatocytes. Environ. Toxicol. Pharmacol. 28, 357–362. doi: 10.1016/j.etap.2009.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Chung, J. H., Seo, D. C., Chung, S. H., Lee, J. Y., and Seung, S. A. (1997). Metabolism and cytotoxicity of menadione and its metabolite in rat platelets. Toxicol. Appl. Pharmacol. 142, 378–385. doi: 10.1006/taap.1996.8048

PubMed Abstract | CrossRef Full Text | Google Scholar

Dieckhaus, C. M., Fernández-Metzler, C. L., King, R., Krolikowski, P. H., and Baillie, T. A. (2005). Negative ion tandem mass spectrometry for the detection of glutathione conjugates. Chem. Res. Toxicol. 18, 630–638. doi: 10.1021/tx049741u

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, H. H., Slain, D., Cheng, J. C., Ma, W. H., and Liang, W. F. (2013). Eighteen cases of liver injury following ingestion of Polygonum multiflorum. Complement. Ther. Med. 22, 70–74. doi: 10.1016/j.ctim.2013.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, P. P., Xia, Q., Lin, G., and Chou, M. W. (2004). Pyrrolizidine alkaloids–genotoxicity, metabolism enzymes, metabolic activation, and mechanisms. Drug Metab. Rev. 36, 1–55. doi: 10.1081/DMR-120028426

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Z. R. (2012). Toxicity risks and drug design. Prog. Pharm. Sci. 36, 1–13.

Google Scholar

He, L. N., Yang, A. H., Cui, T. Y., Zhai, Y. R., Zhang, F. L., Chen, J. X., et al. (2015). Reactive metabolite activation by CYP2C19-mediated rhein hepatotoxicity. Xenobiotica 45, 361–372. doi: 10.3109/00498254.2014.984794

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishihara, Y., Ishii, S., Sakai, Y., Yamamura, N., Onishi, Y., and Shimamoto, N. (2011). Crucial role of cytochrome P450 in hepatotoxicity induced by 2,3-dimethoxy-1,4-naphthoquinone in rats. J. Appl. Toxicol. 31, 173–178. doi: 10.1002/jat.1578

PubMed Abstract | CrossRef Full Text | Google Scholar

James, L. P., Mayeux, P. R., and Hinson, J. A. (2003). Acetaminophen-induced hepatotoxicity. Drug Metab. Dispos. 31, 107–116. doi: 10.1124/dmd.31.12.1499

PubMed Abstract | CrossRef Full Text | Google Scholar

Jung, K. A., Min, H. J., Yoo, S. S., Kim, H. J., Choi, S. N., Ha, C. Y., et al. (2011). Drug-induced liver injury: twenty five cases of acute hepatitis following ingestion of Polygonum multiflorum Thunb. Gut Liver 5, 493–499. doi: 10.5009/gnl.2011.5.4.493

PubMed Abstract | CrossRef Full Text | Google Scholar

Kouzi, S. A., Mcmurtry, R. J., and Nelson, S. D. (1994). Hepatotoxicity of germander (Teucrium chamaedrys L.) and one of its constituent neoclerodane diterpenes teucrin A in the mouse. Chem. Res. Toxicol. 7, 850–856. doi: 10.1021/tx00042a020

PubMed Abstract | CrossRef Full Text | Google Scholar

Laird, A. R., Ramchandani, N., Degoma, E. M., Avula, B., Khan, I. A., and Gesundheit, N. (2008). Acute hepatitis associated with the use of an herbal supplement (Polygonum multiflorum) mimicking iron-overload syndrome. J. Clin. Gastroenterol. 42, 861–862. doi: 10.1097/MCG.0b013e3181492515

PubMed Abstract | CrossRef Full Text | Google Scholar

Lei, X., Chen, J., Ren, J. T., Li, Y., Zhai, J. B., Mu, W., et al. (2015). Liver damage associated with Polygonum multiflorum Thunb.: a systematic review of case reports and case series. Evid. Based Complement. Alternat. Med. 2015:459749. doi: 10.1155/2015/459749

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Niu, M., Bai, Z., Zhang, C., Zhao, Y., Li, R., et al. (2017). Screening for main components associated with the idiosyncratic hepatotoxicity of a tonic herb, Polygonum multiflorum. Front. Med. 11:253–265. doi: 10.1007/s11684-017-0508-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C. L., Ma, J., Zheng, L., Li, H. J., and Li, P. (2012). Determination of emodin in L-02 cells and cell culture media with liquid chromatography-mass spectrometry: application to a cellular toxicokinetic study. J. Pharm. Biomed. Anal. 71, 71–78. doi: 10.1016/j.jpba.2012.07.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, G., Wang, J. Y., Li, N., Li, M., Gao, H., Ji, Y., et al. (2011). Hepatic sinusoidal obstruction syndrome associated with consumption of Gynura segetum. J. Hepatol. 54, 666–673. doi: 10.1016/j.jhep.2010.07.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, L. F., Lin, H. M., Yin, X. B., Zhao, Y., Xia, Z. W., Zhang, M., et al. (2015a). Characterization of the constituents in rat plasma after oral administration of Radix Polygoni Multiflori extracts by ultra-performance liquid chromatography/quadrupole time-of-flight mass spectrometry. Biomed. Chromatogr. 29, 1541–1547. doi: 10.1002/bmc.3456

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, L. F., Lin, H. M., Zhang, M., Ni, B., Yin, X. B., Qu, C. H., et al. (2015b). A novel method to analyze hepatotoxic components in Polygonum multiflorum using ultra-performance liquid chromatography-quadrupole time-of-flight mass spectrometry. J. Hazard. Mater. 299, 249–259. doi: 10.1016/j.jhazmat.2015.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, L. F., Ni, B., Lin, H. M., Zhang, M., Li, X. C., Yin, X. B., et al. (2015c). Traditional usages, botany, phytochemistry, pharmacology and toxicology of Polygonum multiflorum Thunb.: a review. J. Ethnopharmacol. 159, 158–183. doi: 10.1016/j.jep.2014.11.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Lv, G. P., Meng, L. Z., Han, D. Q., Li, H. Y., Zhao, J., and Li, S. P. (2015). Effect of sample preparation on components and liver toxicity of Polygonum multiflorum. J. Pharm. Biomed. Anal. 109, 105–111. doi: 10.1016/j.jpba.2015.02.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, J., Qi, L. W., Li, H. J., and Li, P. (2012). A segmental monitoring strategy based on variable wavelength detection for quality control of three Polygonaceae herbs. J. Pharm. Biomed. Anal. 25, 155–161. doi: 10.1016/j.jpba.2011.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, J., Zheng, L., Deng, T., Li, C. L., He, Y. S., Li, H. J., et al. (2013). Stilbene glucoside inhibits the glucuronidation of emodin in rats through the down-regulation of UDP-glucuronosyltransferases 1A8: application to a drug–drug interaction study in Radix Polygoni Multiflori. J. Ethnopharmacol. 147, 335–340. doi: 10.1016/j.jep.2013.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, J., Zheng, L., He, Y. S., and Li, H. J. (2015). Hepatotoxic assessment of Polygoni Multiflori Radix extract and toxicokinetic study of stilbene glucoside and anthraquinones in rats. J. Ethnopharmacol. 162, 61–68. doi: 10.1016/j.jep.2014.12.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, Z. J. (2013). The Preliminary Studies of Objectivity, Clinical Biomarkers and Mechanism of Liver Injury Induced by Polygoni Multiflori Radix. Chengdu: Chengdu University of Traditional Chinese Medicine.

Google Scholar

Mazzanti, G., Battinelli, L., Daniele, C., Mastroianni, C. M., Lichtner, M., Coletta, S., et al. (2004). New case of acute hepatitis following the consumption of Shou Wu Pian, a Chinese herbal product derived from Polygonum multiflorum. Ann. Intern. Med. 140:W30. doi: 10.7326/0003-4819-140-7-200404060-00042-w3

PubMed Abstract | CrossRef Full Text | Google Scholar

Panigrahi, G. K., Ch, R., Mudiam, M. K., Vashishtha, V. M., Raisuddin, S., and Das, M. (2015). Activity-guided chemo toxic profiling of Cassia occidentalis (CO) seeds: detection of toxic compounds in body fluids of CO-exposed patients and experimental rats. Chem. Res. Toxicol. 28, 1120–1132. doi: 10.1021/acs.chemrestox.5b00056

PubMed Abstract | CrossRef Full Text | Google Scholar

Panis, B., Wong, D. R., Hooymans, P. M., De Smet, P. A., and Rosias, P. P. (2005). Recurrent toxic hepatitis in a Caucasian girl related to the use of Shou-Wu-Pian, a Chinese herbal preparation. J. Pediatr. Gastr. Nutr. 41, 256–258. doi: 10.1097/01.MPG.0000164699.41282.67

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, J. H., Mann, S. P., and Meng, C. N. (2001). Acute hepatitis induced by Shou-Wu-Pian, a herbal product derived from Polygonum multiflorum. J. Gastroen. Hepatol. 16, 115–117. doi: 10.1046/j.1440-1746.2001.02309.x

CrossRef Full Text | Google Scholar

Peterson, L. A. (2013). Reactive metabolites in the biotransformation of molecules containing a furan ring. Chem. Res. Toxicol. 26, 6–25. doi: 10.1021/tx3003824

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, B. Y., Xu, Y., Chen, J. M., Huang, W. L., Peng, Y., and Zheng, J. (2016). Chemical reactivity of emodin and its oxidative metabolites to thiols. Chem. Res. Toxicol. 29, 2114–2124. doi: 10.1021/acs.chemrestox.6b00191

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruan, J., Yang, M., Fu, P., Ye, Y., and Lin, G. (2014). Metabolic activation of pyrrolizidine alkaloids: insights into the structural and enzymatic basis. Chem. Res. Toxicol. 27, 1030–1039. doi: 10.1021/tx500071q

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimizu, S., Atsumi, R., Nakazawa, T., Izumi, T., Sudo, K., Okazaki, O., et al. (2011). Ticlopidine-induced hepatotoxicity in a GSH-depleted rat model. Arch. Toxicol. 85, 347–353. doi: 10.1007/s00204-010-0594-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, R., Lin, H., Zhang, Z. J., Li, Z., Xu, L., Dong, H. J., et al. (2009). Profiling the metabolic differences of anthraquinone derivatives using liquid chromatography/tandem mass spectrometry with data-dependent acquisition. Rapid Commun. Mass Spctrom. 23, 537–547. doi: 10.1002/rcm.3907

PubMed Abstract | CrossRef Full Text | Google Scholar

Stadlbauer, V., Weiss, S., Payer, F., and Stauber, R. E. (2008). Herbal does not at all mean innocuous: the sixth case of hepatotoxicity associated with Morinda citrifolia (noni). Am. J. Gastroenterol. 103, 2406–2407. doi: 10.1111/j.1572-0241.2008.02010_8.x

CrossRef Full Text | Google Scholar

Tang, W., Stearns, R. A., Bandiera, S. M., Zhang, Y., Raab, C., Braun, M. P., et al. (1999). Studies on cytochrome P-450-mediated bioactivation of diclofenac in rats and in human hepatocytes: identification of glutathione conjugated metabolites. Drug Metab. Dispos. 7, 365–372.

PubMed Abstract | Google Scholar

Tu, C., Jiang, B. Q., Zhao, Y. L., Li, C. Y., Li, N., Li, X. F., et al. (2015). Comparison of processed and crude Polygoni Multiflori Radix induced rat liver injury and screening for sensitive indicators. China J. Chin. Mater. Med. 40, k654–660.

PubMed Abstract | Google Scholar

United States Food and Drug Administration (2001). Guidance for Industry: Bioanalytical Method Validation. Available at: http://www.fda.gov/cvm

Google Scholar

Wang, J. B., Zhao, Y. L., Xiao, X. H., Li, H. F., Zhao, H. P., Zhang, P., et al. (2009). Assessment of the renal protection and hepatotoxicity of rhubarb extract in rats. J. Ethnopharmacol. 124, 18–25. doi: 10.1016/j.jep.2009.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilkinson, G. R. (2005). Drug metabolism and variability among patients in drug response. New Engl. J. Med. 352, 2211–2221. doi: 10.1056/NEJMra032424

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, W. J., Yan, R., Yao, M. C., Zhan, Y., and Wang, Y. T. (2014). Pharmacokinetics of anthraquinones in rat plasma after oral administration of a rhubarb extract. Biomed. Chromatogr. 28, 564–572. doi: 10.1002/bmc.3070

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, B. H., Liu, W., Chen, K. X., Wang, Z. T., and Wang, C. H. (2014). Metabolism of diosbulbin B in vitro and in vivo in rats: formation of reactive metabolites and human enzymes involved. Drug Metab. Dispos. 42, 1737–1750. doi: 10.1124/dmd.114.058222

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, J., Xie, J., Ma, X. J., Wang, M. J., Li, N., Wang, J., et al. (2011). Hepatoxicity of major constituents and extractions of Radix Polygoni Multiflori and Radix Polygoni Multiflori Praeparata. J. Ethnopharmacol. 137, 1291–1299. doi: 10.1016/j.jep.2011.07.055

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., Yang, X. H., and Deng, Y. Y. (2009). Evaluation and consideration on safety information abroad of Polygonum multiflorum and its preparations. Chin. J. Chin. Mater. Med. 34, 2414–2418.

PubMed Abstract | Google Scholar

Zhang, R. C., Liu, B., Sun, Z. X., and Xu, D. Y. (2010). Effects of extract of Polygonum multiflorum on cell cycle arrest and apoptosis of human liver cell line L02. J. Chin. Integr. Med. 8, 554–561. doi: 10.3736/jcim20100608

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, D. S., Jiang, L. L., Fan, Y. X., Dong, L. C., Ma, J., Dong, X., et al. (2017). Identification of urine tauro-β-muricholic acid as a promising biomarker in Polygoni Multiflori Radix-induced hepatotoxicity by targeted metabolomics of bile acids. Food Chem. Toxicol. doi: 10.1016/j.fct.2017.02.030 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Polygoni Multiflori Radix, liver injury, glutathione depletion, reactive metabolite, emodin-glutathione adduct

Citation: Jiang L-L, Zhao D-S, Fan Y -X, Yu Q, Li P and Li H-J (2017) Detection of Emodin Derived Glutathione Adduct in Normal Rats Administered with Large Dosage of Polygoni Multiflori Radix. Front. Pharmacol. 8:446. doi: 10.3389/fphar.2017.00446

Received: 24 April 2017; Accepted: 22 June 2017;
Published: 06 July 2017.

Edited by:

Yeong Shik Kim, Seoul National University, South Korea

Reviewed by:

Jia-bo Wang, Beijing 302 Hospital of China, China
Maria Luisa Del Moral, Universidad de Jaén, Spain

Copyright © 2017 Jiang, Zhao, Fan, Yu, Li and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Hui-Jun Li, cpuli@163.com