Quantitative Measurement of Phosphoproteome Response to Osmotic Stress in Arabidopsis Based on Library-Assisted eXtracted Ion Chromatogram (LAXIC)

Global phosphorylation changes in plants in response to environmental stress have been relatively poorly characterized to date. Here we introduce a novel mass spectrometry-based label-free quantitation method that facilitates systematic profiling plant phosphoproteome changes with high efficiency and accuracy. This method employs synthetic peptide libraries tailored specifically as internal standards for complex phosphopeptide samples and accordingly, a local normalization algorithm, LAXIC, which calculates phosphopeptide abundance normalized locally with co-eluting library peptides. Normalization was achieved in a small time frame centered to each phosphopeptide to compensate for the diverse ion suppression effect across retention time. The label-free LAXIC method was further treated with a linear regression function to accurately measure phosphoproteome responses to osmotic stress in Arabidopsis. Among 2027 unique phosphopeptides identified and 1850 quantified phosphopeptides in Arabidopsis samples, 468 regulated phosphopeptides representing 497 phosphosites have shown significant changes. Several known and novel components in the abiotic stress pathway were identified, illustrating the capability of this method to identify critical signaling events among dynamic and complex phosphorylation. Further assessment of those regulated proteins may help shed light on phosphorylation response to osmotic stress in plants.Phosphorylation plays a pivotal role in the regulation of a majority of cellular processes via signaling transduction pathways. During the last decade, quantitative phosphoproteomics has become a powerful and versatile platform to profile signaling pathways at a system-wide scale. Multiple signaling networks in different organisms have been characterized through global phosphorylation profiling (1–3), which has evolved over the years with highly optimized procedures for sample preparation and phosphopeptide enrichment, high resolution mass spectrometry, and data analysis algorithms to identify and quantify thousands of phosphorylation events (4–8).Quantitative phosphoproteomics can be achieved mainly by two major techniques, stable isotope labeling and label-free quantitation. Isotope labeling prior to liquid chromatography-mass spectrometry (LC-MS)¹ has been widely used, including metabolic labeling such as stable isotope labeling by amino acids in cell culture (SILAC), chemical labeling such as multiplexed isobaric tags for relative and absolute quantification (iTRAQ) and isotope-coded affinity tags (ICAT) (9–12). On the other hand, label-free quantitation has gained momentum in recent years (13–15), as no additional chemistry or sample preparation steps are required. Compared with stable isotope labeling, label-free quantitation has higher compatibility with the source of the samples, the number of samples for comparison, and the instrument choice.Many label-free approaches, in particular to phosphoproteomics, are based on ion intensity (16, 17), but they are relatively error-prone because of run-to-run variations in LC/MS performance (18). In theory, such systematic errors can be corrected by spiking an internal standard into every sample to be compared. Several methods based on internal standard spiking have been reported so far. Absolute quantification peptide technology (AQUA) (19), for example, uses synthetic peptides with isotope labeling for absolute quantitation. For a global quantitative proteomics study, it is unrealistic to spike-in all reference peptides. Another labeling reference method, spike-in SILAC appears to be a promising technique to quantify the proteome in vivo with multiplex capability and it can be extended to clinical samples (20). One solution to large-scale quantitation without any isotope labeling is pseudo internal standard approach (21), which selects endogenous house-keeping proteins as the internal standard so that no spike-in reagent is required. However, finding a good pseudo internal standard remains a challenge for phosphoproteome samples. Spike-in experiments are an alternative way to improve normalization profile. Some internal standard peptides such as MassPREP^TM (Waters) were already widely used. Other groups spiked an identical amount of standard protein into samples prior to protein digestion (22–24). There are two major normalization procedures. In one approach, sample peptides were normalized to the total peak intensity of spike-in peptides (25). Alternatively, the digested peptides were compared at first and the normalization factor was determined in different ways such as the median (26) or average of ratios (27). However, spiking an identical amount of standard proteins into phosphoproteomic samples before protein digestion may not be compatible with phosphoproteomic analyses which typically require a phosphopeptide enrichment step. Spectral counting has been extensively applied in large sets of proteomic samples because of its simplicity but the method is often not reliable for the quantitation of phosphoproteins, which are typically identified by single phosphopeptides with few spectra (28–30). Many software packages have been implemented to support the wide variety of those quantitation techniques, including commercial platforms such as Progenesis LC-MS^TM, Mascot Distiller^TM, PEAKS Q^TM, etc., as well as open-source software packages including MaxQuant (31), PEPPeR (32), Skyline (33), etc.In this study, we have devised a novel label-free quantitation strategy termed Library Assisted eXtracted Ion Chromatogram (LAXIC) for plant phosphoproteomic analyses with high accuracy and consistency (Fig. 1). The approach employs synthetic peptide libraries as the internal standard. These peptides were prepared to have proper properties for quality control assessments and mass spectrometric measurements. In particular, peptides were designed to elute sequentially over an entire LC gradient and to have suitable ionization efficiency and m/z values within the normally scanned mass range. Local normalization of peak intensity is performed using Loess Procedure, a data treatment adopted from cDNA microarray data analysis (34). To monitor the diverse ion suppression effect across retention time, the local normalization factors (LNFs) are determined by internal standard pairs in individual time windows. Finally, samples will be quantified with LNFs in order to correct variance of LC-MS conditions. This quantification occurs in a small time frame centered to each target peptide.Open in a separate windowFig. 1.Work flow for the LAXIC strategy to quantify the phosphorylation change in response to osmotic stress. A, Schematic representation of the LAXIC algorithm. First, all the chromatographic peaks were aligned and the ratios were calculated. Second, the normalization factors which equal to ratios of library peptide peaks between MS runs were chosen to construct normalization curve. Third, sample peptide peak ratios were normalized against predicted normalization factor corresponding to certain retention time. B, Schematic representation of quantitative phosphoproteomics. Plants either treated with mannitol or PBS were lysed and mixed proportionally at first. Following peptide digestion and enrichment, phosphopeptides were identified and further quantified through LAXIC incorporated with self-validating process using thelinear regression model to analyze the fold change (fold), linearity (R²) and accuracy (%Acc).Water deficit and salinity causes osmotic stress, which is a major environmental factor limiting plant agricultural productivity. Osmotic stress rapidly changes the metabolism, gene expression and development of plant cells by activating several signaling pathways. Several protein kinases have been characterized as key components in osmotic stress signaling. Arabidopsis histidine kinase AHK1 can complement the histidine kinase mutant yeast, which can act as the osmosensor in yeast (35). Overexpression of AHK1 gene increases the drought tolerance of transgenic plants in Arabidopsis (36). Similar to yeast, the MAPK kinase cascade is also involved in osmotic stress response in plants. It is reported that AtMPK3, AtMPK6, and tobacco SIPK can be activated by NaCl or mannitol, and play positive roles in osmotic signaling (37, 38). MKK7 and MKKK20 may act as the up-stream kinase in the kinase cascade (39). Involvement of some calcium-dependent protein kinases, such as AtCPK21, AtCPK6, and OsCPK7 (CDPK) in osmotic stress signaling has also been reported (40–42). Another kinase family, SNF1-related protein kinase (SnRK) 2, also participates in osmotic stress response. In Arabidopsis, there are ten members in the SnRK2 family. Five from the ten SnRK2s, SnRK2, 3, 6, 7, and 8, can be activated by abscisic acid (ABA) and play central roles in ABA-receptor coupled signaling (43, 44). Furthermore, all SnRK2s except SnRK2.9 can be activated by NaCl or mannitol treatment (43). The decuple mutant of SnRK2 showed a strong osmotic hypersensitive phenotype (45). It is proposed that protein kinases including MAPK and SnRK2s have a critical function in osmotic stress (46), but the detailed mechanism and downstream substrates or target signal components are not yet clarified. We applied, therefore, the LAXIC approach with a self-validating method (47) to profile the osmotic stress-dependent phosphoproteome in Arabidopsis by quantifying phosphorylation events before and after mannitol treatment. Among a total of over 2000 phosphopeptides, more than 400 of them are dependent on osmotic stress. Those phosphoproteins are present on enzymes participating in signaling networks that are involved in many processes such as signal transduction, cytoskeleton development, and apoptosis. Overall, LAXIC represents a powerful tool for label-free quantitative phosphoproteomics.     

Development of rapid and highly sensitive HSPA1A promoter-driven luciferase reporter system for assessing oxidative stress associated with low-dose photodynamic therapy

Photodynamic therapy (PDT) is a regulatory-approved modality for treating a variety of malignant tumors. It induces tumor tissue damage via photosensitizer-mediated oxidative cytotoxicity. The heat shock protein 70 (HSP70-1) is a stress protein encoded by the HSPA1A gene and is significantly induced by oxidative stress associated with PDT. The aim of this study was to identify the functional region of the HSPA1A promoter that responds to PDT-induced oxidative stress and uses the stress responsiveness of HSPA1A expression to establish a rapid and cost-effective photocytotoxic assessment bioassay to evaluate the photodynamic potential of photosensitizers. By constructing luciferase vectors with a variety of hspa1a promoter fractions and examining their relative luciferase activity, we demonstrated that the DNA sequence from −218 to +87 of the HSPA1A gene could be used as a functional promoter to detect the PDT-induced oxidative stress. The maximal relative luciferase activity level of HSPA1A (HSP70-1) induced by hypericin-PDT was nearly nine times that of the control. Our results suggest that the novel reporter gene assay using a functional region of the HSP70A1A promoter has significant advantages for the detection of photoactivity in terms of both speed and sensitivity, when compared with a cell viability test based on ATP quantification and ROS levels. Furthermore, phthalocyanine zinc and methylene blue both induced significantly elevated levels of relative luciferase activity in a dose-dependent manner.     

Inhibition of cholesteryl ester transfer protein by anacetrapib does not impair the anti-inflammatory properties of high density lipoprotein

Cholesteryl ester transfer protein (CETP) is a target of therapeutic intervention for coronary heart disease. Anacetrapib, a potent inhibitor of CETP, has been shown to reduce LDL-cholesterol by 40% and increase HDL-cholesterol by 140% in patients, and is currently being evaluated in a phase III cardiovascular outcomes trial. HDL is known to possess anti-inflammatory properties, however with such large increases in HDL-cholesterol, it is unclear whether CETP inhibition perturbs HDL functionality such as anti-inflammatory effects on endothelial cells. The purpose of the present study was to determine whether CETP inhibition by anacetrapib affects the anti-inflammatory properties of HDL. HDL was isolated from either hamsters treated with vehicle or anacetrapib for 2 weeks, or from normal human subjects treated either placebo, 20 mg, or 150 mg anacetrapib daily for 2 weeks. Anacetrapib treatment increased plasma HDL cholesterol levels by 65% and between 48 and 82% in hamsters and humans, respectively. Pre-incubation of human aortic endothelial cells with HDL isolated from both control and anacetrapib treated hamsters suppressed TNFα induced expression of vascular cell adhesion molecule-1 (VCAM-1), intercellular adhesion molecule-1 (ICAM-1) and E-selectin. Similar results were obtained with human HDL samples pre and post treatment with placebo or anacetrapib. Further, HDL inhibited TNFα-induced MCP-1 secretion, monocyte adhesion and NF-κB activation in endothelial cells, and the inhibition was similar between control and anacetrapib treated groups. These studies demonstrate that anacetrapib treatment does not impair the ability of HDL to suppress an inflammatory response in endothelial cells.     

Development of InDel markers linked to Fusarium wilt resistance in cabbage

Cabbage Fusarium wilt (CFW) is a destructive disease causing great losses to cabbage (Brassica oleracea L. var. capitata L.) production worldwide. At present, there are few reports concerning molecular marker research on cabbage resistance to CFW. In this study, 160 double haploid (DH) lines were obtained from the F1 population of a 99–77 (highly resistant to CFW) × 99–91 (highly susceptible to CFW) cross. Insertion–deletion (InDel) markers were designed according to the reference genome sequence of cabbage and the whole-genome re-sequencing data of the two parents. A genetic map of chromosome C06 including seven InDel markers was constructed based on the DH population. Thus, FOC (resistance gene to Fusarium oxysporum f. sp. conglutinans) was located on chromosome C06 and two InDel markers out of the seven, M10 and A1, flanked the gene at 1.2 and 0.6 cM, respectively. Marker A1 revealed a significant consistency with the phenotype assay in the F2 population as well as in 40 inbred lines (96 and 82 %, respectively). This study lays the foundation for fine mapping and cloning of the FOC gene and for marker-assisted selection in cabbage resistance breeding.     

Physiological mechanisms for plant distribution pattern: responses to flooding and drought in three wetland plants from Dongting Lake, China

Both flooding and drought are important in determining plant distribution in wetlands. However, the roles of plant’s physiological response to flooding and drought in accounting for plant distribution are far from clear. To this end, three typical wetland plants with different distribution patterns (high-elevation species Miscanthus sacchariflorus, low-elevation species Carex brevicuspis and Polygonum hydropiper) in Dongting Lake were treated with three water levels (flooding 25 cm, control 0 cm, drought ?25 cm), and relative growth rate (RGR), malondialdehyde (MDA) content, electrolyte leakage and proline content were investigated. The RGR of the three species decreased significantly in both flooding and drought treatments. Compared to the control, the RGR of M. sacchariflorus decreased more in the flooding treatment but less in the drought treatment compared to the other two species. The contents of MDA in the three species increased in both flooding and drought treatments, except for P. hydropiper in the flooding treatment. MDA contents increased more in M. sacchariflorus in the flooding treatment but less in the drought treatment compared to the other two species. Only M. sacchariflorus had a higher electrolyte leakage in the flooding treatment, and drought led to a higher electrolyte leakage in P. hydropiper and C. brevicuspis. Proline content increased 69.2, 66.7 and 39.6 % in P. hydropiper, C. brevicuspis and M. sacchariflorus in the flooding treatment, and increased 44.2, 13.0 and 45.3 % in the drought treatment, respectively. These results suggest that M. sacchariflorus has a higher tolerance to drought but a lower tolerance to flooding than do the other two species, which might be the intrinsic mechanisms accounting for their different distribution patterns.     

Caenorhabditis elegans as a model for intracellular pathogen infection

The genetically tractable nematode Caenorhabditis elegans is a convenient host for studies of pathogen infection. With the recent identification of two types of natural intracellular pathogens of C. elegans, this host now provides the opportunity to examine interactions and defence against intracellular pathogens in a whole‐animal model for infection. C. elegans is the natural host for a genus of microsporidia, which comprise a phylum of fungal‐related pathogens of widespread importance for agriculture and medicine. More recently, C. elegans has been shown to be a natural host for viruses related to the Nodaviridae family. Both microsporidian and viral pathogens infect the C. elegans intestine, which is composed of cells that share striking similarities to human intestinal epithelial cells. Because C. elegans nematodes are transparent, these infections provide a unique opportunity to visualize differentiated intestinal cells in vivo during the course of intracellular infection. Together, these two natural pathogens of C. elegans provide powerful systems in which to study microbial pathogenesis and host responses to intracellular infection.     

Characterization of water and wildlife strains as a subgroup of Campylobacter jejuni using DNA microarrays

Campylobacter jejuni is the leading cause of human bacterial gastroenteritis worldwide, but source attribution of the organism is difficult. Previously, DNA microarrays were used to investigate isolate source, which suggested a non‐livestock source of infection. In this study we analysed the genome content of 162 clinical, livestock and water and wildlife (WW) associated isolates combined with the previous study. Isolates were grouped by genotypes into nine clusters (C1 to C9). Multilocus sequence typing (MLST) data demonstrated that livestock associated clonal complexes dominated clusters C1–C6. The majority of WW isolates were present in the C9 cluster. Analysis of previously reported genomic variable regions demonstrated that these regions were linked to specific clusters. Two novel variable regions were identified. A six gene multiplex PCR (mPCR) assay, designed to effectively differentiated strains into clusters, was validated with 30 isolates. A further five WW isolates were tested by mPCR and were assigned to the C7‐C9 group of clusters. The predictive mPCR test could be used to indicate if a clinical case has come from domesticated or WW sources. Our findings provide further evidence that WW C. jejuni subtypes show niche adaptation and may be important in causing human infection.     

Roles of Eukaryotic Initiation Factor 5A2 in Human Cancer

Eukaryotic initiation factor 5A (eIF5A), the only known cellular protein containing the amino acid hypusine, is an essential component of translation elongation. eIF5A2, one of the two isoforms in the eIF5A family, is reported to be a novel oncogenic protein in many types of human cancer. Both in vitro and in vivo studies showed that eIF5A2 could initiate tumor formation, enhance cancer cell growth, and increase cancer cell motility and metastasis by inducing epithelial-mesenchymal transition. Accumulatied evidence suggests that eIF5A2 is a useful biomarker in the prediction of cancer prognoses and serves as an anticancer molecular target. In this review, we will focus on updating current knowledge of the EIF5A2 gene in human cancers. The molecular mechanisms of EIF5A2 related to tumorigenesis will also be discussed.