Parallel subfunctionalisation of PsbO protein isoforms in angiosperms revealed by phylogenetic analysis and mapping of sequence variability onto protein structure
© Duchoslav and Fischer. 2015
Received: 12 February 2015
Accepted: 11 May 2015
Published: 9 June 2015
PsbO, the manganese-stabilising protein, is an indispensable extrinsic subunit of photosystem II. It plays a crucial role in the stabilisation of the water-splitting Mn4CaO5 cluster, which catalyses the oxidation of water to molecular oxygen by using light energy. PsbO was also demonstrated to have a weak GTPase activity that could be involved in regulation of D1 protein turnover. Our analysis of psbO sequences showed that many angiosperm species express two psbO paralogs, but the pairs of isoforms in one species were not orthologous to pairs of isoforms in distant species.
Phylogenetic analysis of 91 psbO sequences from 49 land plant species revealed that psbO duplication occurred many times independently, generally at the roots of modern angiosperm families. In spite of this, the level of isoform divergence was similar in different species. Moreover, mapping of the differences on the protein tertiary structure showed that the isoforms in individual species differ from each other on similar positions, mostly on the luminally exposed end of the β-barrel structure. Comparison of these differences with the location of differences between PsbOs from diverse angiosperm families indicated various selection pressures in PsbO evolution and potential interaction surfaces on the PsbO structure.
The analyses suggest that similar subfunctionalisation of PsbO isoforms occurred parallelly in various lineages. We speculate that the presence of two PsbO isoforms helps the plants to finely adjust the photosynthetic apparatus in response to variable conditions. This might be mediated by diverse GTPase activity, since the isoform differences predominate near the predicted GTP-binding site.
KeywordsGene duplication GTPase Homology modelling Manganese-stabilizing protein (MSP) Oxygen evolving complex Parallel evolution Protein structure PsbO
Photosynthetic conversion of light into chemical energy in oxygenic phototrophs is accompanied with evolution of molecular oxygen released from water molecules. This process is realized in the oxygen evolving complex of photosystem II present in thylakoid membranes. Photosystem II (PSII) is a multisubunit protein–cofactor complex that uses light energy to oxidize water and to reduce plastoquinone. PsbO, also known as the manganese-stabilising protein, is one of the extrinsic subunits of photosystem II, located on the luminal side of the thylakoid membrane. PsbO is present in all known oxygenic photosynthetic organisms . Despite the ability of the cyanobacterium Synechocystis sp. PCC 6803 mutant to grow photoautotrophically with deleted psbO gene , PsbO seems to be crucial for PSII function. Neither the mutant of green alga Chlamydomonas reinhardtii lacking PsbO, nor Arabidopsis thaliana (A. thaliana) with silenced expression of both psbO paralogs were able to grow photoautotrophically or even assemble PSII [3, 4].
Three-dimensional structure of PsbO from cyanobacterium Thermosynechococcus was resolved as a part of PSII by X-ray crystallography with a resolution down to 1.9 Å . The crystal structure of PSII or PsbO alone from plants or other eukaryotes is not available. Some information about the structure of the whole PSII dimer surrounded by antenna complexes (the PSII-LHCII supercomplex) from higher plants was obtained by single particle cryo-electron microscopy and cryo-electron tomography [6–8]. Unfortunately, the resolution is insufficient to provide any plant-specific knowledge about the PsbO structure. Still, relatively high pairwise identity between PsbO sequences of Thermosynechococcus and higher plants (around 45 %) allows construction of homologous models for plant PsbOs [9, 10].
The X-ray crystallography of cyanobacterial PSII revealed that PsbO is a β-barrel protein (structural features of PsbOs are discussed in connection with our results and PsbO functions in chapter Discussion). It is located in the vicinity of the water splitting Mn4CaO5 cluster, but it is not directly involved in binding of the cluster . The main function of the PsbO is to stabilise the Mn4CaO5 cluster, in particular to modulate the calcium and chloride requirements for efficient water splitting (for review see ). Besides this “basic” function, PsbO seems to be involved also in other processes (for review see [12, 13]). Spinach PsbO was shown to be able to bind GTP  and also to hydrolyse it, although very slowly . It was proposed that the GTPase activity of PsbO in plants might be involved in D1 repair cycle .
In plants and algae, the PsbO protein is encoded by a nuclear psbO gene . Transport to chloroplasts and thylakoids is ensured by two consecutive N-terminal transit peptides, that are cleaved to produce the mature PsbO . A. thaliana expresses two psbO genes, psbO1 [TAIR:At5g66570] and psbO2 [TAIR:At3g50820], encoding for PsbO1 and PsbO2 proteins [16, 17]. The two isoforms differ in only 11 amino acids ; nevertheless, their function seems to be slightly different. Murakami et al.  reported that A. thaliana PsbO2 recovered oxygen evolution of PsbO-depleted spinach PSII particles less efficiently than PsbO1. The activity with PsbO2 reached only 80 % of that with PsbO1, while the binding efficiency of the isoforms was very similar. In contrast, the oxygen evolution of PSII membranes isolated from A. thaliana mutants lacking PsbO1 or PsbO2 was similar when corrected for the amount of PSII .
The amount of PsbO1 in wild-type A. thaliana plants is higher than that of PsbO2 [18–20]. The expression of the isoforms stays similar during plant development and during various short time stresses . Only after 40 days of cold stress, noticeable change in relative abundance of isoforms was observed in favour of PsbO2 .
In A. thaliana mutants with an impaired psbO1 or psbO2 gene, the compensatory upregulation of the remaining isoform was observed. The expression level of PsbO2 in psbo1 mutant was increased several times, reaching 75 % of the total amount of PsbO in wild-type. The expression level of PsbO1 in psbo2 mutant was 125 % of the total PsbO in wild-type. The amount of other PSII proteins was affected similarly, leading to the same stoichiometry of PsbO per PSII as in wild type .
The psbo1 mutant plants have pale green leaves, reduced rosette size and slower growth rate as compared to wild-type plants [17, 19, 23]. Descriptions of the psbo2 mutant phenotype slightly differ from each other, probably because of different growth conditions and age of used plants . Lundin et al.  observed growth rate slower than in wild-type and the leaf weight was even lower than that of psbo1, while Allahverdiyeva et al.  reported a phenotype very similar to that of wild-type.
Under growth light (120 μmol photons m−2 s−1), the psbo2 mutant had characteristics of electron transport chain very similar to wild-type, whereas investigation of the psbo1 mutant showed malfunction of both the donor and acceptor sides of PSII and high sensitivity of PSII centres to photodamage . Bricker and Frankel  reported that many of the defects of psbo1 photosystems are reverted by higher concentration of CaCl2, but Allahverdiyeva et al.  did not observe similar effect. Nevertheless, the importance of the PsbO2 seems to be exhibited under high light conditions. For example, the maximum quantum efficiency (FV/FM) values of wild-type and mutant plants became similar after 3 weeks of moderate light (500 μmol photons m−2 s−1) . Lundin et al.  reported that after 15 days of high light (1000 μmol photons m−2 s−1), the psbo1 mutant did not have significantly reduced leaf weight, whereas the leaf weight of psbo2 mutant was reduced drastically.
Lundin et al.  also showed that psbo2 mutant has lower level of phosphorylation of D1 and D2 subunits and that the degradation of photo-damaged D1 protein is impaired in this mutant. This, together with a finding that PSII membranes with PsbO2 have higher GTPase activity than PSII membranes with PsbO1 , led to a conclusion, that PsbO1 has a main function in the stabilisation of Mn4CaO5 cluster and the facilitation of the water oxidation reaction, whereas PsbO2 regulates the turnover of D1 subunit [19, 21, 23].
The presence of two PsbO isoforms is not unique for A. thaliana. Our previous study focused on the analysis of a spontaneously tuberising potato mutant revealed that potato plants also express two PsbO isoforms, one of which is missing in the mutant . A comparison of the two characterised A. thaliana and two potato PsbO isoforms showed that sequences of the two paralogs in each species are more related than isoforms coming from different species. It indicated independent duplication of psbO gene in these two species. To understand this unexpected phylogeny and evolution of PsbO isoforms, we did a detailed analysis of psbO sequences from a number of land plant species. Mapping the sequence differences between PsbO proteins from various species and families and between PsbO isoforms in individual species on their tertiary structure, we found that the evolution of the two isoforms was parallel in numerous angiosperm lineages. Based on the location of isoform-specific differences and literature data about A. thaliana and spinach PsbOs, we hypothesise that the pairs of isoforms present in many species differ in GTPase activity and that the presence of proteins diversified in this way helps to improve photosynthetic performance under varying conditions.
Materials and methods
Retrieval and analysis of psbO sequences
Sequences of expressed psbO genes were retrieved as ESTs (expressed sequence tags) and assembled ESTs (PUTs, PlantGDB-assembled unique transcripts) in public sequence databases NCBI GenBank  and PlantGDB (Plant Genome Database) , respectively. The database searches were performed using tBLASTn [28, 29] with potato PsbO protein sequence (sequence “Solanum tuberosum 2”, translation of [PlantGDB:PUT-157a-Solanum_tuberosum-55973153]) as a query. ESTs were aligned into contigs for each species using “De Novo Assemble” tool of Geneious R6 . Formation of consensus sequences from multiple overlapping ESTs strongly increased reliability of analysed sequences compared to individually submitted annotated cDNAs, some of which contain evident errors. All retrieved sequences were aligned using MAFFT v7.017  and incomplete and unreliable sequences were excluded from further analyses (see analysed sequences in Additional file 1). Spinach psbO sequence was retrieved as cDNA [GenBank:X05548.1] because of the lack of ESTs and included in alignment for comparison (Additional file 2). Indexing of isoforms in each family was random and does not reflect relation to A. thaliana isoforms.
Phylogenetic trees were built from psbO coding sequences by maximum likelihood (ML) method using CIPRES Science Gateway . ML analysis was implemented in tool RAxML v7.6.6  using GTRGAMMA approximation with 1000 bootstrap replicates.
The presence and position of introns was analysed by comparing psbO cDNAs (Additional file 1) and corresponding genomic sequences, obtained using BLASTn [28, 29] searches in Phytozome database  for the following representative species with easily available genomic sequence: Arabidopsis lyrata, Arabidopsis thaliana, Brassica rapa, and Thellungiella halophila from Brassicaceae family and Oryza sativa, Physcomitrella patens, Populus trichocarpa, Solanum lycopersicum, and Vitis vinifera from other families.
Evaluation of PsbO sequence variability
The frequency of differences between isoforms, between species and between families were calculated for each position in the alignment independently using scripts written in R language  and partially using SeqinR package . Plant families represented with just a single PsbO sequence were not included in the calculation. Only two most divergent isoforms were considered in case of species expressing more than two isoforms. All sequences excluded from calculation are marked with an asterisk in Additional file 2. To estimate the between-isoform and between-species variability across all angiosperms, both types of differences were first calculated for every family independently and afterwards the values were averaged, in order to avoid bias caused by different numbers of analysed species within each family.
where n is the number of compared species, ∆ species is the non-linear average value of between-species variability (the mean from pair wise comparisons) and ∆ species linear is the linearised value of the between-species variability.
To estimate the between-family variability, the above mentioned method for the calculation of the between-species differences was applied on sets containing sequences from just one species from each family. A mean values obtained from all such combinations of species (53,760 in total) included both between-species and between-family differences, so the values of between-species differences were subtracted from it, giving the net between-family differences.
Homology modelling and mapping of variability on the protein structure
Homology model of potato PsbO (sequence “Solanum tuberosum 2”) was built using Swiss-Model server [37, 38] based on PsbO from cyanobacterium Thermosynechococcus vulcanus [PDB:3ARC] (chain O) . Extra 13 amino acids present on the N-terminus of potato PsbO were pasted to the model manually using Swiss-PdbViewer v4.1.0  without attempt to show any folding.
Homology model of potato PsbO was coloured according to the frequency of the respective type of variability using Swiss-PdbViewer v4.1.0  and scripts written in R language . The images were rendered using POV-Ray v3.6 .
Determination of spatial centres of differences
Spatial centres of the differences were calculated using coordinates of α-carbon atoms of amino acids in the PsbO homology model using scripts written in R language . The arithmetic mean of the coordinates was weighed by frequency of the respective difference on each position. The 13 N-terminal amino acids with unknown folding were excluded from the calculation. Overall spatial centres of the differences between isoforms and the differences between species in angiosperms were calculated as an arithmetic mean of spatial centres calculated for all families. The statistical significance of the divergence in the location of the spatial centres of the between-isoform and the between-species differences was assessed using a randomisation test. Variable positions in the alignment were randomly shuffled and the spatial centres for the between-isoform and the between-species variability were calculated. Difference between means of the two types of spatial centres projected on the axis of highest variability was compared with the value obtained for real alignment. The p-value was calculated from 50,000 randomisations.
The majority of angiosperm species express two psbO genes
Searching public databases for expressed sequences of psbO genes from land plants (Embryophyta) we obtained 91 sequences from 49 species and 36 genera. Analysis of these sequences showed that the majority of the analysed angiosperm species express more than one, in most cases two psbO isoforms (Additional file 3). In contrast, all analysed representatives of gymnosperms (from both Cycadophyta and Coniferophyta groups) seem to express only one psbO isoform.
In monocots, psbO sequences were available from only two families: Zingiberaceae species have two psbO isoforms, whereas most Poaceae species with available ESTs express only one psbO gene. A single psbO gene was found also in the genomic sequence of Oryza sativa. Zea mays, a recent tetraploid, expresses two isoforms with little divergence (Additional file 3).
Among dicots, Malvaceae, Myrtaceae, Phrymaceae and Rutaceae seem to express only one psbO gene. Asteraceae, Euphorbiaceae, Fabaceae, Salicaceae, Solanaceae and Vitaceae seem to express two psbO genes (or four in the case of recent tetraploids such as Glycine max or Nicotiana tabacum). Brassicaceae have various numbers of psbO isoforms; however, most of them can be sorted into two groups. While Arabidospsis thaliana expresses just two isoforms (psbO1, psbO2), each from one group, genus Brassica expresses three to five genes - one gene corresponds to psbO2 of A. thaliana, while the gene orthologous to psbO1 of A. thaliana is present in several very similar sub-isoforms (4 in B. napus, 3 in B. rapa and 2 in B. oleracea; Additional files 3 and 4). Thellungiella halophila expresses three psbO genes, two of which correspond to psbO1 and psbO2 of A. thaliana, the third one is most similar to pseudogenes that can be found in genomic sequences of A. thaliana [TAIR:At4g37230], Arabidopsis lyrata [GenBank:XM_002866937] and Brassica rapa [Phytozome:Bra017790] (data not shown).
Pairs of PsbO isoforms evolved in every angiosperm family independently
The basic topology of the phylogenetic tree does not contain dichotomous branching to two groups of functionally diverged orthologs at the tree base, but it reflects basic phylogeny of land plant families. The branching to two isoforms is also absent at the base of angiosperms. Instead, the branching events are clearly present at the bases of several families (for example Solanaceae, Fabaceae, Brassicaceae, Zingiberaceae; Fig. 1). This unexpected topology indicates that duplications of psbO gene occurred independently in each plant family that contains species with multiple PsbO isoforms. Moreover, these families do not form any cluster in the phylogenetic tree of psbO or in the consensual phylogeny of angiosperms.
To further confirm the independent duplication of psbO genes in ancestor of each angiosperm family, the presence and position of introns was analysed in available genomic sequences of psbO genes. According to this analysis, all land plants have an intron at a conserved site, 12 nucleotides upstream the boundary between sequences encoding the transit peptide and the mature protein. In addition, all psbO genes from Brassicaceae family contain an additional intron, 282 nucleotides downstream the boundary between the transit peptide and the mature protein. The intron is present at a conserved site in all psbO genes in this family, including the most divergent isoform of Thellungiella halophila. This indicates that all these psbO genes evolved from one common Brassicaceae-specific ancestor gene containing the additional intron, absent in psbOs in other families.
Extent of divergence of PsbO isoforms is similar in all species
The extent of differences between protein sequences of PsbO isoforms in every species is in the same range, even though the duplication seems to have occurred in each family independently. The numbers of different amino acid residues range from six in a recent tetraploid Zea mays to 23 in Populus deltoides and Populus x canadensis (2–9 % of total residues; Additional file 3). Interestingly, similar divergence between isoforms can be found also in the moss Physcomitrella patens (24 different amino acid residues between two most divergent isoforms).
The level of differences between PsbO isoforms is kept within this range even if the duplication events of psbO occurred at different times in evolutionary history. For instance, pairwise identity of nucleotide sequences encoding mature PsbOs of V. vinifera (80 %) is much lower than that of Populus trichocarpa (92 %). This indicates that the duplication of the Vitis psbO gene probably occurred earlier compared to that of the Populus gene. However, pairwise identity of the protein sequences of PsbO isoforms of V. vinifera (93 %) is similar to that of P. trichocarpa (92 %).
Three classes of PsbO sequence variability
Considering that many angiosperm species express two isoforms of psbO, we asked whether the differences between the isoforms are similar in multiple families despite the independent duplications of psbO genes. Detailed analysis of the sequence alignment failed to identify any compact region in the primary sequence that would be specific for one or the other isoform across the analysed plant families. Also, single positions with similar differences between isoforms in the majority of species were rare (see the alignment in Additional file 2).
Amino acid residues varying between isoforms differ predominantly in the length of side chains
Exchanges in other amino acid residues were less conserved among various families. But generally, substitutions between residues, which differed only in the length of the side chain and had similar physicochemical properties, predominated over substitutions between residues with more variable character. The three most frequent amino acid substitutions (D-E, I-V and S-T; Fig. 3) match these criteria and comprise together almost 50 % of all exchanges. Though seemingly synonymous, these substitutions are strongly conserved in orthologous isoforms within families and in some cases even shared across more families (see the alignment in Additional file 2).
Residues varying between isoforms cluster together on the tertiary structure of PsbO
The positions with amino acids varying predominantly between isoforms did not cluster together in the primary sequence. As protein function is tightly connected with tertiary structure, we decided to analyse spatial location of amino acid substitutions between PsbO isoforms on the protein structure. Because no crystal structure of eukaryote PsbO is available, we constructed homologous model of PsbO2 from Solanum tuberosum using PsbO structure from Thermosynechococcus vulcanus  as a template (identity of the protein sequences is 47 %). All PsbO sequences of angiosperms are well comparable on a single model of structure thanks to a very high conservation of both the amino acid sequence and the length of the chain. In the alignment of 78 protein sequences of PsbOs from angiosperms, 59 % of positions are fully identical and the length of the chain of mature proteins varies mostly between 247 and 248 amino acid residues (Additional file 2).
Before drawing any conclusions, we had to prove that this spatial location is specific for differences between isoforms and does not reflect a high level of general variability in these regions. We compared the position of isoform differences with between-species differences in all families (Fig. 4c). We found that differences between species (red-coloured in the figure) are more dispersed over the PsbO structure. To allow statistical analysis, we calculated spatial centres of between-isoform differences and between-species differences (green and red spheres in Fig. 4c, respectively) for each family. The spatial centres of isoform differences are shifted towards the luminal end of the β-barrel (with one exception, the Salicaceae family, which has the centre of differences between isoforms shifted towards the β1-β2 loop due to high frequency of differences in this part of the structure). The shift of spatial centres of isoform differences compared with the centres of between-species differences is significant according to a randomization test (p = 0.002).
PSII-exposed surface is conserved, while differences between families are mainly on the luminal side of the β5-β6 loop
Differences between families are the most frequent class of differences among PsbO sequences (Additional files 2 and 5). Fig. 4d depicts differences between families merged with the differences between isoforms and overall spatial centres of the three classes of differences (represented with green, red and blue spheres). The differences between families are more spread over the PsbO structure compared to the differences between isoforms, similarly to the differences between species within families. The highest frequency of differences between families is in the part of β5-β6 loop that is not interacting with PSII core proteins (the amino acid side chains are pointing towards thylakoid lumen) and the adjoining part of the β6 strand.
Mechanism of duplication and subfunctionalisation of psbO
Several studies demonstrated that A. thaliana expresses two psbO paralogs [17–19, 21]. Here we show that A. thaliana is not an exception and that species from 9 out of 15 investigated angiosperm families also express two distinct psbO genes (Additional file 3). Unexpectedly, the phylogenetic analysis revealed that the psbO gene was not duplicated in the common ancestor of angiosperms, but the duplication occurred many times independently in individual families (Fig. 1).
There are various mechanisms by which gene duplication can occur. In terms of its extent, duplication can involve single genes, larger segments, chromosomes or entire genomes . In A. thaliana and Populus trichocarpa we found that chromosomal segments containing psbO paralogs are collinear (i.e. contain homologous genes in a similar order; Duchoslav, Vosolsobě, and Fischer, unpublished results), which suggests that psbO was duplicated within the context of a larger-scale duplication.
The phylogenetic tree topology indicates that the duplication event occurred in ancestors of numerous families prior to extensive species radiation. The radiation that involved many extant plant lineages in Paleogene, was likely facilitated by the whole genome duplications (WGD) dated to the last global extinction period at the Cretaceous–Paleogene boundary about 66 million years ago [42, 43]. Based on this indirect evidence, we suggest that the psbO duplication was not gene specific, but rather that the paralogs were in many cases retained after WGD events that occurred independently in ancestors of many successful angiosperm families.
After WGD, most duplicated genes gradually accumulate deleterious mutations and vanish from the genome (within millions of years). More rarely the duplication leads to neo- or subfunctionalisation of the paralogs if these changes improve fitness . Currently, one of the best models explaining stabilisation of duplicated genes is the EAC model of subfunctionalisation (escape from adaptive conflict)  based on the fact that a single protein can perform multiple catalytic or structural functions. In such case, the selective optimization of one function may lead to a decline in another function, creating an adaptive conflict that preserves the single copy gene/protein in an intermediate state. Casual gene duplication can provide a solution – escape from the adaptive conflict via functional specialisation of the resulting paralogs .
Multiple angiosperm species contain just two psbO paralogs with similar extent of diversification, so we assume that the presence of two different PsbO proteins gives an advantage to these species. Although many plants prosper with a single PsbO gene, in species with two isoforms, the loss of one isoform negatively affects growth and photosynthesis, e.g. in A. thaliana [18, 19, 23, 45] or potato . It indicates that functions of current diversified PsbO isoforms are no more equivalent due to subfunctionalisation after the duplication.
Structural aspects in PsbO diversification
Protein functions are connected with protein structure. Therefore, identification of common structural differences between isoforms in multiple species can indicate common functional adaptation. If various plants used duplicated psbOs to solve the same adaptive conflict, the structural and functional differentiation of PsbO isoforms would be similar or identical irrespective of independent duplication in individual families.
To evaluate the between-isoform differences, we first divided the overall variability of PsbO sequences on each position of the primary structure into three classes. The variability in current PsbO sequences reflects both differences present already in the ancestor species before psbO duplication (between-families variability) and differences obtained after the duplication, including specific diversification of isoforms (between-isoforms variability) and species-specific changes (between-species variability; Fig. 2; see quantification bellow the alignment in Additional file 2). The frequency of each variability class on specific positions was mapped on the homology model of PsbO (Fig. 4). The model corresponded to other published homology models of higher plants’ PsbO [9, 10]. The mapping showed that occupancy of the differences on the PsbO surface was unequal and the locations of the three classes of variability significantly differ.
There are practically no differences between isoforms on the PSII-binding surface of PsbO (Fig. 5). However, Murakami et al.  reported that PsbO2 of A. thaliana is less efficient in reconstitution of oxygen evolution in vitro compared to PsbO1. Our analysis showed that the PSII-binding surface is highly conserved in all angiosperms. It indicates that the differences in water oxidation observed by Murakami et al.  were not caused by direct modulation of water oxidation on Mn4CaO5 cluster, but rather by some indirect effect.
The biggest contrast in localisation of between-isoform differences and other types of differences (between-family and between-species) is in the part of β5-β6 loop that is not interacting with PSII core proteins (the amino acid side chains are pointing towards lumen) and the adjoining part of the β6 strand. In this part of PsbO, there is a very high frequency of between-family differences and a high frequency of between-species differences, whereas between-isoform differences are nearly absent. This suggests that this part of the PsbO surface might be involved in binding of some other protein, whose interaction surface can differ in individual species or families. As isoforms do not differ in this region, it seems that both isoforms need to retain this interaction identical. The presence of a hypothetical interactor is further supported by the fact that an unassigned density was detected in this part of PSII supercomplex structure by cryo-electron tomography .
The between-isoform differences were located mostly at the end of the β-barrel protruding into the lumen and on the β1-β2 loop. This pattern was similar in all analysed families and even in the moss Physcomitrella and the relatively recently duplicated psbO in maize (Fig. 4, Additional file 7). This indicates that the differences between isoforms probably enabled the same or similar functional adaptation of PsbOs in all analysed families. Since the psbO duplications were independent, the functional divergence of PsbO isoforms likely represents a parallel evolution, further supporting the impact of observed diversification of PsbO isoforms.
Functional differences between PsbO isoforms
We found that the location of the differences between PsbO isoforms of A. thaliana fits the pattern found in other angiosperms. Nine out of 11 different amino acids are located at the luminal base of β-barrel and one is located on the β1-β2 loop (Additional file 7). Both PsbO isoforms of A. thaliana are able to stabilise the manganese-calcium cluster and enable water splitting [18, 19]. PsbO1 was demonstrated to provide more efficient water splitting , whereas PsbO2 was reported to have higher GTPase activity and was proposed to participate in D1 repair cycle [19, 21, 23].
The highest frequency of between-isoform differences is located just around the hypothetic GTP-binding site predicted by Lundin et al. , which is situated inside the luminal end of the β-barrel. Lundin et al. found hypothetic non-canonical GTP-binding domains in spinach  and A. thaliana PsbO sequence . G1 domain, binding α-phosphate, was predicted in β1 sheet, G2-G3 domain, binding γ-phosphate, in β2 sheet and G4 domain, binding guanine ring, in β4-β5 loop (marked in the alignment in Additional file 2). Regions surrounding the G2-G3 domain, i.e. β1-β2 and β2-β3 loops, were predicted to be Switches I and II, respectively. These switches could have different conformations in GDP- and GTP-bound state.
The proposed G2-G3 domain is rather conserved in angiosperms, whereas G1 and G4 domains and Switches I and II have high frequency of differences between isoforms. The position with the most frequent differences between isoforms (T46 in spinach PsbO) is located in G1 domain, the position with the second most frequent differences between isoforms (E139 in spinach PsbO) is in G4 domain (see positions 47 and 140 in alignment in Additional file 2).
Since PsbO isoforms in A. thaliana were found to differ in GTPase activity  and the location of differences between them corresponds to the situation in other species, we propose that, in general, the differences between isoforms might modulate the GTPase activity of PsbO.
It is important to mention that the amino acids varying between isoforms have side chains mostly pointing outside of the β-barrel and not inside towards the predicted GTP-binding site. However, these amino acids might not only change the binding of GTP, but also modulate accessibility of the binding site for GTP or binding of some regulatory protein. Consistently with this assumption, we found that amino acid residues with longer side chains often predominated in one isoform, whereas residues with shorter chains in the other (mainly in the case of chains with interaction-competent carboxylic groups). Though the substitution between glutamic and aspartic acid is often regarded to be nearly synonymous, these residues alternate within isoform pairs very regularly on some positions. Their strong conservation in orthologous isoforms within families (and in some cases even across more families) indicates that they can affect the PsbO function. The simple increase in the length of amino acid side chain can allow the protein to keep the structure, isoelectric point and other important properties unchanged, whereas it can strongly facilitate (affect) the interactions with other proteins, e.g. PsbP, another extrinsic protein of PSII. Interestingly, PsbP was reported to have a crystal structure similar to the GTPase Ran regulatory protein [21, 46]. The importance of certain carboxylic amino acid residue in PsbO can be demonstrated on spontaneously tuberising potato mutant, where the absence of one PsbO isoform could not be complemented with another allele with a point mutation substituting glutamic acid for aspartic acid (Fischer and Duchoslav, unpublished results; ).
Based on the situation in A. thaliana and suggested isoform evolution by the “escape from adaptive conflict” scenario, we can hypothesise that the amino acid substitutions adapt the isoform to its specific function, which is obviously connected with losing in the other of the two functions. Murakami et al.  tested efficiency of in vitro oxygen evolution with native PsbO isoforms and their chimers that combined fragments (or selected amino acids) from both A. thaliana isoforms. They found that the weaker performance of PsbO2 was connected mainly with the C-terminal third of the protein, namely with two substitutions, V186S (in β5-β6 loop) and L246I (with side chain pointing towards the hypothetical GTP-binding site). On the contrary, substitution V204I (also with side chain pointing towards the hypothetical GTP-binding site) improved the activity and obviously compensated for the decrease in oxygen evolution caused by substitutions on the positions 186 and 246. Unfortunately, the GTPase activity of PsbO was proposed later, so the activity was not determined in these chimeric proteins to test the trade-off hypothesis. The importance of the C-terminal part (often differing between isoforms) was also demonstrated in spinach PsbO. The last three amino acid residues, including leucin at position 246 (245 in spinach), were shown to be critical for binding to PSII and for restoration of oxygen evolving activity .
Besides the modulation of the GTPase activity, there are some other presumptive interactions of PsbO that can be affected by the observed differences between isoforms. An obvious possibility is modulation of interaction with PsbP (mentioned above) or PsbQ, the other extrinsic proteins of PSII. The extrinsic proteins of PSII were suggested to be involved in the interaction between the thylakoid membranes across the lumen [12, 48]. Thus, the changes on the luminal end of the β-barrel of PsbO might modulate this interaction regardless whether the interaction is direct (PsbO-PsbO) or via other extrinsic proteins. Another possibility is that PsbO is interacting with some protein bound to the same membrane through the luminal end of the β-barrel, which is protruding not only to the lumen but also slightly over the outer rim of PSII towards the strongly bound (S) LHCII trimer . This speculation is supported by the evidence that removal of PsbO has an effect on the position of LHCII in PSII-LHCII supercomplex  and that there is an unassigned density on the luminal side of the LHCII trimers in cryo-electron tomography structure of the supercomplex . The N-terminus of PsbO might also be involved in such interaction, because angiosperm PsbOs are longer by usually 13 amino acids compared to the cyanobacterial protein analysed by Umena et al. . Numerous biochemical studies were focused on site-directed mutagenesis of PsbO (reviewed in ref.  and ). However, these studies cannot be used for interpretation of our results, since the mutated residues were mostly the highly conserved ones.
Speculative model on the role of PsbO isoforms
Our findings together with literature data led us to formulate a speculative model explaining the functional specialisation of PsbO isoforms. Hong et al.  showed that the oxidation state of the Mn4CaO5 cluster influences the structure of PsbO. We suggest that when the D1 protein and Mn4CaO5 cluster get damaged, PsbO shifts in the GTPase cycle. This induces a change in conformation of the β1-β2 loop. As the β1-β2 loop is interacting with the other monomer of PSII in the PSII dimer , such change might induce monomerisation of PSII  and its release from the semi-stabile organization in grana that allows diffusion to stroma lamellae. In stroma lamellae, PsbO released from PSII changes the GTPase state and D1 is degraded . After reassembly of PSII, PsbO binds back.
We propose that the PsbO1 isoform is adapted for low light conditions. It has lower GTPase activity  (or lower affinity to GTP) and thus it does not induce as rapid monomerisation of PSII and exchange of damaged D1 subunit as the PsbO2 isoform. It was shown that under low irradiance, a part of PSII-LHCII supercomplexes form ordered crystalline arrays  and the grana are tightly stacked . Thus, PsbO1 might prevent premature disturbance of the effectively working crystalline arrays of PSII-LHCII supercomplexes and support tight stacking of thylakoid membranes. We also suggest that PsbO2 is, on the other hand, optimized for high light conditions, when efficient light harvesting is not important due to the excess of energy, but the repair of damaged D1 subunits must be very fast to prevent further damage of photosynthetic complexes. Correspondingly, Herbstová et al.  demonstrated that high light treatment leads to an increase of protein mobility in grana thylakoids. Moreover, the proportion of PSII in PSII-LHCII supercomplexes in A. thaliana is high in psbo2 mutant and low in psbo1 mutant compared to wild-type plants . Thus, the fast GTPase activity (or high affinity to GTP) of PsbO2 might facilitate rapid disorganisation of PSII-LHCII supercomplexes and crystalline arrays allowing rapid diffusion of damaged PSII out of the grana and their redeployment after D1 replacement. Modulation of the ratio of PsbO1 and PsbO2 bound to PSII can help to adjust and optimise the photosynthetic performance in response to light conditions. An increase in the amount of PsbO2 in A. thaliana after long-term cold stress , which is in many features similar to high light stress , supports our speculation. It suggests that the PsbO1/2 availability is regulated transcriptionally in the long term. The momentary binding of isoforms might be regulated, for example, by different dynamics of the GTPase cycle.
In the species with single PsbO isoform, the function of the PsbO is likely in the mid-way between the two specialised isoforms present in other species and supports both of the specialised functions, at least at the basic level. Depending on the growth strategy and habitat, the single isoform could be principally closer to A. thaliana PsbO1 (shade plants), to PsbO2 (sunny plants) or intermediate. For example, PsbO protein sequences of Poaceae species are closer to PsbO2 than to PsbO1 of A. thaliana, especially regarding the positions with alternating glutamic and aspartic acid residues (Additional files 2 and 6), which is consistent with the typical sunny habitat of these species. However, the species with single isoform that cluster to A. thaliana PsbO1 cannot be considered as shade plants. Consequently, the isoform sequence cannot be simply related to the actual habitat certain species, but it is obviously also affected by other factors such as evolutionary history or life strategy. Nevertheless, we can conclude from these considerations that A. thaliana PsbO1 and its functional analogs cannot be simply assumed as the main isoforms, because the importance of certain isoform depends on the habitat or even on actual growth conditions.
Our study showed that the pairs of PsbO isoforms evolved in numerous angiosperm lineages independently. Yet, the pairs of PsbO isoforms differ at similar regions in the protein structure, mostly on the luminally exposed end of the β-barrel structure near the predicted GTP-binding site. Mapping of conserved and variable positions on PsbO surface also indicated new potential interaction regions. Observed analogy in divergence between isoforms in various species indicates that structural diversification and subfunctionalisation of PsbO isoforms represents an example of parallel evolution. Features evolved in parallel likely bring a significant advantage, so we assume that diversification of PsbO isoforms improve photosynthetic performance under variable conditions. However, the predicted subfunctionalisation related to diverse GTPase activity will require further experimental confirmation.
Availability of supporting data
The phylogenetic data used in this study have been deposited in TreeBASE database (http://purl.org/phylo/treebase/phylows/study/TB2:S17601).
We thank Göran Samuelsson and Marek Romášek for their helpful comments and language corrections, Petr Janšta and Stanislav Vosolsobě for assistance with phylogenetic analysis, and Martin Weiser for helpful discussion about statistical tests. This work was supported by Charles University in Prague (projects GA UK No. 362211 and No. 1472314) and by Ministry of Education, Youth and Sports of Czech Republic (project No. LO1417).
- De Las RJ, Balsera M, Barber J. Evolution of oxygenic photosynthesis: genome-wide analysis of the OEC extrinsic proteins. Trends Plant Sci. 2004;9:18–25.Google Scholar
- Burnap RL, Sherman LA. Deletion mutagenesis in Synechocystis sp. PCC6803 indicates that the manganese-stabilizing protein of photosystem II is not essential for oxygen evolution. Biochemistry. 1991;30:440–6.PubMedView ArticleGoogle Scholar
- Mayfield SP, Bennoun P, Rochaix JD. Expression of the nuclear encoded OEE1 protein is required for oxygen evolution and stability of photosystem II particles in Chlamydomonas reinhardtii. EMBO J. 1987;6:313–8.PubMed CentralPubMedGoogle Scholar
- Yi X, McChargue M, Laborde S, Frankel LK, Bricker TM. The Manganese-stabilizing Protein Is Required for Photosystem II Assembly/Stability and Photoautotrophy in Higher Plants. J Biol Chem. 2005;280:16170–4.PubMedView ArticleGoogle Scholar
- Umena Y, Kawakami K, Shen J-R, Kamiya N. Crystal structure of oxygen-evolving photosystem II at a resolution of 1.9 Å. Nature. 2011;473:55–60.PubMedView ArticleGoogle Scholar
- Nield J, Balsera M, Rivas JDL, Barber J. Three-dimensional Electron Cryo-microscopy Study of the Extrinsic Domains of the Oxygen-evolving Complex of Spinach ASSIGNMENT OF THE PsbO PROTEIN. J Biol Chem. 2002;277:15006–12.PubMedView ArticleGoogle Scholar
- Caffarri S, Kouřil R, Kereïche S, Boekema EJ, Croce R. Functional architecture of higher plant photosystem II supercomplexes. EMBO J. 2009;28:3052–63.PubMed CentralPubMedView ArticleGoogle Scholar
- Kouřil R, Oostergetel GT, Boekema EJ. Fine structure of granal thylakoid membrane organization using cryo electron tomography. Biochim Biophys Acta BBA - Bioenerg. 1807;2011:368–74.Google Scholar
- De Las RJ, Barber J. Analysis of the Structure of the PsbO Protein and its Implications. Photosynth Res. 2004;81:329–43.View ArticleGoogle Scholar
- Lundin B, Thuswaldner S, Shutova T, Eshaghi S, Samuelsson G, Barber J, et al. Subsequent events to GTP binding by the plant PsbO protein: Structural changes, GTP hydrolysis and dissociation from the photosystem II complex. Biochim Biophys Acta BBA - Bioenerg. 2007;1767:500–8.View ArticleGoogle Scholar
- Bricker TM, Roose JL, Fagerlund RD, Frankel LK, Eaton-Rye JJ. The extrinsic proteins of Photosystem II. Biochim Biophys Acta BBA - Bioenerg. 1817;2012:121–42.Google Scholar
- Suorsa M, Aro E-M. Expression, assembly and auxiliary functions of photosystem II oxygen-evolving proteins in higher plants. Photosynth Res. 2007;93:89–100.PubMedView ArticleGoogle Scholar
- Bricker TM, Frankel LK. Auxiliary functions of the PsbO, PsbP and PsbQ proteins of higher plant Photosystem II: A critical analysis. J Photochem Photobiol B. 2011;104:165–78.PubMedView ArticleGoogle Scholar
- Spetea C, Hundal T, Lundin B, Heddad M, Adamska I, Andersson B. Multiple evidence for nucleotide metabolism in the chloroplast thylakoid lumen. Proc Natl Acad Sci U S A. 2004;101:1409–14.PubMed CentralPubMedView ArticleGoogle Scholar
- Seidler A. The extrinsic polypeptides of Photosystem II. Biochim Biophys Acta BBA - Bioenerg. 1996;1277:35–60.View ArticleGoogle Scholar
- The Arabidopsis Genome Initiative. Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000;408:796.View ArticleGoogle Scholar
- Murakami R, Ifuku K, Takabayashi A, Shikanai T, Endo T, Sato F. Characterization of an Arabidopsis thaliana mutant with impaired psbO, one of two genes encoding extrinsic 33-kDa proteins in photosystem II. FEBS Lett. 2002;523:138–42.PubMedView ArticleGoogle Scholar
- Murakami R, Ifuku K, Takabayashi A, Shikanai T, Endo T, Sato F. Functional dissection of two Arabidopsis PsbO proteins. FEBS J. 2005;272:2165–75.PubMedView ArticleGoogle Scholar
- Lundin B, Hansson M, Schoefs B, Vener AV, Spetea C. The Arabidopsis PsbO2 protein regulates dephosphorylation and turnover of the photosystem II reaction centre D1 protein. Plant J. 2007;49:528–39.PubMedView ArticleGoogle Scholar
- Dwyer SA, Chow WS, Yamori W, Evans JR, Kaines S, Badger MR, et al. Antisense reductions in the PsbO protein of photosystem II leads to decreased quantum yield but similar maximal photosynthetic rates. J Exp Bot. 2012;63:4781–95.PubMed CentralPubMedView ArticleGoogle Scholar
- Lundin B, Nurmi M, Rojas-Stuetz M, Aro E-M, Adamska I, Spetea C. Towards understanding the functional difference between the two PsbO isoforms in Arabidopsis thaliana—insights from phenotypic analyses of psbo knockout mutants. Photosynth Res. 2008;98:405–14.PubMedView ArticleGoogle Scholar
- Goulas E, Schubert M, Kieselbach T, Kleczkowski LA, Gardeström P, Schröder W, et al. The chloroplast lumen and stromal proteomes of Arabidopsis thaliana show differential sensitivity to short- and long-term exposure to low temperature. Plant J. 2006;47:720–34.PubMedView ArticleGoogle Scholar
- Allahverdiyeva Y, Mamedov F, Holmström M, Nurmi M, Lundin B, Styring S, et al. Comparison of the electron transport properties of the psbo1 and psbo2 mutants of Arabidopsis thaliana. Biochim Biophys Acta BBA - Bioenerg. 2009;1787:1230–7.View ArticleGoogle Scholar
- Bricker TM, Frankel LK. The psbo1 Mutant of Arabidopsis Cannot Efficiently Use Calcium in Support of Oxygen Evolution by Photosystem II. J Biol Chem. 2008;283:29022–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Fischer L, Lipavska H, Hausman J-F, Opatrny Z. Morphological and molecular characterization of a spontaneously tuberizing potato mutant: an insight into the regulatory mechanisms of tuber induction. BMC Plant Biol. 2008;8:117.PubMed CentralPubMedView ArticleGoogle Scholar
- Database of Expressed Sequence Tags, NCBI GenBank [http://www.ncbi.nlm.nih.gov/nucest]
- Plant Genome Database [http://www.plantgdb.org/]
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.PubMed CentralPubMedView ArticleGoogle Scholar
- Altschul SF, Wootton JC, Gertz EM, Agarwala R, Morgulis A, Schäffer AA, et al. Protein database searches using compositionally adjusted substitution matrices. FEBS J. 2005;272:5101–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30:3059–66.PubMed CentralPubMedView ArticleGoogle Scholar
- Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In: Gateway Computing Environments Workshop (GCE), 2010. 2010. p. 1–8.View ArticleGoogle Scholar
- Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.PubMedView ArticleGoogle Scholar
- Goodstein DM, Shu S, Howson R, Neupane R, Hayes RD, Fazo J, et al. Phytozome: a comparative platform for green plant genomics. Nucleic Acids Res. 2012;40:D1178–86.PubMed CentralPubMedView ArticleGoogle Scholar
- The R Project for Statistical Computing [http://www.r-project.org/]
- Charif D, Lobry JR: SeqinR 1.0-2: A Contributed Package to the R Project for Statistical Computing Devoted to Biological Sequences Retrieval and Analysis. In Structural Approaches to Sequence Evolution. Edited by Bastolla U, Porto M, Roman HE, Vendruscolo M. Springer Berlin Heidelberg; 2007:207–232. [Biological and Medical Physics, Biomedical Engineering]Google Scholar
- Schwede T, Kopp J, Guex N, Peitsch MC. SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res. 2003;31:3381–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Arnold K, Bordoli L, Kopp J, Schwede T. The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006;22:195–201.PubMedView ArticleGoogle Scholar
- Guex N, Peitsch MC. SWISS-MODEL and the Swiss-Pdb Viewer: An environment for comparative protein modeling. Electrophoresis. 1997;18:2714–23.PubMedView ArticleGoogle Scholar
- POV-Ray - The Persistence of Vision Raytracer [http://www.povray.org/]
- Flagel LE, Wendel JF. Gene duplication and evolutionary novelty in plants. New Phytol. 2009;183:557–64.PubMedView ArticleGoogle Scholar
- Fawcett JA, Maere S, de Peer YV. Plants with double genomes might have had a better chance to survive the Cretaceous–Tertiary extinction event. Proc Natl Acad Sci. 2009;106:5737–42.PubMed CentralPubMedView ArticleGoogle Scholar
- Vanneste K, Baele G, Maere S, Peer YV de: Analysis of 41 plant genomes supports a wave of successful genome duplications in association with the Cretaceous-Paleogene boundary. Genome Res 2014;24:1334–47.Google Scholar
- Hittinger CT, Carroll SB. Gene duplication and the adaptive evolution of a classic genetic switch. Nature. 2007;449:677–81.PubMedView ArticleGoogle Scholar
- Liu H, Frankel LK, Bricker TM. Functional Analysis of Photosystem II in a PsbO-1-Deficient Mutant in Arabidopsis thaliana. Biochemistry. 2007;46:7607–13.PubMedView ArticleGoogle Scholar
- Ifuku K, Nakatsu T, Kato H, Sato F. Crystal structure of the PsbP protein of photosystem II from Nicotiana tabacum. EMBO Rep. 2004;5:362–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Betts SD, Lydakis-Simantiris N, Ross JR, Yocum CF. The Carboxyl-Terminal Tripeptide of the Manganese-Stabilizing Protein Is Required for Quantitative Assembly into Photosystem II and for High Rates of Oxygen Evolution Activity†. Biochemistry. 1998;37:14230–6.PubMedView ArticleGoogle Scholar
- De Las RJ, Heredia P, Roman A. Oxygen-evolving extrinsic proteins (PsbO, P, Q, R): Bioinformatic and functional analysis. Biochim Biophys Acta BBA - Bioenerg. 2007;1767:575–82.View ArticleGoogle Scholar
- Dekker JP, Boekema EJ. Supramolecular organization of thylakoid membrane proteins in green plants. Biochim Biophys Acta BBA - Bioenerg. 2005;1706:12–39.View ArticleGoogle Scholar
- Boekema EJ, van Breemen JFL, van Roon H, Dekker JP. Conformational Changes in Photosystem II Supercomplexes upon Removal of Extrinsic Subunits. Biochemistry. 2000;39:12907–15.PubMedView ArticleGoogle Scholar
- Williamson A. Structural and functional aspects of the MSP (PsbO) and study of its differences in thermophilic versus mesophilic organisms. Photosynth Res. 2008;98:365–89.PubMedView ArticleGoogle Scholar
- Hong SK, Pawlikowski SA, Vander Meulen KA, Yocum CF. The oxidation state of the photosystem II manganese cluster influences the structure of manganese stabilizing protein. Biochim Biophys Acta BBA - Bioenerg. 2001;1504:262–74.View ArticleGoogle Scholar
- Aro E-M, Suorsa M, Rokka A, Allahverdiyeva Y, Paakkarinen V, Saleem A, et al. Dynamics of photosystem II: a proteomic approach to thylakoid protein complexes. J Exp Bot. 2005;56:347–56.PubMedView ArticleGoogle Scholar
- Kouřil R, Wientjes E, Bultema JB, Croce R, Boekema EJ. High-light vs. low-light: Effect of light acclimation on photosystem II composition and organization in Arabidopsis thaliana. Biochim Biophys Acta BBA - Bioenerg. 1827;2013:411–9.Google Scholar
- Kirchhoff H. Structural changes of the thylakoid membrane network induced by high light stress in plant chloroplasts. Philos Trans R Soc B Biol Sci. 2014;369:20130225.View ArticleGoogle Scholar
- Herbstová M, Tietz S, Kinzel C, Turkina MV, Kirchhoff H. Architectural switch in plant photosynthetic membranes induced by light stress. Proc Natl Acad Sci. 2012;109:20130–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Huner NPA, Öquist G, Sarhan F. Energy balance and acclimation to light and cold. Trends Plant Sci. 1998;3:224–30.View ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.