|Home | About | Journals | Submit | Contact Us | Français|
β-Glucosidase is an important component of the cellulase enzyme system. It does not only participate in cellulose degradation, it also plays an important role in hydrolyzing cellulose to fermentable glucose by relieving the inhibition of exoglucanase and endoglucanase from cellobiose. Therefore, the glucose-tolerant β-glucosidase with high specific activity for cellobiose might be a potent candidate for industrial applications.
The β-glucosidase gene bgl that encodes a 443-amino-acid protein was cloned and over-expressed from Thermoanaerobacterium thermosaccharolyticum DSM 571 in Escherichia coli. The phylogenetic trees of β-glucosidases were constructed using Neighbor-Joining (NJ) and Maximum-Parsimony (MP) methods. The phylogeny and amino acid analysis indicated that the BGL was a novel β-glucosidase. By replacing the rare codons for the N-terminal amino acids of the target protein, the expression level of bgl was increased from 6.6 to 11.2 U/mg in LB medium. Recombinant BGL was purified by heat treatment followed by Ni-NTA affinity. The optimal activity was at pH 6.4 and 70°C. The purified enzyme was stable over pH range of 5.2–7.6 and had a 1 h half life at 68°C. The activity of BGL was significantly enhanced by Fe2+ and Mn2+. The Vmax of 64 U/mg and 120 U/mg were found for p-nitrophenyl-β-D-glucopyranoside (Km value of 0.62mM) and cellobiose (Km value of 7.9mM), respectively. It displayed high tolerance to glucose and cellobiose. The Kcat for cellobiose was 67.7s-1 at 60°C and pH 6.4, when the concentration of cellobiose was 290mM. It was activated by glucose at concentrations lower that 200mM. With glucose further increasing, the enzyme activity of BGL was gradually inhibited, but remained 50% of the original value in even as high as 600mM glucose.
The article provides a useful novel β-glucosidase which displayed favorable properties: high glucose and cellobiose tolerance, independence of metal ions, and high hydrolysis activity on cellobiose.
Cellulosic biomass is the most abundant renewable resource on earth, whose natural degradation represents an important part of the carbon cycle within the biosphere . β-Glucosidase (EC 184.108.40.206) is a glucosidase enzyme that acts upon β 1–4 bonds linking two glucose or glucose-substituted molecules. It is an important component of the cellulase enzyme system. The limiting step in the enzymatic saccharification of cellulosic material is the conversion of short-chain oligosaccharides and cellobiose, which was resulted from the synergistic action of endogucanases (EC 220.127.116.11) and cellobiohydrolases (EC 18.104.22.168), to glucose, a reaction catalyzed by β-glucosidases . It is well established that cellobiose inhibits the activities of most cellobiohydrolases and endoglucanses . β-glucosidases reduce cellobiose inhibition by hydrolyzing this disaccharide to glucose, thus allowing the cellulolytic enzymes to function more efficiently [4,5]. Furthermore, β-glucosidase is used as a flavor enzyme to enhance the flavor of wine, tea and fruit juice [6,7]. In fruits and other plant tissues many secondary metabolites, including flavor compounds, are accumulated in their glucosylated form [8,9]. Because β-glucosides constitute the majority of the known glycoconjugated flavor compounds, β-glucosidases play an important role in flavor liberation from these precursors. Therefore, producing high-activity and glucose-tolerant β-glucosidase has become important.
Recently, the search for β-glucosidases insensitive to glucose has increased significantly, for these enzymes would improve the process of saccharification of lignocellulosic materials. A few microbial β-glucosidases have been reported to tolerate glucose [10-14]. For example, β-glucosidases from Aspergillus tubingensis CBS 643.92, A. oryzaeA. niger CCRC 31494, A. foetidus, and marine microbial metagenome displayed high inhibition constant by glucose (Ki) of 600mM, 1390mM, 543mM, 520mM, and 1000mM, respectively. But these β-glucosidases have considerably lower specific activity for cellobiose than for p-nitrophenyl-β-D-glucopyranoside. Therefore, over-expression of thermostable β-glucosidase with high glucose tolerance and specific activity for cellobiose abilities will help shed light on degradation of cellulosic biomass.
Thermostable enzymes have several generic advantages, allowing a decreased amount of enzyme needed because of higher specific activity and elongated hydrolysis time due to higher stability. In addition, thermostable enzymes are generally more tolerant and allow more flexibility in process configurations [15,16]. Although some glucose-tolerant β-glucosidases from fungi and bacteria have been reported [10-14], the glucose-tolerant β-glucosidases genes have not been expressed and characterized from thermophilic bacteria. Bacterium Thermoanaerobacterium thermosaccharolyticum is a strict anaerobe that grows on wide range of hexose and pentose at temperature from 37°C to 75°C, which have attracted considerable interests to hydrogen production and thermostable enzyme production . T. thermosaccharolyticum DSM 571 could utilize cellobiose, but the gene for β-glucosidase, the key enzyme in degradation cellobiose, was not reported in the Genbank (NC_014410.1). Because the optimal growth temperature for T. thermosaccharolyticum DSM 571 was at 60°C, the thermostable β-glucosidase could have a considerable potential for industrial applications. Owing to the inherent difficulty of cultivation of T. thermosaccharolyticum DSM 571, it is difficult to obtain a sufficient amount of cells for large-scale enzyme production. For the production of the recombinant protein, genetic engineering is the first choice because it is easy, fast, and cheap.
In this paper, we report the phylogenesis analysis, cloning, over-expression, and detailed biochemical characterization of the β-glucosidase from T. thermosaccharolyticum DSM 571. The favorable properties make the β-glucosidase a good candidate for utilization in biotechnological applications.
By analysis of the genome sequence of T. thermosaccharolyticum DSM 571, a protein (Tthe_1813), defined as β-galactosidase in Genbank, consists of a 1,329-bp fragment encoding 443 amino acids, which belonged to family 1 of the glycoside hydrolases. It shares the highest sequence similarity of 66% with the β-glucosidses from Thermoanaerobacter mathranii (Genbank No. YP_003676178.1) and Thermoanaerobacter pseudethanolicus ATCC 33223 (Genbank No. YP_001665894.1), which were revealed by whole-genome sequencing but has not been biochemically characterized. Alignment of the BGL cluster with several representative members of GH1 indicated that they share similar blocks. The catalytic proton donor, Glu135 and Glu351 in BGL are well conserved among all GH1 proteins (Figure (Figure1).1). The sequence around Glu351 in BGL is [LYT-NGAA], which is consistent with the consensus pattern of PS00572. The results indicated that the protein (Tthe_1813) could be a novel β-glucoside. Then the DNA fragment of a protein (Tthe_1813) gene was amplified from genomic DNA of T. thermosaccharolyticum DSM 571, and ligated to pET-20b at Nde I and Xho I sites to generate plasmid pET-20-BGL.
In order to increase the expression level of BGL in E. coli, site-directed mutagenesis were designed and performed to optimize condons of BGL for E. coli expression system. pET-20-BGLII was obtained from pET-20-BGL in which the rare condons for the N-terminal amino acid residues were replaced by optimal codons in E. coli without and change of amino acid sequence (Figure (Figure2),2), so pET-20-BGLII encodes the same β-glucosidase as that encoded by the wild-type gene. The β-glucosidase activity expression from pET-20-BGLII was 7.5 U/mL (11.2 U/mg total of cell protein) and was estimated to be about 30% of the total protein, which was about 1.7 times higher than the expressed from pET-20-BGL (Figure (Figure3,3, lane 2 and 3).
The protein in the cell-free extract was purified to gel electro homogeneity after a heat treatment and a Ni-NTA affinity. The final preparation gave a single band on SDS-PAGE gel and the molecular mass of the enzyme was estimated to be 52kDa (Figure (Figure3,3, lane 4).
The biochemical properties of BGL were investigated by using the purified recombinant BGL. The optimal pH of the BGL was determined to be 6.4 (Figure (Figure4a),4a), while the β-glucosidase activity was higher than 50% of the maximum activity at the pH range from 5.6 to 7.2. The enzyme was stable for about 1h at pH 5.6 to 8.0 at 60°C in the absence of the substrate (Figure (Figure4c).4c). The optimal temperature of the BGL was 70°C, which the β-glucosidase activity was higher than 40% of the maximum activity at the temperature range from 45 to 75°C (Figure (Figure4b).4b). Thermostability assays indicated that its residual activity was more than 80% after being incubated at 60°C for 2h (pH 6.4, Figure Figure4d4d).
The effects of metal ions and some chemicals on the enzyme activity were shown in (Table (Table1).1). In various assays, the enzyme activity was significantly enhanced by Fe2+, or Mn2+, and completely inactivated by Zn2+, Cu2+, Ag2+, or Hg2+. The effects of Mg2+, Ca2+, K+, Li2, or EDTA (10mM) on the enzyme activity were not so significant.
The enzyme was able to hydrolyze p-nitrophenyl-β-D-glucopyranoside, cellobiose, and p-nitrophenyl-β-D-galactopyranoside, while no activity was detected upon p-nitrophenyl-α-L-arabinofuranoside, p-nitrophenyl-β-D-xylopyranoside, maltose, CMC, and sucrose. p-nitrophenyl-β-D-Galactopyranoside was hydrolyzed at 40% of that of p-nitrophenyl-β-D-glucopyranoside. The dependence of the rate of the enzymatic reaction on the substrates concentration followed Michaelis-Menten kinetics, with Km and Vmax values of 0.62mM and 64 U/mg for p-nitrophenyl-β-D-glucopyranoside, and for cellobiose 7.9mM and 120 U/mg under optimal conditions. The effects of the substrate, cellobiose (290mM), on the enzyme activity were not significant. The Kcat/Km value for cellobiose 13.3mM-1s-1 was less than the β-glucosidase from A. oryzae, but the activity of β-glucosidase from A. oryzae was inhibited by cellobiose, and rapidly decreased above 50°C (Table (Table2).2). Furthermore, the enzyme activity was enhanced by the concentrations of glucose below 200mM, and the enzyme activity was increased 110% when adding 100mM glucose into reaction mixtures (Figure (Figure5).5). When glucose was increased, the enzyme activity of BGL was gradually inhibited, with a Ki of 600mM glucose (Figure (Figure5).5). The properties of the glucose-tolerant β-glucosidase from other microorganisms are summarized in Table Table2.2. As Table Table22 shows, these enzymes have many distinct features, especially in their catalytic properties [12,13,18-21].
Production of glucose from 290mM cellobiose (10%) by the purified BGL was examined. Even if the final concentration of glucose in reaction reached about 580mM, cellobiose (290mM) was found to be degraded completely (Figure (Figure6a,6a, b). At the beginning of the reaction, the Kcat was 67.7s-1 within one hour at 60°C which was identical to the theoretical value. During the whole degradation process, the Kcat was 28.2s-1.
To gain insights into the evolutionary relationship among β-glucosidases, we constructed the phylogenetic trees of 40 candidate sequences using he NJ method and the MP method respectively, both supporting almost the same topology. The results revealed the presences of five well-supported clades: Clade II was GH1 β-glucosidases from fungi, and Clade III was the GH3 β-glucosidases from bacteria, and Clade IV was the GH3 β-glucosidases from fungi. The GH1 β-glucosidases from bacteria was divided into two clades: Clade I mainly contained mesophilic bacteria; Clade V mainly contained thermophile, which is formed by further divided into two subclades, of which one contains all thermophile, and the other Bacillus GH1 β-glucosidases. Clade II and clade III had a relatively close relationship, and the GH1 β-glucosidases from thermophile were distant from the other clades (Figure (Figure77).
A classification of glycoside hydrolases based on amino acid sequence similarities was proposed a few years ago, wherein β-glucosidases were mainly grouped into two superfamilies of glycoside hydrolases I (GH1), and GH3 . Although, the amino acid sequence analysis indicated that BGL belongs to GH1, it shared the highest sequence similarity of 66% with the β-glucosidses from Thermoanaerobacter mathranii (YP_003676178.1). Moreover, it shared only the 63% with the putative β-glucosidase (YP_004471891.1) the Thermoanaerobacterium xylanolyticum LX-11, both belonging to the genus Thermoanaerobacterium. The Phylogenies analysis showed that the BGL was distant with the glucose-tolerant β-Glucosidases from fungi and ADD96762.1 (Figure (Figure7).7). The results indicated that the BGL could be a novel β-glucoside with some different properties. On the other hand, β-Glucosidases may be divided into three groups on the basis of their substrate specificity. The first group is known as aryl-β-glucosidases owing to strong affinity to aryl-β-glucose. The second group consists of cellobiases that hydrolyze oligosaccharides only. The third group is broad specific β-glucosidases that exhibit activity on a wide range of substrates, and are the most commonly observed form of β-glucosidases . The BGL, which was high affinity to p-nitrophenyl-β-D-glucopyranoside, hydrolyzed cellobiose, p-nitrophenyl-β-D-glucopyranoside, and p-nitrophenyl-β-D-galactopyranoside, but not p-nitrophenyl-α-L-arabinofuranoside, p-nitrophenyl-β-D-xylopyranoside, maltose, sucrose, and CMC. These results indicated that BGL belonged to the first group.
Enzymatic hydrolysis of cellulose is a complex process, the last step being a homogenous catalysis reaction involving the action of β-glucosidase on cellobiose. Cellobiose is a strong inhibitor of both cellobiohydrolases and endocellulases. Therefore, β-glucosidase with high tolerance for glucose has become heated in these fields. Fungi, especially Aspergillus species, are generally considered to be a good producer with high yield of β-glucosidases . But the major β-glucosidases belonging to family 3 of the glycoside hydrolases (GH3) from Aspergillus species were subject to competitive inhibition of glucose to produce glucose, the Ki is generally 1–20mM [10,14]. The minor β-glucosidases, which molecular weights are 40–50kDa, exhibited a tolerance to glucose (Table (Table2).2). The effect of glucose on the BGL activity revealed that the enzyme is not only resistant to end-product inhibition, but is activated by glucose at concentrations from 0 to 0.2M. Only two β-glucosidases, activated by glucose, have been reported from Scytalidium thermophilum and marine microbial (Table (Table2)2) [13,20].
Moreover, high specific activity for cellobiose and tolerance to substrate inhibition are other advantages for β-glucosidase in enzymatic hydrolysis of cellulose. Although, several β-glucosidases from a few fungi and bacteria show high glucose tolerant with Ki values of more than 200mM, the Vmax values of these enzymes for cellobiose were much lower than for p-nitrophenyl-β-D-glucopyranoside. The Vmax value of BGL for cellobiose was 120 U/mg, which was about 2 times higher than the Vmax value of BGL for p-nitrophenyl-β-D-glucopyranoside. To our knowledge, in only one other study have workers described the purification and characterization (from A. oryzae) of a β-glucosidase having such a high tolerance to glucose and high specific activity for cellobiose . But the specific activity of β-glucosidase from A. oryzae for cellobiose was much lower than for p-nitrophenyl-β-D-glucopyranoside (Table (Table2).2). The BGL was only the β-glucosidase been reported that it is not only resistant to glucose, but had higher specific activity for cellobiose than for p-nitrophenyl-β-D-glucopyranoside. In addition, the BGL had high tolerance to substrate inhibition, cellobiose. The Kcat of BGL was 67.7s-1 at 60°C and pH 6.4, when the concentration of cellobiose was 10% (Table (Table22).
The chemical agents had various effects on the activity of BGL. The chelating agent EDTA displayed no influence on the β-glucosidase activity, indicating that the β-glucosidase is not a metalloprotein. However, the β-glucosidase activity was greatly stimulated by Fe2+ or Mn2+, which implied that Fe2+ or Mn2+ is required for the maximal activity of BGL. These results distinguish BGL from the other bacteria β-glucosidases, on which Ca2+ show positive effects . In practical applications, the high thermostability of the enzyme is desired because the longer active life means the less consumption of the enzyme. The BGL residual activity was more than 80% after being incubated at 60°C for 2h, and it in enzymatic hydrolysis of cellulose exhibited high activity in broad temperature, which could keep at high levels at temperatures from 45 to 70°C.
The properties of the BGL demonstrated a great potential of the gene in the genetic modification of strains for biomass degradation. Differences in codon usage preference among organisms lead to a variety of problems concerning heterologous gene expression, which can be overcome by rational gene design and gene synthesis. Protein with multiple repetitive rare codons especially within the first 20 amino acids of the amino terminus of the protein may significantly reduce the protein expression. Sometimes, it shuts down the expression completely. Since the rare codons of bgl from 1–20 amino acids were all changed into optimized codons, the activity of BGL was increased by about 70% (Figure (Figure3).3). More optimization of codons for the other amino acid residues in the ORF of bgl may give further improvement in the gene expression levels.
With this study, we successfully over-expressed the novel β-glucosidase (BGL) gene bgl from T. thermosaccharolyticum DSM 571 by replacing the rare codons with the optimal codons in E. coli. The Phylogenies analysis showed that the BGL had close relationship with the β-Glucosidases from thermophile, and was distant from the other glucose-tolerant β-Glucosidases. As compared on the enzyme properties, the BGL was higher tolerant to glucose and cellobiose, more efficient in hydrolysis of cellibiose, more thermal stability than β-glucosidases from other microorganisms. Thus, this study provides a useful novel β-glucosidase, which may be used to improve the enzymatic conversion of cellulosic to glucose through synergetic action.
Thermoanaerobacterium thermosaccharolyticum DSM 571 was purchased from DSMZ (http://www.dsmz.de). It was grown anaerobically at 60°C as described previously . Escherichia coli JM109 and JM109(DE3) was grown at 37°C in Luria-Bertani medium (LB) and supplemented with ampicillin when required. The expression vectors pET-20b (Novagen) were employed as cloning vector and expression vector.
DNA was manipulated by standard procedures . QIAGEN Plasmid Kit and QIAGEN MinElute Gel Extraction Kit (Qiagen, USA) were employed for the purification of plasmids and PCR products. DNA restriction and modification enzymes were purchased form TaKaRa (Dalian, China). DNA transformation was performed by electroporation using GenePulser (Bio-Rad, USA). Site-directed mutagenesis of genes and the modification of the plasmids were performed by inverse-PCR followed by phosporylation and self-ligation using T4 polynucleotide kinase and T4 DNA ligase.
The β-glucosidase gene bgl was amplified from T. thermosaccharolyticum DSM 571 genomic DNA by PCR using primers bgl-1 and bgl-2 (Table (Table3),3), the PCR products were digested with Nde I and Xho I and inserted into pET-20b at Nde I and Xho I sites, yielding the plasmid pET-20-BGL.
In order to improve the expression level of recombinant BGL, the internal region from 1st to 19th amino acids in open reading frame of bgl was mutated in situ by inverse-PCR to replace the rare codons with the optimal codons of E. coli; the primers for the inverse-PCR were designated as bgl-3 and bgl-4 (Table (Table3).3). Inverse-PCR with primers was carried out using Pyrobest with pET-20-BGL as template, generating the plasmid pET-20-BGLII.
Plasmids pET-20-BG and pET-20-BGLII were transformed into E. coli JM109(DE3), and induced to expressed recombinant BGL by adding isopropyl-β-D-thiogalactopyranoside (IPTG) to final concentration of 0.8mM at OD600 about 0.7, and incubated further at 30°C for about 6h.
One liters of the recombinant cells carrying pET-20-BGLII were harvested by centrifugation at 5,000g for 10min at 4°C, and washed twice with distilled water, resuspended in 50mL of 5mM imidazole, 0.5mM NaCl, and 20mM Tris–HCl buffer (pH 7.9), and French-pressured for three times. The cell extracts were heat treated (60°C, 30min), and then cooled in an ice bath, and centrifuged (20,000g, 4°C, 30min). The resulting supernatants were loaded on to an immobilized metal affinity column (Novagen, USA), and eluded with 1M imidazole, 0.5M NaCl, and 20mM Tris–HCl buffer (pH 7.9). Protein was examined by SDS-PAGE , and the protein bands were analyzed by density scanning with an image analysis system (Bio-Rad, USA). Protein concentration was determined by the Bradford method using BSA as a standard.
The reaction mixture, containing 50mM imidole-potassium buffer (pH 6.4), 1mM p-nitrophenyl-β-D-glucopyranoside, and certain amount of β-glucosidase in 0.2mL, was incubated for 5min at 70°C. The reaction was stopped by adding 1mL of 1M Na2CO3. The absorbance of the mixture was measured at 405nm. One unit of enzyme activity was defined as the amount of enzyme necessary to liberate 1μmol of pNP per min under the assay conditions.
The optimum pH for activity β-glucosidase was determined by incubation at 70°C for 5min in the 50mM imidole-potassium buffer from pH 4.8 to 8.4. The optimum temperature for the enzyme activity was determined by standard assay ranging from 45 to 85°C in the 50mM imidole-potassium buffer, pH 6.0. The results were expressed as percentages of the activity obtained at either the optimum pH or the optimum temperature.
The pH stability of the enzyme was determined by measuring the remaining activity after incubating the enzyme (0.1μg) at 50°C for 1h in the 50mM imidole-potassium buffer from pH 5.2 to 8.0. To determine the effect of temperature on the stability of BGL, the enzyme (0.1μg) in the 50mM imidole-potassium buffer (pH 6.4) was pre-incubated for various times at 50°C, 65°C, 68°C and 70°C in the absence of the substrate. The activity of the enzyme without pre-incubation was defined as 100%.
The effects of metals and chemical agents on β-glucosidase activity of purified enzyme (0.1μg) were determined. Fe2+, Mg2+, Zn2+, Mn2+, Ca2+, K+, Al3+, Li2+, Cu2+, Co2+, and Hg2+ were assayed at concentrations of 1mM in the reaction mixture. The chemical agents EDTA (10mM) were assayed. The enzyme was incubated with each reagent for 10min at 50°C before addition of p-nitrophenyl-β-D-glucopyranoside to initiate the enzyme reaction. Activity was determined as described above and was expressed as a percentage of the activity obtained in the absence of the chemical agents and metal cations.
The substrate specificity of the enzyme (0.1μg) was tested by using following p-nitrophenyl-β-D-glucopyranoside, p-nitrophenyl-β-D-xylopyranoside, p-nitrophenyl-α-L-arabinofuranoside, maltose, sucrose, and cellobiose. Kinetic constant of BGL was determined by measuring the initial rates at various p-nitrophenyl-β-D-glucopyranoside concentrations (0.2, 0.4, 0.6, 0.8, 1, 2, and 4.0mM) or various cellobiose concentration (2, 4, 6, 8, 10, 12, 14, and 16mM) under standard reaction conditions. The Ki value of glucose was defined as amount of glucose required for inhibiting 50% of the β-glucosidase activity and was given as the averages of three separate experiments performed in duplicate.
The condon usage preference of E. coli in translation initiation region of pET-20-BGL was analyzed by using codon usage tool (http://gcua.schoedl.de/). The potential ORF of bgl was searched using the ORF search tool provided by the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov). Database searching was performed with Blast at NCBI and against CAZy (http://www.cazy.org). The active site of the enzyme was analyzed with the prosite tool (http://prosite.expasy.org/scanprosite). The multiple sequence alignment tool Clustal X2.0 was used for multiple protein sequence alignment . Sequences were further edited and aligned manually, when necessary, using the Mega 5 for editing. For phylogenetic analyses of conserved domains, sequences were trimmed so that only the relevant protein domains remained in the alignment . Phylogenetic relationships were inferred using the Neighbor-Joining (NJ) and Maximum-Parsimony (MP) method as implemented in Paup 4.0 for the NJ and MP trees, the results were evaluated with 1000 bootstrap replicates . The generated trees were displayed using TREEVIEW 1.6.6 (http://taxonomy.zoology.gla.ac.uk/rod/treeview.html).
The cellobiose was treated with purified BGL, and the degradation was subjected to analysis of thin-layer chromatography (TLC) and HPLC. The reaction mixture (20 μL) contained 290mM cellobiose, and 1μg of BGL in 50mM imidole-potassium buffer (pH 6.4). The reaction was performed for various times at 60°C, and stopped by heating for 5min in a boiling water bath. After centrifuged for 10min at 10,000g, supernatants of the reaction mixtures were applied on silica gel TLC plates (60F254, Merck Co.). Sugars on the plates were partitioned with a solvent system consisting of n-butanol, acetic acid, and water (2:1:1, by vol/vol), and detected using the orcinol reagent . The concentration of glucose was examined by HPLC on a carbohydrate analysis column (Waters Sugarpak1, USA) with water as a mobile phase.
Jianjun Pei and Qian Pang contributed equally to this work
The authors declare that they have no competing interests.
JP carried out the cloning and over-expression and drafted the manuscript. QP and SF helped to purify and characterize the BGL. LZ directed the over-all study and drafted the manuscript. HS helped to perform phylogenies analysis of β-glucosidases. All authors read and approved the final manuscript.
This work was supported by the National Natural Science Foundation of China (Grant No. 31070515 and 30871990) and A Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).