- Open Access
Quantitative proteomics identifies 38 proteins that are differentially expressed in cucumber in response to cucumber green mottle mosaic virus infection
Virology Journalvolume 12, Article number: 216 (2015)
Since it was first reported in 1935, Cucumber green mottle mosaic virus (CGMMV) has become a serious pathogen in a range of cucurbit crops. The virus is generally transmitted by propagation materials, and to date no effective chemical or cultural methods of control have been developed to combat its spread. The current study presents a preliminary analysis of the pathogenic mechanisms from the perspective of protein expression levels in an infected cucumber host, with the objective of elucidating the infection process and potential strategies to reduce both the economic and yield losses associated with CGMMV.
Isobaric tags for relative and absolute quantitation (iTRAQ) technology coupled with liquid chromatography-tandem mass spectrometric (LC-MS/MS) were used to identify the differentially expressed proteins in cucumber plants infected with CGMMV compared with mock-inoculated plants. The functions of the proteins were deduced by functional annotation and their involvement in metabolic processes explored by KEGG pathway analysis to identify their interactions during CGMMV infection, while their in vivo expression was further verified by qPCR.
Infection by CGMMV altered both the expression level and absolute quantity of 38 proteins (fold change >0.6) in cucumber hosts. Of these, 23 were found to be up-regulated, while 15 were down-regulated. Gene ontology (GO) analysis revealed that 22 of the proteins had a combined function and were associated with molecular function (MF), biological process (BP) and cellular component (CC). Several other proteins had a dual function with 1, 7, and 2 proteins being associated with BP/CC, BP/MF, CC/MF, respectively. The remaining 3 proteins were only involved in MF. In addition, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified 18 proteins that were involved in 13 separate metabolic pathways. These pathways were subsequently merged to generate three network diagrams illustrating the interactions between the different pathways, while qPCR was used to track the changes in expression levels of the proteins identified at 3 time points during CGMMV infection. Taken together these results greatly expand our understanding of the relationships between CGMMV and cucumber hosts.
The results of the study indicate that CGMMV infection significantly changes the physiology of cucumbers, affecting the expression levels of individual proteins as well as entire metabolic pathways. The bioinformatic analysis also identified several pathogenesis-related (PR) proteins that could be useful in the development of disease-resistant plants.
Cucumber green mottle mosaic virus (CGMMV), which belongs to the Tobamovirus genus of the Virgaviridae family, is a damaging pathogen of cucurbit crops. The virus is transmitted mechanically, via seeds, pollen and other propagation materials , and produces typical mosaic patterning on the leaves of infected plants, as well as fruit distortions  that result in a reduced yield and lower market value . Since CGMMV was first reported in the UK in 1935 , it has been detected in many other countries including Israel, China, Greece, America and Canada [5–9]. Although the spread of this disease between different countries and regions has been rapid, to date there are no effective control strategies to reduce the economic losses it caused.
The susceptibility of particular plants to disease infection and the severity of the ensuing disease are closely associated with the defense response of the host plants . The invasion of plant pathogens activates various defense mechanisms intended to prevent colonization. These internal changes, often mediated at the transcriptional level, alter gene expression and can also indirectly affect plant performance [11, 12]. Understanding how viruses invade plants and the post-transcriptional mechanisms mediating the interactions between them is fundamental to developing effective management strategies for the control of viral disease in the field. Previous research has shown that both proteins and metabolites of the host plant are involved in the replication of viruses and the successful infection of plant cells . Proteomics is a useful technique for the study of these kinds of interactions and has been widely applied in the identification of proteins induced by abiotic and biotic stress and the characterization of their function during virus infections. For example, 2-DE and LC-MS/MS analysis identified many proteins produced when rice plants were infected with rice yellow mottle virus (RYMV), including those associated with stress related responses, as well as cellular metabolism and host mRNA translation . A similar study identified 203 proteins produced by sugar beet infected with beet necrotic yellow vein virus (BNYVV) , most of which were associated with photosynthesis and energy metabolism, in addition to host responses to stimulus and metabolism, while a proteomics approach has also been used to demonstrate the importance of the N resistance gene in the response of Nicotiana benthamiana to TMV infection . Many so-called effectors are generated during the interaction of pathogens and their hosts, while the proteins produced by phytopathogens themselves can enter plant cells and trigger host defense mechanisms. For example, it has been found that effector-triggered immunity (ETI) can subvert pathogen-associated molecular patterns (PAMPS) by the recognition of microbial proteins, which can result in the suppression of microbial growth via the activation of protein-mediated resistance . For example, silencing of the host ATP-synthase γ-subunit (AtpC) and Rubisco activase (RCA) genes in N. benthamiana resulted in the accumulation and pathogenicity of TMV in its leaves , while the overexpression of the hexose transporter gene LeHT1 in R tomato plants restricted the development of the tomato yellow leaf curl virus (TYLCV) . Research has shown that there is a great diversity of such effectors, which can be either intracellular or extracellular components, and that differential expression of the proteome and post-transcriptional regulation can play an important role during infection of host cells by viruses. Consequently, proteomics approaches to investigate host pathogen interactions have attracted increasing attention in recent years .
Proteomic analysis can generate an abundance of information on the individual proteins involved in biological responses, while the combination of differential-expression proteomics with novel techniques using isobaric tags for relative and absolute quantitation (iTRAQ) provide an opportunity to investigate proteins at the transcriptional level. High-throughput techniques such as iTRAQ are particularly valuable as they can be used to evaluate gene expression across the entire genome. In contrast to 2D electrophoresis, the iTRAQ method can be used to directly identify and compare the relative quantity of proteins derived from different tissues in a single experiment. The technique utilizes isobaric tags linked to different reporter ions, which attach to the side chains of lysine residues and the N-terminus of peptides resulting from a protein digest. Mass spectrometry can then be used to determine the ratio of labeled peptides from different sources and thus provide quantitative data . Furthermore, iTRAQ has the additional advantage of resolving very large (>200 kD) and very small molecules (<10 kD), which can be lost during electrophoresis, as well as reducing the degree of analytical bias as a result of the large number of isoelectric points . The application of iTRAQ therefore more closely reflects the actual situation during pathogen infection. The iTRAQ technique has been used to examine the differential expression of proteins in many studies, and the combination of iTRAQ and liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) can reveal detailed information regarding virulence factors and pathogenicity . Indeed, a recent study of Sugarcane mosaic virus (SCMV), identified 65 and 59 phosphoproteins in resistant and susceptible genotypes of maize (Zea mays L), respectively. Further analysis revealed that many of the phosphoproteins were differentially expressed during SCMV infection, which altered major cellular processes as a consequence of the phosphorylation status of the proteins. Such studies are extremely useful to expand our understanding of the mechanisms of plant-virus interactions . Another study of the Plum pox virus (PPV) in NahG-expressing N. benthamiana plants identified more than 1000 non-redundant proteins, 23 of which were differentially expressed in response to viruses expressing a modified P1 serine protease relative to the wild-type virus . The current study utilized a similar approach using iTRAQ coupled with LC-MS/MS to compare the proteome of cucumber plants infected with CGMMV with healthy mock-inoculated control plants.
Confirmation of CGMMV infection
At 45 days post inoculation (dpi), typical symptoms of CGMMV were observed in the leaves of the inoculated seedlings, which exhibited chlorosis and mottling (Fig. 1b). The mock-inoculated leaves (using phosphate buffer as inoculum) remained healthy and exhibited no symptoms of disease (Fig. 1a). The presence of CGMMV particles in the inoculated leaves was confirmed by SEM (Fig. 1c), while the CGMMV status of both the inoculated and mock-inoculated samples was confirmed by RT-PCR, in which only the inoculated samples produced PCR products (Fig. 1d).
Protein identification and quantification
The total protein yield from the CGMMV-infected plants was much lower than that of the mock-inoculated control, being 3.8 μg/μL and 7.4 μg/μL, respectively. The results of the iTRAQ analysis revealed that the protein profile of the two treatments also differed greatly and identified 191 unique proteins only expressed in the inoculated plants. Statistical analysis combining the data from the GO  and KEGG  analyses was used to interpret the novel proteomic data in relation to genomic information contained in databases such as SwissPROT (http://www.expasy.org), GenBank (http://www.ncbi.nlm.nih.gov/genbank/), EMBL (http://www.embl.org) and DDBJ (http://www.ddbj.nig.ac.jp/). The results revealed that 38 proteins had significantly different expression levels. Of these, 35 were assigned function by the GO analysis (http://www.geneontology.org). According to the results, 22 of the proteins had a combined function being assigned to all three categories: molecular function (MF), biological process (BP), and cellular component (CC). Several other proteins were found to have dual function with 1, 7 and 2 proteins being assigned to BP/CC, BP/MF, CC/MF, respectively, while 3 proteins were found to be associated with a single function, namely MF (Table 1). The KEGG analysis (http://www.genome.jp/kegg/) identified 18 proteins that had roles in 13 separate pathways. Five of the pathways including ko00360, ko00940, ko00680, ko03013 and ko04146 were found to involve the up-regulation of proteins. In particular, pathways ko00360, ko00940, and ko00680, which corresponded to phenylalanine metabolism, phenylpropanoid biosynthesis and methane metabolism, respectively, were associated with a substantial up-regulation of five peroxidase proteins (Q39650, P19135, Q40559, Q39652 and Q39653). However, the majority of the pathways including ko04745, ko04530, ko04810, ko04510, ko04520, ko00630, ko00195 and ko00190 were associated with down-regulation (Fig. 2, Table 2).
The thirteen pathways highlighted by the KEGG analysis were investigated further to determine their interrelationship during CGMMV infection. The KGML of each KEGG document was obtained from its corresponding KEGG webpage, and using the KEGG Parser software package [MathWorks, Massachusetts, U.S.] the information associated with the pathways and the data for individual proteins was extracted to an excel file [MergeKEGG, Kangyusheng, Shanghai, China]. The merged excel file was then imported into Cytoscape software [http://apps.cytoscape.org/apps/kgmlreader], which generated 3 separate network diagrams (Fig. 4). The nodes in the diagrams correspond to the individual metabolites associated with the pathways identified by the KEGG analysis, while the black lines illustrate the interrelationships between them. The results therefore show how the 13 KEGG pathways interact during CGMMV infection.
Expression levels of 25 proteins identified by iTRAQ during CGMMV infection
The expression levels of 25 proteins selected randomly from the iTRAQ analysis were assessed at three time points post inoculation (1, 10 and 30 dpi). The levels of several proteins including P08216, Q4VZI6, Q4VZH2, C1M2W0, A4GWU2, B1PE19, Q8RW69, P19135, P29602, A1BQL5 and Q5DJS5 were found to decline, while the levels of Q4VZK4 and O04379 were found to increase. Other proteins were found to exhibit more complex patterns of expression with Q2QD76, Q9SXL9, A1BQK8, Q9FRW1, B6V8E7, Q2MK92, Q8LK13 and Q94G09 being found to initially rise before falling, while P62776, Q84V66, A3F9M6 and P83966 were down-regulated before being up-regulated (Fig. 3).
The current study identified thirty-eight proteins with significantly altered expression in cucumber plants infected with CGMMV. Gene Ontology (GO) analysis (Table 1) classified the proteins into three groups based on their function, data that could be useful when investigating the function of proteins in vivo and screening their impact during CGMMV infection and ultimately in the development of disease resistant plants. The individual gene products were not allocated to single cellular processes, but were involved in several biological processes with multiple molecular functions [28, 29]. This effect is not uncommon; for example, the C-function activity associated with organ differentiation in Medicago truncatula is controlled by two euAG genes, MtAGa and MtAGb , while conversely it has also been shown that single genes can control multiple phenotypes and different characteristics . The down-regulation of cytochrome b6 (Q4VZK4, −1.4 fold change) in the current study, which affected the photosynthesis pathway (ko00195), is a further example of this phenomenon. Cytochrome b6 is an essential component of the cytochrome b6-f complex, which mediates electron transfer between photosystem I (PSI) and photosystem II (PSII), cyclic electron flow around PSI, as well as being an integral component of the chloroplast and thylakoid [31, 32].
Both in vitro and in vivo studies have demonstrated that the stresses caused by pathogens can alter the physiological characteristics of host plants . The current study found that cucumber plants infected with CGMMV had significantly altered levels of phenylalanine metabolism (ko00360) and phenylpropanoid biosynthesis (ko00940). Phenylalanine was the most important component in the ko00360 and ko00940 pathways, and its metabolites such as flavonoids and lignin can play key roles in the resistance of plants to pathogen infection. Flavonoids are an important group of plant secondary metabolites that have been studied extensively using both physiological and molecular approaches . It has been found that, in addition to acting as attractants for pollinators and seed dispersal agents as a result of their role in flower and fruit pigmentation, they also affect fertility and disease resistance . For example, it has been shown that high levels of flavonoids in Arabidopsis thaliana can affect auxin transport and repress plant growth . Similarly it has also been documented that the accumulation of high levels of lignin in cotton is linked to resistance to Verticillium dahliae . Such studies indicate that the altered plant growth characteristics resulting from increased levels of flavonoids and lignin would be useful traits to consider when developing disease resistant varieties. Furthermore, the current study also found that the effects on phenylalanine metabolism (ko00360) and phenylpropanoid biosynthesis (ko00940) were linked to the well-characterized phenylalanine ammonia lyase (PAL) pathway (k10775), which is known to promote the metabolism of phenylalanine into monolignol, a process which is thought to play a key role in the cell wall lignification that enhances disease resistance by creating a barrier to the further invasion by pathogens . Indeed, it has been shown that the over-expression of PAL in tobacco can enhance host resistance to infection by the fungus Cercospora nicotianae . Plant peroxidases (PO) are also involved in lignin formation and represent another important factor in the resistance of plants to pathogen attack . The current study found that both the ko00360 and ko00940 pathways were involved in the response of cucumbers to infection by CGMMV, with their expression being up-regulated 4.7 fold. Furthermore, the KEGG analysis revealed that the pathways associated with ko00360 and ko00680, which involved the up-regulation of 5 separate proteins, were the most significantly affected of all the pathways identified. All of the up-regulated proteins in these pathways were peroxidases, which implies that the over-expression of POs could be an important factor in the disease response of cucumber to CGMMV infection.
Several other proteins were also found to be up-regulated in the inoculated cucumber plants including the small ras-related protein (A1BQL2) and elongation factor 1-alpha (A4GWU2), which were both associated with the RNA transport pathway (ko03013), and it is interesting to note that both proteins have been associated with cucumber resistance in previous studies [41, 42]. Although not highlighted by the KEGG analysis, the up-regulation of phospholipase D (A3F9M6) was also of interest given that it has been associated with the defense signaling of Arabidopsis infected with powdery mildew (Blumeria graminis f. sp. hordei) , which indicates that it could also be an important component of the defense response of cucumbers.
However, the majority of the pathways (8 in total) were associated with the down-regulation of proteins. For example, it was found that the down-regulation of both isoforms of actin (B8XX44, B8XX45) impacted many pathways in the KEGG analysis (5 in total). Actin, which is evolutionarily highly conserved, is a key component of the cytoskeleton in all eukaryotes [44, 45] and plays an important role in protein-protein interactions [46, 47], as well as being involved in developmental processes, membrane trafficking and disease resistance [46, 48]. Meanwhile, the actin analogues in prokaryotes function as a cytoskeleton  that stabilizes the shape of the cell . In the current study, both isoforms of (B8XX and BXX5) were found to be involved in the phototransduction pathway, which was significantly up-regulated. Given that previous studies have demonstrated that this process can affect the expression of chloroplastic genes , it is possible that the altered expression of actin could be associated with the chlorotic symptoms observed in CGMMV-infected leaves. It was also noted that both cytochrome b and cytochrome b6 were significantly down-regulated in the CGMMV infected plants. Cytochrome b (cyt b) is a membrane protein, which forms the core of the mitochondrial bc 1 (complex III) in the respiratory chain, and is responsible for transmembrane electron transfer. The cyt b subunit possesses many inhibitor sites that can dramatically affect oxidoreductase activity as well as facilitate the action of quinone antagonists . Consequently, cyt b has frequently attracted the attention of microbiologists , and it has been shown that in addition to its role in respiratory electron transport in the phytopathogenic bacteria Pseudomonas cichorii and Pseudomonas aptata , it is also associated with their pathogenicity . In addition, previous studies have also shown that cyt b is associated with fungicide resistance, with two separate point mutations in the cyt b gene resulting in resistance to Qo inhibitors (QoIs) , which inhibit mitochondrial respiration by binding to the Qo site of complex III . It is also interesting to note that the protein CmPP36, a member of the cytochrome b5 family, is required for the shuttle of Red clover necrotic mosaic virus (RCNMV) particles into the phloem translocation pathway , which has been shown to facilitate the systemic movement of CGMMV in infected cucumber plants . This aspect of cyt b associated metabolism provides an intriguing opportunity for the development of control strategies to prevent CGMMV infection.
Meanwhile, the analogous protein cyt b6 is a component of the cyt b6f complex, which is associated with both PSI and PSII, and functions as a proton translocase in the cyclic electron transport of PSI. The cyt b6f complex is comprised of both nuclear and plastid encoded subunits, and it has been shown that nuclear mutants can result in defects in the cyt b6f complex that impair photosynthesis due to the uncoupling of electron transport . Previous research has also shown that some viruses can affect the structure and function of the photosystems, reducing the rate of photosynthesis by as much as 60 %. Moreover, as an indispensable physiological process in all green plants, the altered photosystems also result in a reduction in the synthesis of cellular components in the affected plants [61–65], especially the synthesis of photosynthetic pigments, which has been associated with the symptom of virus infection in Nicotiana tabacum . It is therefore possible that similar effects could explain the chlorosis and mosaic patterning in the leaves of CGMMV-infected cucumbers, while the alterations to actin and cytochrome b could account for the abnormal shape of the plants and fruits, which results in an indirect loss of yield .
Two of the most significantly down-regulated proteins were histones (P62776, P48557). Histone acetylation is important in the activation of genes at the posttranslational level, and a previous study  has shown that histone acetyltransferase (AtHAC1) can regulate flowering time in Arabidopsis via inhibition of Flowering Locus C (FLC). It is therefore possible that the down-regulation of the histones in cucumbers infected with CGMMV could be responsible for the delayed flowing that can indirectly affect production. It is also interesting that both malate synthase (P08216) and isocitrate lyase (ICL) (P49296) were down regulated in plants inoculated with CGMMV. Previous studies have shown that endogenous malate synthase is exploited by Stagonospora nodorum during its colonization of wheat , while isocitrate lyase is a key component of the glyoxylate cycle , and has been found to be the most important factor when Leptosphaeria maculans successfully colonizes canola . It was also interesting to note that galactinol synthase (Q84V66) was down-regulated in the CGMMV inoculated plants, given that a recent study has shown that the overexpression of cucumber galactinol synthase (CsGolS) improved resistance to Botrytis cinerea in transgenic tobacco plants .
The relationships between host and pathogen have long been recognized to have a great influence on the evolution of the host and its biological diversity  as well as it resistance to pathogens . The current study utilized qPCR to validate the expression of the differentially expressed proteins identified by the iTRAQ LC-MS/MS analysis at the transcriptional level, and assess the expression levels of the corresponding genes at 3 time points during CGMMV infection. All of the genes assessed exhibited altered patterns of expression compared to the mock-inoculated control, indicating that these genes were involved in the disease-response of cucumber plants to CGMMV infection. Of particular interest were the genes corresponding to O04379 (Protein argonaute 1) and Q4VZK4 (Cytochrome b6–f complex subunit 4), which exhibited increasing levels of expression as the CGMMV infection progressed. Previous research  has shown that protein argonaute 1 (AGO1), which is an essential component in a multitude of developmental processes including those of leaves, floral organs and the formation of axillary meristems, also has a role in RNA-mediated post-transcriptional gene silencing (PTGS) and antiviral RNA silencing, which indicates that the elevated levels of O04379 could be associated with CGMMV-resistance in cucumber. In contrast, the increased expression of Q4VZK4, which is a component of the cytochrome b6-f complex (cyt-b6f) that mediates electron transfer between photosystem II(PSII) and photosystem I(PSI) , is more likely to be associated with the disease-response of the cucumber host and the symptoms of CGMMV infection. Cyt-b6f is located in the thylakoid membranes of plant chloroplasts , and previous research has shown that the coat protein of tobacco mosaic virus (TMV-CP), which like CGMMV belongs to the Tobamovirus genus, can become embedded inside or otherwise complex with the thylakoid membranes impairing their function .
While GO annotation can establish the basic function of particular proteins, KEGG analysis can be used as an indication of the proteins role within a network of interacting molecules in the cell. This can be a very useful tool to elucidate the role of genes of interest. The merged network diagrams generated in the current study highlight the relationships between the annotated proteins and the metabolic pathways that they impact. The current study found 38 proteins that were differentially expressed in CGMMV-infected cucumber and highlighted the role of 18 differentially expressed proteins in the context of 13 separate metabolic pathways. Based on the results of the GO annotation and the KEGG analysis these 13 KEGG pathways were selected for the merged analysis performed using Cytoscape software, which provided an insight into how the different KEGG pathways interacted (Fig. 4). The results highlighted several pathways that were significantly affected by CGMMV infection including phenylalanine metabolism (ko00360), phenylpropanoid biosynthesis (ko00940), and methane metabolism (ko00680), which were all characterized by a 4.7 fold up-regulation, as well as glyoxylate and dicarboxylate metabolism (ko00630), which were 2.2 fold down-regulated. The regulatory mechanisms induced during the response of plants to disease infection are complex and operate at physiological, cellular and transcriptome levels . The results of the current study showed that CGMMV infection affected not only regulatory mechanisms specific to disease response but also had a crossover effect with other regulatory mechanisms. For example, L-phenylalanine and L-tyrosine both mediated the interactions between pathways ko00360 and ko00940, while pyruvate was associated with the function of the ko00360 and ko00680 pathways (Fig. 4a). Statistical analysis revealed that pathways ko00360 and ko00940 were the most significantly affected producing a p-value of 0.0026. Furthermore, the network diagrams illustrated in Fig. 4 shown that pathway ko04745 can affect other pathways including ko04510, ko04520 and ko04530 via protein-protein interactions (Fig. 4b). The network of relationships identified by such studies provides valuable data regarding plant pathogen interactions and can highlight potential areas for further study. Furthermore, they also suggest that research should not be focused on single regulatory factors but rather the interaction of a range of regulatory mechanisms via various processes including gene expression as well as enzyme-enzyme, protein-protein, and protein-metabolite interactions. It is likely that such an approach would greatly increase our understanding of the regulatory network operating during the CGMMV infection of cucumber plants, and could provide useful clues that could facilitate the long term control of CGMMV by the development of resistant varieties.
Sample preparation and protein extraction
Cucumber (Cucumis sativus) seedlings (cv. Zhongnong 16) obtained from the Institute of Vegetables and Flowers at the Chinese Academy of Agricultural Sciences, (Beijing, China) were grown in sandy soil fertilized with chicken manure. Plants were inoculated with CGMMV obtained from the Chinese Academy of Inspection and Quarantine (CAIQ) at the three true-leaf stage. The plants were grown in a greenhouse at the China Agricultural University under natural light conditions with the temperature maintained between 21 °C and 32 °C, and watered when necessary. The plants were kept in net cages (Zhiguang wire mesh products, Hebei, China) to prevent insect feeding and cross infection by other pathogens. Mock-inoculated control plants inoculated with phosphate buffer (0.05 M, pH7.0) were prepared and maintained under identical conditions. The plants were assessed after 45 days when the inoculated plants displayed typical symptoms of CGMMV infection including yellowing, as well as mottled and mosaic patterning of the leaves. The presence of CGMMV in leaves exhibiting symptoms of disease was confirmed by SEM and RT-PCR as described in a previous study . The proteins from the leaves of the control and CGMMV-positive plants were extracted according to the protocol of Fan  and quantified using the RC DC™ Protein Assay Kit [Bio-Rad, CA, U.S.] according to the protocol of the manufacturer.
Protein processing and iTRAQ labeling
The extracted protein samples were transferred to fresh 1.5 mL centrifuge tubes and suspended in 250 μL sample buffer (8 M UREA, 4 % CHAPS, 30 mM HEPES, 1 mM PMSF, 2 mM EDTA, 10 mM DTT). The samples were dissolved completely using a KQ-250DE ultrasonic sonicator (Kun Shan, China) set to pulse on 2 s, pulse off 3 s, power 180 W. After centrifugation at 12,000 rpm for 25 min, the supernatant was collected and transferred to fresh tubes before 20 μL of 55 mM iodoacetamide (IAM) was added to each tube. After dark-incubation in a water bath at 56 °C for 1 h, four times the sample volume of cold acetone was added to each tube, and the samples precipitated at −20 °C for 3 h. The samples were then centrifuged at 13,000 rpm for 15 min at 4 °C and the supernatant discarded. The protein pellet produced was dissolved in 300 μL dissolution buffer (0.5 M TEAB, 0.1 % SDS) and the quantity of the protein determined by the Bradford method. The prepared soluble protein samples were stored at −80 °C until required.
Protein digests were prepared by adding 3 μg trypsin Gold [Promega, Madison, U.S.] to 100 μg of extracted proteins using a 0.5 M TEAB, 0.1 % SDS solution to balance the sample volume, and incubated at 37 °C for 24 h. The digested samples were re-dissolved in 0.5 M TEAB and labeled with isobaric tags using the iTRAQ Reagents Multiplex Kit [Applied Biosystems Inc., CA, U.S.] according to the protocol of the manufacturer. The samples were labeled using 119 and 121 tags for the CGMMV-positive and mock-inoculated control samples, respectively.
Separation and identification of peptide fragments
The two samples, which had been labeled separately with either the 119 or 121 tag, were mixed, re-suspended and diluted 10-fold in 25 % acetonitrile (ACN), 10 mM potassium dihydrogen phosphate (KH2PO4), adjusted to pH 3.0 using phosphoric acid, and purified using a C18 reverse-phase LC column (Shimadzu corporation, Japan) to remove salts and the by-products of the iTRAQ reaction. The samples were washed and eluted according to the protocol of the manufacturer, and the eluates dried down and fractionated using a Luna®5 μm, 100 A, 250 × 4.60 mm strong cation exchange (SCX) column [Phenomenex Inc., CA, U.S.]. The peptides were eluted at a flow rate of 1.00 mL min−1, with a linear gradient of 0–100 % solvent B (10 mM KH2PO4, 2 M KCl, 25 % ACN, pH3.0) over 86.1 min, and the peptide fraction collected for use in the subsequent analysis.
The peptide fractions were centrifuged and dried at low temperature before being re-suspended in 0.1 % formic acid, and ionized by matrix-assisted laser desorption/ionization (MALDI). The resulting ions were then detected using the Agilent 6550 iFunnel Q-TOF LC/MS mass spectrometry system (Agilent Technologies Inc., U.S.). The data collected was used to search online databases using the MASCOT software package (Matrix Science, Inc., Boston, U.S.) to identify the component proteins. The data and search results were then imported into Scaffold software (http://www.proteomesoftware.com/) [Proteome Software, Inc., Portland, U.S.] for protein quantification, and the results verified by analysis of variance (ANOVA) and the student’s T-test using the SPSS 19.0 software package [IBM, CA, U.S.]. The variation generated between the CGMMV-infected and mock-inoculated samples was used to estimate the level of differential protein expression. The analysis was made using two biological replicates, with an absolute value for the fold change of greater than 0.6 indicating a differentially expressed protein.
qPCR analysis of 25 differentially expressed proteins identified by iTRAQ analysis
The leaves (100 mg) of CGMMV-infected and mock-inoculated cucumber plants were collected at three time points post inoculation: 1, 15 and 30 dpi. The samples were ground to a fine powder in liquid nitrogen using a pestle and mortar and the total RNA extracted using TRIzol Reagent (Invitrogen Life Technologies) before the cDNA was prepared (cDNA synthesis kit, Invitrogen, U.S.) using 3 μg RNA according to the protocol of the manufacturer. The primers listed in Table 3 were used to compare the expression levels of 25 proteins selected at random from the 38 proteins identified by the iTRAQ analysis using tubulin as the internal standard. The qPCR itself was conducted in 25 μl reaction mixtures containing 1.0 μl cDNA, 0.5 μl each of the forward and reverse primers (10 μM), 12.5 μl 2 × TaqMix, and 1.25 μl SYBR Green (BIO-RAD, U.S.), and processed using a real-time PCR thermocycler (BIO-RAD, U.S.) with the following program: denaturation at 95 °C for 5 mins followed by 40 cycles of denaturation at 95 °C for 15 s, annealing for 20 s, and extension at 72 °C for 20 s. The annealing temperature used was specific to individual primer sets depending on their GC content (Table 3). A total of three biological replicates were performed (SD <0.2).
In conlusion, total of 38 proteins were identified for the first time, and the expression level of 23 proteins were up-regulated and 15 proteins down-regulated in cucumber in response to CGMMV infection. Among of them, 18 proteins were involved in 13 individual metabolic pathways. The analysis results that 13 pathways also have crossing relation between of each other. These results provide the basic of investing the differences in protein expression and bring new insights into the relation of CGMMV and cucumber.
cucumber green mottle mosaic virus
isobaric tags for relative and absolute quantitation
kyoto encyclopedia of genes and genomes
liquid chromatography-tandem mass spectrometric
Mink GI. Pollen- and seed-transmitted viruses and viroids. Annu Rev Phytopathol. 1993;31:375–402.
Liu HW, Luo LX, Li JQ, Liu PF, Chen XY, Hao JJ. Pollen and seed transmission of cucumber green mottle mosaic virus in cucumber. Plant Pathol. 2014;63:72–7.
Shim CK, Han KS, Lee JH, Bae DW, Kim DK, Kim HK. Isolation and characterization of watermelon isolate of Cucumber green mottle mosaic virus (CGMMV-HY1) from watermelon plants with severe mottle mosaic symptoms. J Plant Pathol. 2005;21:167–71.
Ainsworth GC. Mosaic diseases of the cucumber. Ann Appl Biol. 1935;22:55–67.
Antignus Y, Pearlsman M, Ben-Yoseph R, Cohen S. Occurrence of variant of cucumber green mottle mosaic virus in Israel. Phys Chem Chem Phys. 2002;18:50–6.
Varveri C, Vassilakos N, Bem F. Characterization and detection of cucumber green mottle mosaic virus in Greece. Phys Chem Chem Phys. 2002;30:493–501.
Zhang YJ, Li GF, Li MF. Occurrence of cucumber green mottle mosaic virus on cucurbitaceous plants in China. Plant Dis. 2009;93:200.
Baker C. Cucumber green mottle mosaic virus (CGMMV) found in the United States (California) in melon. Pest alert, Florida department of agriculture and consumer services, division of plant industry (DACS-P-01863). 2013. http://www.freshfromflorida.com/content/search?SearchText=CGMMV
Ling KS, Li R, Zhang W. First Report of Cucumber green mottle mosaic virus Infecting Greenhouse Cucumber in Canada. Plant Dis. 2014;98(5):701.
O’Donnell PJ, Schmelz EA, Moussatche P, Lund ST, Jones JB, Klee HJ. Susceptible to intolerance -a range of hormonal actions in a susceptible Arabidopsis pathogen response. Plant J. 2003;33:245–57.
Jones JD, Dangl JL. The plant immune system. Nature. 2006;444:323–9.
Vogt T. Phenylpropanoid biosynthesis. Mol Plant. 2010;1:2–20.
Xu K, Nagy PD. Dissecting virus-plant interactions through proteomics approaches. Curr proteomics. 2010;7(4):316–27.
Ventelon-Debout M, Delalande F, Brizard JP, Diemer H, Van Dorsselaer A, Brugidou C. Proteome analysis of cultivar-specific deregulations of Oryza sativa indica and O. sativa japonica cellular suspensions undergoing rice yellow mottle virus infection. Proteomics. 2004;4:216–25.
Weeb KM, Broccardo CJ, Prenni JE, Wintermantel WM. Proteomic profiling of sugar beet (Beta vulgaris) leaves during rhizomania compatible interactions. Protromes. 2014;2:208–23.
Caplan JL, Zhu X, Mamillapalli P, Marathe R, Anandalakshmi R, Dinesh-Kumar SP. Induced ER chaperones regulate a receptor-like kinase to mediate antiviral innate immune response in plants. Cells host microbe. 2009;6(5):457–69.
Chisholm ST, Coaker G, Day B, Staskawicz BJ. Host-microbe interactions: shaping the evolution of the plant immune. Cell. 2006;124(4):803–14.
Bhat S, Folimonova SY, Cole A, Ballard KD, Lei ZT, Watson B, et al. Influence of host chloroplast proteins on tobacco mosaic virus accumulation and intercellar movement. Plant Physiol. 2013;161:134–47.
Czosnek H, Eybishtz A, Sade D, Gorovits R, Sobol I, Bejarano E, et al. Discovering host genes involved in the infection by the tomato yellow leaf curl virus complex and in the establishment of resistance to the virus using tobacco rattle virus-based post transcriptional gene silencing. Viruses. 2013;5:998–1022.
Jenner RG, Young RA. Insights to host responses against pathogens from transcriptional profiling. Nat Rev Microbiol. 2005;3:281–94.
Unwin RD. Quantification of proteins by iTRAQ. Methods Mol Biol. 2010;658:205–16.
Wu WW, Wang GH, Baek SJ, Shen RF. Comparative study of three proteomic quantitative methods, DIGE, clCAT, and iTRAQ, using 2D gel- or LC-MALDI TOF/TOF. J Proteome Res. 2006;5(3):651–8.
Raquel GF, Elena P, Jesús VJN. Proteomics of plant pathogenic fungi. Journal of biomedicine and biotechnology. 2010;doi:10.1155/2010/932527.
Wu LJ, Wang SX, Wu JY, Han ZP, Wang R, Wu, LC, et al. Phosphoproteomic analysis of the resistant and susceptible genotypes of maize infected with sugarcane mosaic virus. Amino acids. 2014;doi:10.1007/s00726-014-1880-2.
Pasin F, Simón-Mateo C, Garia JA. The hypervariable amino-terminus of P1 protease modulates potyviral replication and host defense responses. PLoS Pathog. 2014;10(3):e1003985. doi:10.1371/journal.ppat.1003985.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Lobo I. Pleiotropy: One gene can affect multiple traits. Nature Education. 2008;1:10.
Kuromori T, Takahashi S, Kondou Y, Shinozaki K, Matsui M. Phenome analysis in plant species using loss-of-function and gain-of-function mutants. Plant Cell Physiol. 2009;50:1215–31.
Serwatowska J, Roque E, Gómez-Mena C, Constantin GD, Wen JQ, Mysore KS, et al. Two euAGAMOUS genes control C-function in Medicago truncatula. PLoS ONE. 2014;doi:10.1371/journal.pone.0103770.
Plader W, Yukawa Y, Sugiura M, Malepszy S. The complete structure of the cucumber (Cucumis sativus L.) chloroplast genome: its composition and comparative analysis. Cell Mol Biol Lett. 2007;12:584–94.
Spetea C. Role of chloroplast thylakoid lumen in photosynthetic regulation and plant cell signaling. Progress in botany. 2012;73:207–30.
Ahuja I, de Vos RC, Bones AM, Hall RD. Plant molecular stress responses face climate change. Trends Plant Sci. 2010;15:664–74.
Shirley BW. Flavonoid biosynthesis: ‘new’ functions for an ‘old’ pathway. Trends Plant Sci. 1996;1:377–82.
Schijlen EGWM, Chr-De V, Aj-Van T, Bovy AG. Modification of flavonoid biosynthesis in crop plants. Phytochemistry. 2004;65:2631–48.
Besseau S, Hoffmanna L, Geoffroya P, Lapierreb C, Polletb B, Legranda M. Flavonoid accumulation in Arabidopsis repressed in lignin synthesis affects auxin transport and plant growth. The Plant Cell January. 2007;19:148–62.
Xu L, Zhu LF, Tu LL, Liu LL, Yuan DJ, Jin L, et al. Lignin metabolism has a central role in the resistance of cotton to the wilt fungus Verticillium dahliae as revealed by RNA-Seq-dependent transcriptional analysis and histochemistry. J Exp Bot. 2011;62:5607–21.
Kao YY, Harding SA, Tsai CJ. Differential expression of two distinct phenylalanine ammonia-lyase genes in condensed tannin-accumulating and lignifying cells of quaking aspen. Plant Physiol. 2002;130:796–807.
Shadle GL, Wesley SV, Korth KL, Chen F, Lamb C, Dixon RA. Phenylpropanoid compounds and disease resistance in transgenic tobacco with altered expression of L-phenylalanine ammonia-lyase. Phytochemistry. 2003;64:153–61.
Gayoso C, Pomar F, Novo-Uzal E, Merino F, Martinez de Ilarduya O. The Ve-mediated resistance response of the tomato to Verticillium dahliae involves H2O2, peroxidase and lignins and drives PAL gene expression. BMC Plant Biol. 2010;10:232–51.
Mao WH, Gong YM, Song XS, Xia XJ, Shi K, Zhou YH, et al. Construction of a cucumber cDNA microarray and its application in the study of response of cucumber plants to magnesium deficiency stress. Acta horticulturae sinica (in Chinese). 2006;33:767–72.
Li JW. Identification of cucumber genes associated with resistance to downy mildew and mechanism for the hypersensitive response. Huazhong agricultural university. Academic thesis. 2010.
Pinosa F, Buhot N, Kwaaitaal M, Fahlberg P, Thordal-Christensen H, Ellerstrom M, et al. Arabidopsis phospholipase d is involved in basal defense and nonhost resistance to powdery mildew fungi. Plant Physiol. 2013;163:896–906.
Egelman EH. Actin’s prokaryotic homologs. Curr Opin Struct Biol. 2003;13:244–8.
Walker SK, Garrill A. Actin microfilaments in fungi. Mycologist. 2006;20:26–31.
Day B, Graham T. The host-pathogen interface cell wall and membrane dynamics of pathogen-induced responses. Ann N Y Acad Sci. 2007;1113:123–34.
Dominguez R, Holmes KC. Actin structure and function. Annu Rev Biophys. 2011;40:169–86.
Day B, Henty JL, Porter KJ, Staiger CJ. The pathogen-actin connection: a platform for defense signaling in plants. Annu Rev Phytopathol. 2011;49:483–506.
van den Ent F, Amos LA, LÖwe J. Prokaryotic origin of the actin cytoskeleton. Nature. 2001;413:39–44.
Pollard TD, Cooper JA. Actin, a center player in cell shape and movement. Science. 2009;326:1208–12.
Millar AJ, McGrath RB, Chua NH. Phytochrome phototransduction pathways. Annu Rev Genet. 1994;28:325–49.
Esposti MD, Vries SD, Crimi M, Ghelli A, Patarnello T, Meyer A. Mitochondrial cytochrome b: evolution and structure of the protein. Biochim Biophys Acta. 1993;1143:243–71.
Brandt U, Trumpower BL. The proton-motive Q cycle in mitochondria and bacteria. Crit Rev Biochem Mol Biol. 1994;29:165–97.
Zannoni D. Cytochrome b type oxidases in the respiratory chains of the phytopathogenic fluorescent bacteria Pseudomonas cichorii and Pseudomonas aptata. Arch microbial. 1982;133:267–73.
Zannoni D. The respiratory chains of pathogenic pseudomonads. Biochim Biophys Acta. 1989;975:299–316.
Semar M, Strobel D, Koch A, Klappach K, Stammler G. Field efficacy of pyraclostrobin against populations of Pyrenophora teres containing the F129L mutation in the cytochrome b gene. Journal of plant diseases and protection. 2007;114:117–9.
Bartlett DW, Clough JM, Godwin JR, Hall AA, Hamer M, Parr-Dobrzanski B. The strobilurin fungicides. Pest Manag Sci. 2002;58:649–62.
Xoconostle-Cázares B, Ruiz-Medrano R, Lucas WJ. Proteolytic processing of CmPP36, a protein from the cytochrome b5 reductase family, is required for entry into the phloem translocation pathway. Plant J. 2000;24(6):735–47.
Simn-Buela L, Garca-Arenal F. Virus particles of cucumber green mottle mosaic tobamovirus move systemically in the phloem of infected cucumber plants. Mol Plant Microbe Interact. 1999;12(2):112–8.
Miles D. The role of high chlorophyll fluorescence photosynthesis mutants in the analysis of chloroplast thylakoid membrane assembly and function. Maydica. 1994;39:35–45.
Hutcheson SW, Buchanan BB. Enzyme regulation in Crassulacean acid metabolism photosynthesis. Plant Physiol. 1979;72:870–6.
Balachandran S, Hurry VM, Kelley SE, Osmond CB, Robinson SA, Rohozinski J, et al. Concepts of plant biotic stress. Some insights into the stress physiology of virus-infected plants, from the perspective of photosynthesis. Physiol Plant. 1997;100:203–13.
Jaber R, Isabel GL, Matilde B. Inhibition of photosynthesis by viral infection: Effect on PSII structure and function. Physiol Plant. 2000;110:286–92.
Shalitin D, Wolf S. Cucumber mosaic virus infection affects sugar transport in melon plants. Plant Physiol. 2000;123:597–604.
Nolte KD, Palumbo JC, Bealmear SR. The relative physiology and phenology of cucurbit yellow stunting disorder virus (CYSDV) in desert grown muskmelon (Cucumis melo, var. reticulatus). Journal of the national association & county agricultural agents (NACAA). 2011;4(2). http://www.nacaa.com/journal/index.php?jid=119
Lu J, Du ZX, Kong J, Chen LN, Qiu YH, Li GF, et al. Transcriptome analysis of Nicotiana tabacum infected by cucumber mosaic virus during systemic symptom development. PLoS One. 2012;7(8):e43447. doi:10.1371/journal.pone.0043447.
Deng W, Liu C, Deng X, Niu L, Cao X. Involvement of the histone acetyltransferase AtHAC1 in the regulation of flowering time via repression of FLOWERING LOCUS C in Arabidopsis. Plant Physiol. 2007;143:1660–8.
Solomon PS, Lee RC, Greer TJW, Oliver RP. Pathogenicity of Stagonospora nodorum requires malate synthase. Mol Microbiol. 2004;53:1065–73.
Wang ZY, Thornton CR, Kershaw MJ, Li DB, Talbot NJ. The glyoxylate cycle is required for temporal regulation of virulence by the plant pathogenic fungus Magnaporthe grisea. Mol Microbiol. 2003;47:1601–12.
Idnurm A, Howlett BJ. Isocitrate lyase is essential for pathogenicity of the fungus Leptosphaeria maculans to canola Brassica napus. Eukaryot Cell. 2002;1:719–24.
Kim MS, Cho SM, Kang EY, Im YJ, Hwangbo H, Kim YC, et al. Galactinol is a signaling component of the induced systemic resistance caused by Pseudomonas chlororaphis O6 root colonization. Mol Plant Microbe Interact. 2008;21:1643–53.
Van VL. A new evolutionary law. Evolutionary Theory. 1973;1:1–30.
De Meaux J, Mitchell-Olds T. Evolution of plant resistance at the molecular level: Ecological context of species interactions. Heredity. 2003;91(4):345–52.
Bohmert K, Camus I, Bellini C, Bouchez D, Caboche M, Benning C. AGO1 defines a novel locus of Arabidopsis controlling leaf development. EMBO J. 1998;17:170–80.
Baniulis D, Yamashita E, Zhang H, Hasan SS, Cramer WA. Structure-function of the cytochrome b6f complex. Photochem Photobiol. 2008;84(6):1349–58.
Berg JM, Tymoczko JL, Biochemistry SL. Biochemistry. New York: W.H. Freeman; 2007. ISBN 978-0-7167-8724-2.
Reinero A, Beachy RN. Association of TMV coat protein with chloroplast membranes in virus-infected leaves. Plant Mol Biol. 1986;6:291–301.
Atkinson NJ, Urwin PE. The interaction of plant biotic and abiotic stresses: from genes to the field. J Exp Bot. 2012;63:3523–43.
Fan YH, Chen J, Zhang CY, Guo H, Lin Y. A two-dimensional electrophoresis protocol suitable for proteomic analysis of cucumber leaves. Journal of Shenyang agricultural university. 2007;39:365–7.
We appreciate the support for this study from the National Science Foundation of China (NSFC) project (Grant No. 31371910), the Program for Changjiang Scholars and Innovative Research Team in University (Grant No. IRT1042) and the Special Fund for Agroscientific Research in the Public Interest of China (Grant No. 201303028).
The authors have declared that no competing interests exist.
HWL, LXL and JQL designed the study. HWL and CQL collected the samples and carried out experiments. HWL and PFL analyzed the data. HWL drafted the manuscript. HWL, LXL and JQL finalized the manuscript. All authors read and approved the final manuscript.