Zhu-Jun Mao1*
Lu-Ping Qin- 1College of Pharmaceutical Sciences, Zhejiang Chinese Medical University, Hangzhou, China
- 2College of Basic Medical Sciences, Zhejiang Chinese Medical University, Hangzhou, China
Insulin resistance (IR) is likely to induce metabolic syndrome and type 2 diabetes mellitus (T2DM). Gluconeogenesis (GNG) is a complex metabolic process that may result in glucose generation from certain non-carbohydrate substrates. Chinese herbal medicine astragalus polysaccharides and berberine have been documented to ameliorate IR, and combined use of astragalus polysaccharide (AP) and berberine (BBR) are reported to synergistically produce an even better effect. However, what change may occur in the GNG signaling pathway of IR-HepG2 cells in this synergistic effect and whether AP-BBR attenuates IR by regulating the GNG signaling pathway remain unclear. For the first time, we discovered in this study that the optimal time of IR-HepG2 cell model formation was 48 h after insulin intervention. AP-BBR attenuated IR in HepG2 cells and the optimal concentration was 10 mg. AP-BBR reduced the intracellular H2O2 content with no significant effect on apoptosis of IR-HepG2 cells. In addition, a rapid change was observed in intracellular calcium current of the IR-HepG2 cell model, and AP-BBR intervention attenuated this change markedly. The gene sequencing results showed that the GNG signaling pathway was one of the signaling pathways of AP-BBR to attenuate IR in IR-Hepg2 cells. The expression of p-FoxO1Ser256 and PEPCK protein was increased, and the expression of GLUT2 protein was decreased significantly in the IR-HepG2 cell model, and both of these effects could be reversed by AP-BBR intervention. AP-BBR attenuated IR in IR-HepG2 cells, probably by regulating the GNG signaling Pathway.
Introduction
Diabetes mellitus is classified as type 1 (T1DM) and type 2 (T2DM). Clinically, more than 90% DM patients belong to T2DM. Insulin resistance (IR) is a pathological condition in which a target cell or a target organ fails to respond normally to insulin, so that insulin is unable to stimulate glucose disposal, which is a serious risk factor leading to the onset of T2DM. Other than T2DM, IR is forerunner for the development of many metabolic diseases, such as hypertension, obesity, and atherosclerosis. Therefore, ameliorating the IR status is an effective approach for the treatment of T2DM (Rodriguez-Gutierrez et al., 2018). IR is a state in which the response of target organs to insulin is decreased, or the normal dose of insulin secreted by pancreatic islet β cells produces less than normal biological effects. The IR theory has been generally accepted as the theoretical foundation for the study of T2DM by the medical community.
Our previous study (Zhu-jun et al., 2017)showed that combined use of astragalus polysaccharide (AP) and berberine (BBR) could synergistically produce a better therapeutic effect than either of them alone by promoting the basic secretion ability of IR-INS-1 cells, and this research achievement is on the way of applying for a patent in China (patent application no.: CN201810266347.9; publication/announcement no.: CN108478591A).
The liver is the main organ of glucose and lipid metabolism and plays an important role in energy metabolism of various tissues and organs including muscles and adipose tissues through regulation of insulin, glucagon, epinephrine, growth hormone, and other hormones (von Horn and Minor, 2019). The construction of an IR cell model is an important way to explore the pathological mechanism and drug development of IR-related diseases.
HepG2 cells originating from human hepatic embryonal tumor cells are a phenotype very similar to hepatocytes. It was reported that the number of insulin receptors on the surface of HepG2 cells decreased by 58% under the high-level insulin condition, and the remaining receptors then exhibited a 50% decrease in insulin internalization and degradation on a per receptor basis (Williams and Olefsky, 1990). The degree of such a decrease was positively correlated with the insulin level and the duration of stimulation, indicating that HepG2 is a good cell model to study IR pathogenesis and the mechanism of hypoglycemic drugs in vitro.
The results of gene sequencing in a previous study (Tsujimoto et al., 2015) showed that there were significant differences in glycolysis/gluconeogenesis pathway between the control and model groups, and between the model and drug administration groups. Gluconeogenesis (GNG) is a metabolic process of converting non-carbohydrates, such as amino acids, pyruvate, and glycerin, into glycogen or glucose, which is important for maintaining blood glucose levels within the physiological range.
Under normal physiological conditions, the liver is the main organ of GNG. It has rapid adaptation to metabolic shift between fed and fasting states through a reciprocal mechanism (Wahren and Ekberg, 2007): Hepatic gluconeogenic activity is suppressed in response to postprandial insulin secretion so as to limit glucose production. Conversely, hepatic GNG is stimulated in response to reduced insulin action and elevates glucagon secretion during fasting. The renal glycogenic ability is only 1/10 of the liver (Edgerton et al., 2009).GNG change is not obvious under normal dietary and metabolic conditions, but can undergo significant change in prolonged fasting or exercise and under low-carbohydrate diet conditions, or when there is a significant decrease in glycogen storage in the liver. GNG is an energy-consuming process that starts with two molecules of pyruvate and ends with the synthesis of one molecule of glucose, which requires six molecules of ATP/GTP. It produces a net 2ATP compared to glycolysis. In T2DM patients, the main source of endogenous glucose production is from GNG rather from glycogende composition. Increased hepatic GNG is a hallmark of T2DM, so maintaining the normal rate of GNG is the key to avoiding T2DM (Matsumoto et al., 2007).
However, what change may occur in the GNG signaling pathway of IR-HepG2 cells, and whether AP-BBR can attenuate IR by regulating the GNG signaling pathway remain unclear.
In the present study, we established a model of IR- HepG2 cells by insulin induction in vitro to determine the concentration and timing of stable insulin induction, and detect glucose uptake, reactive oxygen species (ROS) and apoptosis, CCK8 cell viability, the H2O2 concentration, and intracellular calcium ion in HepG2 cells to see whether AP-BBR attenuated IR by regulating the GNG signaling pathway.
Materials and Methods
Determination of Glucose Content in HepG2 Cells
HepG2 Cell Culture
HepG2 cells (Cell Bank of the Chinese Academy of Sciences, Shanghai, China) were cultured in 1640 medium (Hyclone, Beijing, China) containing 10% fetal bovine serum (FBS, Hyclone, Beijing, China) and 1× streptomycin in a 37°C 5% CO2 saturated humidity incubator. The normally cultured HepG2 cell lines in log phase were centrifuged at 100 grpm for 5 min, and 20000 cells/well were placed in a 96-well plate and incubated at 37°C.
Grouping and Drug Administration
The experiment was performed in 24-h control group, 24-h model group, 36-h control group, 36-h model group, 48-h control group, 48-h model group, 72-h control group and 72-h model group.
Insulin (Gibco, NY, USA) was diluted to a final concentration of 10-6 mol/L in complete medium. 200-µl insulin preparation was added into each well for the model group and an equal amount of complete medium was added into each well for the control group. Culture was performed in a 37°C 5% CO2 and saturated humidity incubator. The supernatant of the corresponding medium was collected according to the time point by centrifugation at 3000 r/min for 5 min and stored at −80°C for use.
Determination of the Glucose Content
The reagent for determining the glucose content (RSBIO, Shanghai, China) was balanced and configured at room temperature. 20-ml R1 reagent and 20-ml R2 reagent were mixed well. The EP tubes were marked as a blank tube, a calibration tube, a quality control tube, and a sample tube accordingly. 1000 µl working fluid was added to each tube; 10 µl distilled water was added to the blank tube; 10 µl calibration product was added to each tube; 10 µl quality control product was added to each tube of quality control tube; and 10 µl sample was added to each sample tube. Six samples were added in each group. After full mixing, the EP tubes were placed in 37°C water bath for 15 min. The 200-µl sample from each tube was transferred to a 96-well plate, and the absorbance was measured at 505 nm.
Detection of HepG2 Cell Viability by CCK8
HepG2 cell culture was conducted in the following groups: a normal control group; an insulin model group; a 1-mg AP-BBR group; a 5-mg AP-BBR group, a 10-mg AP-BBR group; a 20-mg AP-BBR group; and a 40-mg AP-BBR group. Astragalus polysaccharide (AP) used in this study was provided by Shanghai Yuanye Biotechnology Co., Ltd (Shanghai, China; Lot No. B20562; AP is a kind of macromolecular active substance extracted from the dried roots of the leguminous plant Astragalus membranaceus or Astragalus membranaceus. It is mainly composed of 75.19% glucose and a small amount of fructose, galactose, arabinose, and xylose. AP is a neutral polysaccharide that makes the iodine liquid blue and has a melting point above 200°C.), and BBR was provided by Beijing Century Aoko Biotechnology Co., Ltd (Beijing, China; Lot No. BWB50136; Molecular formula: C20H18NO4; Molecular weight: 336.37; Melting point: 145°C; boiling point: 354.2°C at 760 mm Hg). AP-BBR was administered at a 1:1 mass ratio of AP: BBR. The above ratio was set up according to our previous studies on AP and BBR (Zitai et al., 2018; Zhu-jun et al., 2018). After successful establishment of the IR model, different concentrations of AP-BBR were added. The system of each pore was 200 µl. The experiment was repeated three times independently. After different time points of drug action, 10 μl CCK8 detection solution (UNOCI, Hangzhou, China) was added to each pore, and the reaction time was 2 h at 37°C. A blank control group was set up and treated in the same way without anything added. Optical density (OD) was measured at 450 and 650 nm. The results were calculated using the following equation: Cell survival rate (%) = OD value of experimental group/OD value of non-drug group × 100%.
Detection of Apoptosis and ROS of HepG2 Cells by Flow Cytometry
Grouping
HepG2 cells were divided into three groups: a normal control group, a model group, and a 10-mg AP-BBR group. The original culture medium was retained in the wells to which 100 µl of the respective insulin concentration was added and treated for 48 h.
ROS Probe Incubation
The cry preserved tube containing 1.5-ml cells was taken out from the liquid nitrogen tank and quickly placed in water bath at 37°C for about 2 min. The cell suspension in the tube was transferred into a 15-ml centrifuge tube, to which 5 ml complete medium was added, and centrifuged at 300×g for 5 min at room temperature. After removing the supernatant, cells were re-suspended with a moderate amount of complete medium heavy precipitation, inoculated in a 10-cm petri dish, added with complete medium to 10 ml, and cultured again in 37°C and 5% CO2, saturated humidity. 1 μl ROS probe (Beyotime, Shanghai, China) was added to the resuscitated cells in the proportion of 1:1000, mixed, incubated at 37°C for 20 min, and oscillated several times per 5 min. After 5 centrifugations, cell precipitation was collected, 1 ml PBS was resuscitated, and centrifuged at 500 × g for 5 min. 1 ml PBS was re-suspended. A negative control group was set up and treated in the same way without adding the probe.
Apoptosis Probe Incubation
The prepared cells were precipitated; 0.1 ml banding buffer was resuscitated; FITC probe 5 μl and PI probe 10 μl were added, mixed, incubated at room temperature for 5 min, and finally tested after addition of the banding buffer to 0.5 ml. The negative control group was treated in the same way without anything added. The single staining group was treated in the same way with addition of the FITC probe or PI probe only.
Detection of ROS and the Apoptosis Rate by Flow Cytometry
The ROS signal was detected by flow cytometry (BIO-RAD, USA) and FL-1A channel. The proportion of ROS-positive cells and the mean fluorescence intensity of cells were recorded. FITC signal was detected by flow cytometry and FL-1A channel, PI signal was detected by FL-2A channel, and the apoptosis rate was recorded.
Determination of the Intracellular H2O2 Concentration in HepG2 Cells
Cells in each group were cultured according to “Methods 1.1”. The Krebs-Ringer buffer was prepared with double steamed water, and the pH was 7.5. The final concentration of Amplex Red was 100 mg, and the final concentration of HRP was 0.25 U/ml with 1× buffer. H2O2 standard was diluted by gradient, and the standard curve was drawn. After removing the original liquid from the 96-well plate and washing PBS, 50-μl working fluid was added to each cell and incubated away from light at 37°C for 30 min. The fluorescence intensity was measured at the excitation wavelength of 530 nm and the emission wavelength of 590 nm.
Detection of Intracellular Ca2+ in HepG2 Cells by Laser Confocal Technique
Hepg2 Cells Were Cultured and an IR-HepG2 Model Was Established
Fluo-4 probe loading was performed as follows: the prepared cell confocal dish was washed with 1 ml PBS buffer without Ca 2 + and Mg 2 +; after addition of the Fluo-4 probe at 1:1 000 and thorough mixing, cells were incubated at 37°C for 30 min, and then washed with PBS without Ca 2 + and Mg 2 + 3 times before testing.
Administration and Detection of Intracellular Ca2+ by Laser Confocal Microscopy
The working module of living cells was opened in advance, and the conditions were set as follows: 37°C, 5% CO 2, saturated humidity. The incubated cell confocal dish samples were placed in the workbench, the control samples and the model samples were put immediately under 100 times the laser confocal microscope x 100 (OLYMPUS FV1000, Japan). The blue light channel (excitation wavelength 488 nm and emission wavelength 526 nm) was scanned by laser every 5 s, and the fluorescence intensity was recorded. The control group and model samples were taken from the blue light channel. After dripping AP-BBR with the final concentration of 10 mg, the fluorescence intensity was recorded, and the control group and the model sample were taken from the blue light channel (488 nm excitation wavelength and 526 nm emission wavelength). Immediately, the blue light channel (488 nm excitation wavelength and 526 nm emission wavelength) was scanned by laser every 5 s, and the fluorescence intensity was recorded. Cells were photographed continuously for 0 to 1,000 s.
Transcriptome Sequencing and Pathway Enrichment Analysis
RNA isolation, purification, and quantification: Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample were quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0.
cDNA Library Construction: Poly(A) RNA was purified from total RNA (5 µg) using poly-T oligo-attached magnetic beads using two rounds of purification. Then, the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which was subsequently used to synthesize U-labeled second-stranded DNAs with Escherichia coli DNA polymerase I, RNase H, and dUTP. An A-base was then added to the blunt ends of each strand, which were prepared for ligation to the indexed adapters. Each adapter contained a T-base overhang for ligating the adapter to the A-tailed fragmented DNA. Single- or dual-index adapters were ligated to the fragments, and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, the ligated products were amplified by PCR under the following conditions: initial denaturation at 95°C for 3 min, 8 cycles of denaturation at 98°C for 15 s, annealing at 60°C for 15 s, extension at 72°C for 30 s, and then final extension at 72°C for 5 min. The mean insert size for the final cDNA library was 300 bp (±50 bp). Finally, 150-bp paired-end sequencing was performed on an Illumina Hiseq 4000 (LC Bio, China) following the vendor's recommended protocol.
Pathway enrichment analysis: Using the DAVID database and mouse genome as background control, the differentially expressed genes were analyzed by gene ontology under “FunctionalAnnotation Chart” functional module. The differentially expressed genes were divided into three categories according to their functions: the biological process, the cell component, and the molecular function. Pathway analysis was carried out by KEGG analysis function.
Analysis of Protein Expression
HepG2 cells were lysed in NP 40 lysis buffer and centrifuged at 12,000×g for 20 min at 4°C. The protein concentration was quantitated, and the protein lysate (20 μg proteins) was subjected to SDS-PAGE on 10% polyacrylamide gels and then electrotransferred onto PVDF membranes (Millipore, Billerica, MA, USA). The membranes were then blocked by 5% fat-free milk, and incubated respectively with FoxO1 antibody (78 kDa; 1:1000), Phospho-FoxO1 (Ser256) antibody (78 kDa; 1:1000), PEPCK (PCK2) antibody (71 kDa; 1:1000), GLUT2 antibody (53 kDa; 1:1000), and anti-GAPDH (36 kDa; 1:2000) (CST, Cambridge, MA, U.S.) at 4°C overnight. Horseradish peroxidase (HRP)-conjugated secondary antibodies (rabbit or mouse) were used to detect bound antibodies, which were visualized by enhanced chemiluminescence (ECL, Millipore, Billerica, MA, USA) on a VersaDoc 4000 MP (BIO-RAD, USA) workstation. Using the ImageJ software, densitometric analysis was performed to determine the relative expression level of the target protein, which was standardized to GAPDH (Bethesda, Maryland. USA).
Statistical Analysis
Data concerning miRNA expression were analyzed with independent sample t-test. All the other data were analyzed with one-way ANOVA followed by LSD when equal variances assumed, and Dunnett's-T3 when equal variances not assumed. The results are expressed as the means ± SEM. All statistical analyses were performed by using SPSS23. p < 0.05 was considered statistically significant. Figures were drawn by GraphPad Prism 6.02 and Adobe Illustrator CS6.
Results
Determination of Glucose Content in HepG2 Cells
The results of intracellular glucose content determination in HepG2 cells are shown in Figure 1A. The best time for the HepG2 insulin model was at 48 h of insulin intervention.
Figure1 The glucose concentration of HepG2 cells (A). Data are shown by mean ± SEM (n = 3), 48 h control group vs. 48 h model group, ΔΔP < 0.01. Proliferation activity of HepG2 cells (B), AP-BBR group vs. model group, **P < 0.01. The ROS content in HepG2 cells (C), model group vs. control group, ΔΔP < 0.01, AP-BBR group vs. model group, **P < 0.01. Apoptosis rate of HepG2 cells (D), model group vs. control group, P > 0.05, AP-BBR group vs. model group, P > 0.05.
Result of the Cell Survival Rate of HepG2 Cells
CCK8 assay was performed to detect the cell survival rate of HepG2 cells. As shown in Figure 1B, the survival rate of IR-HepG2 cells was increased after AP-BBR administration. The optimal dosage of AP-BBR was 10 and 20 mg. For convenience, 10 mg was used in the subsequent experiments.
Determination of ROS Content in HepG2 Cells
ROS content in HepG2 cells was detected by flow cytometry (Figure 1C). The results showed that the level of ROS in IR-HepG2 cells increased significantly, and the level of ROS decreased significantly in AP-BBR group. The results showed that the level of ROS in IR-HepG2 cells increased significantly, and AP-BBR decreased the level of ROS in the IR-HepG2 cell model.
The Apoptosis Rate of HepG2 Cells
The effect of AP-BBR on apoptosis of IR-HepG2 cells was detected by flow cytometry, and the results are shown in Figure 1D. The results showed that the IR model and AP-BBR treatment had no significant effect on apoptosis of HepG2 cells.
H2O2 Concentration in HepG2 Cells
The effect of AP-BBR on the concentration of H2O2 in IR-HepG2 cells was detected. As shown in Figure 2, the H2O2 content of IR-HepG2 cells was significantly increased, and the H2O2 content of IR-HepG2 cells was significantly decreased in AP-BBR group.
Figure 2 H2O2 standard curve (A). X axis: concentration (mM); Y axis: fluorescence OD value. Four-parameter Logistic curve fitting, equation: y = (A -D)/[1 + (x/C)^B] + D.A = 274.48194; B = -0.92136; C = 0.01225; D = -1.30186; r^2 = 0.99990. H2O2 concentration in HepG2 cells (B). Data are shown by mean ± SEM (n = 3). ΔΔP < 0.01 showing a significant difference as compared with control group; **P < 0.01 showing a critical difference as compared with model group (n = 3).
Change in Intracellular Ca2+ in HepG2 Cells
The photos of each HepG2 cell group under the laser confocal microscope are shown in Figure 3A.
The slope of intracellular Ca2+ fluorescence intensity dynamic (0–150 s) curve fitting is shown in Figure 3B. The results suggest that change in intracellular calcium flow became stronger and faster after modeling, and slower and weaker after AP-BBR administration.
The experimental results of laser confocal detection of intracellular Ca2+ change in HepG2 cells are shown in Figure 3C (the raw data is shown in Supplementary Table 1). These results suggest that the intracellular calcium flow changed rapidly in IR-HepG2 cells, and the intracellular calcium flow change was decreased after 10 mg AP-BBR administration.
Figure 3 The photos of each HepG2 cells group under laser confocal microscope (×400) (A). Calculation of slope by curve fitting of intracellular Ca2+ fluorescence intensity dynamics (0–150 s) (B). Changes of intracellular Ca2+ fluorescence intensity in HepG2 cells(0–1000 s) (C).
Results of Pathway Enrichment Analysis
The statistical enrichment of differentially expressed genes in KEGG pathway was analyzed by local perl script and R software package ggplot2. Red indicates a significant difference, and blue indicates no significant difference. The size of the circle reflects the number of genes. Among the top 20 rich factors, pathways related to IR included the PPAR pathway, glucose metabolism pathway, and glycolysis/gluconeogenesis pathway. The results of significant differential expression analysis showed significant differences in the effect of glycolysis/GNG and other pathways on IR (Figure 4A). After AP-BBR intervention, there was significant difference expression in Glycolysis/GNG pathway (Figure 4B).
Figure 4 Enrichment of differentially expressed genes in KEGG pathway. Model group vs. control group (A). AP-BBR group vs. model group (B).
Protein Expressions
The results of Western blot in each HepG2 cell group are shown in Figure 5. Compared with the control group, the protein expression of p-FoxO1Ser256 and PEPCK in the model group was significantly increased, while the protein expression of FoxO1 and GLUT2 was significantly decreased (n = 3). Compared with the model group, the protein expression of p-FoxO1Ser256 and PEPCK in AP-BBR group was significantly decreased, while the protein expression of FoxO1 and GLUT2 was significantly increased (n = 3).
Figure 5 The protein expression in HepG2 cells. The protein expressions of FoxO1, p-FoxO1Ser256, PEPCK, and GLUT2 were detected by Western blot. The figure represents one of three experiments with similar results (A). In (B), ΔΔP < 0.01 vs. control group, **P < 0.01, vs. model group (n = 3).
Discussion
Deficiency of insulin signaling and IR in the liver has a great impact on energy balance and metabolism. Liver IR can lead to a series of metabolic abnormalities, such as hyperglycemia, dyslipidemia, and increased secretion of inflammatory factors. The inhibitory effect of insulin on hepatic glucose production is impaired during liver IR. When IR progresses to T2DM, the lack of insulin signaling molecules is common as hyperglycemia. Insulin signaling can impact T2DM by regulating glucolipid metabolism and energy homeostasis via action on the liver, skeletal muscle, and adipose tissue. It is reported that disruption of insulin signaling in the moue liver, skeletal muscle or adipose tissue can lead to diabetes (Boucher et al., 2014). Mice lacking insulin receptor substrate-2 (IRS-2) in the liver developed severe IR and hyperglycemia (Williams and Olefsky, 1990). All in all, according to the data of human and animal research, liver IR leads to the development of hyperglycemia and T2DM.
Astragalus membranaceus is one of the most widely used traditional Chinese herbal medicines. It can be used as an antidiabetic agent. The main component of Astragalus membranaceus is AP, which can improve overall IR (Sun et al., 2017). Coptis chinensis has been used as an anti-diabetic drug in Traditional Chinese Medicine (TCM) for centuries, and BBR is the main bioactive alkaloid of Coptis chinensis and can moderate T2DM by alleviating IR (Abd et al., 2013).
AP attenuated IR mainly by decreasing the phosphorylation of protein kinase and inhibiting GLUT4 translocation in KKAy mice (Zhang and Pugliese, 2015). AP also promoted the expression of phosphorylated PI3Kp85 and GLUT4 by inhibiting the expression of JNK and phosphorylating insulin receptor substrate-1 (Jiang et al., 2015). In addition, AP decreased serum TNF-α, MCP-1, and ICAM-1 and inhibited the expression of NF- κ B mediated inflammation gene (Gui et al., 2013).
BBR attenuated IR mainly by inhibiting the Rho GTPase signaling pathway to reduce renal inflammation in diabetic rats (Xie et al., 2013), activating various human cell lines and their promoter-dependent protein kinase C-(PKC-), and increasing insulin receptor mRNA and protein expression (Pang et al., 2015). It also inhibited liver GNG of the LKB1-AMPK-TORC2 signaling pathway and improved LPS-induced β-cell damage through LPS4-dependent JNK/NF-kB signaling pathway (Jiang et al., 2015). In addition, BBR exerted its anti-inflammatory effect through insulin sensitization, reduced cytokine production and serine phosphorylation, and increased insulin-mediated tyrosine phosphorylation (Lou et al., 2011).
The liver is the central site of glucose synthesis and metabolism, and also the main target organ for IR (Reinke and Asher, 2016). Therefore, it is very important to establish a stable and reliable IR hepatocyte model for studying the molecular mechanism of IR in vitro and screening drugs for prevention and treatment of IR. HepG2 cells originating from human hepatic embryonal tumor cells have similar morphological and biological functions to normal human hepatocytes, and therefore, are most commonly used for studying hepatic IR. Insulin has been widely recognized as an inducer to establish IR cell models.
In the present study, we explored the conditions in establishing the IR- HepG2 cell model and the effect of AP-BBR on IR. The glucose content in the supernatant of HepG2 cells treated with insulin for 24, 36, 48, and 72 h was measured. It was found that the optimal time for modeling was at 48 h of insulin intervention. CCK8 was used to detect the effect of AP-BBR on the viability of IR-HepG2 cells after 48 h. The results showed that AP-BBR attenuated IR in HepG2 cells.
Oxidative stress in hepatocytes refers to the pathological process of hepatocyte injury caused by excessive or/and endogenous antioxidant capacity of ·OH, H2O2, and other active ROS, and disturbed balance between the oxidant and antioxidant systems (Cortes-Rojo and Rodriguez-Orozco, 2011). Oxidative stress causes lipid peroxidation damage, DNA oxidative damage, and abnormal protein expression in the liver by producing excessive ROS. It is involved in the pathogenesis of viral hepatitis, non-alcoholic/alcoholic hepatitis, drug-induced liver disease, and other liver diseases (Souli et al., 2015; Yu et al., 2015). When oxides or alkylates enter cells, they induce the production of a large amount of ROS. If ROS is not cleared in time, free radicals, the intermediate products of ROS, can directly act on nucleic acid, resulting in base modification and DNA strand breaks, resulting in oxidative DNA damage (Cadet et al., 2011; Rinna et al., 2015). Alterations in ROS and apoptosis of HepG2 cells were detected by flow cytometry. The results showed that the level of ROS increased significantly in IR-HepG2 cells and decreased after AP-BBR treatment. However, AP-BBR treatment had no significant effect on apoptosis of HepG2 cells.
Ca2+ is an important second messenger in cells. There is a close interaction between the ROS signaling system and the calcium signaling system in the intracellular environment. On the one hand, ROS can change the amplitude and temporal and spatial characteristics of local and whole-cell calcium signals by modifying the key protein components of the calcium signaling system, such as voltage-dependent calcium and intracellular calcium release channels. On the other hand, Ca2+ can accelerate mitochondrial metabolism by regulating the activities of dehydrogenase in tricarboxylic acid cycle (Rizzuto and Pozzan, 2006), ATP synthase, and adenine nucleoside transporter so as to accelerate mitochondrial metabolism and increase ROS production. [Ca2+]i can adjust mitochondrial proton driving force (Δ ψ m and Δ PH) and affect mitochondrial energy conversion and ROS production. Ca2+ is also involved in the regulation of cell cycle, cell differentiation, muscle contraction, cell movement, exocytosis, endocytosis, chemical chemotaxis, synaptic transmission and synaptic plasticity, and apoptosis (Liu and O'Rourke, 2009). [Ca2+]i is both a survival signal and a death signal. [Ca2+]i in cells is harmful to cells when it is too high or too low, so cell [Ca2+]i is under very strict regulation and control. The dynamic change of cell [Ca2+]i is the core link in cell signal transduction and regulatory network (Feske et al., 2001). Physiological levels of ROS and Ca2+ to adapt to the changes in energy requirements continue to adjust adaptively, so as to regulate the dynamic balance of mitochondrial morphology and distribution.
When the regulation of Ca2+ and ROS at physiological level fails to adapt to the change of energy demand, the abnormal level of Ca2+ and ROS leads to mitochondrial fusion and division, thus disrupting the balance of dynamic distribution, causing dysfunction of mitochondrial oxidative phosphorylation function and further affecting the cell function. The body is in a pathological state, and at the same time, mitochondrial oxidative phosphorylation dysfunction further leads to abnormal levels of Ca2+ and ROS, forming a vicious circle. ROS plays an important role in regulating cell proliferation, differentiation, apoptosis and cell senescence. Low concentrations of H2O2 can be used as intracellular signaling molecules to initiate gene transcription, thus inducing cell growth. High concentrations of ROS oxidized lipid, protein, and DNA can cause damage to cell integrity (Scherz-Shouval and Elazar, 2007).
Apoptosis induced by H2O2 is closely related to calcium overload. The molecular mechanism of apoptosis induced by oxidative stress is very complex and involves many signal transduction pathways, including the classical mitochondrial pathway (Chiara et al., 2012), death receptor pathway (Gomez-Quiroz et al., 2008) and endoplasmic reticulum pathway (Ip et al., 2011). Oxidative stress damage of cells leads to the release of a large amount of Ca2+ from the endoplasmic reticulum accompanied by endoplasmic reticulum stress, leading to unfolded protein reaction (UPR), which is used to reconstruct the normal function of endoplasmic reticulum (Greotti et al., 2019). If the intracellular Ca2+ concentration continues to increase and endoplasmic reticulum stress lasts for a prolonged time or severely, Ca2+-dependent kinases and phosphatase, such as calpain, caspase-12, and caspase-3 cascade reactions, will be activated, eventually leading to apoptosis (Das et al., 2006). In our experiment, change in the intracellular H2O2 concentration in IR-HepG2 cells was detected after AP-BBR treatment. The results showed that the H2O2 content of IR-HepG2 cells increased after modeling and decreased after AP-BBR treatment. The change in intracellular calcium flow in IR-HepG2 cells was rapid and reduced after AP-BBR treatment. These results suggest that Ca2+ was activated in the IR HepG2 cell model, and AP-BBR could inhibit the activation of Ca2+ to a certain extent.
The IR-INS-1 cell model in our previous study (Zhu-jun et al., 2017) showed that combined use of AP and BBR could promote basic secretion of IR-INS-1 cells significantly and showed a better therapeutic effect than that when either of them was used alone. The results showed that the compatibility of the two drugs had an optimal synergistic effect, and that AP-BBR could promote the proliferation of IR-HepG2 cells and reduce intracellular calcium flow, though the underlying mechanism remains unclear at present.
Insulin is a key hormone that suppresses GNG. It plays a vital role in the regulation of hepatic glucose output. Insulin regulates hepatic GNG by transcriptional modulation and activation of the relative signaling pathways (Hatting et al., 2018). Even an indirect effect of insulin on extra hepatic tissues can be sufficient to maintain normal glucose metabolism (O-Sullivan et al., 2015). T2DM is characterized by excess glucagon and IR. In T2DM patients, insulin could not suppress glucagon secretion or hepatic glucose production, which underscores the importance of determining the GNG rate in the potential treatment for T2DM.
Forkhead box protein O1 (FoxO1) is a member of the forehead box transcription factor class O (FOXO) family, and a direct transcriptional regulator of GNG. It is regulated by insulin reciprocally and plays a major role in regulating insulin sensitivity (Tikhanovich et al., 2013). The activity of FOXO1 as a regulator of blood glucose is modulated by several mechanisms including phosphorylation, acetylation, and deacetylation balance, as well as some other more novel mechanisms (Zhou et al., 2011).
In the high blood glucose state, insulin is released from the pancreas into the blood stream. Insulin signaling causes the activation of PI3kinase and the subsequent production of PIP3 phosphorylates Akt. Akt then phosphorylates FoxO1 to its nuclear exclusion, and is ubiquitinated and degraded by the proteosome (Matsuzaki et al., 2003). That is to say, FOXO1 is a counter of insulin-induced Akt activation in glycogen synthesis-GNG balance, which can subsequently suppress GNG during the fed state (Daitoku and Fukamizu, 2007).
Phosphoenolpyruvate carboxykinase (PEPCK) is an enzyme in the lyase family used in the metabolic pathway of GNG. It converts oxaloacetate into phosphoenolpyruvate and carbon dioxide (Mendez-Lucas et al., 2014).
The transcriptional induction of gluconeogenic enzyme PEPCK-C is an irreversible step of GNG, which therefore is of great importance in glucose homeostasis, as evidenced by the overexpression of PEPCK-C in T2DM laboratory mice. The role that PEPCK-C plays in GNG may be mediated by the citric acid cycle, the activity of which was found to be directly related to PEPCK-C abundance. Previous studies suggested that PEPCK-C level alone was not highly correlated with GNG in the mouse liver (Marandel et al., 2019). While the mouse liver almost exclusively expresses PEPCK-C, humans equally present a mitochondrial isozyme (PEPCK-M). Therefore, the role of PEPCK-C and PEPCK-M in GNG may be more complex and involve more factors than was previously believed.
Glucose metabolism depends on the uptake of glucose by cells. However, glucose is not free to enter cells through the lipid bilayer structure of the cell membrane. Glucose uptake by cells requires the help of glucose transporter (glucose transporters) on the cell membrane, which is known as the glucose transporter (GLUT) transport function.
GLUTs facilitate the movement of glucose from the plasma membrane across cell membranes. Member 2 (SLC2A2) of Glucose transporter 2 (GLUT2) is a principal transmembrane transporter for transfer of glucose between the liver and blood (Gould et al., 1991). It is of great significance for hepatic glucose uptake. Unlike GLUT4, GLUT2 does not rely on insulin for facilitation of glucose diffusion.
GLUT2 has low affinity (high Km, ca. 15-20 mM) for glucose, but it has high capacity for glucose. What's more, GLUT2 isoform is found abundant in the liver and pancreatic β-cells, so it has been suggested to function as a glucose sensor (Jung et al., 2006). GLUT2 may play a minor role in human β-cells, and can be an efficient glucose carrier so as to improve glucose and insulin metabolism.
In our study, the protein expression of p-FoxO1Ser256 and PEPCK in IR-HepG2 cells was increased, but the protein expression of FoxO1 and GLUT2 was decreased, and AP-BBR could reverse the above protein expression changes. The results suggest that AP-BBR attenuates IR in IR-HepG2 cells probably via the regulation of the GNG signaling Pathway.
Data Availability Statement
The datasets generated for this study can be found in the NCBI GEO https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE139929.
Author Contributions
Z-JM generated the data, wrote the manuscript, and approved its final version. Z-JM and L-PQ designed the study. Z-JM, ML, and XZ carried out the experiments. All authors have read and approved the final manuscript.
Funding
This study was supported by a grant from the National Natural Science Foundation of China (No. 81603351).
Conflict of Interest
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.
Abbreviations
IR, Insulin resistance; T2DM, type 2 diabetes mellitus; GNG, Gluconeogenesis; IR-HepG2, Insulin resistance HepG2 cells; AP, Astragalus polysaccharide; BBR, Berberine; AP-BBR, Astragalus polysaccharide plus berberine; ROS, reactive oxygen species; FoxO1, Forkhead box protein O1; PEPCK, Phosphoenolpyruvate carboxykinase; GLUT2, Glucose transporte; 2PPAR, Peroxisome proliferators-activated receptor; IRS-2, Insulin receptor substrate-2; GLUT4, Glucose transporter 4; UPR, Unfolded protein reaction.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2019.01508/full#supplementary-material
Supplementary Table 1 | The raw data of intracellular Ca2+ fluorescence intensity.
References
Abd, E. A., Ghareeb, D. A., Sarhan, E. E., Abu-Serie, M. M., El, D. M. (2013). In vitro biological assessment of Berberis vulgaris and its active constituent, berberine: antioxidants, anti-acetylcholinesterase, anti-diabetic and anticancer effects. BMC Complement Altern. Med. 13, 218. doi: 10.1186/1472-6882-13-218
Boucher, J., Kleinridders, A., Kahn, C. R. (2014) Insulin receptor signaling in normal and insulin-resistant states. Cold Spring Harb. Perspect. Biol. 6, 1. doi: 10.1101/cshperspect.a009191
Cadet, J., Douki, T., Ravanat, J. L. (2011). Measurement of oxidatively generated base damage in cellular DNA. Mutat. Res. 711(1–2), 3–12. doi: 10.1016/j.mrfmmm.2011.02.004
Chiara, F., Gambalunga, A., Sciacovelli, M., Nicolli, A., Ronconi, L., Fregona, D., et al. (2012). Chemotherapeutic induction of mitochondrial oxidative stress activates GSK-3alpha/beta and Bax, leading to permeability transition pore opening and tumor cell death. Cell Death Dis. 3, e444. doi: 10.1038/cddis.2012.184
Cortes-Rojo, C., Rodriguez-Orozco, A. R. (2011). Importance of oxidative damage on the electron transport chain for the rational use of mitochondria-targeted antioxidants. Mini Rev. Med. Chem. 11, 625–632. doi: 10.2174/138955711795906879
Daitoku, H., Fukamizu, A. (2007). FOXO transcription factors in the regulatory networks of longevity. J. Biochem. 141, 6. doi: 10.1093/jb/mvm104
Das, A., Banik, N. L., Ray, S. K. (2006). Mechanism of apoptosis with the involvement of calpain and caspase cascades in human malignant neuroblastoma SH-SY5Y cells exposed to flavonoids. Int. J. Cancer 119, 11. doi: 10.1002/ijc.22228
Edgerton, D. S., Johnson, K. M., Cherrington, A. D. (2009). Current strategies for the inhibition of hepatic glucose production in type 2 diabetes. Front. Biosci. (Landmark Ed) 14, 1169–1181. doi: 10.2741/3301
Feske, S., Giltnane, J., Dolmetsch, R., Staudt, L. M., Rao, A. (2001). Gene regulation mediated by calcium signals in T lymphocytes. Nat. Immunol. 2, 4. doi: 10.1038/86318
Gomez-Quiroz, L. E., Factor, V. M., Kaposi-Novak, P., Coulouarn, C., Conner, E. A., Thorgeirsson, S. S. (2008). Hepatocyte-specific c-Met deletion disrupts redox homeostasis and sensitizes to Fas-mediated apoptosis. J. Biol. Chem. 283, 21. doi: 10.1074/jbc.M707733200
Gould, G. W., Thomas, H. M., Jess, T. J., Bell, G. I. (1991). Expression of human glucose transporters in Xenopus oocytes: kinetic characterization and substrate specificities of the erythrocyte, liver, and brain isoforms. Biochemistry-US 30, 5139–5145. doi: 10.1021/bi00235a004
Greotti, E., Capitanio, P., Wong, A., Pozzan, T., Pizzo, P., Pendin, D. (2019). Familial Alzheimer's disease-linked presenilin mutants and intracellular Ca(2+) handling: A single-organelle, FRET-based analysis. Cell Calcium. 79, 44–56. doi: 10.1016/j.ceca.2019.02.005
Gui, D., Huang, J., Guo, Y., Chen, J., Chen, Y., Xiao, W., et al. (2013). Astragaloside IV ameliorates renal injury in streptozotocin-induced diabetic rats through inhibiting NF-kappaB-mediated inflammatory genes expression. Cytokine 61, 3. doi: 10.1016/j.cyto.2013.01.008
Hatting, M., Tavares, C. D. J., Sharabi, K., Rines, A. K., Puigserver, P. (2018). Insulin regulation of gluconeogenesis. Ann. N Y Acad. Sci. 1411, 21–35. doi: 10.1111/nyas.13435
Ip, S. W., Wu, S. Y., Yu, C. C., Kuo, C. L., Yu, C. S., Yang, J. S., et al. (2011). Induction of apoptotic death by curcumin in human tongue squamous cell carcinoma SCC-4 cells is mediated through endoplasmic reticulum stress and mitochondria-dependent pathways. Cell Biochem. Funct. 29, 8. doi: 10.1002/cbf.1800
Jiang, S., Wang, Y., Ren, D., Li, J., Yuan, G., An, L., et al. (2015). Antidiabetic mechanism of Coptis chinensis polysaccharide through its antioxidant property involving the JNK pathway. Pharm. Biol. 53, 7. doi: 10.3109/13880209.2014.952838
Jiang, S. J., Dong, H., Li, J. B., Xu, L. J., Zou, X., Wang, K. F., et al. (2015). Berberine inhibits hepatic gluconeogenesis via the LKB1-AMPK-TORC2 signaling pathway in streptozotocin-induced diabetic rats. World J. Gastroenterol. 21, 25. doi: 10.3748/wjg.v21.i25.7777
Jung, U. J., Lee, M. K., Park, Y. B., Kang, M. A., Choi, M. S. (2006). Effect of citrus flavonoids on lipid metabolism and glucose-regulating enzyme mRNA levels in type-2 diabetic mice. Int. J. Biochem. Cell Biol. 38, 1134–1145. doi: 10.1016/j.biocel.2005.12.002
Liu, T., O'Rourke, B. (2009). Regulation of mitochondrial Ca2+ and its effects on energetics and redox balance in normal and failing heart. J. Bioenerg. Biomembr. 41, 2. doi: 10.1007/s10863-009-9216-8
Lou, T., Zhang, Z., Xi, Z., Liu, K., Li, L., Liu, B., et al. (2011). Berberine inhibits inflammatory response and ameliorates insulin resistance in hepatocytes. Inflammation 34, 6. doi: 10.1007/s10753-010-9276-2
Marandel, L., Kostyniuk, D. J., Best, C., Forbes, J., Liu, J., Panserat, S., et al. (2019). Pck-ing up steam: Widening the salmonid gluconeogenic gene duplication trail. Gene 698, 129–140. doi: 10.1016/j.gene.2019.02.079
Matsumoto, M., Pocai, A., Rossetti, L., Depinho, R. A., Accili, D. (2007). Impaired regulation of hepatic glucose production in mice lacking the forkhead transcription factor Foxo1 in liver. Cell Metab. 6, 208–216. doi: 10.1016/j.cmet.2007.08.006
Matsuzaki, H., Daitoku, H., Hatta, M., Tanaka, K., Fukamizu, A. (2003). Insulin-induced phosphorylation of FKHR (Foxo1) targets to proteasomal degradation. Proc. Natl. Acad. Sci. U.S.A. 100, 20. doi: 10.1073/pnas.1934283100
Mendez-Lucas, A., Hyrossova, P., Novellasdemunt, L., Vinals, F., Perales, J. C. (2014). Mitochondrial phosphoenolpyruvate carboxykinase (PEPCK-M) is a pro-survival, endoplasmic reticulum (ER) stress response gene involved in tumor cell adaptation to nutrient availability. J. Biol. Chem. 289, 32. doi: 10.1074/jbc.M114.566927
O-Sullivan, I., Zhang, W., Wasserman, D. H., Liew, C. W., Liu, J., Paik, J., et al. (2015). Erratum: FoxO1 integrates direct and indirect effects of insulin on hepatic glucose production and glucose utilization. Nat. Commun. 6, 7681. doi: 10.1038/ncomms8861
Pang, B., Zhao, L. H., Zhou, Q., Zhao, T. Y., Wang, H., Gu, C. J., et al. (2015). Application of berberine on treating type 2 diabetes mellitus. Int. J. Endocrinol. 2015, 12. doi: 10.1155/2015/905749
Reinke, H., Asher, G. (2016). Circadian clock control of liver metabolic functions. Gastroenterology 150, 3. doi: 10.1053/j.gastro.2015.11.043
Rinna, A., Magdolenova, Z., Hudecova, A., Kruszewski, M., Refsnes, M., Dusinska, M. (2015). Effect of silver nanoparticles on mitogen-activated protein kinases activation: role of reactive oxygen species and implication in DNA damage. Mutagenesis 30, 1. doi: 10.1093/mutage/geu057
Rizzuto, R., Pozzan, T. (2006). Microdomains of intracellular Ca2+: molecular determinants and functional consequences. Physiol. Rev. 86, 1. doi: 10.1152/physrev.00004.2005
Rodriguez-Gutierrez, R., Salcido-Montenegro, A., Gonzalez-Gonzalez, J. G. (2018). Early clinical expressions of insulin resistance: the real enemy to look for. Diabetes Ther. 9, 1. doi: 10.1007/s13300-017-0348-2
Scherz-Shouval, R., Elazar, Z. (2007). ROS, mitochondria and the regulation of autophagy. Trends Cell Biol. 17, 9. doi: 10.1016/j.tcb.2007.07.009
Souli, A., Sebai, H., Chehimi, L., Rtibi, K., Tounsi, H., Boubaker, S., et al. (2015). Hepatoprotective effect of carob against acute ethanol-induced oxidative stress in rat. Toxicol. Ind. Health 31, 9. doi: 10.1177/0748233713475506
Sun, S., Yang, S., Dai, M., Jia, X., Wang, Q., Zhang, Z., et al. (2017). The effect of Astragalus polysaccharides on attenuation of diabetic cardiomyopathy through inhibiting the extrinsic and intrinsic apoptotic pathways in high glucose-stimulated H9C2 cells. BMC Complement Altern. Med. 17, 1. doi: 10.1186/s12906-017-1828-7
Tikhanovich, I., Cox, J., Weinman, S. A. (2013). Forkhead box class O transcription factors in liver function and disease. J. Gastroenterol. Hepatol. 28, 125–131. doi: 10.1111/jgh.12021
Tsujimoto, T., Yamamoto-Honda, R., Kajio, H., Kishimoto, M., Noto, H., et al. (2015). Prediction of 90-day mortality in patients without diabetes by severe hypoglycemia: blood glucose level as a novel marker of severity of underlying disease. Acta Diabetol. 52, 2. doi: 10.1007/s00592-014-0640-9
von Horn, C., Minor, T. (2019). Modern concepts for the dynamic preservation of the liver and kidneys in the context of transplantation. Pathologe 40 (3), 292–298. doi: 10.1007/s00292-019-0595-2
Wahren, J., Ekberg, K. (2007). Splanchnic regulation of glucose production. Annu. Rev. Nutr. 27, 329–345. doi: 10.1146/annurev.nutr.27.061406.093806
Williams, J. F., Olefsky, J. M. (1990). Defective insulin receptor function in down-regulated HepG2 cells. Endocrinology 127, 1706–1717. doi: 10.1210/endo-127-4-1706
Xie, X., Chang, X., Chen, L., Huang, K., Huang, J., Wang, S., et al. (2013). Berberine ameliorates experimental diabetes-induced renal inflammation and fibronectin by inhibiting the activation of RhoA/ROCK signaling. Mol. Cell. Endocrinol. 381, 1–2. doi: 10.1016/j.mce.2013.07.019
Yu, S. J., Bae, S., Kang, J. S., Yoon, J. H., Cho, E. J., Lee, J. H., et al. (2015). Hepatoprotective effect of vitamin C on lithocholic acid-induced cholestatic liver injury in Gulo(-/-) mice. Eur. J. Pharmacol. 762, 247–255. doi: 10.1016/j.ejphar.2015.06.008
Zhang, K., Pugliese, M., Pugliese, A., Passantino, A. (2015). Biological active ingredients of traditional Chinese herb Astragalus membranaceus on treatment of diabetes: a systematic review. Mini Rev. Med. Chem. 15, 315–329. doi: 10.2174/1389557515666150227113431
Zhou, Y., Lee, J., Reno, C. M., Sun, C., Park, S. W., Chung, J., et al. (2011). Regulation of glucose homeostasis through a XBP-1-FoxO1 interaction. Nat. Med. 17, 356–365. doi: 10.1038/nm.2293
Zhu-jun, M., Dan, S., Ke-fu, C. (2017). Effects of astragalus polysaccharides and berberine on attenuating insulin resistance in IR-INS-1 cells by down-regulating miR-126-3p 32, 7, 42–45.
Zhu-jun, M., Si-ying, W., Lei, W., Ke-fu, C. (2018). Mechanism study of berberine in regulating miR-29a-3p and insulin receptor pathway in IR-HepG2 cells. CJTCMP 33 (12), 5559–5562.
Keywords: astragalus polysaccharide, berberine, insulin resistance, IR-HepG2 cell model, intracellular calcium flux, gluconeogenesis signaling pathway
Citation: Mao Z-J, Lin M, Zhang X and Qin L-P (2019) Combined Use of Astragalus Polysaccharide and Berberine Attenuates Insulin Resistance in IR-HepG2 Cells via Regulation of the Gluconeogenesis Signaling Pathway. Front. Pharmacol. 10:1508. doi: 10.3389/fphar.2019.01508
Received: 18 April 2019; Accepted: 20 November 2019;
Published: 23 December 2019.
Edited by:
Adolfo Andrade-Cetto, National Autonomous University of Mexico, MexicoReviewed by:
Marcia Hiriart, National Autonomous University of Mexico, MexicoSubbaRao V. Madhunapantula, JSS Academy of Higher Education and Research, India
Copyright © 2019 Mao, Lin, Zhang and Qin. 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) and the copyright owner(s) 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: Zhu-Jun Mao, maozhujun0107@163.com; Lu-Ping Qin, qinsmmu@126.com