<title>Background</title> <p>ADP-glucose pyrophosphorylase (AGPase; EC 2.7.7.27) catalyses a rate-limiting step in starch synthesis, the formation of ADP-glucose from glucose-1-P and ATP. ADP-glucose is the predominant, if not sole, precursor for starch synthesis. While AGPase is a homotetramer in bacteria (including cyanobacteria), it is a heterotetramer in angiosperms and green algae. This heterotetramer comprises two identical large and two identical small subunits. They exhibit a high degree of identity to each other and to the cyanobacterial AGPase, pointing to an origin by gene duplication early in the evolution of plants and green algae (Figure ##FIG##0##1A##) (Additional file ##SUPPL##0##1##) [##REF##1318389##1##,##REF##17496118##2##]. The two subunits have complementary rather than redundant functions, and knockout mutations in either abolish more than 90% of AGPase activity in some experimental systems [##REF##985379##3##].</p> <p>Although both subunits are necessary for full AGPase activity, the angiosperm small subunit appears more conserved than the large subunit throughout its sequence (small subunits exhibit an average of 91.3% amino acid identity while large subunits an average of 70.8% identity) [##REF##17496118##2##]. However, mean percent identities might be misleading since the genes encoding both subunit genes of AGPase underwent a number of duplications after the initial duplication that generated the two subunits. The potential confusion due to the comparison of paralogs rather than orthologs can be overcome by methods that incorporate phylogeny, such as the use of maximum likelihood (ML) to estimate ω (the ratio of non-synonymous substitutions per non-synonymous site [K<sub>A</sub>] to synonymous substitutions per synonymous site [K<sub>S</sub>]). The ML estimate of ω for the large subunit is ~2.7-fold greater than the estimate for the small subunit [##REF##17496118##2##], suggesting a higher rate of amino acid replacement for the large subunit.</p> <p>Although ω provides a convenient and commonly used method to examine evolutionary constraints, it has typically been used to examine sequences that have diverged relatively recently. The rate of synonymous evolution is relatively high in plant nuclear genes [##REF##3480529##4##, ####REF##8816790##5##, ##UREF##0##6##, ##REF##10545473##7####10545473##7##] and estimates of K<sub>S </sub>appear saturated in analyses of some angiosperm gene families, even for relatively shallow evolutionary divergences [##REF##10471724##8##]. Hence, the accuracy of ω estimates for ancient divergences is unclear. Another potential problem for the use of ω is the assumption that mutations at synonymous sites are neutral. It has been suggested that synonymous sites are subject to both positive and purifying selection [##REF##8754221##9##, ####REF##16418745##10##, ##REF##16221894##11##, ##REF##17522087##12####17522087##12##]. The action of selection on synonymous sites may explain why adding among-sites rate variation for synonymous sites to models of codon evolution improves their fit to empirical data [##REF##16107593##13##,##UREF##1##14##]. Saturation and among-sites rate variation both have the potential to cause K<sub>S </sub>to be underestimated (and ω to be overestimated since ω = K<sub>A</sub>/K<sub>S</sub>); biased estimates of ω will lead to incorrect inferences regarding evolutionary constraints on the proteins being analyzed. Finally, ω cannot detect changes in the evolutionary rate when rates of synonymous and non-synonymous substitution increase or decrease simultaneously [##REF##15014159##15##].</p> <p>The almost 3-fold difference in evolutionary rates for the AGPase subunits is a paradox because random mutagenesis revealed that maize endosperm AGPase subunits expressed in bacteria are equally susceptible to activity-altering amino acid changes [##REF##17496118##2##]. Georgelis et al. [##REF##17496118##2##] proposed that the difference in evolutionary rates between AGPase subunits reflected, at least in part, the differences between the subunits in their tissue-expression patterns and the fact that the small subunit has to interact with multiple large subunits in plants. Here, we establish the pattern and timing of duplications in the AGPase gene family and estimate absolute rates of AGPase sequence evolution. Functional divergence has been observed among AGPase subunits based on biochemical criteria [##REF##12748181##16##, ####REF##14520572##17##, ##REF##10557246##18##, ##REF##16299180##19##, ##REF##18024561##20####18024561##20##]. One of our primary goals was to identify candidate sites for functional divergence. We identify specific AGPase sites apparently subject to either positive selection or branch-specific patterns of rate variation (types-I and -II divergence as defined by Gu [##REF##10605109##21##,##REF##11264396##22##]).</p>
<title>Methods</title> <title>Sequence Retrieval and Alignment</title> <p>Full-length AGPase sequences from plants were retrieved from the NCBI database and the DOE Joint Genome Institute (JGI) web site, and the source and accession numbers of all sequences are presented in Additional file ##SUPPL##9##10##. DNA and protein sequence alignments were obtained using the MEGA software [##UREF##7##68##] with BLOSUM matrix followed by manual inspection. The poorly aligned N-termini (~70–80 amino acids for the large subunit and ~40 amino acids for the small subunit) were excluded from alignment. The large subunit amino acid numbers used correspond to the protein encoded by the maize <italic>Shrunken-2 </italic>(<italic>Sh2</italic>) gene (Accession #: P55241) while the small subunit amino acid numbers used correspond to the protein, encoded by the maize <italic>Brittle-2 </italic>(<italic>Bt2</italic>) gene (Accession #: AAQ14870).</p> <title>Phylogenetic analysis</title> <p>Estimates of AGPase gene trees based upon nucleotide data were obtained using the GARLI (Genetic Algorithm for Rapid Likelihood Inference) software [##UREF##8##69##]. Estimates of phylogeny estimated for protein sequence alignments were obtained either using neighbor joining in the MEGA software or ML in the RAxML software [##REF##16928733##70##], using the JTT model [##REF##1633570##71##] in both cases (to estimate distances or calculate likelihoods). Bootstrap support [##UREF##9##72##] was calculated using 100 replicates.</p> <p>Branch lengths were estimated by ML, with those form amino acid trees based estimated using AAML, using the Dayhoff (PAM) model [##UREF##10##73##] with Γ-distributed rates. Nonsynonymous substitutions per nonsynonymous site (K<sub>A</sub>) and synonymous substitutions per synonymous site (K<sub>S</sub>) and the ratio of these values (ω = K<sub>A</sub>/K<sub>S</sub>) were estimated using CODEML. AAML and CODEML are programs in the PAML (Phylogenetic Analysis by Maximum Likelihood) package [##REF##9367129##74##].</p> <title>Reconciled tree analysis</title> <p>Reconciled tree analyses map a gene tree onto a species tree [##UREF##11##75##,##REF##9918954##76##], and most commonly used procedure for doing this is gene tree parsimony, which minimizes the number of duplication events needed to explain a specific gene tree given the species tree [##UREF##5##32##]. However, some error is typically associated with estimates of phylogeny for individual genes, and accommodating error in gene trees is difficult in reconciled tree analyses [##UREF##11##75##]. Chen et al. [##REF##11108472##77##] suggested an algorithm that rearranged gene trees to increase congruence with the species tree when nodes in the gene tree were poorly supported to limit the impact of error in the gene tree. We implemented this idea manually, by rearranging nodes in the gene tree with limited support (those with bootstrap support &lt; 70%; see [##UREF##12##78##]) to increase the congruence with the species tree. This yielded two estimates of the numbers of duplications, one based on the optimal gene trees and a conservative estimate based on well-supported nodes.</p> <title>Estimation of the absolute rate of evolution</title> <p>Estimates of the absolute rate of amino acid evolution and synonymous site evolution for each subunit were obtained using the "tip procedure", which uses the average number of amino acid substitutions per site along unique paths from each tip (extant sequence) to a dated speciation event (allowing us to avoid pseudoreplication of rate estimates). In addition to this method, absolute rates of amino acid substitution were also obtained by estimating the age of each node after smoothing rates using penalized likelihood in r8s ("PL method") [##REF##12538260##79##]. The PL method was used because the tip method is biased towards recent branches, although the PL method also has the potential to be affected by saturation, especially for synonymous sites.</p> <title>Detection of branch-specific patterns of rate variation among sites and positive selection</title> <p>Type-I and type-II functional divergence among large or small subunit groups was examined using the DIVERGE software [##REF##11934757##48##], which implements the tests suggested by Gu [##REF##11264396##22##,##REF##16864604##47##] that can be used to determine whether the coefficients of divergence (θ<sub>I </sub>and θ<sub>II</sub>) are significantly greater than zero. Amino acid sites likely to have undergone types-I or -II divergence were detected as those with a posterior probability &gt; 0.5–0.6.</p> <p>Sites subject to positive selection were identified using the site, branch and branch-site models implemented in CODEML, using the model comparisons recommended by Yang et al. [##REF##11080365##80##]. The first comparison was between model M1a, which includes two ω values (one for sites subject to purifying selection with ω &lt; 1 and one for neutral sites with ω = 1) and model M2a (which adds positively selected sites [with ω &gt; 1] to model M1a). The second comparison was between model M7, which assumes values of ω at different sites are β-distributed, and model M8 (which adds positively selected sites to model M7). We also searched for positive selection by estimating different values of ω for each branch and using branch-site models. For the first test, we compared a model specified as model = 2 NSsites = 2 to model M1a [##REF##12032247##81##]. For the second test, we compared a model specified as model = 2 NSsites = 2 to a model whose specifications are model = 2 NSsites = 2 fix_omega = 1 and omega = 1. When the likelihood ratio test was significant, the Bayes empirical Bayes method was used to calculate posterior probabilities that sites were subject to positive selection [##REF##15689528##82##].</p>
<title>Results</title> <title>Patterns of AGPase gene duplication</title> <p>It is well known that genes can have three possible fates after duplication [##REF##11073452##23##, ####UREF##2##24##, ##UREF##3##25##, ##REF##15252449##26##, ##REF##17267422##27####17267422##27##]: (1) nonfunctionalization, in which one duplicate is lost, (2) subfunctionalization, in which the functions of the original single-copy gene are partitioned between the duplicates, and (3) neofunctionalization, in which one duplicate gains a novel function. The latter two processes can result in paralogs that persist for a substantial length of time (although a few exceptions have been proposed, such as pseudogene resurrection [##REF##17267422##27##]). Throughout this work, the term duplication will be limited to the description of the latter two processes. Although gene loss can be as important as duplication for shaping genomes [##REF##15130831##28##], we have avoided making major conclusions based on gene loss since many organisms included in these analyses lack complete genome sequences.</p> <p>Inclusion of AGPase sequences from the moss <italic>Physcomitrella patens </italic>[##REF##18079367##29##], which has 7 large subunit and 4 small subunit genes, placed a major constraint on the earliest divergence within the large subunit family since the moss sequences were intermixed with angiosperm sequences. This indicates that the earliest duplication in the large subunit family occurred prior to the divergence of angiosperms and mosses, more than 400 million years (MY) ago [##UREF##4##30##]. Since the rate of synonymous evolution in angiosperms varies from ~2 × 10<sup>-9 </sup>to ~10 × 10<sup>-9 </sup>synonymous substitutions per synonymous site per year [##REF##3480529##4##, ####REF##8816790##5##, ##UREF##0##6##, ##REF##10545473##7####10545473##7##], values of K<sub>S </sub>in excess of 2 are expected for some comparisons, which may make estimates of K<sub>S </sub>problematic [##REF##17975267##31##].</p> <p>Since many divergence times for plants can be constrained to reasonable ranges, it should be possible to estimate absolute rates of AGPase subunit amino acid evolution and establish whether they correlate with estimates of ω. However, this requires differentiating between speciation and gene duplication events in AGPase phylogeny. Gene family phylogenies reflect both speciation and duplication events, and these events can be distinguished by reconciled tree analyses if the gene and species trees are known. Gene tree parsimony [##UREF##5##32##] is the most commonly used reconciled tree method, and the only approach practical for even moderately sized phylogenies at this time. Reconciling the AGPase gene trees with the best available estimate of the land plant species tree (Additional file ##SUPPL##1##2##) revealed 11–14 large subunit duplications (Figure ##FIG##0##1B##) and 5–7 small subunit duplications (Figure ##FIG##0##1C##). It may be appropriate to view the lower estimates, which are based upon well-supported nodes, as the primary results since they are based on modified versions of the gene trees (Additional file ##SUPPL##2##3A, B##) in which the topology was rearranged near poorly supported nodes to increase congruence with the species tree (Methods). In contrast, the higher estimates were based only on reconciling the optimal estimates of the gene trees (Figure ##FIG##0##1B, C##) with the species tree. Regardless, both analyses indicate that the large subunit genes underwent a larger number of duplications than the small subunit genes.</p> <p>AGPase large subunits have narrower tissue-specificity than small subunits [##REF##7700228##33##, ####REF##9469935##34##, ##REF##15821022##35##, ##REF##16275672##36##, ##REF##16957017##37##, ##REF##15598655##38####15598655##38##], and the large subunit phylogeny appears more complex (with four major clades some of which include both monocots and eudicots; Figure ##FIG##0##1B##) than the small subunit phylogeny (Figure ##FIG##0##1C##). Large subunit group 1 genes are predominantly expressed in leaves, group 2 genes are expressed both in source and sink tissues, group 3 genes are expressed sink tissues (these genes are subdivided into group 3a in eudicots and group 3b in monocots), and group 4 corresponds to a clade of two sequences that have not been characterized yet in terms of function and expression patterns (Figure ##FIG##0##1B##). Some of the major large subunit clades arose prior to the divergence of monocots and eudicots, and the optimal placement of the AGPase large subunit sequences from <italic>Physcomitrella </italic>suggests that the first duplication in the large subunit happened around 400 MY ago (Figure ##FIG##0##1A##). In contrast, there is no evidence that angiosperm small subunits underwent a duplication prior to the divergence of monocots and eudicots, and we have divided them into a monocot clade (group 1) genes and a eudicot clade (group 2). These results emphasize that the large subunit underwent a larger number of duplications than did the small subunit and that only large subunit duplications began before the divergence of monocots and eudicots.</p> <title>Absolute rates of AGPase evolution</title> <p>Absolute rates of amino acid evolution for AGPase subunits were estimated by examining terminal branch lengths for divergences that reflect speciation events with known divergence times (these divergence times are presented in Additional file ##SUPPL##1##2##). This approach is called the tip procedure since it involves only terminal branches (Methods), and it revealed that the average rate of evolution for the large subunit was 2.7-fold faster than that of the small subunit (Figure ##FIG##1##2##). This rate difference was both congruent with the difference in ML estimates of ω [##REF##17496118##2##] and highly significant (<italic>P </italic>= 0.0006 by Student's unpaired <italic>t</italic>-test). Our conclusions were unchanged if we limited consideration to strongly supported duplication events (those retained when bootstrap support was considered; see Additional file ##SUPPL##2##3##).</p> <p>Estimates of the absolute rate of amino acid substitution for AGPase subunits obtained by the penalized likelihood (PL) method (Figure ##FIG##2##3##) were very similar to those obtained using the tip procedure (Figure ##FIG##1##2##). Using gene trees in which poorly supported nodes were rearranged to minimize number of duplications yielded similar results (Additional file ##SUPPL##3##4##, Figure ##FIG##1##2##). Thus, very similar estimates of the absolute rate of amino acid substitution were obtained despite the different assumptions made by the tip procedure and the PL method.</p> <p>The absolute rate of synonymous evolution was estimated using ML estimates of K<sub>S </sub>(Additional file ##SUPPL##4##5A, B##). The tip procedure resulted in virtually identical rates for both large and small subunit genes (6.5 × 10<sup>-9 </sup>synonymous substitutions per synonymous site per year; Additional file ##SUPPL##4##5C##). PL rate estimates were 5.5 × 10<sup>-9</sup>± 0.3 × 10<sup>-9 </sup>and 6.3 × 10<sup>-9 </sup>± 0.2 × 10<sup>-9 </sup>(mean ± standard error) synonymous substitutions per synonymous site per year for the large and small subunit, respectively (data not shown). The slightly lower estimates based upon PL are consistent with saturation being a problem, but presumably only for the deepest branches in the tree. All of these values are well within the range of previous estimates for a variety of angiosperm genes (which range from approximately 2 × 10<sup>-9 </sup>to 10 × 10<sup>-9 </sup>synonymous substitutions per synonymous site per year and exhibit some variation among lineages [##REF##3480529##4##, ####REF##8816790##5##, ##UREF##0##6##, ##REF##10545473##7####10545473##7##,##REF##16973872##39##,##REF##16702410##40##]). This suggests that there are little to no constraints on the synonymous sites of angiosperm AGPase genes and, when combined with estimate of the absolute rate of sequence evolution, that there was minimal bias in our estimates of ω.</p> <title>Does the large subunit show temporary or permanent elevation of ω?</title> <p>Estimates of the mean rate of evolution for AGPase, whether based upon ω [##REF##17496118##2##] or the absolute rate of amino acid substitution (Figure ##FIG##1##2##), show a substantially higher rate for the large subunit. The existence of these rate differences despite identical sensitivities to mutations when expressed in bacteria suggests that there are important differences <italic>in planta</italic>. Transient increases in the evolutionary rate might explain the observed rate differences if they are more common for the large subunit. Large-scale analyses provide evidence for transient increase in the rate of non-synonymous evolution. Indeed, paralogs with a recent origin (defined as those with K<sub>S </sub>&lt; 0.1) accumulate more non-synonymous mutations per non-synonymous site (and thus have a higher K<sub>A</sub>) relative to the number of the synonymous mutation per synonymous site than older duplicates [##REF##11073452##23##]. Thus, ω is expected to be elevated for paralogs with K<sub>S </sub>&lt; 0.1 relative to those with K<sub>S </sub>&gt; 0.1. Lynch and Conery [##REF##11073452##23##] interpreted this phenomenon as reflecting a temporary relaxation of constraints, positive selection, or a combination of both phenomena. Thus, it is important to consider the potential impact of the elevation of ω on our analyses of the evolutionary processes that shaped the AGPase gene family.</p> <p>The large subunit underwent more gene duplications than did the small subunit (Figure ##FIG##0##1##). Thus, the higher mean estimates of ω for the large subunit might reflect a larger number of periods during which ω is elevated (due to temporary relaxation of constraints and/or positive selection) rather than permanent relaxation of constraints for the large subunit. To distinguish between transient and permanent relaxation of constraint we tested whether the non-synonymous rate was increased after duplication and if the increase is sufficient to explain the observed differences in the mean rate. Branches in both large and small subunit gene trees (Figure ##FIG##0##1##) were placed into two groups, the first of which (class 1) contained branches that follow a duplication event with K<sub>S </sub>= ~0.1 (these branches are shown in Figure ##FIG##0##1B, C## and Additional file ##SUPPL##2##3##). The second group (class 2) contained all other branches (branches that follow either a speciation event or a duplication and have K<sub>S </sub>&gt; 0.1). We examined two nested models using the likelihood ratio test (LRT) [##REF##7288891##41##, ####REF##7679448##42##, ##UREF##6##43##, ##REF##9092465##44##, ##REF##12679544##45####12679544##45##] to determine whether ω for class 1 branches is higher than ω for class 2 branches for either subunit using CODEML (included in PAML software). The more complex model, which assumes two different ω values (one ω for class 1 and one ω for class 2), was favored over the null hypothesis model, which assumes a single ω for both classes, for both the small subunit (2δ<italic>ln</italic>L = 18.4; <italic>P </italic>&lt; 0.001) and the large subunit (2δ<italic>ln</italic>L = 7.48; <italic>P </italic>= 0.005) (for details see Table ##TAB##0##1##). The ω estimate for short branches following duplications was 1.3 to 1.5-fold greater than the ω estimate for all the other branches when the large subunit was examined and 2 to 2.8-fold greater for the small subunit (Table ##TAB##0##1##). These results support periods of increased ω after duplications in the AGPase gene family, due to the temporary relaxation of constraints, positive selection, or both. However, these results also indicate that this phenomenon cannot explain the differential rates of amino acid sequence divergence of the two AGPase subunits, since estimates of ω for the large subunit are 2.6-fold greater than estimates for the small subunit (Table ##TAB##0##1##). Instead these results suggest that the small subunit is permanently subject to greater purifying selection than is the large subunit.</p> <title>What is the role of positive selection in AGPase evolution?</title> <p>To examine the potential role of positive selection in AGPase evolution, we used ML to compare two distinct sets of models. The first model set contains a neutral (null) model M1a allowing two categories of sites, one with ω = 0 and one with ω = 1, and model M2a that adds an extra category of sites with ω &gt; 1. The second includes a neutral (null) model M7 assuming that ω is β-distributed among sites and model M8 that adds an extra category of sites with ω &gt; 1 [##REF##15514074##46##]. Neither of the models with positive selection was significantly better than the null model based upon the LRT for either of the subunits when the tests were applied to the complete trees for the small subunit (Figure ##FIG##0##1C##) or large subunit (Figure ##FIG##0##1B##) (data not shown). Likewise, neither of the models that include positive selection was significantly better when the tests were applied to the individual groups (Figure ##FIG##0##1B, C##) within each subunit (groups 1, 2, 3A, and 3B for the large subunit as well as groups 1, and 2 for the small subunit) (data not shown). These results indicate either that positive selection has not played a role in the evolution of the large and small subunit of AGPase in the angiosperms or that these tests have insufficient power.</p> <p>To increase the power of our tests for positive selection, we used branch-site models to examine the potential for positive selection at specific sites in all tree branches separately (Methods). There are branches in the large or small subunit tree on which specific sites may be subject to positive selection (Figure ##FIG##3##4##), with a total of 0.8 and 0.2 sites/branch potentially affected by positive selection for the large and the small subunit respectively (Additional file ##SUPPL##5##6##). However, the limited number of sites potentially affected by positive selection suggests that purifying selection is the major force in the evolution of both AGPase subunits in angiosperms. Thus, positive selection cannot explain the different rates of amino acid evolution for the AGPase subunits.</p> <title>Functional divergence of AGPase subunits</title> <p>Gu [##REF##10605109##21##,##REF##11264396##22##,##REF##16864604##47##] proposed that specific sites in proteins have the potential to undergo two distinct types of divergences after gene duplication (e.g., divergence among the different groups within the large and small subunits), and these types were designated type-I and type-II divergence. Sites that have undergone type-I divergence are conserved in one group but variable in another group while type-II sites are fixed in both groups but differ between groups [##REF##10605109##21##,##REF##11264396##22##].</p> <p>Tests for types-I and -II divergence based upon the relevant coefficients of divergence (θ), which correspond to the probability that a specific site has undergone type-I or -II divergence in a pairwise comparison, were proposed by Gu [##REF##10605109##21##,##REF##11264396##22##,##REF##11934757##48##]. These tests for types-I and -II establish whether the relevant θ value is significantly greater than zero. We conducted all possible pairwise comparisons among small (groups 1 and 2) and large (groups 1,2, 3a and 3b) subunits with the exception of group 4 of large subunits (that were excluded because the group included only two sequences). All of type-I coefficients (θ<sub>I</sub>) of functional divergence were significantly greater than zero while none of the type-II coefficients (θ<sub>II</sub>) were significantly greater than zero (Table ##TAB##1##2##). The estimates of θ<sub>I </sub>are also much larger than the estimates of θ<sub>II</sub>, suggesting that type-I divergence is the dominant pattern of sequence evolution for AGPase large and small subunit groups.</p> <title>Sites contributing to functional divergence among AGPase groups</title> <p>We identified the sites likely to be involved in the changes of functional constraints between groups revealed by the significant values of θ<sub>I </sub>using a posterior probability analysis (Additional file ##SUPPL##6##7##). A greater number of large subunit sites appear to have undergone type-I divergence than the number of small subunit sites. Specifically, for an alignment with 453 amino acids we found evidence that 78 large subunit sites and 13 small subunit sites show evidence for type-I divergence (Additional file ##SUPPL##7##8##). These sites appear randomly distributed with respect to secondary structure (data not shown), but pairwise comparisons among groups in the large and the small subunit reveal some non-random patterns in the distribution of sites that are conserved in one group but not another. For instance, large subunit group 1 proteins (leaf isoforms) had more conserved sites in the N-terminus while group 3a proteins (sink isoforms) exhibited more conserved in the C-terminus. The large subunit crystal structure has not been elucidated yet. However, the high degree of amino acid sequence identity (~43%) and similarity (~61%) [##REF##17496118##2##] (Additional file ##SUPPL##0##1##) to the small subunit along with structure modelling (unpublished data) strongly suggest that the 3D structure of the large subunit is almost identical to the known structure of the small subunit. The N-terminal domain includes the active site that resembles a Rossman fold while the C-terminal domain is a β-helix that extensively interacts with the N-terminal domain based on the structure of the potato tuber small subunit elucidated by Jin et al. [##REF##15692569##49##]. Both domains are important for stability and catalytic/allosteric properties of the enzyme [##REF##15692569##49##, ####REF##11339804##50##, ##REF##11237727##51##, ##REF##9628013##52####9628013##52##], but the non-random spatial distribution of type I sites clearly suggests these sites should be targeted in mutational studies focused on the analysis of AGPase structure/function relationships.</p> <p>In contrast to the significant estimates of θ<sub>I</sub>, estimates of θ<sub>II </sub>(the type-II coefficient) were not significantly greater than zero. However, the power of the test for type-II divergence is unclear. To determine whether there was any evidence for type-II divergence we examined the sequences to determine whether any specific sites show evidence for this type of divergence. The likelihood that these sites reflect <italic>bona fide </italic>instances of type-II divergence would be increased if they correspond to sites that have been identified in mutagenesis studies. The posterior ratio test, using a cut-off 2 (i.e., a posterior probability of 0.5), identified a few potential type-II sites (Additional file ##SUPPL##8##9##). Some sites were also identified by analyses of positive selection and one has an important function revealed by mutagenesis (see Discussion), so it is reasonable to postulate that at least some of these sites reflect genuine instances of type-II divergence that have contributed to the specialization among the large and small subunits.</p>
<title>Discussion</title> <p>The large and the small subunit of AGPase in plants exhibit considerable sequence identity [##REF##1318389##1##,##REF##17496118##2##] and they reflect a gene duplication that occurred prior to the divergence of land plants and green algae. Both subunits are equally sensitive to activity-altering amino acids, at least when expressed in a bacterial system [##REF##17496118##2##]. However, the small subunit of angiosperms is more conserved than large subunit based upon estimates of the rate of evolution (both estimates of ω [##REF##17496118##2##] and estimates of the absolute rate of amino acid evolution). These results suggest saturation has not had a major impact upon the estimation of ω for angiosperm AGPases, and this may reflect, at least in part, the limited codon bias of these genes (data not shown) since codon bias can have a major impact on the estimation of K<sub>S </sub>[##REF##10471724##8##]. They also suggest that estimation of the absolute rate of amino acid evolution provides a valid method that can be used as an alternative to ω analysis when the use of ω is not appropriate (e.g., ancient divergences, especially for gene families with strong codon bias, and for gene families where there is evidence for strong selection on synonymous sites). Taken as a whole, these results confirm that plant AGPases represent a genuine paradox: the large and small subunits exhibit similar sensitivities to activity-altering changes but differ almost 3-fold in their rates of non-synonymous evolution.</p> <p>Although a temporary elevation of ω after duplication was observed for both large and small subunit gene families, this transient elevation cannot account for the overall difference in ω value (and the related difference in the overall rate of amino acid substitutions) between the two subunits. Additionally, although both subunits appear to have been subject to positive selection the observed rate differences are too large to be explained by postulating that they reflect greater differences in positive selection. Based upon the falsification of these hypotheses, we conclude that the small subunit has been evolving under stronger purifying selection than the large subunit.</p> <p>Consistent with the numbers of large and small subunit genes found in the sequenced plant genomes [##REF##17496118##2##,##REF##15821022##35##,##REF##15598655##38##], reconciled tree analyses indicated that large subunit genes underwent more duplications than small subunit genes. Both phylogeny and molecular clock analyses indicate that the initial duplication in the large subunit family of angiosperms occurred <italic>ca</italic>. 400 MY ago, close to the divergence of angiosperms and bryophytes [##UREF##4##30##] (Figure ##FIG##2##3A## and Additional file ##SUPPL##3##4A##). In contrast, the oldest retained small subunit duplicates date back to 120–140 MY ago, after the divergence of monocots and eudicots (Figure ##FIG##2##3B## and Additional file ##SUPPL##3##4B##). The reason why the large subunit had more ancient and duplications than did the small subunit remains enigmatic. Macroevolutionary models that can be applied to gene families are poorly developed, so it remains possible that the observed difference is coincidental. Alternatively, the large subunit might have a greater ability to undergo subfunctionalization after duplication. The fact that 7 large subunit genes and 4 small subunit genes can be identified in the <italic>Physcomitrella </italic>genome suggests similar patterns in both mosses and angiosperms. However, more rigorous tests to distinguish among these scenarios must await both the acquisition of additional data from additional deep-branching land-plant lineages (e.g., liverworts and hornworts) and the development of better models of gene family macroevolution.</p> <p>Studies of expression patterns of AGPase genes in several species, including rice, <italic>Arabidopsis</italic>, potato, tomato and barley, have shown that the large subunit is tissue specific while the small subunit is more broadly expressed [##REF##7700228##33##, ####REF##9469935##34##, ##REF##15821022##35##, ##REF##16275672##36##, ##REF##16957017##37##, ##REF##15598655##38####15598655##38##]. Based on these studies, the major large subunit groups (Figure ##FIG##0##1B##) are likely to be expressed in different tissues in most or all plants. The tissue-specificity of large subunit genes suggests the expression patterns of these genes might undergo subfunctionalization after duplication, as predicted by the "DDC model" of gene duplication [##REF##10101175##53##]. The DDC model predicts that duplicated genes are preserved by complementary changes in their expression pattern (e.g., a broadly expressed gene might undergo duplication and have one duplicate expressed in a specific tissue like leaves while the other duplicate is expressed elsewhere). Although the potential for subfunctionalization due to changes in gene expression to preserve duplicated genes is generally accepted [##REF##16140417##54##,##REF##16280546##55##], it also remains possible that distinct AGPase genes have specialized in terms of protein function (e.g., their pH optima might have shifted based upon the specific tissues in which the paralogs are expressed). The relative contributions of subfunctionalization and specialization or neofunctionalization to gene family evolution are open questions [##REF##17988397##56##,##REF##15654095##57##], and it is unclear that there is any reason why the large subunit would be more likely to undergo specialization at either the protein or gene expression level. Such a model, which postulates that subfunctionalization of gene expression was followed by specialization, is similar to a combined model called subneofunctionalization [##REF##15654095##57##]. The subneofunctionalization model postulates that subfunctionalization occurs shortly after duplication while neofunctionalization is a more prolonged process [##REF##17988397##56##,##REF##15654095##57##]. If the combined model were applied to AGPases, the initial preservation of paralogous AGPase genes immediately after duplication might reflect subfunctionalization but this process would be followed by adaptation to the more specialized domains of expression in which each of the paralogs are expressed. Either a neofunctionalization or a subneofunctionalization model would be consistent with the evidence that different large subunit groups have functionally diverged from each other at the protein level (Figure ##FIG##3##4##, Table ##TAB##1##2##). However, subneofunctionalization provides a means to directly link the divergence of expression patterns and divergence at the protein level. Corroborating the subneofunctionalization will require correlating gene expression and sequence divergence for a large number of plants.</p> <p>Differences in their patterns of expression represent the major difference between the large and small subunits that could explain the differences in their rates of evolution, since broadly expressed genes are more conserved than tissue-specific genes [##REF##11682310##58##, ####REF##10666707##59##, ##REF##15201397##60####15201397##60##]. However, this raises the question of why broadly expressed genes, like AGPase small subunit genes, exhibit slower rates of evolution. Although it may simply be that mutations in broadly expressed genes have a greater impact on fitness or because these genes have to function in multiple cellular environments [##REF##11682310##58##], a simpler explanation might be that small subunit genes have to function with multiple large subunit genes. Georgelis et al. [##REF##17496118##2##] presented data consistent with this possibility, since they showed that the effects of several amino acid changes in the maize endosperm small subunit on enzyme activity depended on the identity of the large subunit [maize endosperm large subunit (SHRUNKEN-2)(SH2) and maize embryo large subunit (AGPLEMZM) were used]. Both SH2 and AGPLEMZM are members of group 3b (Figure ##FIG##0##1B##), so these results suggest that even fairly similar large subunit genes can interact differently with small subunits.</p> <p>In addition to the potential for subfunctionalization due to changes in gene expression, the observation that estimates of θ<sub>I </sub>for large subunit groups were significantly greater than zero suggests that it will be possible to attribute differences among AGPase genes to specific amino acid changes. We found 99 candidate sites of the large subunit likely to have been involved in rate shifts (either type-I or -II divergence; Additional files ##SUPPL##7##8## and ##SUPPL##8##9##). At least some of these putative rate shift residues are likely to have contributed to functional changes among the different large subunit groups. The estimate of θ<sub>I </sub>for the small subunit groups 1 (monocot) and 2 (eudicot) is also significantly greater than zero. It was possible to find evidence for 13 type-I candidate sites in the small subunit alignment. Like the large subunit, the estimate of θ<sub>II </sub>for the small subunit was not significantly greater than zero. Nonetheless, there were two potential type-II sites could be identified (Additional file ##SUPPL##8##9##).</p> <p>A total of 21 candidate sites for positive selection could be identified in the large subunit branches following duplications that led to different groups (Branch numbers: 4,6,9,11,12,13,16 shown in Additional file ##SUPPL##5##6##), and six of these sites overlapped with the set of sites that appear to have undergone either type-I (sites 341, 364, 445) or type-II (sites 106, 114, 382) divergence. Biochemical and genetic studies confirm that at least one of the sites (sites 106) is important for AGPase activity. This site is a threonine (T) in large subunit groups 3a and 3b but a lysine (K) in groups 1 and 2 and in all small subunits. The potato tuber large subunit, which falls into group 3a, has a T at site 106 and it forms an inactive complex if it is combined with an inactivated potato tuber small subunit [##REF##15632142##61##]. Changing this T in the potato tuber large subunit to a K actually results (the T106K mutant) in a complex with some activity with the same inactivated potato tuber small subunit [##REF##15632142##61##]. These results were interpreted as evidence that the large subunit lost its catalytic ability partly because of the K to T change. Although the K residue at site 106 may be necessary for catalysis if the small subunit is inactive, another model that explains the data would be one in which the wild-type large subunits (which have a T at site 106) require prior catalysis by the small subunit before they perform catalysis by themselves. Such a catalytic mechanism has been proposed for the <italic>Escherichia coli </italic>AGPase [##REF##363717##62##]. The T residue at site 106 is absolutely conserved in large subunit groups 3 and the branch-site model suggests positive selection for the T immediately following the duplication that generated group 3 (Additional file ##SUPPL##5##6##). This suggests that this change was important for enzymatic activity and beneficial for the plants. Indeed, the overall activity of a complex that includes the T106K mutant of the potato tuber large subunit and the wild-type potato tuber small subunit showed significantly reduced activity relative to wild-type potato tuber AGPase [##REF##17137579##63##].</p> <p>One potential type-II site (site 507) and four sites that represent candidates for positive selection on branches that immediately follow duplications (sites 104, 230, 441, 445) have been shown to be important for the allosteric properties of AGPase [[##REF##15692569##49##,##REF##9733546##64##,##REF##11524424##65##], Hannah personal communication]. However, most of the candidate sites for either rate shifts or positive selection do not have a known function.</p> <p>The existence of type-I and type-II divergence among AGPase subunit groups along with the detection of positively selected sites after duplications that led to different groups in the large subunit provide evidence for functional divergence especially among the large subunit groups. Our data are consistent with biochemical studies showing that the four possible AGPase complexes in <italic>Arabidopsis</italic>, which have a single functional small subunit gene and four distinct large subunit genes (belonging to different groups in Figure ##FIG##0##1B##), have different kinetic and allosteric properties [##REF##12748181##16##]. There is further evidence for functional divergence among plant AGPases, since the maize and barley endosperm AGPases are less dependent than potato tuber AGPase on the allosteric activator 3-PGA for activity and the maize endosperm AGPase is more heat labile than potato tuber AGPase [##REF##14520572##17##, ####REF##10557246##18##, ##REF##16299180##19##, ##REF##18024561##20####18024561##20##]. Functional divergence among the different subunit groups was also suggested by Georgelis et al. [##REF##17496118##2##], who showed that all groups of large subunit genes have ω values (which range from 0.073 to 0.132) that exceed those for small subunit genes (which range from 0.027 to 0.054). These results are consistent with the rate shifts within the large and the small subunit families that were observed in the present study and they further imply that the various groups of plant AGPases have undergone functional divergence. The present study also identifies specific residues that are likely to have contributed towards that divergence. Site directed mutagenesis of these candidate sites is likely to shed some light on their functions and reveal the proportion of candidate sites that reflect type II error for tests to identify sites subject to rate shifts and positive selection.</p> <p>A number of angiosperm AGPases have been successfully expressed in <italic>E. coli </italic>and purified, including maize endosperm AGPase [##REF##8650177##66##], potato tuber AGPase [##REF##8380404##67##], and all possible <italic>Arabidopsis </italic>AGPase complexes [##REF##12748181##16##]. The most straightforward candidate sites to test are the type-II sites and those subject to positive selection, where it is possible to change the relevant residue either to the amino acid present in the other group (for type-II sites) or the ancestral amino acid (for positively selected sites). Testing type-I sites may be more challenging, since there is not a clear way to swap the residues present in a member of the focal group of proteins with that in a different group. However, it is reasonable to predict that any of the residues present in the group of proteins for which the site is variable should alter the biochemical activity of a protein in which the site is invariant. Regardless of the specific strategies for generating mutants, mutant large or small subunits could be expressed in <italic>E. coli </italic>along with a wild-type version of the other subunit, the relevant complex purified, and the properties of the enzyme determined to allow the impact of the mutations to be studied. This will allow elucidation of the importance of the sites in enzyme activity. Ultimately, it will be interesting to determine whether these sites are important for the kinetic and allosteric properties of AGPase, for enzyme stability for the pH optimum, or for multiple properties.</p>
<title>Conclusion</title> <p>Herein, we validated and extended the observation, initially based upon estimates of ω [##REF##17496118##2##], that the AGPase large subunit accumulated non-synonymous substitutions more rapidly than the small subunit in angiosperms by estimating absolute rates of amino acid change. The earliest duplication in the large subunit family of angiosperms was close to the time that angiosperms and mosses diverged (~400 MY ago;[##UREF##4##30##]). The large subunit underwent a larger number of duplications than the small subunit, which only began to duplicate after the divergence of monocots and eudicots. We suggest that the large subunit evolved faster due to permanently relaxed constraints since positive selection and sporadic episodes of relaxed constraints cannot account for the different rates of evolution between the large and the small subunit.</p> <p>Large subunit genes exhibit narrower tissue specificity than small subunit genes in terms of their gene expression patterns [##REF##7700228##33##, ####REF##9469935##34##, ##REF##15821022##35##, ##REF##16275672##36##, ##REF##16957017##37##, ##REF##15598655##38####15598655##38##], and they are likely to have experienced subfunctionalization in terms of expression patterns. However, we use analyses of rate shifts and positive selection to demonstrate that different groups of both large and small subunits are likely to have diverged at the protein level. We have identified candidate amino acid sites with the potential to account for the functional divergence and described strategies for site-directed mutagenesis experiments that could shed light into the specific roles of these sites.</p>
<title>Background</title> <p>ADP-glucose pyrophosphorylase (AGPase), which catalyses a rate limiting step in starch synthesis, is a heterotetramer comprised of two identical large and two identical small subunits in plants. Although the large and small subunits are equally sensitive to activity-altering amino acid changes when expressed in a bacterial system, the overall rate of non-synonymous evolution is ~2.7-fold greater for the large subunit than for the small subunit. Herein, we examine the basis for their different rates of evolution, the number of duplications in both large and small subunit genes and document changes in the patterns of AGPase evolution over time.</p> <title>Results</title> <p>We found that the first duplication in the AGPase large subunit family occurred early in the history of land plants, while the earliest small subunit duplication occurred after the divergence of monocots and eudicots. The large subunit also had a larger number of gene duplications than did the small subunit. The ancient duplications in the large subunit family raise concern about the saturation of synonymous substitutions, but estimates of the absolute rate of AGPase evolution were highly correlated with estimates of ω (the non-synonymous to synonymous rate ratio). Both subunits showed evidence for positive selection and relaxation of purifying selection after duplication, but these phenomena could not explain the different evolutionary rates of the two subunits. Instead, evolutionary constraints appear to be permanently relaxed for the large subunit relative to the small subunit. Both subunits exhibit branch-specific patterns of rate variation among sites.</p> <title>Conclusion</title> <p>These analyses indicate that the higher evolutionary rate of the plant AGPase large subunit reflects permanent relaxation of constraints relative to the small subunit and they show that the large subunit genes have undergone more gene duplications than small subunit genes. Candidate sites potentially responsible for functional divergence within each of the AGPase subunits were investigated by examining branch-specific patterns of rate variation. We discuss the phenotypes of mutants that alter some candidate sites and strategies for examining candidate sites of presently unknown function.</p>
<title>Authors' contributions</title> <p>NG conducted the experiments and carried out the analyses. NG, ELB and LCH conceived and designed the experiments. NG, ELB and LCH wrote the manuscript. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank Rebecca Kimball and members of the Kimball, Braun and Hannah laboratories for many useful comments and discussions. This research was supported by National Science Foundation Grants IOB-0444031, DBI-0077676, DBI-0606607, and IOB-9982626 and USDA Grant 2006-35100-17220.</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>Reconciled large and small subunit trees</bold>. A) Amino acid tree of the large and small subunits from angiosperms, <italic>Physcomitrella patens </italic>and <italic>Chlamydomonas reinhardtii</italic>. The topology of the tree was determined by ML using aligned amino acid sequences using PhyML. Branch lengths reflect numbers of amino acid substitutions per site. The branches within groups have been replaced by the grey triangles. ML bootstrap values are indicated above branches, and the bar shows the number of amino acid substitutions per site. B) Angiosperm large subunit reconciled tree. C) Angiosperm small subunit reconciled tree. The topology of the trees shown in B) and C) was determined by ML using aligned cDNA sequences using GARLI. Branch lengths reflect numbers of amino acid substitutions per site as estimated by AAML and the scale bar shows the number of amino acid substitutions per site. ML bootstrap values &gt; 50% are indicated above branches. Reconciled tree analyses (using the gene trees shown and the species tree shown in Additional file ##SUPPL##0##1##) were conducted using GENETREE. Black boxes at nodes indicate duplication events. The arrow in B) indicates the divergence of <italic>Physcomitrella patens </italic>from angiosperms. The trees in B) and C) were rooted with the AGPase large and small subunit from <italic>Chlamydomonas reinhardtii </italic>respectively. Thicker lines indicate branches that follow duplication events and have K<sub>S </sub>&lt; 0.1 (based upon ML estimates of synonymous branch lengths).</p></caption></fig> <fig id="F2" position="float"><label>Figure 2</label><caption><p><bold>Absolute rate of evolution of the large and the small subunit of AGPase from angiosperms (measured in aass MY<sup>-1</sup>)</bold>. The blue bars indicate the average rates (amino acid substitutions per site per million years; aass MY<sup>-1</sup>) estimated from the most recent dated speciation events to present sequences of the trees shown in Figure 1B and 1C. The red bars indicate the average aass MY<sup>-1 </sup>estimated from the most recent dated speciation events to present sequences of the trees shown in Additional file ##SUPPL##2##3A## and ##SUPPL##1##2B##. The yellow bars indicate the average aass MY<sup>-1 </sup>estimated from all branches in Figure 3A and 3B. The green bars indicate the average aass MY<sup>-1 </sup>estimated from all branches in Additional file ##SUPPL##3##4A## and ##SUPPL##2##3B##. The error bars indicate 2× standard error.</p></caption></fig> <fig id="F3" position="float"><label>Figure 3</label><caption><p><bold>Phylogenetic trees of the large and the small subunits from angiosperms after rate-smoothing</bold>. The trees in A) and B) are the trees shown in Figure 1B and 1C respectively after rate-smoothing. Rate-smoothing was done by using penalized likelihood as implemented in the r8s software [##REF##12538260##79##]. Black boxes indicate duplication events.</p></caption></fig> <fig id="F4" position="float"><label>Figure 4</label><caption><p><bold>Positive selection in the large and the small subunit of the angiosperms</bold>. The trees shown in A) and B) have the same topology as the trees shown in Figure 1B and 1C, although the shown here are unrooted. White circles indicate branches where positive selection was detected only by Test 1 but not Test 2 (described in Methods). Black circles indicate branches where positive selection was detected by both Test 1 and Test 2. The branches where positive selection is detected are numbered.</p></caption></fig>
<table-wrap id="T1" position="float"><label>Table 1</label><caption><p>Temporary relaxation of purifying selection, positive selection, or both after duplications in the large and small subunit of AGPase from angiosperms.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">Branches</th><th align="center">Large subunit</th><th align="center">p-value</th><th align="center">Small subunit</th><th align="center">p-value</th></tr></thead><tbody><tr><td align="center">Class 1 (Figure 1)</td><td align="center">0.114</td><td align="center">0.005</td><td align="center">0.058</td><td align="center">&lt; 0.001</td></tr><tr><td align="center">Class 2 (Figure 1)</td><td align="center">0.086</td><td/><td align="center">0.029</td><td/></tr><tr><td align="center">All branches (Figure 1)</td><td align="center">0.090</td><td/><td align="center">0.033</td><td/></tr><tr><td/><td/><td/><td/><td/></tr><tr><td align="center">Class 1 (Additional file ##SUPPL##2##3##)</td><td align="center">0.125</td><td align="center">&lt; 0.001</td><td align="center">0.081</td><td align="center">&lt; 0.001</td></tr><tr><td align="center">Class 2 (Additional file ##SUPPL##2##3##)</td><td align="center">0.086</td><td/><td align="center">0.029</td><td/></tr><tr><td align="center">All branches (Additional file ##SUPPL##2##3##)</td><td align="center">0.089</td><td/><td align="center">0.034</td><td/></tr></tbody></table></table-wrap> <table-wrap id="T2" position="float"><label>Table 2</label><caption><p>Type I and II functional divergence between large and small subunit groups.</p></caption><table frame="hsides" rules="groups"><thead><tr><th/><th align="center" colspan="6">Large subunit</th><th align="center">Small subunit</th></tr><tr><th/><th colspan="7"><hr/></th></tr><tr><th/><th align="center">Group 3b/</th><th align="center">Group 3b/</th><th align="center">Group 3b/</th><th align="center">Group 3a/</th><th align="center">Group 3a/</th><th align="center">Group 2/</th><th align="center">Group 1/</th></tr><tr><th/><th align="center">Group 3a</th><th align="center">Group 2</th><th align="center">Group 1</th><th align="center">Group 2</th><th align="center">Group 1</th><th align="center">Group 1</th><th align="center">Group 2</th></tr></thead><tbody><tr><td/><td/><td/><td align="center" colspan="2">Type I divergence</td><td/><td/><td/></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="center">Theta-I</td><td align="center">0.261*</td><td align="center">0.383*</td><td align="center">0.448*</td><td align="center">0.220*</td><td align="center">0.331*</td><td align="center">0.273*</td><td align="center">0.263*</td></tr><tr><td align="center">SE</td><td align="center">0.069</td><td align="center">0.095</td><td align="center">0.097</td><td align="center">0.060</td><td align="center">0.072</td><td align="center">0.076</td><td align="center">0.084</td></tr><tr><td colspan="8"><hr/></td></tr><tr><td/><td/><td/><td align="center" colspan="2">TypeII divergence</td><td/><td/><td/></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="center">Theta-II</td><td align="center">0.016</td><td align="center">0.072</td><td align="center">0.080</td><td align="center">0.016</td><td align="center">0.058</td><td align="center">0.001</td><td align="center">0.028</td></tr><tr><td align="center">SE</td><td align="center">0.061</td><td align="center">0.052</td><td align="center">0.054</td><td align="center">0.059</td><td align="center">0.061</td><td align="center">0.053</td><td align="center">0.033</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>Alignment of large and small subunits AGPases from angiosperms with protein domains highlighted</bold>. The blue domain indicates the hypervariable N terminus of the large and the small subunit. The pink and green domains indicate the catalytic domain and the β-helix domain respectively. The yellow domain indicates the loop that connects the catalytic to the β-helix domain.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p><bold>Species tree</bold>. Times of divergence are indicated in million of years (MY) at nodes [##REF##8896372##83##, ####REF##10759555##84##, ##REF##10908680##85##, ##REF##11674868##86##, ##UREF##13##87##, ##REF##15489274##88##, ##REF##15499401##89##, ##UREF##14##90##, ##REF##16677401##91##, ##REF##16632644##92####16632644##92##]. All divergence times were examined for consistency with the fossil record [##REF##21652316##93##].</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p><bold>Reconciled large and small subunit trees</bold>. A) Angiosperm large subunit reconciled tree. B) Angiosperm small subunit reconciled tree. The topology of the trees shown in A) and B) was determined by ML using aligned cDNA sequences analyzed by GARLI. Nodes with bootstrap values &lt; 70% (Figure ##FIG##0##1##) were then rearranged to minimize the number of duplications (to increase congruence with the species tree). Branch lengths reflect numbers of amino acid substitutions per site, estimated AAML (with the scale bar showing the number of amino acid substitutions per site). Reconciled tree analyses were conducted using GENETREE and the species tree in Additional file 1. Black boxes indicate duplication events. The trees in A) and B) were rooted with the AGPase large and small subunit from <italic>Chlamydomonas reinhardtii </italic>respectively. Thicker lines indicate branches with K<sub>S </sub>&lt; 0.1 following duplication events (using ML estimates of synonymous branch lengths).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p><bold>Phylogenetic trees of the large and the small subunits from angiosperms after rate-smoothing</bold>. The trees in parts A) and B) of this figure are rate-smoothed versions of the gene trees shown in Additional file 3A and 3B that were rearranged to increase congruence with the species tree. Rate-smoothing was done by using the PL method implemented in the r8s software. Black boxes indicate duplication events.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p><bold>Average number of synonymous substitutions per site per year</bold>. The trees in A) and B) have the topology of the trees shown in Figure ##FIG##0##1B## and ##FIG##0##1C## respectively. The length of the branches represents the number of synonymous substitutions per site as estimated by the free model of CODEML. The bars correspond to the number of synonymous substitutions per site. The numbers of synonymous substitutions per site per year, shown in C), were estimated from the most recent dated speciation events to present sequences of the trees shown in A) and B). The error bars represent 2× standard error.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p><bold>Amino acid sites in the large and the small subunit of AGPase from angiosperms under positive selection</bold>. Large subunit site numbers correspond to the amino acid sequence encoded by <italic>Shrunken-2 </italic>(NCBI accession number: P55241). Small subunit site numbers correspond to the amino acid sequence encoded by <italic>Brittle-2 </italic>(NCBI accession number: AAQ14870).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p><bold>Distribution of type I sites along the large (A) and the small (B) subunit</bold>. The cut-off value of posterior probability is empirical and it was set to 0.5 for all group comparisons except for group 1-group 3b and group 2-group 3b where the cut-off value was set to 0.6, since theta was greater for these pairs. The Y-axis corresponds to posterior probability. The X-axis corresponds to the number of the amino acid site based on the subunits encoded by <italic>Shrunken-2 </italic>(A) (NCBI accession number: P55241) and <italic>Brittle-2 </italic>(B) (NCBI accession number: AAQ14870).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S8"><caption><title>Additional file 8</title><p><bold>Type-I sites in the large and the small subunit of AGPase from angiosperms</bold>. Type-I functional divergence between large and small subunit groups was estimated by DIVERGE. Large subunit site numbers correspond to the amino acid sequence encoded by <italic>Shrunken-2 </italic>(NCBI accession number: P55241). Small subunit site numbers correspond to the amino acid sequence encoded by <italic>Brittle-2 </italic>(NCBI accession number: AAQ14870).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S9"><caption><title>Additional file 9</title><p><bold>Type-II sites in the large and the small subunit of AGPase from angiosperms</bold>. Type-II functional divergence between large and small subunit groups was estimated by DIVERGE. Large subunit site numbers correspond to the amino acid sequence encoded by <italic>Shrunken-2 </italic>(NCBI accession number: P55241). Small subunit site numbers correspond to the amino acid sequence encoded by <italic>Brittle-2 </italic>(NCBI accession number: AAQ14870).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S10"><caption><title>Additional file 10</title><p>AGPase subunit accession numbers</p></caption></supplementary-material>
<table-wrap-foot><p>The branches of the trees shown in Figure 1B and Additional file ##SUPPL##2##3A## (large subunit) or Figure 1C and Additional file ##SUPPL##2##3B## (small subunit) were classified into two groups. One group includes the short branches (K<sub>S </sub>&lt; 0.1) following a duplication event (post-duplication; class 1) and the other group includes all the other branches (class 2). The overall ω value of each group and the overall ω value of all branches were estimated by CODEML. The LRT was used to compare the model with two different ω values (one for class1 branches and one for class 2 branches) to the null hypothesis model with a single ω value.</p></table-wrap-foot> <table-wrap-foot><p>Coefficients of type I (θ<sub>I</sub>) and II (θ<sub>II</sub>) functional divergence were estimated using DIVERGE. * denotes statistical significance at 5% level of confidence. SE: standard error.</p></table-wrap-foot>
<graphic xlink:href="1471-2148-8-232-1"/> <graphic xlink:href="1471-2148-8-232-2"/> <graphic xlink:href="1471-2148-8-232-3"/> <graphic xlink:href="1471-2148-8-232-4"/>
<media xlink:href="1471-2148-8-232-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-232-S2.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-232-S3.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-232-S4.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-232-S5.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-232-S6.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-232-S7.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-232-S8.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-232-S9.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-232-S10.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Gaut"], "given-names": ["BS"], "person-group": ["Hecht MK"], "article-title": ["Molecular clocks and nucleotide substitution rates in higher plants"], "source": ["Evolutionary Biology"], "year": ["1998"], "publisher-name": ["New York, Plenum Press"], "fpage": ["93"], "lpage": ["120"]}, {"surname": ["Mayrose", "Doron-Faigenboim", "Bacharach", "Pupko"], "given-names": ["I", "A", "E", "T"], "article-title": ["Towards realistic codon models: among site variability and dependency of synonymous and non-synonymous rates"], "source": ["Bioinformatics"], "year": ["2007"], "volume": ["23"], "fpage": ["319"], "lpage": ["327"], "pub-id": ["10.1093/bioinformatics/btm176"]}, {"surname": ["Ohno"], "given-names": ["S"], "source": ["Evolution by Gene Duplication"], "year": ["1970"], "publisher-name": ["New York, Springer"]}, {"surname": ["Zhang"], "given-names": ["J"], "article-title": ["Evolution by gene duplication: an update"], "source": ["Trends Ecol Evol"], "year": ["2003"], "volume": ["18"], "fpage": ["292"], "lpage": ["298"], "pub-id": ["10.1016/S0169-5347(03)00033-8"]}, {"surname": ["Kenrick", "Crane"], "given-names": ["P", "PR"], "article-title": ["The origin and early evolution of plants on land"], "source": ["Nature"], "year": ["1997"], "volume": ["389"], "fpage": ["33"], "lpage": ["39"], "pub-id": ["10.1038/37918"]}, {"surname": ["Goodman", "Czelusniak", "Moore", "Romero-Herrera", "Matsuda"], "given-names": ["M", "J", "GW", "AE", "G"], "article-title": ["Fitting the gene lineage into its species lineage: a parsimony strategy illustrated by cladograms constructed from globin sequences"], "source": ["Syst Zool"], "year": ["1979"], "volume": ["28"], "fpage": ["132"], "lpage": ["168"], "pub-id": ["10.2307/2412519"]}, {"surname": ["Yang", "Goldman", "Friday"], "given-names": ["Z", "N", "A"], "article-title": ["Maximum likelihood trees from DNA sequences: a peculiar statistical estimation problem"], "source": ["Syst Biol"], "year": ["1995"], "volume": ["44"], "fpage": ["384"], "lpage": ["399"], "pub-id": ["10.2307/2413599"]}, {"surname": ["Kumar", "Tamura", "Nei"], "given-names": ["S", "K", "M"], "article-title": ["MEGA 3: Integrated Software for Molecular Evolutionary Genetics Analysis and Sequence Alignment"], "source": ["Brief Bioinf"], "year": ["2004"], "volume": ["5"], "fpage": ["150"], "lpage": ["163"], "pub-id": ["10.1093/bib/5.2.150"]}, {"surname": ["Zwickl"], "given-names": ["DJ"], "article-title": ["Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion"], "source": ["Ph.D. dissertation"], "year": ["2006"], "publisher-name": ["The University of Texas at Austin"]}, {"surname": ["Felsenstein"], "given-names": ["J"], "article-title": ["Confidence-limits on phylogenies: An approach using the bootstrap"], "source": ["Evolution Int J Org Evolution"], "year": ["1985"], "volume": ["39"], "fpage": ["783"], "lpage": ["791"]}, {"surname": ["Dayhoff", "Schwartz", "Orcutt"], "given-names": ["MO", "RM", "BC"], "person-group": ["Dayhoff MO"], "article-title": ["A model of evolutionary change in proteins"], "source": ["Atlas of Protein Sequence and Structure"], "year": ["1978"], "volume": ["5"], "publisher-name": ["Silver Springs, MD, National Biomedical Research Foundation"], "fpage": ["345"], "lpage": ["352"]}, {"surname": ["Page", "Cotton"], "given-names": ["RDM", "JA"], "person-group": ["Sankoff D, Nadeau J"], "article-title": ["GeneTree: a tool for exploring gene family evolution"], "source": ["Comparative Genomics: Empirical and Analytical Approaches to Gene Order Dynamics, Map Alignment, and the Evolution of Gene Families"], "year": ["2000"], "publisher-name": ["Dordrecht, Kluwer Academic Press"], "fpage": ["525"], "lpage": ["536"]}, {"surname": ["Hillis", "Bull"], "given-names": ["DM", "JJ"], "article-title": ["An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis"], "source": ["Syst Biol"], "year": ["1993"], "volume": ["42"], "fpage": ["182"], "lpage": ["192"], "pub-id": ["10.2307/2992540"]}, {"surname": ["Gaut"], "given-names": ["BS"], "article-title": ["Evolutionary dynamics of grass genomes"], "source": ["New phytologist"], "year": ["2002"], "volume": ["154"], "fpage": ["15"], "lpage": ["28"], "pub-id": ["10.1046/j.1469-8137.2002.00352.x"]}, {"surname": ["Soltis", "Soltis", "Edwards"], "given-names": ["P", "D", "C"], "article-title": ["Angiosperms, Flowering Plants"], "ext-link": ["http://tolweb.org/Angiosperms/20646/2005.06.03", "http://tolweb.org/"], "comment": ["in The Tree of Life Web Project, "]}]
BMC Evol Biol. 2008 Aug 12; 8:232