- Open Access
The quantitative proteomic response of Synechocystis sp. PCC6803 to phosphate acclimation
© Fuszard et al.; licensee BioMed Central Ltd. 2013
- Received: 8 June 2012
- Accepted: 28 January 2013
- Published: 26 February 2013
Inorganic phosphate (Pi) is a critical nutrient for all life and is periodically limiting in marine and freshwater provinces, yet little is understood how organisms acclimate to fluctuations in Pi within their environment. To investigate whole cell adaptation, we grew Synechocystis sp. PCC6803, a model freshwater cyanobacterium, in 3%, and 0.3% inorganic phosphate (Pi) media. The cells were allowed to acclimate over 60 days, and cells were harvested for quantitative high throughput mass spectrometry-based proteomics using the iTRAQ™ labelling technology.
In total, 120 proteins were identified, and 52 proteins were considered differentially abundant compared to the control. Alkaline phosphatase (APase) activities correlated significantly (p < 0.05) with observed relative PhoA abundances. PstS1 and PstS2 were both observed, yet PstS1 was not differentially more abundant than the control. Phycobilisome protein abundances appeared to be coordinated, and are significantly less abundant in 0.3% Pi than 3% Pi cultures. Also, the central metabolic cell function appears to have shifted towards the production of (NADPH) reducing energy and nucleotide sugars.
This acclimation response bears strong similarity to the previously reported response to nitrogen deprivation within Synechocystis sp. PCC 6803. However, it also demonstrates some characteristics of desiccation stress, such as the regulation of fatty acids and increased abundance of rehydrin in the 3% Pi culture.
- Phosphate acclimation
- Pentose phosphate pathway
- Ribose sugars
Cyanobacteria comprise one of the largest groups of prokaryotes. They are found in a wide variety of niches, including terrestrial and hypersaline environments, and in oceans, lakes and rivers. Synechocystis sp. PCC 6803 (hereafter Synechocystis), a model freshwater cyanobacterium, has strong potential biotechnology applications in areas as diverse as biofuels , biopolymers  and other secondary metabolites . Due to this biotechnological interest, its general applicability as a model organism for CO2 fixation and environmental stress responses has been the subject of post-genomic investigations [4–8], including proteomics investigations using both 2-DE  and shotgun proteomics workflows .
For most microorganisms, the preferred source of the essential nutrient phosphorus is inorganic phosphate (Pi). However, phosphate is the planet’s least abundant naturally occurring mineral . Pi limitation occurs for sustained periods in many geographically dispersed freshwater lakes and marine sites. To date, however, there have been relatively few investigations into the cell-wide response of microorganisms to this form of nutrient stress, particularly over longer periods that reflect adaptation rather than short-term response. An understanding of these processes will be of use not only in the basic understanding of microbes and the niches they occupy, but could also be applied to develop biomarkers to follow stages in the phosphorus biogeochemical cycle within aquatic environments. In addition, the identification of novel proteins associated with Pi uptake could be exploited in the removal of Pi for water quality improvement, whilst enzymes involved in the catabolism of organophosphorus molecules are of interest in the pharmaceutical and agribiotech industrial sectors.
There has been a range of phosphate starvation and limitation studies in bacteria carried out through cDNA microarrays to investigate global transcriptome changes [10–13] and qualitative overview analysis using 2-DE workflows to study broad proteome changes [14–17]. To date, there have only been two high-throughput quantitative proteomic assays of Pi limitation in cyanobacteria, focusing on the marine species Prochlorococcus[18, 19]. While it is tempting to compare the observations obtained from these studies, comparison among datasets has also been further complicated by the different types of experiments conducted, which range from evaluation of the short-term responses to starvation conditions (shift to P-free medium) to studies of adaptation to Pi limitation. However, one of the most common observations in all of these studies is the up-regulation of pho regulon genes, as seen in studies conducted on Bacillus subtilis, Escherichia coli, Prochlorococcus spp. [13, 18, 19], Synechococcus WH8102  and Corynebacterium glutamicum.
Unlike other prokaryotic systems, reports on Pi limitation responses in cyanobacteria are relatively few. When available, these studies consist more generally of short-term Pi shock experiments with a shift from replete to Pi free medium [10, 13, 20, 22–26]. However, these studies have revealed that the phosphate-sensing system in Synechocystis is modulated by SphS (sll0337; histidine kinase; analogous to PhoR of E. coli) and SphR (slr0081; response regulator; analogous to PhoB of E. coli) , which form a two-component signal transduction system. Also, an in-depth investigation of Pi-binding proteins within Synechocystis showed differential regulation, functional roles and Pi transport characteristics for its two ABC-type phosphate transporter proteins, PstS1 and PstS2 .
Inactivation of phoU, one of the genes encoding phosphate-specific transport system regulatory proteins, caused an increase in both inorganic phosphate (Pi) uptake and polyphosphate accumulation rates . It has been shown that under phosphate starvation, the Synechocystis pho regulatory system responds by increasing its transcript abundance by more than 7-fold compared to normal Pi replete conditions . It has been suggested that high light levels also triggered the pho regulon, since light accelerated phosphate assimilation despite limited cell capacity to perform the task . Suzuki et al. used DNA microarrays to investigate the Synechocystis phosphate starvation response  and found that the genes encoding the ABC-type phosphate transporter (sphX, PstS1 and PstS2 clusters) and extracellular nuclease (nucH) were significantly up-regulated in addition to the alkaline phosphatase (phoA) gene. The gene encoding a periplasmic binding protein (urtA) was down-regulated. In a separate study under similar conditions, significant up-regulation of the genes for inorganic pyrophosphatase (ppa) and exopolyphosphatase (ppx) was observed .
The goal of this work was to study the adaptation of Synechocystis to phosphate limitation using a quantitative shotgun proteomics workflow [30, 31]. The proteomes of Synechocystis grown in different phosphate concentrations (3% and 0.3% of the non-limiting level in BG-11 medium) were quantitatively measured using iTRAQ reagents. Coupled with enzyme assays, this investigation forms the basis for understanding global cellular responses, such as the central metabolism and biosynthetic pathways, of Synechocystis in phosphate-limited environments. This study is both novel and timely as it specifically deals with an in-depth analysis of Pi-depletion over different levels of phosphate, and over a longer timescale, than has been explored by others. The results of this study provide insights into the ways in which cellular processes are directly influenced by ambient phosphorus concentrations.
In phosphate-depleted conditions, the growth rate of Synechocystis was significantly impaired (Additional file 1: Figure S1), with noticeable reduction of colour within the flasks (data from observation, image not available). The control culture reached an OD730 of 7 at the end of the 60-day cultivation, whereas the phosphate-limited cultures only reached a maximum OD730 of 0.9 (for 3% Pi level) and 0.5 (for 0.3% Pi level). Their specific growth rates were calculated as 0.20 day-1 for the control culture, 0.012 day-1 for 3% Pi culture and 0.002 day-1 for 0.3% Pi culture. This impairment in growth with sub-optimal Pi levels was expected, and agrees with the trend observed for other organisms [21, 32, 33].
Differentially abundant proteins related to control cultures
Nutrient assimilation and regulation
Periplasmic Pi-binding protein
Periplasmic Pi-binding protein
Light harvesting and photosystems
C-phycocyanin alpha chain
C-phycocyanin beta chain
Phycobilisome 32.1 kDa linker polypeptide, phycocyanin-associated, rod 1
Phycobilisome 32.1 kDa linker polypeptide, phycocyanin-associated, rod 2
Phycobilisome rod-core linker polypeptide
Allophycocyanin alpha chain
Allophycocyanin beta chain
Phycobilisome LCM core-membrane linker polypeptide
Ferredoxin NADP reductase
Thylakoid-associated protein slr0729
Uncharacterized thylakoid-associated protein ssl2501
6-phosphogluconate dehydrogenase, decarboxylating
RuBisCO large subunit
RuBisCO small subunit
3-oxoacyl-[acyl-carrier-protein] reductase 2
Nucleoside diphosphate kinase
Transcription, translation and stress
ATP-dependent Clp protease regulatory subunit
groEL protein 1
Elongation factor Tu
DNA-binding protein HU
Superoxide dismutase [Fe]
It is instructive to observe that cells acclimated to a lower initial phosphate concentration demonstrate a significantly higher fraction of differentially more abundant proteins from almost all categories (Figure 1), with the exception of proteins associated with the phycobilisomes. These differences suggest a substantial physiological shift within cells in order to maximise Pi usage throughout. Indeed, when observing relative abundance levels, the functional categories most broadly affected are metabolism, nutrient acquisition and the proteins of the phycobilisomes (Figure 1, Table 1).
Nutrient assimilation and regulation
Alkaline phosphatase activity within samples
Phosphate concentration (mg/L)
APA (nmoles of Pi released per min per mg of protein)
PhoA relative abundance
PhoA standard deviation
Ten proteins involved with nutrient acquisition and regulation (Figure 2A) were quantified. Of these, six of the eight proteins associated with phosphate assimilation and storage were significantly more abundant than in the control (Figure 2A, Table 1). The alkaline phosphatase enzyme, PhoA, and the periplasmic phosphate-binding protein, PstS2, were both significantly more abundant in both 3% and 0.3% Pi cultures, whereas inorganic pyrophosphatase (ppa) was significantly more abundant in 3% Pi culture. We have further confidence with PstS2 abundances, as both data points were significant using the Bonferroni correction. As would be expected, there appears to be no significant observations regarding proteins involved with carbon sources (DctP), iron (bfr) and nitrogen metabolism (GlnA, GlnB, and UreA).
Considering both the relative protein abundance of PhoA in the 100%, 3% and 0.3% Pi cultures, and observed APase activity (Figure 2B, Table 2), there is a direct correlation (p < 0.05, where polynomial regressions fitted both datasets with R2 values > 0.99, and APase activity is approximately 22.5 x PhoA abundance). To the authors’ knowledge, this is the first observed direct correlation between protein abundances and APase activity within Synechocystis.
An interesting observation in relative levels of protein abundances are the significant differences between abundances of proteins in the 3% Pi and 0.3% Pi cultures. PhoA, and PstS12 are both significantly more abundant in the 0.3% Pi than in 3% Pi (Figure 2A), whereas there is no significant difference in abundance of Ppa proteins between the two cultures. This suggests the possibility of a two-tier strategy in cell acclimation: (i) a phased response, whereby protein abundance is putatively associated with/ proportional to internal or external phosphate concentrations, and (ii) a switched response, whereby the cell has a possible biochemical ‘marker’ for activation, for example an as yet undefined external phosphate concentration, and which below this abundance levels of that protein are constitutively maintained within the cell. However, more research is required to investigate this.
It has been previously demonstrated that the two periplasmic phosphate-binding proteins PstS12 have unique enzymatic and regulatory properties . pstS1 encodes a low-affinity, high velocity Pi transporter protein that was expressed to a significantly greater extent than the pstS2 transcript when cells were grown in Pi-replete conditions. However, pstS2, a high-affinity, low-velocity transporter, was expressed at significantly higher levels under Pi-stress conditions . In addition, the disruption of pstS1 caused constitutive expression of the pho regulon, indicating a role in external Pi sensing . The results here, in part, support those findings. However, our results also infer that elevated expression of pstS2 is not limited to a stress response, but is also an acclimation response. PstS2 is significantly higher than PstS1 within 3% and 0.3% Pi cultures, and the abundance of PstS2 within the 0.3% Pi culture is significantly higher than in the 3% Pi culture (Figure 2A).
Light harvesting and photosystems
It is interesting to note the significantly linear relationship within the phycocyanin group and allophycocyanin group of proteins. A statistical analysis of the phycobiliprotein abundances between the 3% and 0.3% Pi cultures revealed that there is no significant difference between the phycocyanins (CpcA-G, Figure 3A) and the allophycocyanins (ApcA-E, Figure 3A) at the 3% Pi level (p = 0.09, df = 32). However, there is a significant difference at the 0.3% Pi level (p = 0.006, df = 32). This difference indicates that the light harvesting antennae of 0.3% Pi culture cultures have significantly less of the phycocyanins in proportion to their allophycocyanin content than do the 3% Pi culture cultures. This potentially suggests independent regulatory systems between the different PBP groups, or differential regulation of group specific proteases, whose activity is putatively dependent upon external Pi concentrations. This is a reasonable hypothesis, as phycocyanins are attached externally to allophycocyanins through linker peptides, and light energy is funnelled from phycocyanins through allophycocyanins to PSII . This regulation may be due to the organism attempting to reduce photosynthetically induced oxidative damage.
Four proteins within photosystem II (PSII) and the photosynthetic electron transport chain (PETC) are differentially less abundant (Figure 3B). Both PsbU and PsbW appear to be less abundant than in the control culture under both conditions, whereas PsbV is significantly less abundant only in the 0.3% Pi culture. The extrinsic protein PsbV, alongside PsbU, is a critical protein in the oxygen evolving complex of cyanobacterial PSII complexes [45, 46]. Indeed, PsbU appears to promote stability within the water-splitting PSII structure . PsbW (characterized recently, and more commonly referred to as Psb28 within Synechocystis) appears to be closely associated with PSII complexes, and more specifically to the biogenesis of the chlorophyll-binding CP47 protein . As a consequence, it appears that when cells acclimate to low levels of ambient Pi, light harvesting and structural integrity of the PSII assembly is compromised.
Ferridoxin NAPDH reductase (PetH, Figure 3B) is commonly associated with PSI and the terminal end of the PETC, nominally with regards to the production of reducing energy. However, in Synechocystis it has also been demonstrated to be closely associated with the phycobilisomes . Within our study, the relative abundance of PetH is significantly less within the 0.3% Pi culture, but not the 3% Pi. This then, coupled with the PsbV observation, implies a strategic down regulation of PETC and PSII activity as ambient Pi concentrations decrease. However, no other proteins associated with PSI were identified, and the suggestion that PSI function was decreased requires further validation in the future.
It is clear that ambient Pi concentrations have a significant effect upon the photosynthetic machinery within Synechocystis. However, it is also clear that this did not mean a complete loss of function, as occurs during prolonged nitrogen deprivation , and represents a proportional loss of function, dependent upon Pi concentrations. Our results present evidence of a response comparable to Synechococcus sp. PCC 7942 .
In contrast, three glycolysis and carbon fixation proteins, phosphoglycerate kinase (Pgk) and the two subunits of rubisco (RbcL and RbcS) were significantly less abundant in 3% Pi and 0.3% Pi vs the 100% Pi controls (Figure 4A, Table 1). These proteins are the central enzymes within the first two phases of the Calvin-Benson cycle (being the fixation of CO2, the reduction of 3-phosphoglycerate, and the regeneration of ribulose 1,5-bisphosphate). This cycle is critical in all photosynthetic organisms, and is fundamentally limiting to growth and reproduction. This result is unsurprising, given the significant reduction in photosynthetic light harvesting earlier observed in both 3% and 0.3% Pi cultures. It is also interesting to note that the enzyme Transketolase (TktA), which functions in both the Calvin-Benson cycle and PPP was not differentially abundant.
Interestingly, at the transcript level, gnd, tal and zwf were all strongly induced under nitrogen deprivation, as well as one of the two glycogen phosphorylases glpP (sll1367). The Calvin-Benson cycle genes rbcL and rbcS were also down regulated . Our results demonstrate that the other glycogen phosphorylase (sll1356), rather than sll1367 is more abundant (see below), however the similarity in results is noteworthy. Osani and co-workers  concluded that nitrogen deprivation induces sugar catabolism, and that the PPP is the primary method of reducing power production. Our results confirm that Pi deprivation also leads to a similar biochemical response, and thus indicate that the central metabolic response of the cell is similar for both N and Pi deprivation. However this response has a significantly different level of severity depending upon N or Pi deprivation .
When considering energy metabolism (defined in this instance as the synthesis and degradation of the glycolytic/PPP precursor glucose 6-phosphate), Figure 4B shows a significant difference between the relative abundance of two glycogen phosphorylation proteins (GlgP(1) & GlgP(2)). Both proteins are assumed to function to degrade glycogen, thereby generating D-glucose-1-phosphate in the preparatory, rate-limiting step before entry into glycolysis. However, the differential relative abundance between both proteins suggests differential functional roles within the cell. The enhancement of central PPP protein abundances highlighted earlier suggests that glycogen is being processed by GlgP(2) to support the production of NADPH and pentose sugars. As mentioned earlier, GlgP(1) rather than GlgP(2) was found to be differentially upregulated at the transcript level when subjected to N deprivation , which then suggests differing regulatory controls for each enzyme and condition. Also, Wegener and associates observed that sll1367 was more abundant than sll1356 under P deplete conditions , and so we must conclude that over short P starvation periods, sll1367 is preferred, and over longer, chronic, periods sll1356 predominates.
Transcription, translation, stress and hypothetical proteins
Six proteins out of 24 relating to transcription, translation and stress were significantly differentially abundant within both reduced Pi cultures. An interesting group of proteins are the ATP-dependant Clp proteases, ClpC, ClpP1, ClpP2 and ClpB. These proteins are vital for growth and stress acclimation, and act to degrade denatured polypeptides and maintain cellular homeostasis . Within Escherichia coli there is evidence that these proteins act in a multi-subunit complex, with the exclusion of ClpB (which acts independently) . However, little is known about cyanobacterial Clp proteases, except that cyanobacterial ClpC is 90% similar to the ClpC in plant chloroplasts, is transcriptionally constitutively expressed at a level which appears to be unaffected by stress conditions, and that loss-of-function manipulation is generally fatal . This last point emphasises the protein’s role as an essential chaperone.
Translation appears to be affected as the elongation factor, EF-Tu (Tuf, Figure 5), which facilitates the entry of aminoacyl tRNA into the free site of the ribosome, was significantly lower in abundance within the 0.3% culture. EF-Tu is critical in successful translation, and acts as a powerful indicator that translation is lower within cells acclimated to long term Pi deprivation. It is instructive to note that the other three elongation factors identified in our study (Gst, Tsf, and FusA) were unaffected, which may indicate a regulatory role for EF-Tu.
Two known stress response proteins were more abundant in low Pi: rehydrin and SodB (Figure 5). The functional role of rehydrin, a protein associated with desiccation stress in mosses , is unknown, and this represents the first recording of this protein being intrinsic to Pi acclimation. In contrast, the function of the iron superoxide dismutase protein, SodB, is well characterised , and is primarily responsible for the cell’s response to oxidative stress.
One biosynthetic protein was significantly differentially abundant within the Pi-deplete cultures than the control, 3-oxoacyl-[acyl-carrier protein] reductase (FabG2, Figure 4B, Table 1), which participates in fatty acid biosynthesis. This is interesting as fatty acids play two main functions within the cell: (i) they are precursors of membrane lipids, and (ii) they accumulate for energy storage. However given their differential abundance, it would be expected that another critical protein required in this process, the acyl-carrier protein (acp) demonstrates no differential abundance in either the 3% Pi, or the 0.3% Pi culture. The reason for this is unclear, however the abundance of FabG2 is a strong indication that the cell over produced fatty acids when acclimating to low Pi environments. Indeed, fatty acids have been speculated to have a role in stress tolerance within cyanobacteria, notably to desiccation, salt and temperature stress (see Singh et al.  for review). As well as energy storage, the production of lipids may help to stabilize nutrient acquisition machinery within the external membrane , and replace lipids damaged by reactive oxygen species, which have concurrent effects upon membrane fluidity and other cellular metabolic functions . The role of fatty acids in desiccation stress is intriguing, specifically due to the increased abundance of the protein rehydrin within the 3% culture (Figure 5, Table 1).
Twenty-one hypothetical proteins were identified within this study (Additional file 2: Table S1). Although at the time of writing, the functions of particular hypothetical proteins are unknown, it is an important observation that three proteins are differentially abundant within the cultures within this study. They are slr0645, slr1619 and slr2025 (Figure 5, Table 1). Of the three, slr2025 was recently observed within a proteome-wide survey of Synechocystis under CO2 limitation. However, it was not considered differentially abundant , and so appears to be unimportant to acclimation to that environment. In our study, it was more abundant within the 3% Pi culture.
Since Synechocystis typically inhabit rivers and lakes where the ambient Pi concentrations range from nanomolar to micromolar ranges , our control cultures were acclimated to “luxury” Pi concentrations, and hence the 3% and the 0.3% Pi cultures, once re-acclimated, are more indicative of natural conditions. This study then emphasises the physiological mechanisms important for normal cellular reproduction and growth. The observed proteome acclimation response to environmental phosphate concentrations involved shifts in cellular abundances of Pi acquisition and storage mechanisms, coinciding with changes in its photosynthetic machinery and a shift in central metabolism towards the production of reducing energy. While the coordination and control of this response was not uncovered, we and others have previously speculated that it is likely to be a complex interaction of internal regulatory systems initiated through the phoBR homologue of Synechocystis – the SphS-SphR two component system , and the activity of PstS1 .
We have shown that the observed abundance response of alkaline phosphatase is proportional to measured APase activity, and interestingly that there appears to be a coordinated and specific reduction of individual phycobiliprotein groups during acclimation. As has been previously reported for higher plants, these results provide some evidence that intracellular Pi concentrations control photosynthetic activity in Synechocystis through the activation or deactivation of RuBisCO at different concentrations . The similarity of the Synechocystis response to Pi deprivation appears to be remarkably similar to prolonged nitrogen deprivation, although the magnitude of the response is greater in low nitrogen than in low Pi. In both cases, the photosynthetic machinery was reduced significantly, and central metabolic processes shift to produce not just reducing energy but nucleotide sugars. We hypothesise that this cellular shift is designed to function in two ways: 1) that an increase in reducing energy provides the ‘power’ for the cell to stabilize and maintain the cell membranes, and 2) nucleotide sugars can be readily used as precursors for both ATP and GTP, as well as enabling the cell to quickly respond to changes in the environment through rapid production of RNA. The cell is therefore able to respond, grow and replicate rapidly.
Synechocystis was principally grown in BG-11 medium  in 1 L conical flasks with a 0.5-L working volume shaken at 100 rpm. The Pi replete BG-11 medium contained 40 mg/L of K2HPO4. For phosphate depletion experiments, the concentration of K2HPO4 within BG-11 was prepared as 3% w/v (1.2 mg/L) and 0.3% w/v of the control Pi replete BG-11 medium. The experimental design utilised biological duplicates of Synechocystis PCC 6803 cultured in each condition: (A) 100% Pi (control), (B) 3% Pi, and (C) 0.3% Pi. Inocula for each culture were prepared by initial centrifugation of stock cultures followed by washing cells twice with sterilized deionised water and once with the respective Pi reduced medium, prior to resuspension into the respective medium and experimental inception. Cultures were incubated in a Sanyo MLR-350H environmental chamber at 25°C on a 12-h light–dark cycle for 60 days with an initial optical density of 0.2. Irradiance was provided by cool-white fluorescent lamps at an incident intensity of 70 μE m-2 s-1 as measured by a three-dimensional QSL-2100 Light Sensor (Biospherical Instrument Inc, San Diego, CA). The growth rate was monitored and recorded using an Ultraspec 2100-Pro spectrophotometer (Biochrom, Cambridge, UK) at an OD of 730 nm.
Protein extraction and quantitation
Cells were harvested at day 60 by centrifugation at 5,000 × g for 5 min at 4°C. The cell pellets were re-suspended in 500 mM triethylammonium bicarbonate buffer (TEAB) at pH 8.5 and proteins were extracted through mechanical cracking with liquid nitrogen . Soluble proteins were recovered from the supernatant by centrifugation at 21,000 × g for 30 min at 4°C. The total soluble protein concentrate was then measured using the RC DC Protein Quantification Assay (Bio-Rad, Hertfordshire, UK) according to the manufacturer’s protocol.
Acid/alkaline phosphatase assays
A subsample of cells broken by mechanical grinding with liquid nitrogen as described above was subsequently re-suspended in Tris–HCl buffer (pH 8.0) containing 5 mM MgCl2. Alkaline phosphatase assays were performed in duplicate as described by O’Loughlin et al.  at 30°C. The assay solution contained 50 mM Tris–HCl (pH 8.0), 20 mM MgCl2, and 5 mM p-nitrophenyl phosphate with 0.01-0.1 mg cell-extract protein in a final reaction volume of 1 mL. For acid phosphatase assays, the Tris–HCl buffer was replaced with 50 mM succinate-NaOH buffer (pH 5.0).
Due to the presence of large quantities of chlorophylls in the cell-extract, detection of p-nitrophenol production by means of the widely used absorbance measurement at 412 nm was not possible. Reactions were therefore quenched with 0.2 vol of 3 M trichloroacetic acid solution following addition of cell extract protein and a duplicate reaction tube quenched 10 min later. The resultant reaction supernatants were cooled on ice for 10 min prior to centrifugation at 20,000 × g for 5 min to remove precipitated protein, following which the phosphate concentration in the samples was determined using the acid molybdate method of Fiske and SubbaRow . ‘No substrate’ and ‘no cell-extract’ controls were included and activity was expressed as nmol of Pi released per min per mg of cell-extract protein, taking the control reactions into account.
SCX prefractionation and iTRAQ
Approximately 100 μg of protein from each sample was placed in PCR clean LoBind microcentrifuge tubes (Eppendorf, Cambridgeshire, UK), then reduced, alkylated and digested as previously reported . The resultant peptide digests were then labelled using 6 labels from an 8-plex iTRAQ chemical reagents kit (AB Sciex UK Limited, UK), with phosphate replete control samples (nominated as 100% Pi) tagged with the 113 and 114 labels, cells grown in 3% Pi tagged with 115 and 116, and the 0.3% phosphate-grown cells tagged with the 117 and 118 labels. Primary peptide fractionation was carried out through strong cation exchange (SCX) on a BioLC HPLC unit (Dionex, Surrey, UK) with a Poly SULFOETHYL A™ column (PolyLC, Columbia, MD, USA). The column used 5 μm particle size, and had dimensions of 200 mm length with an internal diameter of 2.1 mm and a 200 Å pore size. The buffers used were Buffer A (10 mM KH2PO4 with 20% HPLC grade acetonitrile (ACN) at pH 2.5) and Buffer B (10 mM KH2PO4 with 20% HPLC grade acetonitrile (ACN) and 500 mM KCl at pH3) with a gradient of 0% B for 5 min, 0-40% B for 30 min, 40-100% B for 10 min, 100% B for 5 min, and 0% B for 5 min. Before fractionation, the digests were dried in a vacuum centrifuge and resuspended in 200 μL Buffer A. In order to remove salts and unbound iTRAQ reagents, samples were injected onto the column, and manual acquisition was initiated. This allowed Buffer A to flow through the column, leaving bound peptides attached and purging all non-bound material. Two injections were completed before initiation of the gradient. The flow rate was maintained at 200 μL min-1. A UV detector (UVD170U) was used to monitor the fractionation process, and the entire operation was controlled and through Chromeleon software (Dionex/LC packings, the Netherlands). Fractions were collected each minute for the duration of the run with a Foxy Jr. Fraction Collector (Dionex). Fractions were selected for subsequent MS/MS analysis according to SCX chromatogram intensity.
Mass spectrometric analysis was performed on a QStar XL Hybrid ESI Quadrupole time-of-flight tandem mass spectrometer, ESI-qQ-TOF MS/MS (Applied Biosystems, Framingham, MA; MDS-Sciex, Concord, Ontario, Canada), coupled with an Ultimate 3000 online capillary liquid chromatography (LC) system (Applied Biosystems, Framingham, MA), with a PepMap C-18 RP column (LC Packings) set to a constant flow rate of 0.3 μL min-1. Selected fractions were dried and resuspended in 10 μL of Buffer A prior to analysis. To prevent salts, and other contaminants entering the mass spectrometer, a C-18 column clean up (The Nest Group, Inc. Southborough, MA) of each fraction was performed. Cleaned samples were then injected into the online nanoflow LC-MS/MS unit where the buffers were Buffer A (3% (v/v) ACN with 0.1% (v/v) formic acid (FA)), and Buffer B (97% (v/v) ACN with 0.1% (v/v) FA). The gradient applied for the second dimension of chromatography was 3% Buffer B for 3 min, 3 – 40% Buffer B for 70 min, 90% Buffer B for 5 min, and 3% Buffer B for 7 min. The data acquisition protocol in the mass spectrometer was set to positive ion mode, and the addition of formic acid to the buffers was used in order to promote protonation. Peptides were selected by the mass spectrometer for further tandem MS/MS due to observed charge states of +2 and +3.
Mass spectrometry data analysis, peak list conversion, peptide identification and quantitation was carried out using the combination of Analyst QS ver. 1.1 (ABSciex, Foster City, CA) and Phenyx 2.6 (Geneva Bioinformatics, Geneva, Switzerland) software against the UniProt Synechocystis (accessed May 2010). The databases used for searching were Synechocystis PCC6803 (including plasmids), and the NCBI non-redundant database (http://www.ncbi.nlm.nih.gov/, Accessed 01/07/10). Searches of the Synechocystis database were conducted using a target-decoy approach (reversed sequence decoy). All searches were performed with the parameter settings: 0.1 Da MS scan tolerance, 0.15 Da MS/MS scan tolerance, strict trypsin cleavage with one allowable missed cleavage and a minimum of 20% of combined y, y2+, b, b2+ ions. Acceptance threshold for peptide identification were set to p-value ≤ 1 × 10-4 with a peptide score ≥ 5. Protein identification was then set to include only matches with at least a score of 7 (AC score ≥ 7).
Statistical inference was obtained using a t-test comparison between the reporter ions’ intensities, as previously detailed in Mukherjee et al.  and Pham et al. . Within this study, we include observations made with and without a Bonferonni correction.
A principal component analysis (PCA) and a hierarchical clustering were performed on the proteomic data using Mathematica (Wolfram Research) so as to provide an appreciation of biological/stress variation (see Additional file 1: Figures S3 and Additional file 1: Figure S4). These two analyses were carried out on proteins quantified with three or more MS/MS scans. For the PCA, the first two principal components were used; the clustering analysis used Ward’s linkage method.
The authors wish to acknowledge the provision of an EPSRC studentship, Advanced Research Fellowship for CAB (EP/E053556/01) and further EPSRC funding (GR/S84347/01 and EP/E036252/1). Professor Kenneth Reardon and Dr Josselin Noirel also acknowledge the Royal Academy of Engineering Distinguished Visiting Fellowship Programme.
- Schutz K, Happe T, Troshina O, Lindblad P, Leitao E, Oliveira P, Tamagnini P: Cyanobacterial H(2) production - a comparative analysis. Planta. 2004, 218: 350-359. 10.1007/s00425-003-1113-5.View ArticleGoogle Scholar
- Frey KM, Oppermann-Sanio FB, Schmidt H, Steinbuchel A: Technical-scale production of cyanophycin with recombinant strains of Escherichia coli. Appl Environ Microbiol. 2002, 68: 3377-3384. 10.1128/AEM.68.7.3377-3384.2002.View ArticleGoogle Scholar
- Lagarde D, Beuf L, Vermaas W: Increased production of zeaxanthin and other pigments by application of genetic engineering techniques to Synechocystis sp. strain PCC 6803. Appl Environ Microbiol. 2000, 66: 64-72. 10.1128/AEM.66.1.64-72.2000.View ArticleGoogle Scholar
- Burja AM, Banaigsb B, Abou-Mansourc E, Burgess JG, Wright PC: Marine cyanobacteria-a prolific source of natural products. Tetrahedron. 2001, 57: 9347-9377. 10.1016/S0040-4020(01)00931-0.View ArticleGoogle Scholar
- Gill RT, Katsoulakis E, Schmitt W, Taroncher-Oldenburg G, Misra J, Stephanopoulos G: Genome-wide dynamic transcriptional profiling of the light-to-dark transition in Synechocystis sp. strain PCC 6803. J Bacteriol. 2002, 184: 3671-3681. 10.1128/JB.184.13.3671-3681.2002.View ArticleGoogle Scholar
- Burja AM, Dhamwichukorn S, Wright PC: Cyanobacterial postgenomic research and systems biology. Trends Biotechnol. 2003, 21: 504-511. 10.1016/j.tibtech.2003.08.008.View ArticleGoogle Scholar
- Gan CS, Reardon KF, Wright PC: Comparison of protein and peptide prefractionation methods for the shotgun proteomic analysis of Synechocystis sp. PCC 6803. Proteomics. 2005, 5: 2468-2478. 10.1002/pmic.200401266.View ArticleGoogle Scholar
- Fulda S, Mikkat S, Huang F, Huckauf J, Marin K, Norling B, Hagemann M: Proteome analysis of salt stress response in the cyanobacterium Synechocystis sp. strain PCC 6803. Proteomics. 2006, 6: 2733-2745. 10.1002/pmic.200500538.View ArticleGoogle Scholar
- Wassen MJ, Venterink HO, Lapshina ED, Tanneberger F: Endangered plants persist under phosphorus limitation. Nature. 2005, 437: 547-550. 10.1038/nature03950.View ArticleGoogle Scholar
- Suzuki S, Ferjani A, Suzuki I, Murata N: The SphS-SphR two component system is the exclusive sensor for the induction of gene expression in response to phosphate limitation in Synechocystis. J Biol Chem. 2004, 279: 13234-13240.View ArticleGoogle Scholar
- Allenby NE, O’Connor N, Pragai Z, Ward AC, Wipat A, Harwood CR: Genome-wide transcriptional analysis of the phosphate starvation stimulon of Bacillus subtilis. J Bacteriol. 2005, 187: 8063-8080. 10.1128/JB.187.23.8063-8080.2005.View ArticleGoogle Scholar
- Kocan M, Schaffer S, Ishige T, Sorger-Herrmann U, Wendisch VF, Bott M: Two-component systems of Corynebacterium glutamicum: deletion analysis and involvement of the PhoS-PhoR system in the phosphate starvation response. J Bacteriol. 2006, 188: 724-732. 10.1128/JB.188.2.724-732.2006.View ArticleGoogle Scholar
- Martiny AC, Coleman ML, Chisholm SW: Phosphate acquisition genes in Prochlorococcus ecotypes: evidence for genome-wide adaptation. Proc Natl Acad Sci U S A. 2006, 103: 12552-12557. 10.1073/pnas.0601301103.View ArticleGoogle Scholar
- Eymann C, Mach H, Harwood CR, Hecker M: Phosphate-starvation-inducible proteins in Bacillus subtilis: a two-dimensional gel electrophoresis study. Microbiology. 1996, 142: 3163-3170. 10.1099/13500872-142-11-3163.View ArticleGoogle Scholar
- Antelmann H, Scharf C, Hecker M: Phosphate starvation-inducible proteins of Bacillus subtilis: proteomics and transcriptional analysis. J Bacteriol. 2000, 182: 4478-4490. 10.1128/JB.182.16.4478-4490.2000.View ArticleGoogle Scholar
- Madhusudhan KT, McLaughlin R, Komori N, Matsumoto H: Identification of a major protein upon phosphate starvation of Pseudomonas aeruginosa PAO1. J Basic Microbiol. 2003, 43: 36-46. 10.1002/jobm.200390002.View ArticleGoogle Scholar
- He Z, Zhong H, Hu Y, Xiao S, Liu J, Xu J, Li G: Analysis of differential-expressed proteins of Acidithiobacillus ferrooxidans grown under phosphate starvation. J Biochem Mol Biol. 2005, 38: 545-549. 10.5483/BMBRep.2005.38.5.545.View ArticleGoogle Scholar
- Fuszard MA, Wright PC, Biggs CA: Cellular acclimation strategies of a minimal picocyanobacterium to phosphate stress. FEMS Microbiol Lett. 2010, 306: 127-134. 10.1111/j.1574-6968.2010.01942.x.View ArticleGoogle Scholar
- Fuszard M, Wright PC, Biggs CA: Comparative quantitative proteomics of Prochlorococcus ecotypes to a decrease in environmental phosphate concentrations. Aquat Biosystems. 2012, 8: 7-10.1186/2046-9063-8-7.View ArticleGoogle Scholar
- Tetu SG, Brahamsha B, Johnson DA, Tai V, Phillippy K, Palenik B, Paulsen IT: Microarray analysis of phosphate regulation in the marine cyanobacterium Synechococcus sp. WH8102. ISME J. 2009, 3: 835-849. 10.1038/ismej.2009.31.View ArticleGoogle Scholar
- Ishige T, Krause M, Bott M, Wendisch VF, Sahm H: The phosphate starvation stimulon of Corynebacterium glutamicum determined by DNA microarray analyses. J Bacteriol. 2003, 185: 4519-4529. 10.1128/JB.185.15.4519-4529.2003.View ArticleGoogle Scholar
- Pitt FD, Mazard S, Humphreys L, Scanlan DJ: Functional characterization of synechocystis sp. Strain PCC 6803 pst1 and pst2 gene clusters reveals a novel strategy for phosphate uptake in a freshwater cyanobacterium. J Bacteriol. 2010, 192: 3512-3523. 10.1128/JB.00258-10.View ArticleGoogle Scholar
- Adams MM, Gomez-Garcia MR, Grossman AR, Bhaya D: Phosphorus deprivation responses and phosphonate utilization in a thermophilic synechococcus sp. From microbial mats. J Bacteriol. 2008, 190: 8171-8184. 10.1128/JB.01011-08.View ArticleGoogle Scholar
- Orchard ED, Webb EA, Dyhrman ST: Molecular analysis of the phosphorus starvation response in Trichodesmium spp. Environ Microbiol. 2009, 11: 2400-2411. 10.1111/j.1462-2920.2009.01968.x.View ArticleGoogle Scholar
- Juntarajumnong W, Hirani T, Simpson J, Incharoensakdi A, Eaton-Rye J: Phosphate sensing in Synechocystis sp. PCC 6803: SphU and the SphS–SphR two-component regulatory system. Arch Microbiol. 2007, 188: 389-402. 10.1007/s00203-007-0259-0.View ArticleGoogle Scholar
- Hirani TA, Suzuki I, Murata N, Hayashi H, Eaton-Rye JJ: Characterization of a two-component signal transduction system involved in the induction of alkaline phosphatase under phosphate-limiting conditions in Synechocystis sp. PCC 6803. Plant Mol Biol. 2001, 45: 133-144. 10.1023/A:1006425214168.View ArticleGoogle Scholar
- Morohoshi T, Maruo T, Shirai Y, Kato J, Ikeda T, Takiguchi N, Ohtake H, Kuroda A: Accumulation of inorganic polyphosphate in phoU mutants of Escherichia coli and Synechocystis sp. strain PCC6803. Appl Environ Microbiol. 2002, 68: 4107-4110. 10.1128/AEM.68.8.4107-4110.2002.View ArticleGoogle Scholar
- Bhaya D, Vaulot D, Amin P, Takahashi AW, Grossman AR: Isolation of regulated genes of the cyanobacterium Synechocystis sp. strain PCC 6803 by differential display. J Bacteriol. 2000, 182: 5692-5699. 10.1128/JB.182.20.5692-5699.2000.View ArticleGoogle Scholar
- Gomez-Garcia MR, Losada M, Serrano A: Concurrent transcriptional activation of ppa and ppx genes by phosphate deprivation in the cyanobacterium Synechocystis sp. strain PCC 6803. Biochem Biophys Res Commun. 2003, 302: 601-609. 10.1016/S0006-291X(03)00162-1.View ArticleGoogle Scholar
- Chong PK, Gan CS, Pham TK, Wright PC: Isobaric tags for relative and absolute quantitation (iTRAQ) reproducibility: implication of multiple injections. J Proteome Res. 2006, 5: 1232-1240. 10.1021/pr060018u.View ArticleGoogle Scholar
- Gan CS, Chong PK, Pham TK, Wright PC: Technical, experimental and biological variations in isobaric tags for relative and absolute quantitation (iTRAQ). J Proteome Res. 2007, 6: 821-827. 10.1021/pr060474i.View ArticleGoogle Scholar
- Osorio G, Jerez CA: Adaptive response of the archaeon Sulfolobus acidocaldarius BC65 to phosphate starvation. Microbiology. 1996, 142: 1531-1536. 10.1099/13500872-142-6-1531.View ArticleGoogle Scholar
- Ticconi CA, Delatorre CA, Abel S: Attenuation of phosphate starvation responses by phosphite in Arabidopsis. Plant Physiol. 2001, 127: 963-972. 10.1104/pp.010396.View ArticleGoogle Scholar
- Ray JM, Bhaya D, Block MA, Grossman AR: Isolation, transcription, and inactivation of the gene for an atypical alkaline phosphatase of Synechococcus sp. strain PCC 7942. J Bacteriol. 1991, 173: 4297-4309.Google Scholar
- Doonan BB, Jensen TE: Physiological aspects of alkaline phosphatase in selected cyanobacteria. Microbios. 1980, 29: 185-207.Google Scholar
- Bacsi I, Vasas G, Suranyi G, M-Hamvas M, Mathe C, Toth E, Grigorszky I, Gaspar A, Toth S, Borbely G: Alteration of cylindrospermopsin production in sulfate- or phosphate-starved cyanobacterium Aphanizomenon ovalisporum. FEMS Microbiol Lett. 2006, 259: 303-310. 10.1111/j.1574-6968.2006.00282.x.View ArticleGoogle Scholar
- MacColl R: Cyanobacterial phycobilisomes. J Struct Biol. 1998, 124: 311-334. 10.1006/jsbi.1998.4062.View ArticleGoogle Scholar
- Allen MM, Smith AJ: Nitrogen chlorosis in blue-green algae. Arch Microbiol. 1969, 69: 114-120.Google Scholar
- Elmorjani K, Herdman M: Metabolic control of phycocyanin degradation in the cyanobacterium synechocystis PCC 6803: a glucose effect. J Gen Microbiol. 1987, 133: 1685-1694.Google Scholar
- Battchikova N, Vainonen JP, Vorontsova N, Keranen M, Carmel D, Aro E-M: Dynamic changes in the proteome of synechocystis 6803 in response to CO2 limitation revealed by quantitative proteomics. J Proteome Res. 2010, 9: 5896-5912. 10.1021/pr100651w.View ArticleGoogle Scholar
- Wang H-L, Postier BL, Burnap RL: Alterations in global patterns of gene expression in synechocystis sp. PCC 6803 in response to inorganic carbon limitation and the inactivation of ndhR, a LysR family regulator. J Biol Chem. 2004, 279: 5739-5751.View ArticleGoogle Scholar
- Eisenhut M, von Wobeser EA, Jonas L, Schubert H, Ibelings BW, Bauwe H, Matthijs HCP, Hagemann M: Long-term response toward inorganic carbon limitation in wild type and glycolate turnover mutants of the cyanobacterium synechocystis sp. strain PCC 6803. Plant Physiol. 2007, 144: 1946-1959. 10.1104/pp.107.103341.View ArticleGoogle Scholar
- Collier JL, Grossman AR: Chlorosis induced by nutrient deprivation in Synechococcus sp. strain PCC 7942: not all bleaching is the same. J Bacteriol. 1992, 174: 4718-4726.Google Scholar
- Wegener KM, Singh AK, Jacobs JM, Elvitigala T, Welsh EA, Keren N, Gritsenko MA, Ghosh BK, Campll DG, Smith RD, Pakrasi HB: Global proteomics reveal an atypical strategy for carbon/nitrogen assimilation by a cyanobacterium under diverse environmental peturbations. Mol Cell Proteomics. 2010, 9: 12-View ArticleGoogle Scholar
- Seidler A: The extrinsic polypeptides of photosystem II. Biochim Biophys Acta. 1996, 1277: 35-60. 10.1016/S0005-2728(96)00102-8.View ArticleGoogle Scholar
- Shen JR, Inoue Y: Binding and functional properties of two new extrinsic components, cytochrome c-550 and a 12-kDa protein, in cyanobacterial photosystem II. Biochemistry. 1993, 32: 1825-1832. 10.1021/bi00058a017.View ArticleGoogle Scholar
- Inoue-Kashino N, Kashino Y, Satoh K, Terashima I, Pakrasi HB: PsbU provides a stable architecture for the oxygen-evolving system in cyanobacterial photosystem II. Biochemistry. 2005, 44: 12214-12228. 10.1021/bi047539k.View ArticleGoogle Scholar
- Kashino Y, Lauber WM, Carroll JA, Wang Q, Whitmarsh J, Satoh K, Pakrasi HB: Proteomic analysis of a highly active photosystem II preparation from the cyanobacterium synechocystis sp. PCC 6803 reveals the presence of novel polypeptides. Biochemistry. 2002, 41: 8004-8012. 10.1021/bi026012+.View ArticleGoogle Scholar
- Dobakova M, Sobotka R, Tichy M, Komenda J: Psb28 Protein is involved in the biogenesis of the photosystem II inner antenna CP47 (PsbB) in the cyanobacterium synechocystis sp. PCC 6803. Plant Physiol. 2009, 149: 1076-1086.View ArticleGoogle Scholar
- van Thor JJ, Gruters OWM, Matthijs HCP, Hellingwerf KJ: Localization and function of ferredoxin:NADP + reductase bound to the phycobilisomes of Synechocystis. EMBO J. 1999, 18: 4128-4136. 10.1093/emboj/18.15.4128.View ArticleGoogle Scholar
- Stryer L: Biochemistry. 1995, New York: W.H. Freeman Co, 4Google Scholar
- Osanai T, Imamura S, Asayama M, Shirai M, Suzuki I, Murata N, Tanaka K: Nitrogen induction of sugar catabolic gene expression in synechocystis sp. PCC 6803. DNA Res. 2006, 13: 185-195. 10.1093/dnares/dsl010.View ArticleGoogle Scholar
- Clarke AK: ATP-dependent Clp proteases in photosynthetic organisms- a Cut above the rest!. Ann Bot. 1999, 83: 593-599. 10.1006/anbo.1999.0878.View ArticleGoogle Scholar
- Clarke A: Variations on a theme: combined molecular chaperone and proteolysis functions in Clp/HSP100 proteins. J Biosci. 1996, 21: 161-177. 10.1007/BF02703106.View ArticleGoogle Scholar
- Clarke AK, Eriksson M-J: The cyanobacterium Synechococcus sp. PCC 7942 possesses a close homologue to the chloroplast ClpC protein of higher plants. Plant Mol Biol. 1996, 31: 721-730. 10.1007/BF00019460.View ArticleGoogle Scholar
- Stanne TM, Pojidaeva E, Andersson FI, Clarke AK: Distinctive types of ATP-dependent Clp proteases in cyanobacteria. J Biol Chem. 2007, 282: 14394-14402. 10.1074/jbc.M700275200.View ArticleGoogle Scholar
- Pandhal J, Wright PC, Biggs CA: A quantitative proteomic analysis of light adaptation in a globally significant marine cyanobacterium Prochlorococcus marinus MED4. J Proteome Res. 2007, 6: 996-1005. 10.1021/pr060460c.View ArticleGoogle Scholar
- Oliver MJ: Desiccation tolerance in vegetative plant cells. Physiol Plant. 1996, 97: 779-787. 10.1111/j.1399-3054.1996.tb00544.x.View ArticleGoogle Scholar
- Bannister JV, Bannister WH, Rotilio G: Aspects of the structure, function, and applications of superoxide dismutase. CRC Crit Rev Biochem. 1987, 22: 111-180. 10.3109/10409238709083738.View ArticleGoogle Scholar
- Singh SC, Sinha RP, Hader D-P: Role of lipids and fatty acids in stress tolerance in cyanobacteria. Acta Protozool. 2002, 41: 297-308.Google Scholar
- Bandyopadhyay U, Das D, Banerjee RK: Reactive oxygen species: oxidative damage and pathogenesis. Curr Sci. 1999, 77: 658-666.Google Scholar
- Hudson JJ, Taylor WD, Schindler DW: Phosphate concentrations in lakes. Nature. 2000, 406: 54-56. 10.1038/35017531.View ArticleGoogle Scholar
- Marcus Y, Gurevitz M: Activation of cyanobacterial RuBP-carboxylase/oxygenase is facilitated by inorganic phosphate via two independent mechanisms. Eur J Biochem. 2000, 267: 5995-6003. 10.1046/j.1432-1327.2000.01674.x.View ArticleGoogle Scholar
- Allen MM: Simple conditions for growth of unicellular blue-green algae on plates. J Phycol. 1968, 4: 1-4.View ArticleGoogle Scholar
- Axmann IM, Kensche P, Vogel J, Kohl S, Herzel H, Hess WR: Identification of cyanobacterial non-coding RNAs by comparative genome analysis. Genome Biol. 2005, 6: R73-10.1186/gb-2005-6-9-r73.View ArticleGoogle Scholar
- O’Loughlin SN, Graham RLJ, McMullan G, Ternan NG: A role for carbon catabolite repression in the metabolism of phosphonoacetate by Agromyces fucosus Vs2. FEMS. 2006Google Scholar
- Fiske CH, SubbaRow Y: The colourimetric determination of phosphorus. J Biol Chem. 1925, 66: 375-400.Google Scholar
- Mukherjee J, Ow SY, Noirel J, Biggs CA: Quantitative protein expression and cell surface characteristics of Escherichia coli MG1655 biofilms. Proteomics. 2011, 11: 339-351. 10.1002/pmic.201000386.View ArticleGoogle Scholar
- Pham TK, Roy S, Noirel J, Douglas I, Wright PC, Stafford GP: A quantitative proteomic analysis of biofilm adaptation by the periodontal pathogen Tannerella forsythia. Proteomics. 2010, 10: 3130-3141. 10.1002/pmic.200900448.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.