Original Research ARTICLE
Quantitative Proteomic Analysis of Germination of Nosema bombycis Spores under Extremely Alkaline Conditions
- Laboratory of Invertebrate Pathology, College of Animal Sciences, Zhejiang University, Hangzhou, China
The microsporidian Nosema bombycis is an obligate intracellular pathogen of the silkworm Bombyx mori, causing the epidemic disease Pebrine and extensive economic losses in sericulture. Although N. bombycis forms spores with rigid spore walls that protect against various environmental pressures, ingested spores germinate immediately under the extremely alkaline host gut condition (Lepidoptera gut pH > 10.5), which is a key developmental turning point from dormant state to infected state. However, to date this process remains poorly understood due to the complexity of the animal digestive tract and the lack of genetic tools for microsporidia. Here we show, using an in vitro spore germination model, how the proteome of N. bombycis changes during germination, analyse specific metabolic pathways employed in detail, and validate key functional proteins in vivo in silkworms. By a label-free quantitative proteomics approach that is directly based on high-resolution mass spectrometry (MS) data, a total of 1136 proteins were identified with high confidence, with 127 proteins being significantly changed in comparison to non-germinated spores. Among them, structural proteins including polar tube protein 1 and 3 and spore wall protein (SWP) 4 and 30 were found to be significantly down-regulated, but SWP9 significantly up-regulated. Some nucleases like polynucleotide kinase/phosphatase and flap endonucleases 1, together with a panel of hydrolases involved in protein degradation and RNA cleavage were overrepresented too upon germination, which implied that they might play important roles during spore germination. The differentially regulated trends of these genes were validated, respectively, by quantitative RT-PCR and 3 proteins of interest were confirmed by Western blotting analyses in vitro and in vivo. Furthermore, the pathway analysis showed that abundant up- and down-regulations appear involved in the glycolysis, pentose phosphate pathway, purine, and pyrimidine metabolism, suggesting preparations of energy generation and substance synthesis for the following invasion and proliferation inside the host. This report, to our knowledge, provides the first proteomic landscape of N. bombycis spores, and also a stepping stone on the way to further study of the unique infection mode of this economically important pathogen and other microsporidia in general.
Microsporidia are a type of obligate intracellular parasite related to fungi that include over 1500 species in ~187 genera (Vavra and Lukes, 2013). Microsporidia are important pathogens that infect a wide variety of animal groups, including vertebrates and invertebrates (Bhat et al., 2009). For instance, Nosema bombycis, the first reported microsporidian species, can cause a highly mortal disease referred to as pebrine in the silkworm Bombyx mori. Infecting through horizontal and vertical transmissions, N. bombycis proliferates continuously inside the host, which leads to chronic damage to all of larval tissues and organs including the intestines, silk glands, muscles, and malpighian tubules, and ultimately death. Thereby, N. bombycis represents a significant threat to sericulture worldwide (Becnel and Andreadis, 1999).
When surviving outside of the host, microsporidia form environmentally resistant spores that are protected by a thick two-layered wall, and are able to maintain metabolic activity (Vavra and Larsson, 1999). These mature spores, which are not germinated, are considered dormant spores. Upon coming into contact with host cells or being triggered by appropriate stimuli, dormant spores rapidly start the germination process. They fire a special invasion organelle, known as the polar tube, to penetrate the host cell membrane, and transfer the infective sporoplasm into the host cytoplasm, where proliferation and the next round of spore production occurs (Schottelius et al., 2000). As the first step in the microsporidia invasion, biological changes during the spore germination process have been suggested to play important roles (Troemel and Becnel, 2015). However, to date this process still remains poorly understood due to complex host-pathogen interactions and the lack of genetic tools for microsporidia.
By establishing a method for the in vitro spore germination, the polar tube firing process and invasion of the host cell have been observed in several microsporidian species (Troemel and Becnel, 2015). It has been broadly accepted that an appropriate stimulus such as pH, ionic species, osmotic pressure, or radiation can make microsporidian spores germinate, though the exact mechanism has not been fully elucidated.
A number of spore wall proteins (SWPs) have been identified and confirmed to be related to germination and invasion (Zhao W. et al., 2015). SWP5 has been confirmed to have the ability to protect spores from phagocytic uptake by cultured insect cells to modulate host cell invasion and interact with the polar tube (Cai et al., 2011; Li et al., 2012). An IG-like protein denoted BmTLP from B. mori has been proven to interact with SWP26, which is involved in endospore formation, host cell adherence, and infection in vitro (Li et al., 2009; Zhu et al., 2013). SWP16 was identified as a new exospore protein and is probably involved in N. bombycis spore adherence (Wang et al., 2015). However, the protein regulation mechanism of spore germination process is still largely unknown.
With the rapid development of proteomics technologies in combination with increased genome sequence information, proteomic analysis has become one of the most important strategies for characterizing expression profiles and pathway changes in diverse biological systems. Using traditional two-dimensional electrophoresis combined with mass spectrometry (MS), 80 proteins were found to be abundantly changed upon horsetail (Equisetum arvense L.) spore germination (Zhao Q. et al., 2015). Similarly, enzymes that are crucial to stress response were identified by this approach in thermophilic Archaea Thermococcus kodakarensis KOD1 (Jia et al., 2015). Recently, a promising label-free quantitative (LFQ) proteomic approach has been developed, in which peptide intensities measured during individual liquid chromatography (LC) runs are compared across runs (Tao et al., 2013). The advantage of this method is that it can be applied to any proteomic sample directly.
In the current study, we provides the first example of using label-free quantitation to investigate and compare the global proteomic changes during N. bombycis spore germination, which is regarded as the first step in host invasion. A total of 1136 proteins were identified with high confidence, with 127 proteins being significantly changed. Based on these data, proteins involved in germination and invasion were identified, as well as a number of important changes in metabolic pathways during the germination process. An in-depth understanding of spore germination could reveal useful insights into the unique mode of infection of this economically important pathogen and other microsporidia in general.
Materials and Methods
Insect Sample Collection and PH Measurements
The silkworms were provided by the Silkworm Germplasm Bank at the College of Animal Sciences, Zhejiang University, and reared under standard conditions (25° ± 1°C, 70 ± 5% humidity with a photoperiod of 14 h of light and 10 h of dark) using fresh mulberry leaves. The 5th instar larvae were anesthetized on ice for 30 min and dissected carefully. A pH microelectrode (Unisense, Denmark) with a tip diameter of 20–30 μm was calibrated with standard buffers at pH 4.0, 7.0, and 10.0 and directly applied to the freshly extracted guts as previously described (Sudakaran et al., 2012). The pH was recorded in three replicates in different gut regions.
Spore Purification and In vitro Germination Model
N. bombycis spores were propagated and purified from infected silkworms as previously described (Mei and Jin, 1989; Zhang et al., 2007). Firstly, 3rd instar molted larvae were fed mulberry leaves that were artificially contaminated by N. bombycis spores (~3 × 104 spores per larva). Heavily infected silk glands and midguts were collected from 5th instar larvae at the fifth day and gently homogenized in physiological saline with a glass homogenizer. Large silkworm debris was removed by filtrating twice through four layers of cheesecloth. After centrifugation (3000 g, 5 min), the sediment was resuspended in 20 mL of a solution containing 0.05% saponin and 0.05% (v/v) Triton X-100 in phosphate buffered saline (PBS). The cell suspension was passed through a plastic syringe with an inserted glass wool wad. A crude spore suspension was obtained by differential centrifugation after most of the contaminant was removed. The spores were further purified through continuous Percoll (pH = 7.2, GE Healthcare, Chicago, UK) gradient centrifugation (46,000 g, 60 min). After washes with distilled water, the spore purity was examined by phase-contrast microscopy. A total of 1 × 109 purified spores were divided into two parts on average and washed three times with PBS. A non-germinated spore (NGS) control sample was then treated with 200 μL of physiological saline. The germinated spores (GS) sample was treated with 200 μL of GKK germination buffer (0.05 M Glycine, 0.05 M KOH, and 0.375 M KCl, pH = 10.5), and both the NGS and GS samples were incubated at 30°C for 50 min. The experiment was performed for five independent biological replicates with the control and treatment groups.
Protein Isolation and Trypsin Digestion via the FASP Method
After germination, all 10 samples were then separately lysed with 1 mL of SDT-lysis buffer (4% SDS, 0.1 mol/L DTT, and 0.1 mol/L Tris-HCl, pH = 7.6) and 10 μL Protease Inhibitor Cocktail (Sangon, Shanghai, China) using acid-washed glass beads (diameter: 425–600 μm; Sigma, St Louis, USA) in a Precellys-24 (Bertin Technologies, Aix en Provence, France). The mixture was broken for 20 s on ice with 30 s intervals for 5 times and centrifuged at 14,000 g at 4°C. The protein concentration was detected by a Bicinchoninic Acid (BCA) Protein Assay Kit (Beyotime, Haimen, China). Fifty micrograms of each sample was digested by the FASP procedure (Wisniewski et al., 2010). DTT was added to the sample to a final concentration of 0.1 mol/L, and the mixture was incubated at 95°C for 5 min. Then, after adding 200 μL of UA buffer (8 M urea and 0.1 M Tris-HCl, pH = 8.5), each sample was concentrated in 30k Microcon filtration devices (Millipore, Boston, USA) and centrifuged at 14,000 g for 15 min at 20°C. Then, another 200 μL of UA buffer was added and the samples were centrifuged again at 14,000 g for 15 min. This step was repeated once. Next, the concentrate was mixed with 100 μL of 50 mM/L iodoacetamide (IAA) in UA buffer and incubated for an additional 20 min at room temperature in the dark. After that, the IAA was removed by centrifugation at 14,000 g for 15 min. One-hundred microliters of 50 mM/L NH4HCO3 was added and the samples were centrifuged at 14,000 g for 15 min; this step was repeated twice. Finally, 50 μL of 50 mM/L NH4HCO3, and trypsin (1:50, enzyme to protein) was added and incubated with the mixture at 37°C for 16 h. The tryptic peptide mixtures were collected for further analysis.
Reverse phase-high performance liquid chromatography (RP-HPLC) separation was achieved on an Easy-nLC nano-LC system (Thermo Fisher Scientific, Bremen, Germany) equipped with a self-packed column (75 μm × 150 mm; 3 μm ReproSil-Pur C18 beads, 120 Å, Dr. Maisch GmbH, Ammerbuch, Germany) at a flow rate of 250 nL/min. The RP-HPLC mobile phase A was 0.1% formic acid in water and mobile phase B was 0.1% formic acid in acetonitrile. The peptides were eluted using a gradient (2–80% mobile phase B) over a 3 h period. A nano-ESI Q-Exactive mass spectrometer (Thermo Fisher Scientific, Bremen, Germany) operated in data-dependent mode. The mass spectrometer was set such that each full MS scan followed by MS/MS for the 10 highest-intensity ions with the following parameters: precursor ion charges ≥ +2, a precursor ion isolation window of 2 Da, and a normalized collision energy of 30 for HCD. Dynamic Exclusion™ was set at 25 s. The full mass and subsequent MS/MS analyses were performed in an Orbitrap analyzer with R = 70,000 and R = 17,500 (at m/z 200), respectively. Five biological replicates were performed with NGS and GS.
The MS data were analyzed using the MaxQuant software (http://maxquant.org/, version 18.104.22.168). Carbamidomethyl (C) was set as a fixed modification, while oxidation (M, +15.99492 Da) was set as a variable modification. Proteins were identified by searching the MS and MS/MS data for peptides against the N. bombycis protein database from NCBI (http://www.ncbi.nlm.nih.gov/). Trypsin/P was selected as the digestive enzyme with two potential missed cleavages. Label-free quantification of NGS and GS was performed in MaxQuant as previously described (Cox et al., 2009). Protein abundance was calculated on the basis of the normalized spectral protein intensity (label-free quantitation intensity, LFQ intensity). The p-value of Log2LFQ intensity of each protein between the NGS and GS by two-tailed Student's t-test using the Perseus program, and p < 0.05 was set as a cutoff criterion for a significant change. The functional annotations of the identified proteins were initially assigned using the Blast2Go program (Conesa et al., 2005). BLASTp searches were first performed against the NCBI protein database and further analyses included Gene Ontology (GO) and Enzyme Code (EC) annotations. To investigate the significantly changing protein patterns in specific pathways, the Blast2Go program was used to locate the proteins in the Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway database (http://www.genome.jp/kegg; Conesa et al., 2005; Gotz et al., 2008). Briefly, all the amino acid sequences of 1136 identified proteins and 127 significantly changed proteins were firstly input into the Blast2Go program, respectively, then the mapping from KEGG database (online) was performed and the metabolic pathways of identified proteins were obtained.
RNA Extraction and Validation by Quantitative RT-PCR Analysis
After in vitro germination, the NGS and GS were separately lysed with 1 ml of TRIzol reagent (Invitrogen, Carlsbad, USA) using acid-washed glass beads in a Precellys-24. The total RNA from the broken spores was extracted according to the TRIzol Reagent manufacturer's protocol. cDNA synthesis for NGS and GS was performed using the PrimeScript™ RT Reagent Kit with gDNA Eraser (TaKaRa, Dalian, China). To validate the results of the label-free quantitation, quantitative real-time PCR (qRT-PCR) was performed for 16 randomly selected genes with significant changes in protein levels. Sixteen primers were designed using the Primer5 software (Premier Biosoft International) and all of the primer sequences were listed in Table S1. The internal control is β-tubulin. Quantitative real-time PCR (qRT-PCR) experiments were performed using the Roche LightCyler 480 system (Roche, Basel, Switzerland). Each reaction was conducted with a volume of 20 μL using LightCycler 480 SYBR Green 1 Master (Roche, Basel, Switzerland). All of the assays for a particular gene were performed with triplicate biological replicates and technical replicates. The expression level of each target gene was normalized to the level of β-tubulin. The fold change of each gene was calculated by the 2−ΔΔCT method (Schmittgen and Livak, 2008) and determined by three independent quantitative PCR amplifications of RNA that had been independently extracted.
Western Blot Analysis and Enzyme Activity Assay
Mulberry leaves (2.5 × 2.5 cm) were firstly disinfected by bleaching powder and artificially contaminated with N. bombycis spores. 8 × 107 spores were fed to each 5th instar larva (similar body weight). The larvae of “0 h” treatment group were collected when the larvae ingested mulberry leave completely, the larvae of “1” and “2 h” treatment group were collected at 1 and 2 h post-mulberry leaves were ingested, respectively. Then the spore protein of each treatment group was prepared. Briefly, the midgut was pulled out from the larvae to the ice cold PBS. The midgut and spores inside were lysed with 1 mL of SDT-lysis buffer and 10 μL Protease Inhibitor Cocktail (Sangon, Shanghai, China) using acid-washed glass beads in a Precellys-24 tissue homogenizer. The mixture was broken for 25 s on ice with 30 s intervals for 9 times and centrifuged at 14,000 g at 4°C. The total protein of “0,” “1,” and “2 h” were collected and protein concentration was determined using the BCA method.
Total Proteins of NGS, GS, “0,”, “1,” and “2 h” were resolved in 12% SDS-PAGE and transferred onto nitrocellulose membrane, blocked with 5% non-fat milk (in PBS), and incubated with the respective primary antibodies and beta-Actin (as a control for equal protein loading, HuaAn Biotechnology Company, Hangzhou, China) overnight at 4°C. The rabbit anti-SWP4 polyclonal antibody (diluted in 1:2000), rabbit anti-Flap endonucleases 1 polyclonal antibody (diluted in 1:2000) and rabbit anti-Pnkp polyclonal antibody (diluted in 1:2000) were produced by Genscript Biotechnology Company and the antigen sequence for respective antibody preparation were “CYVMKKKLKKSARKK,” “KENSKKGIVERPLSC,” and “CDKIPFMFDDSNFN,” respectively. The immunoblot detection was carried out using an ECL Western blotting analysis system using peroxidase conjugated secondary antibodies (Sigma, St Louis, USA). Every gel was loaded with molecular weight markers including proteins with molecular weight from 10 to 170 kDa (Bio-Rad, Hercules, USA).
Protease activity was measured using a commercial protease assay kit (Zhang et al., 2015). For NGS, purified spores (1 × 1010) were treated with 250 μL of physiological saline, while GKK germination buffer was used for GS. Both NGS and GS were incubated at 30°C for 30 min. Then, the supernatants were collected and the protein concentration was determined using the BCA method. Protease activities in two groups were measured by a Protease Assay™ Kit (G-Biosciences, St. Louis, USA), which uses a dye-labeled protein substrate. Briefly, 20 μL sample supernatant, 2.5 μL protease substrate solution, and 27.5 μL incubation buffer were mixed in a micro tube and incubated at 37°C for 24 h. After incubation, 50 μL precipitation agent was added in the tube and incubated at 37°C for 10 min. Then the samples were centrifuged at 12,000 g for 5 min and 80 μL supernatant was transferred to clean tubes with 120 μL assay buffer. The reaction was recorded at 570 nm against the blank. Chemically stabilized trypsin (MSG-Trypsin™) was supplied in the kit for preparing a general protease activity standard curve. The protease activity is proportional to the intensity of the reaction color, which can be quantified using a standard calibration curve.
The Extremely Alkaline Gut Environment of the Silkworm Stimulates the Germination of N. bombycis Spores
N. bombycis spores generally enter the silkworm host through the esophagus and are carried to the midgut where they germinate rapidly to deliver the germ cell (sporoplasm) into the host cell to initiate intracellular infection, which cause the pebrine disease in silkworm. After infections, larvae are inactive and slow in development. Later brown or black spots, a characteristic disease symptom, appear throughout their bodies and lead to death (Figure 1).
Figure 1. Pebrine disease of silkworms caused by N. bombycis. (A) A normal 5th instar silkworm. (B) Infected larvae show typical symptoms of disease: they are usually covered in brown or black dots and are unable to spin silkworm thread.
The midgut of silkworm, just like that of other herbivorous lepidopteran larvae, is a simple tube comprising the largest part of the whole body (Figure 2A). And midgut conditions are highly complex with various digestive enzymes, ion, and pH shifts, which might impact the N. bombycis spore germination. Among those factors, the extremely high pH is particularly noticeable. The silkworm midgut lumen provides a maximum pH of 11 in situ (Figure 2B), among the highest pH-values known to be generated by any biological system (Dow, 1992). The longitudinal pH profile observed in the gut region from anterior to posterior revealed that the pH rises to a peak in the middle region and declines again in the posterior region (Figure 2B).
Figure 2. The spore germination model of N. bombycis. (A) Dissected alimentary canal from a normal 5th instar larva of B. mori indicating the gut segmentation. (B) Extremely high pH in the gut lumen. The pH was measured in each gut region and data are shown as mean ± S.E.M. of at least three biological replicates. (C) Comparison between non-germinated (left) and germinated spores (right) in the N. bombycis spore germination model. N. bombycis spores were treated with GKK buffer (pH = 10.5) to initiate germination (polar tube firing) and their morphology were analyzed using phase-contrast microscopy. Black arrows indicate non-germinated spores (NGS); red arrows indicate germinated spores (GS). The scale bar represents 10 μM.
To test the contribution of pH on spore germination, N. bombycis spores were treated with a GKK germination buffer with a pH of 10.5 to induce germination (polar tube firing). Under phase-contrast microscopy, the non-germinated spores (NGS) and germinated spores (GS) were observed according to their unique developmental morphology. The bright and smooth spores (Figure 2C, black arrow) were non-germinated, whereas the dark spores (Figure 2C, red arrow) were germinated. Therefore, this simplified in vitro spore germination model successfully initiates the spore germination process and was employed to obtain the total protein from NGS and GS of N. bombycis for further proteomic investigation.
Label-Free Quantitative (LFQ)Proteomic Analysis of NGS and GS of N. bombycis
The total protein extracted from the NGS group and GS group of N. bombycis were analyzed in parallel using a LFQ proteomic approach. The relative LFQ intensities of the proteins across each individual samples were calculated using MaxQuant algorithms. The mass tolerance for peptide precursor ions and product ions was set at 6 and 20 ppm, respectively, and the false discovery rates at both peptide and protein levels were controlled to be <1%. These stringent filtration criteria provided highly confident protein identification. The p-value of Log2LFQ intensity of each protein between the NGS and GS groups were calculated using Student's t-test. As shown in the Figure 3A, the LFQ was highly reproducible between the five biological replicates in each treatment group and between the LC-MS/MS runs from the different treatment groups. The correlation between the normalized LFQ intensities was higher than 0.9 for each comparison of samples. The density plot of the log2 ratio between NGS group and GS group are similar to normal distribution (Figure 3B), which indicated that the experimental procedure was performed without systematic bias toward different samples. It thus allowed an accurate quantitative comparison of the protein expression in N. bombycis before and after the spore germination. After searching a non-redundant predicted protein dataset of N. bombycis, 1136 proteins were identified from analysis of the entire 10 LC-MS/MS raw data files of two treatment groups. The summary of all 1136 proteins is shown in Table S2.
Figure 3. Label-free quantitative proteomic analysis of N. bombycis NGS and GS. (A) Scatter plots and correlations between samples. (B) Density plot of the ratio between NGS and GS. (C) Volcano plot of the 1136 identified proteins between NGS and GS. A total of 127 proteins were identified as significantly changed proteins induced by germination, all of which exhibited a p < 0.05. Red represents the significantly up-regulated proteins, green the significantly down-regulated proteins, and blue the proteins that did not change significantly during spore germination. ***p < 0.01.
To investigate the exact protein content changes in N. bombycis, we performed a quantitative proteomic comparison between NGS and GS to better understand the important biological processes involved in this developmental stage. Using a t-test (p ≤ 0.05) as the variance for an up- or down-shift in abundance, the normalized ratio of protein abundance changes revealed significant and dynamic changes in protein profiles. Following this criteria of selection, 127 proteins were identified as significantly changed, with 60 proteins being up-regulated, and 67 proteins being down-regulated during the germination process. Volcano plot of all the identified proteins from NGS and GS is shown in Figure 3C. The molecular weight (MW), isoelectric point (PI), and grand average of hydropathicity (GRAVY) of the 127 significantly changed proteins were also investigated (Table S3). The top 20 significantly changed proteins with p ≤ 0.05 were listed in Table 1. The significantly up-regulated proteins were mainly infection-associated proteins (such as bifunctional polynucleotide phosphatase/kinase), spore proliferation-related proteins (transcriptional activator, Flap endonuclease 1-A, splicing factor and so on), spore wall proteins (putative SWP9), energy metabolism-related proteins (glyceraldehyde-3-phosphate dehydrogenase 2), proteasome-associated proteins (proteasome subunit alpha type-7, Proteasome component Y13), and so on. In contrast, the significantly down-regulated proteins mainly belonged to ribosomal proteins, SWP30 and putative SWP4, polar tube proteins, histone-binding proteins and so on (Table S3). To further validate the results from LFQ proteomic analysis, a set of 16 significantly changed proteins were selected randomly for quantitative RT-PCR analysis (Figure 4A). The expression patterns for all of these genes were consistent to the changes observed in related proteomic data. Scatter plots for LFQ intensity of 3 interesting proteins including SWP4, flap endonucleases 1 (FEN1) and polynucleotide kinase/phosphatase (PNKP; P < 0.01) were shown in Figure 4B. In addition, Western blotting also confirmed that two selected candidates, FEN1, and PNKP, were highly expressed in the GS group, while expression of SWP4 was higher in NGS (Figure 4C). Most importantly, FEN1 and PNKP were also highly expressed inside the gut of silkworms post infection, in contrast SWP4 was down regulated (Figure 4D).
Figure 4. Validation of significantly changed proteins by quantitative RT-PCR and Western blotting. (A) The fold change of 16 randomly selected genes in the two groups (NGS and GS) measured by quantitative RT-PCR. The internal control is β-tubulin. Hbp, histone-binding protein RBBP4; Pnkp, bifunctional polynucleotide phosphatase/kinase; Fen1, flap endonuclease 1-A; Ep1, Exportin-1; G3pd, glyceraldehyde-3-phosphate dehydrogenase 2; Pm3ep, Pre-mRNA 3-end-processing factor FIP1; Sfa, splicing factor; Ece, exosome complex exonuclease RRP40; Pmsh, PRE-mRNA splicing helicase; Swp30, spore wall protein 30; Ptp3, polar tube protein 3; S3, septin 3; Ef1, elongation factor 1-alpha; Swp4, spore wall protein 4; Ptps, protein transport protein SEC23; R5pi, ribose-5-phosphate isomerase A. The fold change for each gene was determined from biological triplicates and were presented as the mean ± SD. (B) Comparison of the signal intensities of 3 selected differentially (P < 0.01) expressed proteins (SWP4, FEN1, and PNKP) obtained from the LFQ proteomic analysis. (C) Western blot analysis was used for validation of abundance changes for the 3 selected differentially expressed proteins during N. bombycis spore germination in vitro. (D) In vivo validation of 3 proteins of interest inside the gut of B. mori at different time points post infection. β-actin was used as control.
Altogether the quantitative RT-PCR and western blot analyses confirmed the results of the quantitative proteomic analysis, which revealed important players during the spore germination of N. bombycis.
Functional Categorization and Enzyme Code Analysis
The identified protein sequences were assigned to Gene Ontology (GO) terms for functional classification. Three main categories of GO classification including biological process, molecular function, and cellular component, were analyzed separately to obtain the overall functional distribution of the identified proteins from NGS and GS. The resulting GO distributions for all 1136 proteins and for the 127 significantly changed proteins during germination were described in Figure 5, respectively. For all identified proteins, the most representative GO terms were organic substance metabolic process (GO:0071704), organic cyclic compound binding (GO:0097159), heterocyclic compound binding (GO:1901363), and cell part (GO:0044464). Whereas, the most representative GO terms for the 127 significantly changed proteins included: primary metabolic process (GO:0044238), cellular metabolic process (GO:0044237), and organic substance metabolic process (GO:0071704) within biological process; organic cyclic compound binding (GO:0097159), and heterocyclic compound binding (GO:1901363) within molecular function; and cell part (GO:0044464) within cellular component. Taken together, most of the significantly changed proteins were assigned to primary metabolic process, cellular metabolic process, and organic substance metabolic process, which indicates that the metabolic activity changed significantly after spore germination. Besides this, two proteins (transcriptional activator and transcription initiation factor) were identified in the transport (GO:0006810) and molecular function (GO:0003674) categories. Additionally, several ribosomal proteins were identified in the translation (GO:0006412) category, suggesting that the transcriptional and translation activities should be regulated after spore germination.
Figure 5. Functional categorization of identified proteins upon N. bombycis spore germination. A total of 1136 proteins were classified into 44 functional categories by Gene Ontology (GO) annotation (right panel). GO analysis of the 127 significantly changed proteins was showed in the left. GO terms are colored at level 2.
The 127 significantly changed proteins were further analyzed with enzyme codes (Figure 6A). Totally, 18 enzyme proteins were identified including 11 hydrolases, four transferases, one lyase, one isomerase, and one ligase. Interestingly, according to the LFQ analysis, most of enzyme proteins including nine hydrolases, three transferases, one lyase, and one ligase were significantly up-regulated after N. bombycis spore germination. For hydrolases, we paid attention to five proteinases and cleavage and polyadenylation specificity factor (CPSF) subunit 3. The proteasome degrades unneeded or damaged proteins by proteolysis to control protein stability in the organism. CPSF is involved in the cleavage of the 30027-signaling region from a newly synthesized pre-messenger RNA (pre-mRNA) molecule during gene transcription. Leucyl tRNA synthetase, which promotes the ligation of leucine to tRNA, was significantly up-regulated too after spore germination (Table S3). We further compared the specific activity of enzymes important for protein degradation from supernatants of GS and NGS. GS exhibited significant protease activity (Figure 6B).
Figure 6. Characterization of enzymes functioning in N. bombycis spore germination. (A) Enzyme classification (EC) codes for the significantly changed proteins upon N. bombycis spore germination. Note that one sequence can be associated with more than one EC term. (B) The protease activities of GS supernatant and NGS supernatant. The absorbance of the dye-labeled peptide was measured at 570 nm for determination of the protease activity, equivalent to that of trypsin. Triplicate measurements were performed. Data are shown as the mean ± SD. GS supernatant showed proteolytic activity. N. D., not detected.
Metabolic Pathway Analysis Mapping the Identified Proteins to KEGG Pathways
To identify the related biological pathways that are active during N. bombycis spore germination, all the 1136 identified proteins were searched against pathway collections in the KEGG database. All proteins were mapped into 37 metabolic pathways (Table S4). And the enrichment of the 127 significantly changed proteins in each pathway was obtained by calculating the probability using hypergeometric distribution. As shown in Table S4, most identified proteins were mapped to the purine metabolism and the number of mapped proteins was 96. Another 70 identified proteins were mapped to the thiamine metabolism. For the 127 significantly changed proteins, nine pathways (including Glycolysis/Gluconeogenesis, Pentose phosphate pathway, Purine metabolism, Pyrimidine metabolism, Biosynthesis of antibiotics, Fructose and mannose metabolism, Thiamine metabolism) were enriched (Table S5 and Dataset S6). Glycolysis/Gluconeogenesis and pentose phosphate pathways should participate in the energy metabolism process during spore germination and some key relevant proteins were identified (Figure 7). Almost all of the enzymes involved in glycolysis and pentose phosphate pathway (except ribulose-5P epimerase and hexokinase) and two subunits from the downstream pyruvate dehydrogenase enzyme were observed, with the majority of them showing an up-regulation. More importantly, three of these enzymes showed significantly differential regulation (P < 0.05) in protein expression levels after germination (Figure 7). Six enzymes were mapped to the purine metabolism, which are the most abundant enzymes for any of the pathways. Two enzymes were also mapped to the pyrimidine metabolism (Table 2). The two pathways are important for nucleotide synthesis.
Figure 7. Major metabolic pathways anticipated to participate in spore germination process in N. bombycis. The proteomic data was projected to the KEGG pathway database. Significantly changed pathways involve: glycolysis (blue), pentose-phosphate (orange), purine metabolism (green), pyrimidine metabolism (purple), and metabolism in mitosome (gray). The metabolites are in black. The up-regulated enzymes are highlighted in red and down-regulated proteins in green. Red and green arrows represent significantly changed proteins (p < 0.05).
As a turning point from the dormant state to the infected state and the first step in host cell invasion, spore germination is believed to be crucial for the success of microsporidian pathogens. A better understanding of the germination process therefore could make important contributions to the studies on this specific mode of infection and to the development of biological controlling strategies of these economically important and human-related pathogens. The entire protein coding sequences annotated in the N. bombycis genome is 4468 (Pan et al., 2013), and 1136 proteins were identified and quantified in our proteomic analysis of N. bombycis spores. In contrast, only 141 of the total 2573 predicted proteins were obtained in a spore germination study of the microsporidia Spraguea lophii using complex-mix proteomics (Campbell et al., 2013). Therefore, the LFQ proteomic approach employed in our study was sufficient for the microsporidia proteomic study.
Some major structural proteins involved in host-pathogen interactions have been characterized in the current work. We found that polar tube proteins (PTPs) 1 and 3 were significantly down-regulated during N. bombycis spore germination. It has been shown in another microsporidian Encephalitozoon cuniculi that PTPs accumulate in the sporogony stage but play roles in the germination process (Grisdale et al., 2013). Spore wall proteins, such as SWP26, SWP30, and SWP32, changed significantly too. It is known that SWP30 could not interact with the spore coat after alkali treatment (Yang et al., 2014). The N. bombycis germination process is generally triggered by alkali conditions inside the intestine of silkworms. Therefore, the decreased SWP30 level observed in our study is consistent with previous results. Similarly, SWP4 also showed a significant down-regulation, but its biological function has not been revealed until now. These proteins probably play important roles in the spore germination process and deserve further studies. A recent report showed that there was an interaction between SWP7 and SWP9: SWP9 in the chitin layer of the spore wall could act as a scaffolding protein that supports SWP7 and the two proteins were confirmed to mediate the infectious process by enabling spore adherence to host cells (Yang et al., 2015). Both SWP9 and SWP7 were confirmed to accumulate after germination at the proteomics level in our study, which might promote adherence to host cells. Currently many detection methods for N. bombycis SWPs have been established, such as monoclonal antibody (Li et al., 2015), and electrochemical immunosensor (Wang et al., 2014). The enrichment and purification methods of N. bombycis SWPs were also optimized (Yang et al., 2014; Zhao W. et al., 2015). However, the pertinence of detection is limited due to multiple types of SWPs. Spore germination is the important step for invasion. The significantly changed proteins obtained in this study might be related with spore germination and invasion. According to the analysis of SWPs here, the pertinence and efficiency of detection could be enhanced.
Other interesting proteins that showed an significant change in our study included a polynucleotide kinase/phosphatase (PNKP) and a RNA ligase 1 (RNL 1), both of which are the prototypal RNA healing and sealing enzymes involved in the RNA repair pathway (Novogrodsky et al., 1966; Novogrodsky and Hurwitz, 1966; Cameron and Uhlenbeck, 1977; Sugino et al., 1977). It is reported in phage that PNKP provides phages a way to evade RNA-damaging antiviral host response (Zhu et al., 2004). The increase of PNKP indicates that it might help N. bombycis evade the RNA-damaging mediated by host immune response and contribute to the normal growth of the pathogen inside. Flap endonucleases 1 (FEN1), found in diverse organisms, is a key nuclease that plays a role in DNA replication, repair, and recombination (Sayers, 1994; Ceska and Sayers, 1998; Liu et al., 2004). FEN1 was identified as a significantly up-regulated protein too upon spore germination. This increase might be related to the replication of double-stranded DNA during the pathogen proliferative phase. FEN1 could also take part in DNA repair pathway via the base excision repair pathway (Dianov and Lindahl, 1994), suggesting that it probably also protects N. bombycis from DNA-damage. These new findings provide additional perspectives for the complex B. mori-N. bombycis interactions.
According to the enzyme classification of LFQ data, 18 enzymes including five proteinases and the CPSF subunit 3 were found to be significantly up-regulated after spore germination. The protease activity was confirmed in the supernatant of GS. These released proteases may participate in the digestion of the host cell or the resistance of the host immune defense, and then make contributions to the spore infection and proliferation. As an important functional protein, the proteasome could play a role in cell cycle control (Chesnel et al., 2006), apoptosis (Haas et al., 1995; He et al., 2015), response to cellular stress (Garrido et al., 2006), and immunoreaction (Wang and Maldonado, 2006). CPSF is an important protein which could play a role in mRNA transcription. The significant up-regulations of these proteins indicate the levels of protein regulation are enhanced during the germination process. Furthermore, based on the significant up-regulation of leucyl tRNA synthetase, we inferred that there were strong synthetic demands for leucyl-related proteins in the spore germination. As a metabolic substrate, leucyl takes part in protein synthesis and could regulate and control gene transcription, translation, and post-translational modifications. Up to now, the role of leucyl in spore germination has been scarcely reported.
It is well-accepted that the microsporidian proliferative phase begins along with the activation of pathways related to energy generation and substance synthesis after spore germination. Comparative genomics research revealed that N. bombycis lacked genes that are required for oxidative phosphorylation, the tricarboxylic acid cycle, and fatty acid β-oxidation (Pan et al., 2013). Glycolysis, the pentose phosphate pathway, and a trehalose metabolism were thought to be the backbone for microsporidian energy metabolism because their mitochondria are massively reduced (Williams and Keeling, 2003; Van der Giezen et al., 2005). Mitosomes, which are a residual mitochondrion-derived organelle, have been regarded as the main organelle involved in energy metabolism in Encephalitozoon cuniculi (Katinka et al., 2001). In our study, almost all of the enzymes involved in glycolysis and the pentose phosphate pathway and two downstream pyruvate dehydrogenase enzymes were identified, furthermore three of these enzymes showed significantly differential regulation in protein expression levels after germination. Analysis of the two pathways provided the evidence to better understand the mode of energy generation in N. bombycis. We also identified 15 significantly up-regulated proteins involved in DNA transcription, translation, and transport after GO classification (Table S3). In contrast, histone-binding protein (RBBP4) was significantly down-regulated, which was connected to transcriptional silencing (Wolffe et al., 2000). With the help of the KEGG pathway analysis, most of the enzymes were mapped to the purine metabolism and two enzymes to the pyrimidine metabolism (Table 2). Nucleotide metabolism, which supports the DNA and RNA synthesis, is the very important metabolic pathway in free-living bacteria, Archaea, as well as microsporidia (Kanehisa et al., 2002). The nucleotides participate in a wide range of cellular functions, such as energy transduction, biomolecules syntheses, and signaling pathway (Welin and Nordlund, 2010). Purine and pyrimidine can be provided from de novo biosynthetic pathways or supplied via salvage pathways. Some microsporidia acquired genes involved in nucleic acid synthesis by horizontal gene transfer. The utilization of nucleotides and adaptability in host of microsporidia were enhanced by these genes (Alexander et al., 2016). In addition, some enzymes for the transformation of purine and pyrimidine during different phosphorylation and oxidation states were identified in microsporidia. These enzymes play key roles in the nucleotide metabolism too and help to better use of nucleotides for microsporidia (Heinz et al., 2014). Altogether, these results indicate that metabolic activities of N. bombycis are significantly enhanced during spore germination, which provide striking support for the following pathogen invasion and proliferation inside the host.
The spore germination is an intricate and elaborate process. Our study presents the first global quantitative proteomic analysis of N. bombycis spores during the germination process. Despite the proteome is considerably complex, a set of proteins related to cytoarchitecture, energy metabolism, and host immunity resistance have been identified as key players upon spore germination. Not only help to gain a better understanding of the microsporidia germination, this research is also a stepping stone on the way to further study of unique infection mode of this economically important pathogen and other microsporidia in general.
YS and XLu contributed to the initial design of the research. HL conducted the experiment and BC performed bioinformatics analyses with guidance from YS. SH and XLi collected samples. YS, XLu, and HL prepared the draft of this publication and all authors contributed to the subsequent stages of manuscript preparation.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors would like to thank Changwei Liu and Jing Zhao for providing advice in this research, and Jan Palmer for editorial assistance. This study has been supported by the National Natural Science Foundation of China (grant number 31601906 to YS), the “Hundred Talents Program to YS” from Zhejiang University, the Modern Agricultural Industrial Technology System (grant number CARS-22-ZJ0202), the Zhejiang Sericultural Sci-Tech Innovation Team (grant number 2011R50028) and the National Natural Science Foundation of China (grant number 31302033 to XLu).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.01459
Table S1. PCR primers used for the qRT-PCR validation.
Table S2. Summary of the 1136 identified proteins in NGS and GS using LFQ proteomic analysis.
Table S3. Summary of the 127 significantly changed proteins with p ≤ 0.05 (NGS/GS).
Table S4. KEGG mapping analysis for the 1136 proteins between NGS and GS.
Table S5. KEGG mapping analysis for the 127 significantly changed proteins between NGS and GS.
Dataset S6. KEGG pathway diagrams for the 127 significantly changed proteins.
Alexander, W. G., Wisecaver, J. H., Rokas, A., and Hittinger, C. T. (2016). Horizontally acquired genes in early-diverging pathogenic fungi enable the use of host nucleosides and nucleotides. Proc. Natl. Acad. Sci. U.S.A. 113, 4116–4121. doi: 10.1073/pnas.1517242113
Cai, S., Lu, X., Qiu, H., Li, M., and Feng, Z. (2011). Identification of a Nosema bombycis (Microsporidia) spore wall protein corresponding to spore phagocytosis. Parasitology 138, 1102–1109. doi: 10.1017/S0031182011000801
Campbell, S. E., Williams, T. A., Yousuf, A., Soanes, D. M., Paszkiewicz, K. H., and Williams, B. A. P. (2013). The genome of Spraguea lophii and the basis of host-microsporidian interactions. PLoS Genet. 9:e1003676. doi: 10.1371/journal.pgen.1003676
Chesnel, F., Bazile, F., Pascal, A., and Kubiak, J. Z. (2006). Cyclin B dissociation from CDK1 precedes its degradation upon MPF inactivation in mitotic extracts of Xenopus laevis embryos. Cell Cycle 5, 1687–1698. doi: 10.4161/cc.5.15.3123
Conesa, A., Gotz, S., Garcia-Gomez, J. M., Terol, J., Talon, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610
Cox, J., Matic, I., Hilger, M., Nagaraj, N., Selbach, M., Olsen, J. V., et al. (2009). A practical guide to the MaxQuant computational platform for SILAC-based quantitative proteomics. Nat. Protoc. 4, 698–705. doi: 10.1038/nprot.2009.36
Garrido, C., Brunet, M., Didelot, C., Zermati, Y., Schmitt, E., and Kroemer, G. (2006). Heat shock proteins 27 and 70: anti-apoptotic proteins with tumorigenic properties. Cell Cycle 5, 2592–2601. doi: 10.4161/cc.5.22.3448
Gotz, S., Garcia-Gomez, J. M., Terol, J., Williams, T. D., Nagaraj, S. H., Nueda, M. J., et al. (2008). High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 36, 3420–3435. doi: 10.1093/nar/gkn176
Grisdale, C. J., Bowers, L. C., Didier, E. S., and Fast, N. M. (2013). Transcriptome analysis of the parasite Encephalitozoon cuniculi: an in-depth examination of pre-mRNA splicing in a reduced eukaryote. BMC Genomics 14:207. doi: 10.1186/1471-2164-14-207
Haas, A. L., Baboshina, O., Williams, B., and Schwartz, L. M. (1995). Coordinated induction of the ubiquitin conjugation pathway accompanies the developmentally programmed death of insect skeletal muscle. J. Biol. Chem. 270, 9407–9412. doi: 10.1074/jbc.270.16.9407
He, X., Fu, Z., Li, M., Liu, H., Cai, S., Man, N., et al. (2015). Nosema bombycis (Microsporidia) suppresses apoptosis in BmN cells (Bombyx mori). Acta Biochim. Biophys. Sin. (Shanghai) 47, 696–702. doi: 10.1093/abbs/gmv062
Heinz, E., Hacker, C., Dean, P., Mifsud, J., Goldberg, A. V., Williams, T. A., et al. (2014). Plasma membrane-located purine nucleotide transport proteins are key components for host exploitation by microsporidian intracellular parasites. PLoS Pathog. 10:e1004547. doi: 10.1371/journal.ppat.1004547
Jia, B., Liu, J., Van Duyet, L., Sun, Y., Xuan, Y. H., and Cheong, G. W. (2015). Proteome profiling of heat, oxidative, and salt stress responses in Thermococcus kodakarensis KOD1. Front. Microbiol. 6:605. doi: 10.3389/fmicb.2015.00605
Katinka, M. D., Duprat, S., Cornillot, E., Metenier, G., Thomarat, F., Prensier, G., et al. (2001). Genome sequence and gene compaction of the eukaryote parasite Encephalitozoon cuniculi. Nature 414, 450–453. doi: 10.1038/35106579
Li, Y., Tao, M., Ma, F., Pan, G., Zhou, Z., and Wu, Z. (2015). A monoclonal antibody that tracks endospore formation in the microsporidium Nosema bombycis. PLoS ONE 10:e121884. doi: 10.1371/journal.pone.0121884
Li, Y., Wu, Z., Pan, G., He, W., Zhang, R., Hu, J., et al. (2009). Identification of a novel spore wall protein (SWP26) from microsporidia Nosema bombycis. Int. J. Parasitol. 39, 391–398. doi: 10.1016/j.ijpara.2008.08.011
Li, Z., Pan, G., Li, T., Huang, W., Chen, J., Geng, L., et al. (2012). SWP5, a spore wall protein, interacts with polar tube proteins in the parasitic microsporidian Nosema bombycis. Eukaryotic Cell 11, 229–237. doi: 10.1128/EC.05127-11
Novogrodsky, A., Tal, M., Traub, A., and Hurwitz, J. (1966). The enzymatic phosphorylation of ribonucleic acid and deoxyribonucleic acid. II. Further properties of the 5′-hydroxyl polynucleotide kinase. J. Biol. Chem. 241, 2933–2943.
Pan, G., Xu, J., Li, T., Xia, Q., Liu, S., Zhang, G., et al. (2013). Comparative genomics of parasitic silkworm microsporidia reveal an association between genome expansion and host adaptation. BMC Genomics 14:186. doi: 10.1186/1471-2164-14-186
Schottelius, J., Schmetz, C., Kock, N. P., Schuler, T., Sobottka, I., and Fleischer, B. (2000). Presentation by scanning electron microscopy of the life cycle of microsporidia of the genus Encephalitozoon. Microbes Infect. 2, 1401–1406. doi: 10.1016/S1286-4579(00)01293-4
Sudakaran, S., Salem, H., Kost, C., and Kaltenpoth, M. (2012). Geographical and ecological stability of the symbiotic mid-gut microbiota in European firebugs, Pyrrhocoris apterus (Hemiptera, Pyrrhocoridae). Mol. Ecol. 21, 6134–6151. doi: 10.1111/mec.12027
Tao, Y., Fang, L., Yang, Y., Jiang, H., Yang, H., Zhang, H., et al. (2013). Quantitative proteomic analysis reveals the neuroprotective effects of huperzine A for amyloid beta treated neuroblastoma N2a cells. Proteomics 13, 1314–1324. doi: 10.1002/pmic.201200437
Wang, Q., Song, Y., Chai, Y., Pan, G., Li, T., Yuan, Y., et al. (2014). Electrochemical immunosensor for detecting the spore wall protein of Nosema bombycis based on the amplification of hemin/G-quadruplex DNAzyme concatamers functionalized Pt@Pd nanowires. Biosens. Bioelectron. 60, 118–123. doi: 10.1016/j.bios.2014.03.075
Wang, Y., Dang, X., Ma, Q., Liu, F., Pan, G., Li, T., et al. (2015). Characterization of a novel spore wall protein NbSWP16 with proline-rich tandem repeats from Nosema bombycis (microsporidia). Parasitology 142, 534–542. doi: 10.1017/S0031182014001565
Welin, M., and Nordlund, P. (2010). Understanding specificity in metabolic pathways–structural biology of human nucleotide metabolism. Biochem. Biophys. Res. Commun. 396, 157–163. doi: 10.1016/j.bbrc.2010.04.054
Wisniewski, J. R., Nagaraj, N., Zougman, A., Gnad, F., and Mann, M. (2010). Brain phosphoproteome obtained by a FASP-based method reveals plasma membrane protein topology. J. Proteome Res. 9, 3280–3289. doi: 10.1021/pr1002214
Yang, D., Dang, X., Tian, R., Long, M., Li, C., Li, T., et al. (2014). Development of an approach to analyze the interaction between Nosema bombycis (microsporidia) deproteinated chitin spore coats and spore wall proteins. J. Invertebr. Pathol. 115, 1–7. doi: 10.1016/j.jip.2013.10.004
Yang, D., Pan, G., Dang, X., Shi, Y., Li, C., Peng, P., et al. (2015). Interaction and assembly of two novel proteins in the spore wall of the microsporidian species Nosema bombycis and their roles in adherence to and infection of host cells. Infect. Immun. 83, 1715–1731. doi: 10.1128/IAI.03155-14
Zhang, D., Li, S., Hu, L., Sheng, L., and Chen, L. (2015). Modulation of protease-activated receptor expression by Porphyromonas gingivalis in human gingival epithelial cells. BMC Oral Health 15:128. doi: 10.1186/s12903-015-0105-8
Zhang, F., Lu, X., Kumar, V. S., Zhu, H., Chen, H., Chen, Z., et al. (2007). Effects of a novel anti-exospore monoclonal antibody on microsporidial Nosema bombycis germination and reproduction in vitro. Parasitology 134, 1551–1558. doi: 10.1017/S0031182007002934
Zhao, Q., Gao, J., Suo, J., Chen, S., Wang, T., and Dai, S. (2015). Cytological and proteomic analyses of horsetail (Equisetum arvense L.) spore germination. Front. Plant Sci. 6:441. doi: 10.3389/fpls.2015.00441
Zhao, W., Hao, Y., Wang, L., Zhou, Z., and Li, Z. (2015). Development of a strategy for the identification of surface proteins in the pathogenic microsporidian Nosema bombycis. Parasitology 142, 865–878. doi: 10.1017/S0031182015000190
Zhu, F., Shen, Z., Hou, J., Zhang, J., Geng, T., Tang, X., et al. (2013). Identification of a protein interacting with the spore wall protein SWP26 of Nosema bombycis in a cultured BmN cell line of silkworm. Infect. Genet. Evol. 17, 38–45. doi: 10.1016/j.meegid.2013.03.02
Keywords: proteome, germination, Nosema bombycis, microsporidia, metabolic pathway
Citation: Liu H, Chen B, Hu S, Liang X, Lu X and Shao Y (2016) Quantitative Proteomic Analysis of Germination of Nosema bombycis Spores under Extremely Alkaline Conditions. Front. Microbiol. 7:1459. doi: 10.3389/fmicb.2016.01459
Received: 04 July 2016; Accepted: 31 August 2016;
Published: 21 September 2016.
Edited by:Baolei Jia, Chung-Ang University, South Korea
Reviewed by:Abu Hena Mostafa Kamal, University of Texas at Arlington, USA
Shaolei Teng, Howard University, USA
Copyright © 2016 Liu, Chen, Hu, Liang, Lu and Shao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.