Abstract

Physical activity has emerged as a powerful intervention that promotes healthy aging by maintaining the functional capacity of critical organ systems. Here, by combining functional and proteomics analyses, we examined how hepatic phenotypes might respond to exercise treatment in aged rats. 16 male aged (20 months old) SD rats were divided into exercise and parallel control groups at random; the exercise group had 8 weeks of treadmill training with medium intensity. Whole protein samples of the liver were extracted from both groups and separated by two-dimensional gel electrophoresis. Alternatively objective protein spots with >2-fold difference in expression were selected for enzymological extraction and MS/MS identification. Results show increased activity of the manganese superoxide dismutase and elevated glutathione levels in the livers of exercise-treated animals, but malondialdehyde contents obviously decreased in the liver of the exercise group. Proteomics-based identification of differentially expressed proteins provided an integrated view of the metabolic adaptations occurring in the liver proteome during exercise, which significantly altered the expression of several proteins involved in key liver metabolic pathways including mitochondrial sulfur, glycolysis, methionine, and protein metabolism. These findings indicate that exercise may be beneficial to aged rats through modulation of hepatic protein expression profiles.

1. Introduction

Sedentary habits are associated with low levels of physical fitness and a high risk of mortality [1], whereas physical activity has emerged as a powerful intervention that promotes successful aging by maintaining general physical fitness [2]. However, the physiological mechanisms underlying physical activity-induced physical fitness are only partly understood.

Exercise-induced oxygen consumption increases associated with energy metabolism adaptation might be involved in physical fitness [3]. The first direct evidence of this was provided by Davies et al. [4] through the demonstration that high-intensity exercise enhances reactive oxygen species (ROS) production. Although high concentrations of ROS can damage proteins, lipids, and nucleic acids, further evidence suggests that exercise promotes the formation of ROS, which participate in redox regulation. In particular, moderate concentrations of ROS function as regulatory mediators to reestablish “redox homeostasis” in signaling processes [5, 6]. In this context, exercise-induced ROS may be involved in the establishment of physical fitness [7].

Although the potential effects of muscular ROS have been explored in depth, the role and related mechanisms of ROS in the hepatic adaptation induced by exercise have not been extensively studied [8]. Several recent findings have resulted in physical activity being prescribed as part of the overall strategy (or therapy) to reduce the risk of senescence-associated liver metabolic syndrome. The potential mechanisms of this protective role might involve the enhancement of hepatic insulin sensitivity and the activities of antioxidant enzymes, the attenuation of lipid accumulation in the liver, and subsequent prevention of hyperglycemia, hypercholesterolemia, and liver diseases (such as nonalcoholic steatohepatitis or cirrhosis) [9]. In addition, ROS are direct triggers of liver insulin resistance [10] and inflammation [11]. In contrast, however, evidence from rats indicates that regular exercise can attenuate hepatic oxidative stress [12]. Nonetheless, further research is needed to evaluate the detailed mechanism of ROS in the hepatic adaptation induced by exercise.

Given the multitude of changes that occur in each organ during exercise, genetic and epigenetic studies might be able to elucidate the biology underlying the individual response to exercise [13]. 2D gel tandem mass spectrometry-based proteomic methods are valuable tools for quantifying differences in protein abundance between different physical conditions or treatments [14]. Current proteomic methods allow the identification of several thousand proteins from microgram quantities of proteins in skeletal muscle, cardiomyocytes [15], adipose tissue [16], and liver [7, 17, 18]. In addition, proteomics analysis has previously been performed by several groups to determine how liver proteome profiles differ between exercise-trained and untrained individuals [18], which might provide information that could be of benefit in combating the effects of postmenopause or of nonalcoholic fatty liver disease-related loss of liver function including metabolic syndrome by promoting redox homeostasis. However, insight into the adaptive response of liver tissue proteome signatures during aging under a variety of regular exercise challenges is still required [19]. In the current study, we focused our attention specifically on the putative phenotypic/metabolic adaptations in the liver of aged rats that were induced by eight weeks of medium-intensity exercise.

2. Material and Methods

2.1. Subjects

In total, 16 male adult Sprague-Dawley (SD) rats (Rattus norvegicus, Guangdong Medical Laboratory Animal Center (GDMLAC), Guangdong, China) with similar body weights (range of 340–260 g), aged 18 months at the beginning of the training period, were randomly distributed into two groups: treadmill exercise () (EXE) and sedentary control () (CON). The animals were fed standard laboratory chow and water ad libitum while housed (3–5 per cage) in a temperature-controlled room (22°C) with 12 h dark/light cycles. All animal protocols were approved by the local Experimental Animal Use Committee.

2.2. Endurance Training

The animal exercise protocol was as follows: the training intensity was progressively increased. For the first 4 weeks, the speed of the treadmill and duration of the training sessions were gradually increased from a speed of 10 m/min for 10 min to a running speed of 15 m/min for 60 min by the end of week 4. For the next 4 weeks, a 5 min warm-up session at a speed of 10 m/min was followed by the 60 min training session at a speed of 15 m/min [20]. This training protocol was developed based on previous experience training in rodents [21], but minor modifications were applied according to another training protocol [22]. The intensity and duration were chosen based on previous reports from the groups of Song and Bejma in which a similar speed of 15 m/min and a similar parameter of 60 min/bout were applied, which resulted in a series of positive changes such as increased activity of antioxidant enzymes in the aging liver [21, 23]. The 8-week exercise program of treadmill running at 0° incline was chosen based on previous reports from the Siu group [22]. Training animals and control animals were sacrificed 6 h after the last training session. The liver tissues were quickly removed and frozen immediately in liquid nitrogen and stored at −80°C until further analysis.

2.3. Assay of Malondialdehyde (MDA) and Glutathione (GSH) Contents as well as Superoxide Dismutase (SOD) Activity in the Livers of Old Rats

For biochemical analyses of in vivo MDA, GSH, and SOD activity, 20 mg of liver tissue was homogenized on ice in phosphate-buffered saline and centrifuged for 10 min at 12,000 ×g at 4°C. The supernatant was collected and analyzed. MDA level, GSH content, and SOD activity were determined using commercial assay kits produced by Nanjing Jiancheng Bioengineering Institute, China. Protein concentrations were determined using a BCA protein assay kit (Nanjing Jiancheng Biochemistry Co.) as described by the manufacturer.

2.4. Proteomics Analysis
2.4.1. Sample Preparation for Two-Dimensional Gel Electrophoresis (2DGE)

We selected 8 “EXE” and 8 “CON” rats for proteomics analysis. From these animals, 50 mg of liver was homogenized in 500 μL of sample buffer [(7 M urea, 2 M thiourea, 4% (w/v) dimethyl[3 propyl] azaniumyl propane-1-sulfonate (CHAPS), 1% (v/v) immobilized pH gradient (IPG) buffer pH 3–10 NL) (Amersham Biosciences, Piscataway, NJ, USA), 40 mM dithiothreitol (DTT), and Protease Inhibitor Cocktail Complete (Roche, Mannheim, Germany)] and allowed to stand at room temperature for 1 h, after which the homogenate was then centrifuged at 40,000 ×g for 60 min at 4°C. Subsequently, the supernatant was removed and the protein concentration was measured using the Bradford protein assay. Homogenates were stored at −70°C for later use. 2-DE was performed using a 24 cm (nonlinear, pH 3.0–10.0) IPG gel strip, according to Sun et al. [24]. Spot detection, gel matching, and interclass analysis were performed using PDQuest 8.0 software (Amersham Biosciences, USA) [24].

2.4.2. In-Gel Tryptic Digestion and Matrix-Assisted Laser Desorption/Ionization-Time of Flight Mass Spectrometry (MALDI-TOF MS/MS)

In-gel digestion and protein identification were performed as described by Perkins et al. [25]. Candidate spots were excised from the stained gel, destained with 0.1 M ammonium bicarbonate in 50% acetonitrile (Sigma-Aldrich, St. Louis, MO, USA), and dried using a SpeedVac SC110 (Savant Instruments, Holbrook, NY, USA). The gels were then rehydrated and incubated in trypsin solution (Promega, Madison, WI, USA) overnight at 37°C, and the peptide mixtures were then analyzed using MALDI-TOF MS/MS (Applied Biosystems, Foster City, CA, USA) [25].

2.4.3. Database Searching

Proteins were identified by peptide mass fingerprinting with the search engine programs MASCOT [25]. The criteria for positive identification of proteins were according to Perkins et al. [25] and were set as follows: (i) at least four matching peptide masses, (ii) 50 ppm or better mass accuracy, and (iii) matching of the molecular weights (Mr) and isoelectric points (pI) of the identified proteins to the estimated values obtained from image analysis.

2.5. Statistical Analysis

Statistical analysis was performed using predictive analytics software statistics 16.0 (SPSS Inc., Chicago, IL, USA). Comparisons across the experimental groups were performed using Student’s -test. For the analyses, the significance level was set at 5% (). Data are presented in the figures as the mean ± standard error of the mean (SEM).

3. Results

3.1. Liver Weights and Liver Weight/Body Weight Ratios between Sedentary and Exercised Aged Rats

The ratio of the weight of the liver to the whole animal body weight was determined. As shown in Figure 1, no difference in body weight was observed between exercised (EXE) and sedentary control (CON) groups. However, liver weight and liver weight/body weight ratio were higher for the livers of EXE animals than for those of the CON group ().

3.2. Effects of Treadmill Training on the Oxidative Status of the Livers of Aged Rats

Changes in the levels of hepatic MDA and GSH contents as well as in MnSOD activity in response to 8-week moderate-intensity treadmill training were analyzed. As shown in Figure 2, after 8 weeks of treadmill training, MnSOD activity, hepatic GSH level, and GSH/oxidized glutathione (GSSH) ratio in the EXE group were significantly higher in comparison to the CON group (, Figure 2), whereas the hepatic MDA and GSSH levels were decreased (, Figure 2).

3.3. Proteomics Comparison of Hepatic Expression Profiles between Sedentary and Exercised Aged Rats

The hepatic protein profiles are shown in Figure 3. Following gel image analysis, differentially expressed spots were selected and identified using MALDI-TOF MS when the normalized spot intensities showed differences of at least 2-fold between the groups. The 12 proteins identified as differentially expressed are listed in Table 1. The experimental ratios of the Mr and pIs matched those of the theoretical data, suggesting that the identification of proteins by our proteomic method was reliable.

4. Discussion

The liver is the primary organ required for processing of nutrients, hormones, and drugs; therefore, to evaluate exercise-dependent antiaging effects and adaptive alterations in the proteome of the livers of aged rats after 8 weeks of moderate-intensity training exercise, total liver tissue extracts from sedentary aged rats (22 months, control) and aged rats that were provided exercise training were resolved by 2DGE. 2DGE with medium-range IPGs (pH 3–10) revealed 3839 protein spots in the liver that were matched among all rats. Among these protein spots, 35 were found to be differentially expressed in the liver between sedentary and exercised rats, of which 12 were identified by MALDI-TOF MS/MS (Figure 3). From these proteins, we elucidated that age-related hepatic proteome changes were effectively reversed through exercise, which might be beneficial to combat health conditions related to age-related abnormal liver morphology through regulation of mitochondrial energy metabolism, proteolytic systems, detoxification, and quality control of cellular and cytoskeletal proteins.

Among the identified proteins, FAS activity is believed to be a determinant of the maximal capacity of the liver to synthesize fatty acids by de novo lipogenesis. Advanced age has been associated with upregulation of FAS expression and enzymatic activity, which in turn has been found to be an underlying contributor to many age-related conditions including hepatic steatosis [26] and necroinflammation [27]. Aging is also associated with variable degrees of fibrosis in the liver. Our results showed that moderate-intensity training exercise decreased the abundance of FAS protein, a marker of hepatic fibrosis.

SAFB1 has been shown to be important in numerous cellular processes including cell growth, stress response, and apoptosis [28]. Studies performed by the Gelman group have shown that SAFB1 can repress the transcriptional activity of various other nuclear receptors including estrogen receptor [29], heat shock factor 1 [30], and p53 [31]. Additional studies identified that SAFB1 regulates the expression of xanthine oxidoreductase (XOR) [32], peroxisome proliferator-activator receptor (PPAR) gamma [33], and sterol regulatory element-binding protein- (SREBP-) 1c [33]. SREBP-1c and PPAR gamma are transcription factors that control the expression of hepatic genes involved in fatty acid and triglyceride synthesis during overnutrition to facilitate the conversion of glucose to fatty acids and triglycerides for the storage of the excess energy [34]. Additionally, in this context, it should be noted that SAFB1 was found to be downregulated by exercise, which is probably related to a decrease in lipogenesis as suggested by decreased abundance of the FAS protein. Furthermore, Lee et al. [35] showed that SAFB1 is also involved in the later stages of apoptosis, in which it is possibly associated with endonuclease–mediated DNA cleavage. These data clearly suggest that downregulation of SAFB1 in the liver might play an important role in age-related fatty acid metabolism, apoptosis, and oxidative stress, which are the primary causes of fatty liver disease.

In a previous study, perturbations in the energy metabolism of aged livers were indicated by the downregulation of several glycolytic enzymes [36]. Moreover, Rothwell and Sayce found that enhancement of glycolysis via PGM can protect cells from oxidative stress, thus preventing senescence [37]. The classical function of PGM is to link glycolysis and glycogen metabolism by interconverting glucose-1-phosphate and glucose-6-phosphate. PGM has been shown to be expressed ubiquitously; however, high levels of expression are observed in the heart, skeletal muscle, kidney, liver, and lungs [37]. The activity of PGM is under hormonal control and may play a physiological role in glucose metabolism and physical activity [38]. Notably, proteomics analysis has revealed decreased expression of PGM in elderly women [39]. Recent evidence indicates that downregulation of PGM is associated with hepatocyte failure in cirrhotic hepatocytes, especially at the late stages of cirrhosis [40]. In addition, Preisler et al. [41] implied that the liver might be affected by PGM deficiency and that reduced liver glucose levels might contribute to the low blood glucose levels observed during exercise in this condition. In this regard, polymorphism in PGM has been reported to be useful as a genetic biomarker for alcoholic cirrhosis [42], ischemic heart disease [43], and obesity disorders [44, 45]. The results of the current study show that training exercise increased the protein content of PGM in the livers of old rats, which likely contributed to the improved glycolytic metabolism of these animals, despite their advanced cirrhosis.

Notably, α-enolase1 (ENO1), a glycolytic enzyme that catalyzes the conversion of 2-phosphoglycerate to phosphoenolpyruvate, was found to be downregulated in the liver of old rats after moderate-intensity training exercise, as shown in Figure 3. ENO1 is highly expressed in liver tissue and in the cytoplasm of hyperplastic bile ducts [46]. In addition to its key contribution in the regulation of glycolysis and the generation of ATP, ENO1 also plays a crucial role in growth control, hypoxia tolerance, and allergic responses [47]. Notably, a previous study had used an immunoproteomics approach to identify and characterize ENO1 as an autoantigen in liver fibrosis [48]. Although it is unclear whether ENO1 plays an active role during the formation of liver fibrosis, it has also been reported that ENO1 is transcriptionally regulated by the c-Jun NH2-terminal kinase (JNK) signaling pathway, which is frequently activated in liver fibrosis [49]. Furthermore, ENO1 might act as an oxidative stress-related mediator to modulate the tissue and cellular events responsible for the progression of liver fibrosis [50]. Our results show a decreased protein content of ENO1 at 8 weeks after moderate-intensity training exercise in aged rats, in parallel with an increase in MnSOD activity and in the hepatic GSH level and GSH/GSSH ratio (Figure 2), which likely contributes to improved liver fibrosis and liver function. However, whether or not the improvement of liver function mediated by training exercise is associated with modulation of ENO1 in the livers of aged rats requires further study.

A previous study demonstrated that 11β-HSD1 catalyzed the intracellular regeneration of active glucocorticoids, promoting insulin resistance in the liver [51], and elevated mRNA levels for 11β-HSD1 in the liver were associated with aging [52]. Fatty liver occurs in mice that overexpress 11β-HSD1 as a result of increases in gluconeogenesis, fatty acid synthesis, and plasma insulin levels [53]. Thus, the age-dependent increase in the mRNA level of 11β-HSD1 is likely to induce insulin secretion and cause hyperinsulinemia. Recent evidence indicates that exercise protects against the elevations in glucocorticoids associated with the development of insulin resistance by decreased adrenal sensitivity to the 11β-HSD1 content in the livers of diabetic fatty rats [54]. Consistent with these findings, aged rats that exhibited higher hippocampal corticosterone levels also showed increased 11β-HSD1 protein expression as well as increased levels of the enzyme that catalyzes glucocorticoid formation and greater hippocampal glucocorticoid receptor (GR) activation, whereas exercise as a therapeutic intervention significantly reduced total hippocampal GR and 11β-HSD1 expression in these animals [55]. Similarly, the results of the current study show that 11β-HSD1 protein expression significantly decreased after 8 weeks of aerobic exercise; as discussed above, this decrease in 11β-HSD1 expression is believed to play an important role in mitigating age-related fatty acid metabolism, gluconeogenesis, and hyperinsulinemia, which is a main cause of insulin resistance in the liver.

The liver is an important organ for H2S production and clearance, and hepatic H2S metabolism affects multiple metabolic processes. In contrast, malfunction of hepatic H2S metabolism might be involved in the pathogenesis of many liver diseases such as hepatic fibrosis and cirrhosis [56]. The mitochondrial dioxygenase gene ETHE1 encodes a mitochondrial enzyme involved in sulfide detoxification, which uses molecular O2 and water to oxidize the mobile persulfide from quinone oxidoreductase to form sulfite, which in turn is involved in sulfur metabolism and takes part in the pathway in which H2S is oxidized to thiosulfate in three steps [57]. Although the basal endogenous concentration of H2S has been found to have a hepatoprotective effect [58], inhibition of this pathway prevents binding of H2S to quinone oxidoreductase and enables the accumulation of H2S and thiosulfate, which induces oxidative stress and hepatotoxicity (e.g., by inhibiting the respiratory chain and decreasing cellular NAD+/NADH) in the liver during aging [59]. Thus, the decreased levels of ETHE1 in the rats fed a high-fat diet might allow the occurrence of higher concentrations of H2S, which might mitigate some of the risks of high blood pressure and oxidative stress associated with obesity and aging [60]. However, ETHE1 is also an antiapoptotic protein that increases the deacetylase activity of p53 in association with histone deacetylase 1, leading to the suppression of apoptosis [61].

S-Adenosylmethionine (SAMe) is an important methyl donor in various methyltransferase reactions that occur in many lipid metabolic pathways. SAMs catalyzes the synthesis of SAMe from methionine in the presence of ATP, which is the major pathway in mammals for the metabolism of methionine [62]. Previous studies have indicated that a marked loss of SAMs was accompanied by a reduction of the concentration of SAMe in cirrhosis [62] and that treatment with SAMe increased hepatic GSH levels and increased survival in patients with less advanced alcoholic liver cirrhosis [63]. Our data therefore support that SAMs upregulates SAMe via exercise, which in turn increases the levels of hepatic GSH (Figure 2) with a concomitant decline in FAS [64] and in K8 [65] (discussed below), especially in the liver of aged rats.

RNH1 belongs to the leucine-rich repeats (LRR) motif protein superfamily [66] and is a tight-binding cytosolic inhibitor of the RNase A family with both ribonucleolytic and angiogenic activities that are essential for cell survival [67]. In addition, RNH1 is a known regulator of vascularization and a mediator of oxidative stress. Knockdown of RNH1 has been shown to significantly decrease cell survival under stress conditions, accompanied with an increase in cell apoptosis. The reported antioxidant [68] and redox homeostatic [69] effects might also contribute to the regulatory function of RNH1 in cell survival. These results suggest that overexpression of RNH1 has antioxidant effects that protect against hepatic damage and apoptosis induced by H2O2 in mice [68]. In this context, it should be noted that RNH1 was observed to be upregulated by exercise in the current study, which is likely related to the increase in antioxidant activity and antiapoptosis effects as suggested by the increase in MnSOD activity and the levels of hepatic GSH (Figure 2).

Tight regulation of proteostasis is essential for maintaining cellular homeostasis in postmitotic tissues [70]. Advanced age is an intrinsic process that results in alterations of tissue protein homoeostasis. These degenerative changes are intimately associated with abnormal protein aggregation and impaired protein degradation pathways. Direct interventions to modulate protein quality control mechanisms might therefore improve health and potentially increase lifespan. Proteasomes can be regarded as the major proteolytic enzymes responsible for the degradation of the majority of cellular proteins [71]. PSMB6 is a member of the 20S proteasomal subunit family, which forms the proteolytic core of the 26S proteasome [71]. Beta subunits are characterized as containing the active sites of proteases, whereas the subunits serve as docking domains for the regulatory particles and form a gate that blocks unregulated access of substrates to the central pore. The age-associated decrease in the activity of the 26S proteasomes, which are responsible for the degradation of ubiquitinated proteins, is in accordance with previous findings that high-molecular-weight ubiquitin-protein conjugates accumulate in various tissues over time [70].

In this context, it is interesting to note that, in addition to the increase in the activity of antioxidant enzymes known to occur in the liver of aged rats as a result of training exercise, the upregulation of proteasomes might also play a significant role in improved metabolism by reducing the amount of damaged proteins in the liver. Radák et al. [72] previously showed that regular moderate exercise initiated at advanced age and continued for two months in BDF1 mice resulted in significant upregulation of the activity of antioxidant enzymes and 26S proteasomes in the brain. Our findings are also in accordance with previous reports that exercise training mitigated to a substantial degree of oxidative stress caused by advanced age, with a concomitant increase in proteasomal peptidase activities in the brain.

K8, a cytoprotective stress protein, belongs to the keratin family of cytoskeletal intermediate filament proteins that play a central role in protecting hepatocytes from apoptosis; accordingly, K8 mutations predispose their carriers to liver cirrhosis [73]. Emerging evidence also indicates that K8 protects cells from nonmechanical injury by cellular survival kinases and cell cycle progression [74]. Subsequent studies have suggested that keratins modulate apoptosis induced by FAS [73] and tumor necrosis factor [75]. In cirrhotic livers, K8 hyperphosphorylation was previously shown to occur preferentially in hepatic nodule cells adjacent to bridging fibrosis and was associated with increased stress kinase activation, oxidative injury, and susceptibility to undergo apoptosis [65]. Furthermore, according to Kucukoglu et al. [76], a high-fat diet triggers aggregate formation and the development of liver oxidative injury in susceptible individuals through the misfolding and cross-linking of excess K8. The possibility that the improvement in liver fibrosis and fatty liver mediated by exercise is associated with modulation of K8 in the livers of aged rats requires further study.

HP is produced as an acute-phase reactant during liver fibrosis (LF), and it possesses anti-inflammatory and antioxidative properties [77], preventing liver damage. Thus, HP might be a potential molecular target for early LF diagnostics and therapeutic applications [78]. In contrast, lack of HP attenuates the hepatomegaly/steatosis often associated with obesity. When animals were stressed by acute exercise or by a high-fat diet in a previous study, oxidative stress, muscle atrophy, and force drop were exacerbated in HP-knockout mice [79]. Conversely, recent data demonstrated that the increased HP levels in the serum in response to traumatic brain injury were derived from the liver and that IL-6 is an important mediator of this induction [80]. T-kininogen has also been demonstrated to be a cysteine proteinase inhibitor present in rat plasma [81]. In a previous report, preliminary evidence suggested that the degree of reduction in plasma T-kininogen levels was directly related to the severity of impaired liver function [81]. Recent results showed that aged rats present reduced expression of T-kininogen and HP and increased hepatic injury in response to an inflammatory insult [82]. Therefore, the upregulation of HP and T-kininogen following training exercise observed in this study might be the consequence of decreased levels of markers of inflammation and oxidative stress (Figure 2) in the livers of aged rats.

This study has several limitations. First, the changes in protein levels were not confirmed by western blot. Second, comprehensive data were not obtained regarding the phenotype of the aging liver, such as the presence of liver tissue histology, hepatocyte apoptosis, intrahepatic lipid accumulation, or cirrhosis. Third, we only measured levels of primary antioxidative enzymes (MnSOD and GSH), secondary antioxidative enzymes such as glutathione-S-transferase (GST) and glucose 6-phosphate dehydrogenase, and molecules such as glutathione.

In summary, we succeeded in objectively identifying 12 proteins that exhibited adaptive changes in their expression associated with extended exercise in the livers of aged rats. Identification and detailed analysis of the impacted proteins indicated that exercise training significantly altered the expression of several proteins involved in key liver metabolic pathways including sulfur, glycolysis, and methionine metabolism as well as general protein metabolism in conjunction with the modulation of overall hepatic protein expression profiles. These proteins are potential disease state- and intervention-specific biomarkers. Further research on these specific proteins is needed to uncover possible pathways involved in the beneficial effects of exercise on metabolic adaptations in the liver of aged rats.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The current project was supported by the National Natural Science Foundation of China (Grant 31500961) and the Guangdong Scientific Project (Grant 2014A020220015 and Grant 2015A020219015) and partially supported by the Distinguished Young Talents in Higher Education of Guangdong, China (Grant 2014KQNCX225).