Skip to main content


Front. Microbiol., 08 March 2017
Sec. Aquatic Microbiology
Volume 8 - 2017 |

Glucose Uptake in Prochlorococcus: Diversity of Kinetics and Effects on the Metabolism

  • 1Departamento de Bioquímica y Biología Molecular, Campus de Excelencia Internacional Agroalimentario, Universidad de Córdoba, Córdoba, Spain
  • 2Centre for Proteome Research, Institute of Integrative Biology, University of Liverpool, Liverpool, UK
  • 3National Oceanography Centre, Southampton, UK

We have previously shown that Prochlorococcus sp. SS120 strain takes up glucose by using a multiphasic transporter encoded by the Pro1404 gene. Here, we studied the glucose uptake kinetics in multiple Prochlorococcus strains from different ecotypes, observing diverse values for the Ks constants (15–126.60 nM) and the uptake rates (0.48–6.36 pmol min-1 mg prot-1). Multiphasic kinetics was observed in all studied strains, except for TAK9803-2. Pro1404 gene expression studies during the 21st Atlantic Meridional Transect cruise showed positive correlation with glucose concentrations in the ocean. This suggests that the Pro1404 transporter has been subjected to diversification along the Prochlorococcus evolution, in a process probably driven by the glucose availabilities at the different niches it inhabits. The glucose uptake mechanism seems to be a primary transporter. Glucose addition induced detectable transcriptomic and proteomic changes in Prochlorococcus SS120, but photosynthetic efficiency was unaffected. Our studies indicate that glucose is actively taken up by Prochlorococcus, but its uptake does not significantly alter the trophic ways of this cyanobacterium, which continues performing photosynthesis. Therefore Prochlorococcus seems to remain acting as a fundamentally phototrophic organism, capable of using glucose as an extra resource of carbon and energy when available in the environment.


Prochlorococcus is one of the main picophytoplankters in the oligotrophic oceans of our planet, and its outstanding ecological importance has been exposed in a wide series of studies published since its discovery (Chisholm et al., 1988; Partensky and Garczarek, 2010; Biller et al., 2015). Given its capability to perform oxygenic photosynthesis, it was considered as a phototrophic organism for a long time (Partensky and Garczarek, 2003) and its ability to utilize organic compounds has often been overlooked. However, different groups have demonstrated that this cyanobacterium can utilize organic compounds, such as amino acids (Zubkov et al., 2003; Michelou et al., 2007; Talarmin et al., 2011; Bjorkman et al., 2015), or dimethylsulfoniopropionate (Vila-Costa et al., 2006). These molecules can be utilized for catabolism and/or anabolism, while providing limiting elements in some cases (as nitrogen or sulfur).

We have previously shown that Prochlorococcus can take up glucose at high rates in cultures, inducing specific changes in the expression of genes related to glucose metabolism (Gómez-Baena et al., 2008). More recently, we described that the Pro1404 gene product in Prochlorococcus is a glucose transporter with biphasic kinetics and a high affinity constant (Ks) in the nanomolar range (Muñoz-Marín et al., 2013). Furthermore, we observed that glucose is taken up by natural Prochlorococcus populations in the Atlantic Ocean, in environments where glucose is available at nanomolar concentrations (Muñoz-Marín et al., 2013). In good agreement with our findings, some studies have shown that glucose supplementation induces an increase in Prochlorococcus populations in the oligotrophic Pacific Ocean (Moisander et al., 2012). However, some essential aspects of this process have not yet been explored. In particular, the mode of glucose transport, how kinetic parameters of the transporter vary across the Prochlorococcus radiation, and how glucose uptake affects the physiology of Prochlorococcus are not known and could inform the evolutionary history of this process.

In this work, we analyzed the glucose transport kinetics in several strains of Prochlorococcus representative of the different ecotypes of this cyanobacterium. We characterized the mode of transport performed by Pro1404, by using several transport and photosynthesis inhibitors. We measured the expression of several genes related to glucose metabolism under low glucose concentrations; in addition, we studied the Pro1404 expression in field samples obtained during the AMT21 cruise. We have also used proteomics to assess how addition of glucose affects Prochlorococcus sp. SS120 cultures.

Materials and Methods

Seawater Samples Collection

Seawater samples were collected along the Atlantic Meridional Transect 21 cruise (AMT-21), from September 29 to November 14, 2011, at 11 stations (Supplementary Table S1). The seawater was prefiltered through 20 μM mesh and vacuum filtered through 47 mm diameter, 0.2 μm filters (Millipore). Filters were immersed into 4 mL of RNA resuspension buffer (20 mM EDTA; 400 mM NaCl; 0.75 mM sucrose; 50 mM Tris-HCl pH 7) and immediately frozen in liquid nitrogen (Holtzendorff et al., 2002).

Cyanobacterial Laboratory Strains and Growth Conditions

Prochlorococcus spp. strains PCC 9511 [axenic, high-irradiance (HL)-adapted], TAK9803-2 [high-irradiance (HL)-adapted], and SS120, MIT9303, NATL1-A, and NATL2-A [low-irradiance (LL)-adapted] were cultured and collected as previously described (El Alaoui et al., 2001). Cultures were grown in PCR-S11 medium in a culture room set at 24°C under continuous blue irradiances (40 μE/m2/s for HL adapted ecotypes and 4 μE/m2/s for LL adapted ecotypes) using neon Sylvania F18W/154-ST Daylight, covered with a filter Moonlight blue L183 from Lee Filters.

Cell Extracts

For proteomic studies, the cells were broken using a French press cell (SLM/Aminco model FA-079) at 1,000 psi; the obtained extract was centrifuged for 10 min at 16,000 g and 4°C and the supernatant was transferred to a clean tube prior analysis. The cell suspensions stored in Tris buffer were placed in sterile 1.5 mL eppendorf tubes containing glass beads 0.10–0.11 mm of diameter (B. Braun Melsungen AG, Melsungen, Germany). Cell disruption was performed by five cycles of 2 min vortex and 2 min in ice. Tubes were then centrifuged at 16,000 g for 5 min at 4°C and the supernant transferred to a new sterile 1.5 mL tube prior analysis.

Protein concentration was determined using the Bio-Rad Protein Assay kit (Bradford, 1976).

Determination of Glucose Uptake by Cultured Cyanobacterial Strains

[U-14C]-glucose (281 mCi/mmol, American Radiolabeled Chemicals, Inc.) was added to 0.7 mL culture samples to reach final concentrations ranging from 0.01 to 20 μM, depending on the experiment. Mixtures of radiolabeled and unlabeled glucose were added to achieve the desired concentration. Glucose uptake, maximum uptake rates and Ks constants were determined by the method of Eadie–Hofstee plots (Azam and Hodson, 1981), as previously described (Muñoz-Marín et al., 2013).

Studies with Inhibitors

[U-14C]-glucose was added at a final concentration of 1 μM to a 4.5 mL sample of Prochlorococcus SS120 cultures supplemented with different inhibitors (Sigma-Aldrich). Carbonyl cyanide m-chlorophenylhydrazone (CCCP), monensin and valinomicin were disolved in DMSO. N,N-dicyclohexylcarbodiimide (DCCD) and (3,4-dichlorophenyl) 1,1 dimethylurea (DCMU) were dissolved in 95% ethanol. 2,5-dibromo-3-methyl-6-isopropyl-p-benzoquinone (DBMIB) was prepared in 100% methanol. The final concentrations in the cultures were: 50 μM CCCP, 10 μM monensin, 15 μM valinomycin, 200 μM DCCD, 0.3 μM DCMU, and 0.06 μM DBMIB. All inhibitors were prepared as previously described (El Alaoui et al., 2001; Soontharapirakkul et al., 2011). The effect of ethanol, methanol and DMSO were examined by adding the same concentration to the cultures as above: no substantial alteration in the glucose uptake was noted (not shown).

Determination of Quantum Yield of Photosystem II

Culture samples (250 mL) were centrifuged at 26,000 g for 8 min at 4°C. The pellet was resuspended in 2 mL of PCR-S11 medium and placed in a 24 well culture plate (Biofil). Chlorophyll fluorescence was measured using an IMAGING-PAM chlorophyll fluorometer (M-Series, Maxi version, Heinz Walz GmbH). Cells were dark adapted for 30 min prior to determination of Fo and Fm.

RNA Isolation from Laboratory Prochlorococcus Cultures

RNA was isolated from 500 mL cultures aliquots subjected to either control (standard growth conditions, no glucose added) or glucose conditions (supplemented with 0.1 μM glucose), as described (Domínguez-Martín et al., 2014).

RNA Isolation from Natural Samples

RNA from natural oceanic samples was isolated with a method modified from that described by Lindell and Post (2001). Filtered seawater samples were mixed with 200 μL of 10% SDS (final concentration 1%) and 18 U of proteinase K (from Tritirachium album, Sigma-Aldrich) and incubated at 60°C for 20 min. Then samples were transferred to sterilized Corex tubes and 1 mL of phenol was added. Samples were vortexed and centrifuged for 15 min at 15,344 g. After two extractions with phenol:chloroform:isoamilic alcohol (25:24:1) and one extraction with chlorophorm:isoamilic alcohol (24:1), in all cases with 15 min centrifugation at 15,344 g, 2 mL of 70% isopropanol and 800 μL of ammonium acetate, pH 7.5 were added, and the mixture was kept at -20°C overnight. Samples were centrifuged again for 30 min at 15,344 g at 4°C, washed with 70% ethanol, and resuspended in 17 μL of nucleases-free water, 2 μL DNAsa 10x buffer (Ambion) and 2 μL DNAsa (4 U), and incubated for 30–40 min at 37°C. Finally, DNA was removed by addition of 2 μL DNAse Inactivation Reagent. The mixture was incubated for 5 min at room temperature, centrifuged for 5 min at 15,344 g, and the supernatant used for further RNA analysis.

Real-Time Quantitative RT-PCR Analysis of Gene Expression

The synthesis of the cDNA by the reverse transcriptase (RT) reaction was carried out using the iScriptTM cDNA Synthesis kit (Quanta) as recommended by the manufacturers. RNA concentration was quantified using a NanoDrop instrument (Thermo Scientific). For 20 μL of total volume reaction, 0.4 μg of RNA was reverse transcribed for Prochlorococcus cultures and 0.1 μg of RNA for field samples. Specific primers to amplify fragments of the the genes of interest from the Prochlorococcus strain SS120 were designed using the software Oligo 4.05 (Molecular Biology Insights, Inc.), on the basis of the corresponding Prochlorococcus SS120 genome (Dufresne et al., 2003). The primers were tested for Prochlorococcus specificity using Primer-BLAST from NCBI. The sequences of the primers used are compiled in Table 1.


TABLE 1. List of primers utilized in quantification of gene expression by semi-quantitative RT-PCR.

As a quality control step, all samples were tested by PCR to check the presence of residual DNA before and after the DNAse treatment. Furthermore, negative control reactions without RT were used to verify the absence of DNA template contamination.

Real time quantitative PCR reactions were performed in triplicate. The reaction mixtures contained 1x of SsoFastTM EvaGreen® Supermix (Bio-Rad), 0.128–0.384 μM forward and reverse primers, depending on the efficiency calculations (not shown) and the corresponding cDNA. The efficiency of the reactions was calculated and optimized following the method described previously (Pfaffl, 2001).

The products were checked for single amplification of DNA fragments of the expected size by agarose gel electrophoresis and by sequencing in the case of seawater samples (not shown). Control reactions were performed in triplicates for each pair of primers in the absence of template.

An iCycler IQ multicolor real time PCR detection system (Bio-Rad) was used for quantitative detection of amplified PCR products using the following thermal cycling conditions: 95°C for 2 min, and 50 cycles of 95°C for 15 s, followed by 58°C for 30 s and 72°C for 30 s. At the end, reactions were checked to discard false amplifications by verifying the melting point of PCR products, determining the fluorescence between 65 and 100°C, with increases of 0.5°C measured each 10 s.

Measurements were carried out in triplicate from at least three independent biological samples subjected to identical culture conditions. The relative change in gene expression was endogenously normalized to that of the gene rnpB (Table 1), encoding RNase P, calculated using the 2-ΔΔCt method (Pfaffl, 2001). The primers rnpRTfwd/rnpRTrev were used for the endogenous standard (rnpB) in the natural samples (Holtzendorff et al., 2002). No change in the expression of rnpB was detected under our experimental conditions.

Prochlorococcus Genomic Sequences

Cyanobacterial genomic data were obtained from the Joint Genome Institute1,2 .

Phylogenetic Analysis

Proteins similar to the glucose transporter Pro1404 were screened within published sequences of Synechococcus and Prochlorococcus by using the blastp algorithm from the NCBI3 imposing a minimum length of 339 amino acids for alignment (>75% of the sequence) and a minimum of 35% identity for Synechococcus and 65% identity for Prochlorococcus sequences.

GOS sequences were collected by the blastp algorithm against the Pro1404/glcH protein from the NCBI in the CAMERA portal4 imposing the same restrictions as above. GOS sequences giving a maximum score with Prochlorococcus and Synechococcus were aligned, using the software SeaView5, with the sequences of Prochlorococcus and Synechococcus detected in the first analysis. Ambiguous areas of alignment were removed using Gblocks (Thermo Scientific), yielding 198 common amino acids. The phylogenetic tree was based on proteins sequences similarity using the maximum likelihood method in the server RaxML BlackBox6 under default parameters, using 100 bootstrap replicates and edited by the web-based tool Interactive Tree of Life7. The sequences corresponding to Synechococcus sp. PCC 7002, Microcystis aeruginosa, Synechocystis sp. PCC 6803, and Thermosynechococcus elongates BP-1 acted as outgroups in this tree.

Proteomic Analysis

Cell cultures were grown under standard conditions or in medium supplemented with 0.1 μM glucose for 24 h, before harvesting cells to carry out proteomic analysis. Cell extracts containing 100 μg of protein were incubated with 0.05% (final concentration w/v) of RapiGest (Waters Corporation) for 10 min at 80°C. Samples were then reduced with 3 mM DTT for 10 min at 60°C, followed by carbamidomethylation with 9 mM iodoacetamide for 30 min in the dark at room temperature. Finally, trypsin (ratio 50:1) was added and samples incubated overnight at 37°C. Digestions were stopped by addition of trifluoroacetic acid (TFA) at final concentration of 0.5% (v/v) followed by incubation for 30 min at 37°C. To remove acid inactivated detergent and all insoluble material, sample digests were centrifuged at 13,000 g for 15 min. Completeness of digestion was checked by SDS-PAGE (not shown).

Three biological replicates of each condition (control vs. glucose) were used in the analysis, and tryptic digests were analyzed on a QExactive mass spectrometer (Thermo Scientific). Relative and absolute abundances were calculated by a label free approach using Progenesis QI (Waters Corporation). Specifically, the Hi3 method (Silva et al., 2006) was used to gain absolute quantification of the detected proteins. This method is based on the principle that the intensity of the three most intense peptides obtained after trypsin digestion is proportional to the abundance of the protein, and also constant in all proteins, with a 10% error. In our experiments, 50 fmol of rabbit glycogen phosphorylase (Waters Corporation) were added to the sample as internal standard for absolute quantification.

LC–MS/MS analysis was performed using a QExactive instrument (Thermo Scientific) coupled to an Ultimate 3000 chromatography system (Thermo Scientific). Protein digests were trapped on a PepMap 100 C18 precolumn, 300 μm id, 5 mm long, 5 μm particles, 100 Å (Thermo Scientific) for 7 min, at a flow rate of 4 μL/min in 2% (v/v) ACN/0.1% (v/v) TFA, and peptides were separated through a analytical Easy-spray PepMap RSLC C18 column, 75 μm id, 50 cm long, 2 μm particles, 100 Å (Thermo Scientific) at 300 nL/min over a 120 min linear gradient from 3 to 40% (v/v) ACN in 0.1% (v/v) formic acid. The QExactive instrument was operated in data dependent acquisition mode. Full scan MS spectra (m/z 300–2000) and fragment ions were acquired in the orbitrap at 70,000 and 35,000 resolution, respectively. Dynamic exclusion was set to 20 s. The 10 most intense multiply charged ions (z ≥ 2) were sequentially isolated and fragmented by high energy collisional dissociation (HCD) at 30% standardized HCD collision energy.

Proteome Discoverer version 1.4 (Thermo Scientific) was used to generate peaklists using default parameters and MASCOT search engine (version 2.4.0, Matrix Science) to identify peptides and proteins, using a database containing all the entries annotated for Prochlorococcus marinus sp. SS120 (downloaded on 290415, 1,881 sequences) in Uniprot8 and applying a FDR < 1% using the Decoy database tool in MASCOT. Trypsin was selected as specific enzyme allowing one missed cleavage. Carbamidomethylation of cysteines was selected as fixed modification and methionine oxidation as variable, mass tolerance of 10 ppm for precursors and 0.01 Da for fragment ions.

Quantitative data were obtained using Progenesis QI software (Waters Corporation). Features with positive charge states between +2 and +7, and three or more isotopic peaks were taken to further identification. A merged peak list generated by Progenesis QI and searched against the above described database in MASCOT search engine (version 2.4.0, Matrix Science) using the search parameters described above. MASCOT results were imported in Progenesis QI and a peptide cut off score of 20 was applied after manually evaluating the quality of the lowest scored peptides. Hi3 expression data was obtained using the Progenesis QI tool. For this purpose, 50 fmol of glycogen phosphorylase (SwissProt P00489, MassPREP digestion standard, Waters Corporation) were used as internal standard. Similar proteins were grouped and only non-conflictive features (unique peptides) were used for quantification. For a protein to be considered as significantly differentially expressed, it has to be identified and quantified using at least two unique peptides and an ANOVA p-value ≤ 0.05.

Statistical Analysis

Experiments were carried out at least with three independent biological replicates. The results are shown with error bars corresponding to the standard deviation. Significance of data was assessed by using the Student’s t-test (p ≤ 0.05; ∗∗p ≤ 0.01).

Results and Discussion

Diversity of Glucose Uptake Kinetics in Prochlorococcus Strains

In previous studies, we compared the glucose uptake kinetics of Prochlorococcus marinus SS120 with that of Synechococcus WH7803, a potential competitor in oligotrophic areas of the oceans (Muñoz-Marín et al., 2013). By analyzing the Ks constant and corresponding Vmax, we observed that Prochlorococcus was almost seven times more efficient than Synechococcus (Muñoz-Marín et al., 2013). In the present work, we assessed possible differences within the Prochlorococcus genus, by comparing the glucose uptake kinetics among different strains representative of different evolutive clades: Prochlorococcus spp. strains PCC 9511, MIT9303, NATL1A, NATL2A, and TAK9803-2. The strain TAK9803-2 belongs to the recently evolved HL II clade (West et al., 2001). The strains PCC 9511 [genetically identical to MED4 (Rippka et al., 2000)] belongs to the HL I clade (West et al., 2001; Martiny et al., 2009), also late-branching in the Prochlorococcus radiation. The strains NATL1A and NATL2A belong to the LL I clade, and SS120 to the LL II/III: both clades occupy an intermediate position in the phylogeny of Prochlorococus. Finally, the strain MIT9303 is one of early-branching members of the Prochlorococcus radiation and belongs to the LL IV clade (Partensky et al., 1999; West et al., 2001; Martiny et al., 2009). The selected strains cover all Prochlorococcus clades with cultured representatives (Biller et al., 2015). We used one axenic strain (PCC 9511, Palinska et al., 2000), while the rest were non-axenic cultures. The fact that the cultures are non-axenic might raise some concern about the results; however, the level of contamination in our cultures is very low (Gómez-Baena et al., 2008); furthermore, we showed that the kinetics of glucose uptake was the same (including both Ks value and multiphasic kinetics) when comparing the non-axenic Prochlorococcus sp. strain SS120 with an axenic, recombinant strain of Synechococcus sp. strain PCC 7942 (naturally incapable of glucose uptake) expressing the melB gene from Prochloroccus SS120 (Muñoz-Marín et al., 2013). Therefore, the low level of heterotrophic bacteria present in the culture does not interfere in the results and they reflect the actual values of glucose uptake in Prochlorococcus.

Glucose uptake kinetics (maximum uptake rate and Ks constant) were determined in the above mentioned strains using laboratory cultures. Table 2 summarizes the results and Figures 1A,B show the kinetics observed for each strain, including in both cases the previously described results for SS120 and WH7803 (Muñoz-Marín et al., 2013) as a reference. Figure 1C depicts a simplified phylogenetic tree of the Prochlorococcus radiation [adapted from the review by Chisholm and coworkers (Biller et al., 2015)], including also marine Synechococcus, showing the Ks values and the efficiency for each strain.


TABLE 2. Kinetics parameters in the Prochlorococcus sp. strains PCC 9511, MIT9303, NATL1A, NATL2A, and TAK9803-2.


FIGURE 1. Kinetics studies of glucose uptake in Prochlorococcus sp. strains PCC 9511, MIT9303, NATL1A, NATL2A, and TAK9803-2. (A) 14C-Glucose concentration vs. 14C-glucose uptake. (B) 14C-Glucose concentration vs. its uptake with a different scale in the axis to allow a proper representation of the first phase of the kinetics. Previously described results for SS120 and WH7803 (Muñoz-Marín et al., 2013) are plotted in both graphs as a reference. Data presented are the average of at least three independent biological replicates. Bars show standard deviation for each sample. (C) Simplified phylogenetic tree of marine picocyanobacteria based on the rRNA internal transcribed spacer [ITS; adapted from that published by S. W. Chisholm and coworkers (Biller et al., 2015)] showing the evolution of glucose uptake. Numbers next to the name of each strain show the GlcH transporter Ks value (in blue), followed by the efficiency of glucose uptake in pink (as shown in Table 2).

We observed some diversity for both the Ks (15 nM in PCC 9511 vs. 124.6 nM in SS120) and the Vmax values (0.16 pmol min-1 mg prot-1 in WH7803 vs. 6.36 pmol min-1 mg prot-1 in SS120). When both parameters were taken together in order to calculate the efficiency in the glucose uptake (Vmax /Ks), the results ranged from 0.0074 to 0.12 (for WH7803 and PCC 9511, respectively); in other words, Prochlorococcus PCC 9511 is 16 times more efficient for glucose uptake than Synechococcus WH7803. Interestingly, when we correlated the kinetics parameters of glucose transport to the marine picocyanobacteria phylogeny (Figure 1C), we observed an increase in the efficiency of the process, which was highest in late-branching HL I strain PCC 9511, and lowest in Synechococcus WH7803. It is worth mentioning that our results only include one Synechococcus strain, and that another late-branching strain (TAK9803-2, belonging to the HLII clade) shows a value of efficiency similar to earlier branching groups as LL I and LLII/III (i.e., NATL1A, NATL2A, SS120), and thus this is a proposal which should be considered with caution. Further studies addressing a wider range of strains are required to confirm our hypothesis of the trends toward increased efficiency of glucose uptake in marine picocyanobacteria.

The energetic cost of glucose transport is lower than that of its biosynthesis (Muñoz-Marín et al., 2013). These levels of diversity in the glucose uptake suggest that this process might be important for marine cyanobacterial strains, and therefore an increased competence for the uptake of glucose could be a selective advantage under specific conditions (i.e., extended darkness; Coe et al., 2016), thus in part driving the diversification described above. The trend to efficiency increasing for glucose uptake described above fits nicely with this hypothesis.

The maximum glucose uptake rate was much higher in Prochlorococcus sp. SS120 than in other strains. Although SS120 is adapted to live at depth where light availability is low, these environments have higher nutrient concentrations than the surface ocean (Cullen, 2015; Yelton et al., 2016). These results are in good agreement with previous studies, showing that the mixotrophy genes are more abundant in low light-adapted ecotypes of Prochlorococcus (Yelton et al., 2016) and with higher amino acid uptake rates than populations close to the surface (Zubkov et al., 2004). This might suggest that glucose uptake at depth might be more advantageous for Prochlorococcus, although this requires experimental confirmation in the field.

The results showed multiphasic kinetics in all studied strains except in TAK9803-2 (Figures 1A,B). The third phase was not fully characterized due to the high concentrations of glucose required to reach saturation, which make them meaningless from an ecological point of view, since the glucose concentrations observed in the oceans are usually in the nanomolar range (Azam and Hodson, 1981; Rich et al., 1996; Muñoz-Marín et al., 2013). It is not clear yet why TAK9803-2 lacks the biphasic kinetics, despite being quite closely related to PCC 9511 (Figure 1C). The loss of multiphasic kinetics might be another step in the evolution of glucose uptake in Prochlorococcus, taking into account that TAK9803-2 belongs to one of the late-branching clades (HL II). Future studies will address this possibility.

All strains showed a high affinity Ks constant (in the nanomolar range) in the first phase of glucose uptake, suggesting that this is a feature of the Pro1404 transporter from marine cyanobacteria (Figure 1B and Table 2). Therefore we propose to change the current name of this gene (melB, for melibiose transporter) to a more appropriate glcH, since we have shown it is a very high affinity (H) glucose transporter (glc), with higher affinity for glucose than for other sugars (Muñoz-Marín et al., 2013).

Expression of glcH in Natural Populations of Prochlorococcus

We have shown that Prochlorococcus takes up glucose in the Atlantic Ocean (Muñoz-Marín et al., 2013). In order to check whether glcH expression (previously annotated as melB) was also regulated in the natural habitat of Prochlorococcus, as it seems to be the case in laboratory cultures (Gómez-Baena et al., 2008), we determined glcH expression in samples obtained during the AMT21 cruise.

RNA samples were collected at 13 sites9 (Supplementary Table S1; Muñoz-Marín et al., 2013), and the glcH expression was determined using the same glcH/melB primers previously utilized for the SS120 strain (Figure 2A) (Gómez-Baena et al., 2008). Moreover, we compared the glcH expression with the glucose bioavailability and the glucose uptake by Prochlorococcus measured during the transect (Muñoz-Marín et al., 2013). Unfortunately, we could only compare the glucose bioavailability to glcH expression in five stations, since it was not possible to sample for RNA expression in all stations due to time constraints (Figure 2B). We observed that the areas where we found the highest glcH expression matched those where Prochlorococcus population contributed to maximum glucose transport 2.6–3.7% of the total bacterioplankton glucose uptake (Muñoz-Marín et al., 2013).


FIGURE 2. Expression of glcH in natural populations of Prochlorococcus. (A) Relative expression of glcH in natural samples of Prochlorococcus obtained during the AMT-21 cruise. (B) Relationship between the glucose concentration found in the sea (orange) and the glcH expression in natural samples of Prochlorococcus (green). Information on the specific stations of the AMT21 cruise are available in a previous paper (Muñoz-Marín et al., 2013) and in the AMT21 Cruise Report (

The highest glcH expression was observed at station 35, 39 and 50 of AMT21, where glucose concentration was higher as well (Figure 2B). Our results suggest a positive correlation between glucose concentration and glcH expression in Prochlorococcus, but this hypothesis requires further studies in the field for confirmation. This is in good agreement with previous results of our team, showing increased glcH expression after glucose addition to Prochlorococcus sp. SS120 cultures (Gómez-Baena et al., 2008).

Phylogeny of the Pro1404 Gene

We screened for genes encoding proteins similar to the glucose transporter (melB/glcH) from Prochlorococcus SS120 and used them to build a phylogenetic tree with similar proteins in 74 sequenced strains of Prochlorococcus and Synechococcus (Supplementary Figure S1). To further understand glcH transcription in natural samples, we searched in the metagenomic libraries for the glucose transporter gene Pro1404/glcH from Prochlorococcus to check whether similar genes are often detected in the ocean. The samples used were obtained at the Global Ocean Survey (GOS; covering sites in the Atlantic, Pacific, and Indian Oceans; August 2003–January 2006) (Rusch et al., 2007; Yooseph et al., 2007; Williamson et al., 2008).

The phylogenetic analysis between Prochlorococcus and Synechococcus sequences revealed three clades: the entire marine Synechococcus sequences, which represent clades a and b, and another distinct clade c belonging one sub-clade to the Synechococcus genus and the big sub-clade of the Prochlorococcus genus. The glucose transporter (melB/glcH) is highly conserved in the Prochlorococcus strains, being present in all of them (Dufresne et al., 2003; Rocap et al., 2003; Biller et al., 2014). This fact, taken together with the variation in the kinetic parameters (Figure 1A) and the diversity observed in its phylogeny (Supplementary Figure S1), suggests that the glcH gene is important for Prochlorococcus and is selectively maintained. The large diversity of GOS sequences observed in the HLII clade provide additional support to this hypothesis. Therefore the efficiency in glucose transport, and the mixotrophic capabilities it confers, might be an important factor in the diversification of marine picocyanobacterial ecotypes.

Recent studies addressed the global genetic capacity for mixotrophy in marine cyanobacteria (Yelton et al., 2016), showing a trend to loss genes encoding transporters of organic compounds from the LLIV to the high HLII clades in Prochlorococcus. Moreover, the distribution of glcH genes in the ocean (Supplementary Figure S1, Yelton et al., 2016) is consistent with the presence of the different Prochlorococcus and Synechococcus ecotypes, in good agreement with our hypothesis. This fact is even more remarkable if we take into account that coastal isolates of the highly abundant SAR11 heterotrophic marine bacteria could utilize glucose as a carbon source, but some open ocean isolates could not (Schwalbach et al., 2010). Furthermore, all Prochlorococcus clades have a single glucose transporter gene (glcH), except clade LL IV where the presence of glcP [a specific glucose permease found in free-living and symbiotic cyanobacteria (Picossi et al., 2013)] has also been reported (Yelton et al., 2016). This reinforces the importance of the high affinity GlcH glucose transporter for Prochlorococcus; interestingly, all the strains studied in the sub-cluster 5.1 B of Synechococcus possess more than one gene encoding sugar transporter (Yelton et al., 2016).

Another interesting aspect described by Yelton et al. (2016) is that the average number of sugar transporter genes per genome (both in Prochlorococcus and Synechococcus) is roughly similar in coastal vs. open ocean regions. Therefore, the kinetic diversity described in our work (Figures 1A,B and Table 2) suggests that selective adaptation might be acting on the efficiency of the sugar uptake (Figure 1C), rather than on the number of sugar transporter genes.

Characterization of the GlcH Transporter with Specific Inhibitors

Transport Inhibitors

The energetics of the glucose transport system of Prochlorococcus SS120 was also studied. Our initial hypothesis considered that glucose uptake could be most useful for Prochlorococcus under energy-limited conditions; however, we have shown that the glucose uptake rate was higher under light, with a circa 40% decrease when cells were subjected to darkness (Gómez-Baena et al., 2008), in accordance with what has been described in previous studies on uptake of organic compounds (Rubin et al., 1977; Church et al., 2004, 2006; Michelou et al., 2007; Mary et al., 2008; Talarmin et al., 2011; Gómez-Pereira et al., 2013; Bjorkman et al., 2015). To evaluate if GlcH is a primary or secondary active transporter, we used different inhibitors of active transport to evaluate if GlcH is a primary or secondary active transporter (i.e., uses energy provided by ATP or from ionic gradients in order to take up glucose, respectively), namely monensin, valinomycin, DCCD, and CCCP.

Monensin is an antibiotic with high selectivity for sodium over other metal cations while valinomycin shows a high selectivity for potassium (Hilgenfeld and Saenger, 1982). The inhibitor DCCD reacts with the carboxyl group of the conserved amino acid residue in the c-subunit present in all ATP synthases to inactivate the enzyme (Shoshan and Selman, 1980; Yoshida et al., 1981). CCCP is classified as an uncoupler increasing the proton permeability of the membrane, thereby preventing the formation of the proton gradients needed for ATP synthesis (Mitchell, 1979). Therefore, monensyn and valinomycin are inhibitors of secondary transport, while CCCP and DCCD are inhibitors of primary transport.

Figure 3A shows the results obtained after addition of these inhibitors to Prochlorococcus sp. SS120 cultures. Monensin and valinomycin had no effect on the transport, compared to the control values. This suggests that the transport mediated by GlcH is not a symport using either sodium or potassium as ion species. On the other hand, both CCCP and DCCD induced a complete inhibition of glucose uptake, indicating that ATP biosynthesis is essential for this transport, and therefore blocking that process does prevent the uptake of glucose. Therefore, we propose GlcH is an active primary transporter, depending on the generation of ATP to carry out glucose transport. However, the protein encoded by glcH has two conserved domains in the sequence (MFS super family and MFS_2 domains), which are typical features of secondary transporters of the Major Facilitator Superfamily (MFS), using Na+ as ion (Paulsen et al., 1998; Ethayathulla et al., 2014). The reasons of this discrepancy will be explored in future studies.


FIGURE 3. Characterization of the GlcH transporter of Prochlorococcus sp. SS120 with specific inhibitors. (A) Effect of monensin, CCCP, valinomycin and DCCD on glucose uptake. (B) Effect of DCMU and DBMIB on glucose uptake. In both graphs the data are the average of three independent biological replicates. Bars show standard deviation for each sample.

The active nature of glucose transport in Prochlorococcus (Figure 3A) and the fact that glucose uptake is lower in the dark in this organism (Gómez-Baena et al., 2008) fit nicely with the hypothesis proposing that Prochlorococcus might invest a significant amount of energy directly from photophosphorylation into organic nutrient acquisition in very oligotrophic environments, probably at the expense of CO2 reduction (Zubkov, 2009).

Electron Transport Inhibitors

We wanted to explore the participation of photosynthetic electron transport in the glucose uptake regulation. DCMU blocks the transfer of electrons between the PSII and plastoquinone pool (Trebst, 1980) whereas DBMIB prevents electron transfer between the plastoquinone pool and cytochrome b6f (Rich et al., 1991). Hence, DCMU allows cyclic electron flow around photosystem I (Bailey et al., 2008; Gao and Wang, 2012).

DBMIB almost completely inhibited the glucose uptake, while DCMU decreased it ca. 50% (Figure 3B). The different inhibitory strength of DBMIB and DCMU on the glucose uptake could indicate that the ATP generated by the cyclic electron flow of PS I can maintain up to 50% of the glucose uptake. This is in agreement with previous studies with inhibitors, including marine Synechococcus (Bailey et al., 2008), Prochlorococcus (Zorz et al., 2015), and other organisms (Gao and Wang, 2012).

Metabolic Effects of Glucose Uptake in Prochlorococcus

To analyze the effects of glucose uptake on Prochlorococcus metabolism, we utilized three approaches: first we studied the expression of a number of genes involved in the glucose metabolism in Prochlorococus sp. SS120; then we studied the effects on photosynthetic efficiency after glucose addition to cultures of several Prochlorococcus strains, and finally we assessed proteomic changes in Prochlorococcus sp. SS120. We utilized this strain because it belongs to the low light ecotype, which we expected to be the most benefitted from glucose uptake. For this reason, a significant part of our previous studies on glucose uptake has been carried out in Prochlorococcus sp. SS120 (Gómez-Baena et al., 2008; Muñoz-Marín et al., 2013).

Effects of Glucose Addition on Gene Expression in Prochlorococcus SS120

We have previously described some changes on the gene expression in Prochlorococcus sp. SS120 after 1 μM glucose addition (Gómez-Baena et al., 2008). In the present study, we performed 24 h time-course experiments to determine expression of a number of additional genes after addition of 0.1 μM glucose, which is closer to the concentration observed in oligotrophic oceans (Wright and Hobbie, 1966; Azam and Hodson, 1981; Rich et al., 1996; Zubkov et al., 2008; Muñoz-Marín et al., 2013). In order to assess whether glucose availability modifies the main carbon incorporation pathways in Prochlorococcus, genes involved in either glucose metabolization or photosynthesis were selected. These included phosphoribulokinase (frb), fructose 1,6-biphosphate aldolase (fbpA), glucose 1-phosphate-adenylyltransferase (g1pA), transaldolase (tal), the small and large rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase) subunits (rbsSL), pyruvate kinase (prk) and glyceraldehide 3-phosphate dehydrogenase (g3pA).

The results were heterogeneous and no major changes in the expression levels were found (Figure 4). However, for several genes, an increase in expression was detected after 8 and 24 h (i.e., frb and tal, respectively); in other cases, a small increase was observed at the beginning of the experiments (i.e., g1pa, rbcS).


FIGURE 4. Effect of glucose addition on the gene expression of Prochlorococcus sp. SS120. Glucose (0.1 μM final concentration) was added to Prochlorococcus cultures. Samples were taken at the indicated times and the gene expression measured by real time RT-PCR. Expression level = 1 indicates no changes with respect to the control (no addition). Bars indicate the average of six independent determinations. Error bars show standard deviation for each sample. Means significant with respect to the samples obtained at T = 0.

Interestingly, the Rubisco genes rbcL and rbcS decreased their expression at 5 and 8 h, while the gene encoding transaldolase (tal) increased its expression at those times. This is consistent with a redirection of the flux from the Calvin cycle to the pentoses phosphate pathway; similar changes have been reported in phage-infected Prochlorococcus strains (Thompson et al., 2011). Furthermore, this is in agreement with several studies in cyanobacteria and plants which showed that the transcript level of rbcLS and other photosynthetic genes decreases under mixotrophic conditions (Quick et al., 1991; Nieva and Fernández Valiente, 1996; Smeekens and Rook, 1997; Yang et al., 2002).

Overall, our data on gene expression suggest that Prochlorococcus cells were reacting to the availability of glucose, and probably redirecting to some extent the metabolic fluxes; but the transcriptional response was not being profoundly affected, thus indicating that glucose does not significantly change the trophic strategy of Prochlorococcus. This is coherent with the fact that glucose is a central metabolite, and therefore a higher availability of this sugar would not require deep metabolic changes for its utilization.

Effects of Glucose Addition on the Photosynthetic Efficiency in Different Prochlorococcus Strains

Although expression studies suggested that the transcription of some photosynthetic genes was not affected by glucose addition (Figure 4), we addressed whether it might affect the photosynthetic efficiency in different Prochlorococcus strains. To this goal, we measured the effective photochemical quantum yield of PSII (Y(II)) of Prochlorococcus sp. strains SS120, PCC 9511, GP2, NATL1A, NATL2A, MIT9313, and TAK9803-2 cultures 24 h after the supplementation with glucose, and compared to the respective control samples. Y(II) was measured in a range of irradiancies (from 0 to 400 μmol quanta m-2 s-1) (not shown). These studies did not show any significant change in Y(II) when the cultures were grown in the presence of glucose, suggesting that Prochlorococcus retains the functional integrity of its photosynthetic apparatus in the presence of glucose. Nevertheless, in several cyanobacterial species, the heterotrophic growth is known to elicit a number of significant changes in both the structural and functional organization of the photosynthetic apparatus (Evans and Carr, 1975; Smart and Mcintosh, 1991; Vernotte et al., 1992; Mannan and Pakrasi, 1993; Kurian et al., 2006). However, the present study is not the first showing that a cyanobacterium maintains its photosynthetic efficiency during the heterotrophic growth (Raboy et al., 1976).

The presence of glucose might be helpful for Prochlorococcus under specific conditions in natural environments. For instance, a recent report has shown that addition of both glucose and pyruvate to Prochlorococcus cultures extends their ability to survive prolonged darkness (Coe et al., 2016), although glucose uptake is lower in darkness than in the light (Gómez-Baena et al., 2008). This suggests that other metabolites, such as amino acids (Zubkov et al., 2004), should be considered in addition to glucose to fully understand its physiological effects in Prochlorococcus.

Effects of Glucose on the Proteome of Prochlorococcus sp. Strain SS120

The effect in the proteome profile of 0.1 μM glucose addition to Prochlorococcus cultures was measured after 24 h. This time was selected in order to allow the proteome to adapt to the new conditions of glucose availability.

A total of 722 proteins were identified: 601 common to glucose and control conditions, 50 proteins were identified uniquely in the control, and a total of 71 proteins uniquely in the presence of glucose [applying false discovery rate (FDR) < 1% and a least two peptides per identification; Figure 5 and Supplementary Table S2]. Quantification data was confidently obtained for 702 proteins (Supplementary Table S3).


FIGURE 5. Effect of glucose addition on the proteome of Prochlorococcus sp. SS120: summary of the proteomic analysis. (A) Overlapping average quantification data obtained for each protein in control and glucose treated cultures. (B) Venn diagram showing the number of proteins identified in control and glucose treated cultures. (C) Correlation between the average quantification obtained per protein for control and glucose treated cultures. (D) Heat map of protein expression (log 10) and hierarchical clustering analysis: C1, control condition from biological replicate 1; G1, glucose condition from biological replicate 1; C2, control condition from biological replicate 2; G2, glucose condition from biological replicate 2; C3, control condition from biological replicate 3; G3, glucose condition from biological replicate 3.

The addition of 0.1 μM glucose did not promote big changes in the proteome of Prochlorococcus (Figure 6 and Supplementary Figure S2), in good agreement with our gene expression results. Most of the proteins remained unchanged under glucose addition. When clustering analysis was performed it was evident that samples from the same biological replicate are grouped together, supporting the fact that no major changes are promoted after glucose addition (Figure 5D).


FIGURE 6. Changes in the concentration of enzymes involved in glucose utilization in Prochlorococcus SS120 upon addition of 0.1 μM glucose. The scheme depicts the main pathways of glucose utilization in Prochlorococcus SS120 using the KEGG database, showing the detected changes in concentration for the different enzymes. Colors show the fold change (glucose/control) obtained from relative quantitation. Red corresponds to upregulation, green to downregulation.

Figure 6 outlines the pathways of glucose metabolization in Prochlorococcus SS120 with the main observed changes in enzyme concentrations. Glycolysis is not active due to the absence of phosphofructokinase. Therefore, glucose can be metabolized either by the Entner-Doudoroff or by the pentoses phosphate pathways. We observed an increase in a number of enzymes belonging to both pathways: in the Entner-Doudoroff pathway, several enzymes showed increases of around 20% while in the pentoses phosphate pathway, three of them showed increases of 30–40%. Namely, 2-keto-3-deoxy-6-phosphogluconate aldolase (which is the key enzyme of the Entner-Doudoroff pathway) increased 20%; on the other hand 6-phosphogluconate dehydrogenase (catalyzing the first specific step in the pentose phosphate pathway) increased 32%. Although these changes might not seem dramatic, it has been described in Prochlorococcus that small changes in protein level can drive significant redirection in metabolic fluxes (Waldbauer et al., 2012).

A recent study proposed that Prochlorococcus might use the Entner-Doudoroff pathway to metabolize glucose under mixotrophic conditions (Chen et al., 2016). Our results provide support for this hypothesis (Figure 6), but also indicate an increase in some enzymes related to the pentose phosphate pathway, in good agreement with our results on gene expression (Figure 4). Our results suggest that availability of glucose induces a metabolic redirection of fluxes in Prochlorococcus, which would involve obtaining reducing power via the pentoses phosphate pathway, and the synthesis of ATP via the Entner-Doudoroff pathway.

Network analysis was performed using STRING10 on other proteins significantly changed after glucose addition (Supplementary Figure S3). The results emphasized the connection between three of them related to ribosomal machinery. There was no evidence for a decrease in the concentration of proteins related to CO2 fixation (Supplementary Table S3).

These results are consistent with the view that Prochlorococcus uses glucose as an extra source of energy and carbon but does not repress photosynthetic CO2 fixation. This is in contrast with the results shown in Synechocystis PCC 6803 (Kurian et al., 2006) in which the photosynthetic capacity in the presence of glucose was markedly reduced, suggesting a significant down-regulation of CO2-fixation during heterotrophic growth. However, it should be noted that the glucose concentration used in the latter study (5 mM) was much higher than that used in our experiments, and cells were kept in darkness.

Several photosynthetic proteins (as Psb27) showed an increase in concentration after 24 h of glucose addition (Supplementary Figure S4 and Table S3). The increase in proteins involved in photosynthesis might be related to the energy demand for the functioning of the transporter GlcH, which is ultimately provided by photosynthesis. This supports the hypothesis that Prochlorococcus is essentially a phototroph organism that uses glucose as an additional source of energy and carbon. This is consistent with the the normal functioning of the photosynthetic apparatus. The regulatory proteins NtcA and PII did not show a significant increase after 24 h in the presence of glucose, in agreement with the expression results previously described (Gómez-Baena et al., 2008). This might be due to a rapid metabolization of glucose taken up in the cell, preventing internal shifts in C/N ratios.

There is in general a good correlation if the data derived from proteomics and mRNA quantification are compared (Figures 4, 7). For example, in the case of the enzymes phosphoribulokinase and transaldolase, the corresponding gene expressions increased after 5 and 8 h after glucose addition, going back to basal levels after 24 h. Proteomics results are consistent with these data: higher average concentration in glucose treated cultures at 24 h, although this difference was not significant (p-value > 0.05).


FIGURE 7. Effect of glucose addition on the proteome of Prochlorococcus sp. SS120: effect on the abundance of proteins related to glucose metabolism. Absolute quantification data obtained for proteins related to glucose metabolism: orange, control condition, green, glucose condition. Corresponding p-value is shown in each panel.


In this work we have explored different aspects of glucose uptake in Prochlorococccus, such as its kinetic and phylogenetic diversity, the transport mechanism and the effects on the cell metabolism. Our results showed that a significant level of variation exists in the kinetics of glucose uptake in marine cyanobacteria, therefore suggesting that the capability for efficient glucose uptake might be an important feature conferring selective advantages in the oligotrophic oceans. This advantage might be possessed by select members of broadly defined ecotypic groups. The features of the transport mechanism indicate a primary active transport, thus involving the consumption of metabolic energy in order to incorporate glucose. Finally, gene expression and proteomic studies point out to a redirection in some metabolic fluxes, while photosynthetic efficiency remained mostly unchanged. This suggests that Prochlorococcus is actively taking up glucose, when available in the environment, as a way of saving metabolic resources to synthesize sugars; but it does not alter its main role as an organism fundamentally autotrophic.

Author Contributions

MM-M and GG-B performed research; MM-M, GG-B, JD, RB, DG-B, MZ, and JG-F designed research, analyzed data, wrote and approved the final version of the manuscript.


This work was supported by grants BFU-2009-08008/BMC and BFU-20103-44767-P (Spanish Ministerio de Economía y Competitividad, cofunded by the European Social Fund from the European Union), P12-BIO-2141 (Proyectos de Excelencia, Junta de Andalucía, Spain) and Universidad de Córdoba (Programa Propio de Investigación). MM-M received a grant from projects BFU-2009-08008/BMC and P07-CVI-3055 (Proyectos de Excelencia, Junta de Andalucía, Spain). GG-B received a post-doctoral Marie Curie grant from the European Union. DG-B received a Ramón y Cajal grant (RYC-2011-07671) from the Spanish Ministerio de Economía y Competitividad.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


We thank the Roscoff Culture Collection (Station Biologique, Roscoff, France; ASSEMBLE project, grant agreement no. 227799, from the “Capacities” program, 7th Framework Program, European Union), the Pasteur Culture Collection (Institut Pasteur, Paris, France) and the MIT Culture Collection (Masachussets Instute of Technology, USA) for providing cyanobacterial cultures. We also thank Dr. G. Tarran (principal scientist of the cruise D371-AMT21), as well as the crew of the R.R.S. Discovery for their assistance during the AMT21 cruise, and Prof. S. García Navarro (Departamento Biología Celular, Universidad de Córdoba) for kindly providing access to radioactivity facilities. We acknowledge the kind collaboration of the Instituto Español de Oceanografía, the Odón de Buen ship crew, for supplying the seawater. We are grateful to Dr. Philip Brownridge for excellent instrumentation support within the Centre for Proteome Research (University of Liverpool, UK) and Dr. Charles Bachy (Monterey Bay Aquarium Research Institute, USA) for his help with the phylogenetic analysis. This study is a contribution to the international IMBER project and was supported by the UK Natural Environment Research Council National Capability funding to Plymouth Marine Laboratory and the National Oceanography Centre, Southampton. This is contribution number 310 of the AMT programme. Finally, we thank Dr. Ignacio Luque (Instituto de Bioquímica Vegetal y Fotosíntesis, Sevilla, Spain) for critical reading of the manuscript and Dr. Enrique Flores (Instituto de Bioquímica Vegetal y Fotosíntesis, Sevilla, Spain) for insightful discussions.

Supplementary Material

The Supplementary Material for this article can be found online at:

FIGURE S1 | Maximum likelihood phylogenetic tree of glcH based on protein sequence similarity with the Prochlorococcus, Synechococcus and GOS sequences. The tree includes 55 sequences (shown by their NCBI accesion numbers in the Supplementary Table S1) of Prochlorococcus and Synechococcus strains and 30 GOS sequences. Bootstrap values (100 replicates) above 75% are indicated on the branches. The sequences corresponding to Synechococcus sp. PCC 7002, Microcystis aeruginosa, Synechocistis sp. PCC 6803 and Thermosynechococcus elongates BP-1 acted as outgroups in this tree.

FIGURE S2 | Correlation plot between the relative protein abundance (log 10) obtained from individual samples. Bottom panels show the correlation scatter plot and top panel show the correlation coefficient.

FIGURE S3 | Network illustrating the relationship between proteins significantly changed after glucose addition. Proteomics was used to evaluate the effect of glucose addition on the proteome of Prochlorococcus SS120. Analysis performed using String ( requiring low confidence (0.150) and adding two (white) nodes. Different colors represent different types of evidence for association: (green) Neighborhood, (red) Gene fusion, (navy blue) Coocurrence, (black) Coexpression, (pink) Experiments, (light blue) Database.

FIGURE S4 | Effect of glucose addition on the proteome of Prochlorococcus sp. SS120: proteins significantly changed after glucose addition. Each panel shows the absolute quantification data obtained for proteins significantly changed (p-value < 0.05): orange, control condition, green, glucose condition. Corresponding p-value is shown in each panel.


  1. ^
  2. ^
  3. ^
  4. ^
  5. ^
  6. ^
  7. ^
  8. ^
  9. ^
  10. ^


Azam, F., and Hodson, R. (1981). Multiphasic kinetics for D-glucose uptake by assemblages of natural marine-bacteria. Mar. Ecol. Prog. Ser. 6, 213–222. doi: 10.3354/meps006213

CrossRef Full Text | Google Scholar

Bailey, S., Melis, A., Mackey, K. R., Cardol, P., Finazzi, G., Van Dijken, G., et al. (2008). Alternative photosynthetic electron flow to oxygen in marine Synechococcus. Biochim. Biophys. Acta 1777, 269–276. doi: 10.1016/j.bbabio.2008.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Biller, S., Berube, P., Berta-Thompson, J., Kelly, L., Roggensack, S., Awad, L., et al. (2014). Genomes of diverse isolates of the marine cyanobacterium Prochlorococcus. Sci. Data 1:140034. doi: 10.1038/sdata.2014.34

PubMed Abstract | CrossRef Full Text | Google Scholar

Biller, S., Berube, P., Lindell, D., and Chisholm, S. (2015). Prochlorococcus: the structure and function of collective diversity. Nat. Rev. Microbiol. 13, 13–27. doi: 10.1038/nrmicro3378

PubMed Abstract | CrossRef Full Text | Google Scholar

Bjorkman, K. M., Church, M. J., Doggett, J. K., and Karl, D. M. (2015). Differential assimilation of inorganic carbon and leucine by Prochlorococcus in the oligotrophic North Pacific subtropical gyre. Front. Microbiol. 6:1401. doi: 10.3389/fmicb.2015.01401

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradford, M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254. doi: 10.1016/0003-2697(76)90527-3

CrossRef Full Text | Google Scholar

Chen, X., Schreiber, K., Appel, J., Makowka, A., Fahnrich, B., Roettger, M., et al. (2016). The Entner-Doudoroff pathway is an overlooked glycolytic route in cyanobacteria and plants. Proc. Natl. Acad. Sci. U.S.A. 113, 5441–5446. doi: 10.1073/pnas.1521916113

PubMed Abstract | CrossRef Full Text | Google Scholar

Chisholm, S., Olson, R., Zettler, E., Goericke, R., Waterbury, J., and Welschmeyer, N. (1988). A novel free living prochlorophyte abundant in the oceanic euphotic zone. Nature 334, 340–343. doi: 10.1038/334340a0

CrossRef Full Text | Google Scholar

Church, M., Ducklow, H., and Karl, D. M. (2004). Light dependence of [3H]- Leucine incorporation in the oligotrophic North Pacific Ocean. Appl. Environ. Microbiol. 70, 4079–4087. doi: 10.1128/AEM.70.7.4079-4087.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Church, M., Ducklow, H., Letelier, R., and Karl, D. M. (2006). Temporal and vertical dynamics in picoplankton photoheterotrophic production in the subtropical North Pacific ocean. Aquat. Microb. Ecol. 45, 41–53. doi: 10.3354/ame045041

CrossRef Full Text | Google Scholar

Coe, A., Ghizzoni, J., Legault, K., Biller, S., Roggensack, S. E., and Chisholm, S. W. (2016). Survival of Prochlorococcus in extended darkness. Limnol. Oceanogr. 71, 1375–1388. doi: 10.1002/lno.10302

CrossRef Full Text | Google Scholar

Cullen, J. (2015). Subsurface chlorophyll maximum layers: enduring enigma or mystery solved? Ann. Rev. Mar. Sci. 7, 207–239. doi: 10.1146/annurev-marine-010213-135111

PubMed Abstract | CrossRef Full Text | Google Scholar

Domínguez-Martín, M. A., López-Lozano, A., Diez, J., Gómez-Baena, G., Rangel-Zúñiga, O., and García-Fernández, J. M. (2014). Physiological regulation of isocitrate dehydrogenase and the role of 2-oxoglutarate in Prochlorococcus sp. strain PCC 9511. PLOS ONE 9:e103380. doi: 10.1371/journal.pone.0103380

PubMed Abstract | CrossRef Full Text | Google Scholar

Dufresne, A., Salanoubat, M., Partensky, F., Artiguenave, F., Axmann, I., Barbe, V., et al. (2003). Genome sequence of the cyanobacterium Prochlorococcus marinus SS120, a nearly minimal oxyphototrophic genome. Proc. Natl. Acad. Sci. U.S.A. 100, 10020–10025. doi: 10.1073/pnas.1733211100

PubMed Abstract | CrossRef Full Text | Google Scholar

El Alaoui, S., Diez, J., Humanes, L., Toribio, F., Partensky, F., and García-Fernández, J. (2001). In vivo regulation of glutamine synthetase activity in the marine chlorophyll b-containing cyanobacterium Prochlorococcus sp. strain PCC 9511 (Oxyphotobacteria). Appl. Environ. Microbiol. 67, 2202–2207. doi: 10.1128/AEM.67.5.2202-2207.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Ethayathulla, A. S., Yousef, M. S., Amin, A., Leblanc, G., Kaback, H. R., and Guan, L. (2014). Structure-based mechanism for Na(+)/melibiose symport by MelB. Nat. Commun. 5:3009. doi: 10.1038/ncomms4009

PubMed Abstract | CrossRef Full Text | Google Scholar

Evans, E. H., and Carr, N. G. (1975). Dark-light transitions with a heterotrophic culture of a blue-green alga. Biochem. Soc. Trans. 3, 373–376. doi: 10.1042/bst0030373

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, S., and Wang, G. (2012). The enhancement of cyclic electron flow around photosystem I improves the recovery of severely desiccated Porphyra yezoensis (Bangiales, Rhodophyta). J. Exp. Bot. 63, 4349–4358. doi: 10.1093/jxb/ers082

PubMed Abstract | CrossRef Full Text | Google Scholar

Gómez-Baena, G., López-Lozano, A., Gil-Martínez, J., Lucena, J., Diez, J., Candau, P., et al. (2008). Glucose uptake and its effect on gene expression in Prochlorococcus. PLoS ONE 3:e3416. doi: 10.1371/journal.pone.0003416

PubMed Abstract | CrossRef Full Text | Google Scholar

Gómez-Pereira, P. R., Hartmann, M., Grob, C., Tarran, G. A., Martin, A. P., Fuchs, B. M., et al. (2013). Comparable light stimulation of organic nutrient uptake by SAR11 and Prochlorococcus in the North Atlantic subtropical gyre. ISME J. 7, 603–614. doi: 10.1038/ismej.2012.126

PubMed Abstract | CrossRef Full Text | Google Scholar

Hilgenfeld, R., and Saenger, W. (1982). Structural chemistry of natural and synthetic ionophores and their complexes with cations. Top. Curr. Chem. 101, 1–82. doi: 10.1007/3-540-11103-4_6

PubMed Abstract | CrossRef Full Text | Google Scholar

Holtzendorff, J., Marie, D., Post, A., Partensky, F., Rivlin, A., and Hess, W. (2002). Synchronized expression of ftsZ in natural Prochlorococcus populations of the Red Sea. Environ. Microbiol. 4, 644–653. doi: 10.1046/j.1462-2920.2002.00347.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurian, D., Jansen, T., and Maenpaa, P. (2006). Proteomic analysis of heterotrophy in Synechocystis sp. PCC 6803. Proteomics 6, 1483–1494.

PubMed Abstract | Google Scholar

Lindell, D., and Post, A. (2001). Ecological aspects of ntcA gene expression and its use as an indicator of the nitrogen status of marine Synechococcus spp. Appl. Environ. Microbiol. 67, 3340–3349. doi: 10.1128/AEM.67.8.3340-3349.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Mannan, R. M., and Pakrasi, H. B. (1993). Dark heterotrophic growth conditions result in an increase in the content of photosystem II units in the filamentous cyanobacterium Anabaena variabilis ATCC 29413. Plant Physiol. 103, 971–977. doi: 10.1104/pp.103.3.971

PubMed Abstract | CrossRef Full Text | Google Scholar

Martiny, A., Tai, A., Veneziano, D., Primeau, F., and Chisholm, S. (2009). Taxonomic resolution, ecotypes and the biogeography of Prochlorococcus. Environ. Microbiol. 11, 823–832. doi: 10.1111/j.1462-2920.2008.01803.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mary, I., Tarran, G., Warwick, P., Terry, M., Scanlan, D., Burkill, P., et al. (2008). Light enhanced amino acid uptake by dominant bacterioplankton groups in surface waters of the Atlantic Ocean. FEMS Microbiol. Ecol. 63, 36–45. doi: 10.1111/j.1574-6941.2007.00414.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Michelou, V., Cottrell, M., and Kirchman, D. (2007). Light-stimulated bacterial production and amino acid assimilation by cyanobacteria and other microbes in the North Atlantic Ocean. Appl. Environ. Microbiol. 73, 5539–5546. doi: 10.1128/AEM.00212-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, P. (1979). Keilin’s respiratory chain concept and its chemiosmotic consequences. Science 206, 1148–1159. doi: 10.1126/science.388618

CrossRef Full Text | Google Scholar

Moisander, P. H., Zhang, R., Boyle, E. A., Hewson, I., Montoya, J. P., and Zehr, J. P. (2012). Analogous nutrient limitations in unicellular diazotrophs and Prochlorococcus in the South Pacific Ocean. ISME J. 6, 733–744. doi: 10.1038/ismej.2011.152

PubMed Abstract | CrossRef Full Text | Google Scholar

Muñoz-Marín, M. C., Luque, I., Zubkov, M. V., Hill, P. G., Diez, J., and García-Fernández, J. M. (2013). Prochlorococcus can use the Pro1404 transporter to take up glucose at nanomolar concentrations in the Atlantic Ocean. Proc. Natl. Acad. Sci. U.S.A. 110, 8597–8602. doi: 10.1073/pnas.1221775110

PubMed Abstract | CrossRef Full Text | Google Scholar

Nieva, M., and Fernández Valiente, E. (1996). Inorganic carbon transport and fixation in cells of Anabaena variabilis adapted to mixotrophic conditions. Plant Cell Physiol. 37, 1–7. doi: 10.1093/oxfordjournals.pcp.a028905

CrossRef Full Text | Google Scholar

Palinska, K., Jahns, T., Rippka, R., and Tandeau De Marsac, N. (2000). Prochlorococcus marinus strain PCC 9511, a picoplanktonic cyanobacterium, synthesizes the smallest urease. Microbiology 146, 3099–3107. doi: 10.1099/00221287-146-12-3099

PubMed Abstract | CrossRef Full Text | Google Scholar

Partensky, F., and Garczarek, L. (2003). “The photosynthetic apparatus of chlorophyll b- and d-containing oxyphotobacteria,” in Photosynthesis in Algae, eds A. Larkum, S. Douglas, and J. Raven (Dordrecht: Kluwer Academic Publishers), 29–62.

Google Scholar

Partensky, F., and Garczarek, L. (2010). Prochlorococcus: advantages and limits of minimalism. Ann. Rev. Mar. Sci. 2, 305–331. doi: 10.1146/annurev-marine-120308-081034

PubMed Abstract | CrossRef Full Text | Google Scholar

Partensky, F., Hess, W., and Vaulot, D. (1999). Prochlorococcus, a marine photosynthetic prokaryote of global significance. Microbiol. Mol. Biol. Rev. 63, 106–127.

PubMed Abstract | Google Scholar

Paulsen, I., Sliwinski, M., and Saier, M. (1998). Microbial genome analyses: global comparisons of transport capabilities based on phylogenies, bioenergetics and substrate specificities. J. Mol. Biol. 277, 573–592. doi: 10.1006/jmbi.1998.1609

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfaffl, M. (2001). A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29:e45. doi: 10.1093/nar/29.9.e45

CrossRef Full Text | Google Scholar

Picossi, S., Flores, E., and Ekman, M. (2013). Diverse roles of the GlcP glucose permease in free-living and symbiotic cyanobacteria. Plant Signal. Behav. 8:e27416. doi: 10.4161/psb.27416

PubMed Abstract | CrossRef Full Text | Google Scholar

Quick, W., Schurr, U., Fichtner, K., Schulze, E.-D., Rodermel, S., Bogorad, L., et al. (1991). The impact of decreased Rubisco on photosynthesis, growth, allocation and storage in tobacco plants which have been transformed with antisense rbcS. Plant J. 1, 51–58. doi: 10.1111/j.1365-313X.1991.00051.x

CrossRef Full Text | Google Scholar

Raboy, B., Padan, E., and Shilo, M. (1976). Heterotrophic capacities of Plectonema boryanum. Arch. Microbiol. 110, 77–85. doi: 10.1007/BF00416971

PubMed Abstract | CrossRef Full Text | Google Scholar

Rich, J., Ducklow, H., and Kirchman, D. (1996). Concentrations and uptake of neutral monosaccharides along 140°W in the equatorial Pacific: contribution of glucose to heterotrophic bacterial activity and the DOM flux. Limnol. Oceanogr. 41, 595–604. doi: 10.4319/lo.1996.41.4.0595

CrossRef Full Text | Google Scholar

Rich, P., Madgwick, S., and Moss, D. (1991). The interactions of duroquinol, DBMIB and NQNO with the chloroplast cytochrome b6f complex. Biochim. Biophys. Acta 108, 1188–1195.

Google Scholar

Rippka, R., Coursin, T., Hess, W., Lichtlé, C., Scanlan, D., Palinska, K., et al. (2000). Prochlorococcus marinus Chisholm et al. 1992 subsp. pastoris subsp. nov. strain PCC 9511, the first axenic chlorophyll a2/b2-containing cyanobacterium (Oxyphotobacteria). Int. J. Syst. Evol. Microbiol. 50, 1833–1847. doi: 10.1099/00207713-50-5-1833

PubMed Abstract | CrossRef Full Text | Google Scholar

Rocap, G., Larimer, F., Lamerdin, J., Malfatti, S., Chain, P., Ahlgren, N., et al. (2003). Genome divergence in two Prochlorococcus ecotypes reflects oceanic niche differentiation. Nature 424, 1042–1047. doi: 10.1038/nature01947

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubin, P. M., Zetooney, E., and Mcgowan, R. E. (1977). Uptake and utilization of sugar phosphates by Anabaena-Flos-Aquae. Plant Physiol. 60, 407–411. doi: 10.1104/pp.60.3.407

PubMed Abstract | CrossRef Full Text | Google Scholar

Rusch, D., Halpern, A., Sutton, G., Heidelberg, K., Williamson, S., Yooseph, S., et al. (2007). The Sorcerer II Global Ocean Sampling Expedition: northwest atlantic through eastern tropical pacific. PLoS Biol. 5:e77. doi: 10.1371/journal.pbio.0050077

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwalbach, M., Tripp, H., Steindler, L., Smith, D., and Giovannoni, S. (2010). The presence of the glycolysis operon in SAR11 genomes is positively correlated with ocean productivity. Environ. Microbiol. 12, 490–500. doi: 10.1111/j.1462-2920.2009.02092.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shoshan, V., and Selman, B. R. (1980). The relationship between light-induced adenine nucleotide exchange and ATPase activity in chloroplast thylakoid membranes. J. Biol. Chem. 255, 384–389.

PubMed Abstract | Google Scholar

Silva, J. C., Gorenstein, M. V., Li, G. Z., Vissers, J. P., and Geromanos, S. J. (2006). Absolute quantification of proteins by LCMSE: a virtue of parallel MS acquisition. Mol. Cell. Proteomics 5, 144–156. doi: 10.1074/mcp.M500230-MCP200

PubMed Abstract | CrossRef Full Text | Google Scholar

Smart, L., and Mcintosh, L. (1991). Expression of photosynthesis genes in the cyanobacterium Synechocystis sp. PCC 6803: psaA-psaB and psbA transcripts accumulate in dark-grown cells. Plant Mol. Biol. 17, 959–971. doi: 10.1007/BF00037136

PubMed Abstract | CrossRef Full Text | Google Scholar

Smeekens, S., and Rook, F. (1997). Sugar sensing and sugar-mediated signal transduction in plants. Plant Physiol. 115, 7–13. doi: 10.1104/pp.115.1.7

CrossRef Full Text | Google Scholar

Soontharapirakkul, K., Promden, W., Yamada, N., Kageyama, H., Incharoensakdi, A., Iwamoto-Kihara, A., et al. (2011). Halotolerant cyanobacterium Aphanothece halophytica contains an Na+-dependent F1F0-ATP synthase with a potential role in salt-stress tolerance. J. Biol. Chem. 286, 10169–10176. doi: 10.1074/jbc.M110.208892

PubMed Abstract | CrossRef Full Text | Google Scholar

Talarmin, A., Van Wambeke, F., Catala, P., Courties, C., and Lebaron, P. (2011). Flow cytometric assessment of specific leucine incorporation in the open Mediterranean. Biogeosciences 8, 253–265. doi: 10.5194/bg-8-253-2011

CrossRef Full Text | Google Scholar

Thompson, L. R., Zeng, Q., Kelly, L., Huang, K. H., Singer, A. U., Stubbe, J., et al. (2011). Phage auxiliary metabolic genes and the redirection of cyanobacterial host carbon metabolism. Proc. Natl. Acad. Sci. U.S.A. 108, E757–E764. doi: 10.1073/pnas.1102164108

PubMed Abstract | CrossRef Full Text | Google Scholar

Trebst, A. (1980). Inhibitors in the electron flow. Methods Enzymol. 69, 675–715. doi: 10.1016/S0076-6879(80)69067-3

CrossRef Full Text

Vernotte, C., Picaud, M., Kirilovsky, D., Olive, J., Ajlani, G., and Astier, C. (1992). Changes in the photosynthetic apparatus in the cyanobacterium Synechocystis sp. PCC 6714 following light-to-dark and dark-to-light transitions. Photosynth. Res. 32, 45–57. doi: 10.1007/BF00028797

PubMed Abstract | CrossRef Full Text | Google Scholar

Vila-Costa, M., Simo, R., Harada, H., Gasol, J., Slezak, D., and Kiene, R. (2006). Dimethylsulfoniopropionate uptake by marine phytoplankton. Science 314, 652–654. doi: 10.1126/science.1131043

PubMed Abstract | CrossRef Full Text | Google Scholar

Waldbauer, J. R., Rodrigue, S., Coleman, M. L., and Chisholm, S. W. (2012). Transcriptome and proteome dynamics of a light-dark synchronized bacterial cell cycle. PLoS ONE 7:e43432. doi: 10.1371/journal.pone.0043432

PubMed Abstract | CrossRef Full Text | Google Scholar

West, N., Schonhuber, W., Fuller, N., Amann, R., Rippka, R., Post, A., et al. (2001). Closely related Prochlorococcus genotypes show remarkably different depth distributions in two oceanic regions as revealed by in situ hybridization using 16S rRNA-targeted oligonucleotides. Microbiology 147, 1731–1744. doi: 10.1099/00221287-147-7-1731

PubMed Abstract | CrossRef Full Text | Google Scholar

Williamson, S. J., Rusch, D. B., Yooseph, S., Halpern, A. L., Heidelberg, K. B., Glass, J. I., et al. (2008). The Sorcerer II Global Ocean Sampling Expedition: metagenomic characterization of viruses within aquatic microbial samples. PLoS ONE 3:e1456. doi: 10.1371/journal.pone.0001456

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, R. T., and Hobbie, J. E. (1966). Use of glucose and acetate by bacteria in aquatic ecosystems. Ecology 47, 447–464. doi: 10.2307/1932984

CrossRef Full Text | Google Scholar

Yang, C., Hua, Q., and Shimizu, K. (2002). Integration of the information from gene expression and metabolic fluxes for the analysis of the regulatory mechanisms in Synechocystis. Appl. Microbiol. Biotechnol. 58, 813–822. doi: 10.1007/s00253-002-0949-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Yelton, A. P., Acinas, S. G., Sunagawa, S., Bork, P., Pedros-Alio, C., and Chisholm, S. W. (2016). Global genetic capacity for mixotrophy in marine picocyanobacteria. ISME J. 10, 2946–2957. doi: 10.1038/ismej.2016.64

PubMed Abstract | CrossRef Full Text | Google Scholar

Yooseph, S., Sutton, G., Rusch, D., Halpern, A., Williamson, S., Remington, K., et al. (2007). The Sorcerer II Global Ocean Sampling Expedition: expanding the universe of protein families. PLoS Biol. 5:e16. doi: 10.1371/journal.pbio.0050016

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshida, M., Poser, J., Allison, W., and Esch, F. (1981). Identification of an essential glutamic acid residue in the b subunit of the adenosine triphosphatase from the thermophilic bacterium PS3. J. Biol. Chem. 256, 148–153.

Google Scholar

Zorz, J. K., Allanach, J. R., Murphy, C. D., Roodvoets, M. S., Campbell, D. A., and Cockshutt, A. M. (2015). The RUBISCO to photosystem II ratio limits the maximum photosynthetic rate in picocyanobacteria. Life 5, 403–417. doi: 10.3390/life5010403

PubMed Abstract | CrossRef Full Text | Google Scholar

Zubkov, M., Fuchs, B., Tarran, G., Burkill, P., and Amann, R. (2003). High rate of uptake of organic nitrogen compounds by Prochlorococcus cyanobacteria as a key to their dominance in oligotrophic oceanic waters. Appl. Environ. Microbiol. 69, 1299–1304. doi: 10.1128/AEM.69.2.1299-1304.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zubkov, M., Tarran, G., and Fuchs, B. (2004). Depth related amino acid uptake by Prochlorococcus cyanobacteria in the Southern Athlantic tropical gyre. FEMS Microbiol. Ecol. 50, 153–161. doi: 10.1016/j.femsec.2004.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Zubkov, M., Tarran, G., Mary, I., and Fuchs, B. (2008). Differential microbial uptake of dissolved amino acids and amino sugars in surface waters of the Atlantic Ocean. J. Plankton Res. 30, 211–220. doi: 10.1093/plankt/fbm091

CrossRef Full Text | Google Scholar

Zubkov, M. V. (2009). Photoheterotrophy in marine prokaryotes. J. Plankton Res. 31, 933–938. doi: 10.1093/plankt/fbp043

CrossRef Full Text | Google Scholar

Keywords: marine cyanobacteria, mixotrophy, transporter kinetics, proteomics, gene expression

Citation: Muñoz-Marín MC, Gómez-Baena G, Díez J, Beynon RJ, González-Ballester D, Zubkov MV and García-Fernández JM (2017) Glucose Uptake in Prochlorococcus: Diversity of Kinetics and Effects on the Metabolism. Front. Microbiol. 8:327. doi: 10.3389/fmicb.2017.00327

Received: 16 December 2016; Accepted: 16 February 2017;
Published: 08 March 2017.

Edited by:

Hongbin Liu, The Hong Kong University of Science and Technology, Hong Kong

Reviewed by:

Ryan Paerl, University of Copenhagen, Denmark
Douglas Andrew Campbell, Mount Allison University, Canada
Karin M. Björkman, Hawaii University, USA

Copyright © 2017 Muñoz-Marín, Gómez-Baena, Díez, Beynon, González-Ballester, Zubkov and García-Fernández. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: José M. García-Fernández,

These authors have contributed equally to this work.