Volume 14 Supplement 15
The dynamics of functional classes of plant genes in rediploidized ancient polyploids
© Chen et al.; licensee BioMed Central Ltd. 2013
Published: 15 October 2013
To understand the particular evolutionary patterns of plant genomes, there is a need to systematically survey the fate of the subgenomes of polyploids fixed as whole genome duplicates, including patterns of retention of duplicate, triplicate, etc. genes.
We measure the simultaneous dynamics of duplicate orthologous gene loss in rosids, in asterids, and in monocots, as influenced by biological functional class. This pan-angiosperm view confirms common tendencies and consistency through time for both ancient and more recent whole genome polyploidization events.
The gene loss analysis represents an assessment of post-polyploidization evolution, at the level of individual gene families within and across sister genomes. Functional analysis confirms universal trends previously reported for more recent plant polyploidy events: genes involved with regulation and responses were retained in multiple copies, while genes involved with metabolic and catalytic processes tended to lose copies, across all three groups of plants.
In whole genome duplication (WGD) or triplication, the entire gene and chromosome structure of a genome undergoes polyploidization followed by rediploidization of the new larger genome and fractionation, or homeologous gene loss, of many or most of the duplicate genes. The doubling or tripling of all the gene contents of a genome has been hypothesized as an important source of gene innovations and the radiation of species. The inherent variability in these processes may create a variety of beneficial phenotypes similar to those seen due to heterosis , and the increase in the diversity of genetic elements may help drive long-term morphological complexity and adaptation to new environments .
Studies in a variety of organisms have provided evidence that gene retention during fractionation may differ among functional categories. This has given rise to a number of explanatory models, particularly a reformulation of the classical genetic Gene Balance Hypothesis , all of which attempt to explain functional bias in gene retention after WGD.
These results, such as the relatively rapid fractionation of genes labelled "metabolic", and the resistance to fractionation of those labelled "biological regulation", confirmed various accounts in the literature [8–11].
This work raised new questions. The first is to what extent the patterns in Figure 1 are specific to rosids. Which functional categories have the same tendencies in other flowering plant groups? The second concern has to do with the severe restriction on sets of orthologs and paralogs we included in our data - we required that each rosid manifest at least one copy of the gene in a common syntenic context when compared to the other rosids, and the same for paralogs within each genome. This requirement attempted to isolate the process of fractionation from other processes of gene family dynamics, such as gene movement and gene family expansion, by focusing on the "natural experiment" created by the core eudicot triplication, whereby each species started with three copies of the same gene in the three identical syntenic contexts, and each then retained either one, two or three of these copies over time. A working assumption was that fractionation would be the dominant process affecting these particular sets of genes. The requirement, however, reduced the total number of "homology sets" considered to about a third of the number of genes in the extant genomes. Do our results therefore pertain only to genes that are syntenically conserved and comprise the "stable genome" for functional reasons, or are they also valid for some additional class of paralog sets, some of whose elements are not detectible using syntenic criteria? The third and final problem is whether the results we traced on high-level functional categories indicate general tendencies within these categories, or simply reflect the preponderance of some subcategories among many others with diverse, and perhaps contrary fractionation and positional stability behaviour.
In this paper we address these three concerns. First, we replicate our original studies on two other groups of plants, one consisting of three lamiid asterid species: tomato (Solanum lycopersicum) , humped bladderwort (Utricularia gibba) , and monkey flower (Mimulus guttatus) , the other, four Poaceae monocot species: rice (Oryza saliva) , foxtail millet (Setaria italica) , sorghum (Sorghum bicolor) , and purple false brome (Brachypodium distachyon) , thus compiling the first pan-angiosperm study focused on gene fractionation patterns. The asterids have all individually undergone further WGD since the core eudicot whole genome triplication they share with the rosids. Nevertheless, the principles underlying functionally-influenced fractionation can be assumed to hold just as well after several WGD as after one. The divergence of the asterids from the rosids seems to have occurred within five or ten million years after the triplication event they shared approximately 120 million years ago , so that 90 - 95% of their evolution, including much of their fractionation, took place independently in the two groups. The monocot species in our study all share two or three WGD independent of the all polyploidy events in the eudicots. We will show that the independent fractionation patterns of all three diverse groups of angiosperms are highly parallel, following the same patterns of fractionation based on gene functional classes.
Second, we investigate the hypothesis that our results only apply to a non-mobile core of genes that are detected by the syntenic context they share with orthologs in all the genomes in the group. We relax somewhat the requirement that all the extant genomes must contain at least one syntenically validated descendant of the ancestral gene. Even though the additional cohort of genes, part of the more "mobile" genome that is prone to translocation in a genome, is less numerous and thus less conducive to statistical significance of the results, we find that they continue to support the tendencies found for the functionally "stabilized" part of the genome. We demonstrate that the homology sets we study are enriched for some categories and depleted for others, in comparison with random samples of genes from the entire extant genomes, but hypothesize that this has more to do with lineage-specific expansion of gene families, rather than fractionation dynamics.
Finally, we search for tendencies within subcategories of some of the large functional categories. We find that the negative slope within the catalytic activity category reflects in large measure the consistent behavior across certain classes of enzymes. We resolve the apparent conflict between the trends for the high-level terms "organelle" and "organelle part" seen in Figure 1 by de-convoluting lower level terms contained therein.
Methods and data
We compared the retention of homologs in six core eudicot species, namely three rosids: peach, cacao and grape, and three asterids: tomato, Utricularia and Mimulus, and four monocot species: rice, Setaria, sorghum and Brachypodium, forming three independent data sets.
The homology sets were first examined to see whether they contained at least one gene from each species in the group. In the first analysis, all sets with no gene in any of these species were excluded from the analysis. (In a later analysis, described below, other homology sets were used.) The remaining sets were classified according to the number of species in which there was more than one copy, so that in the three-species comparisons, the sets could be classified as 0, 1, 2, or 3, and in the four-species set a score of 4 was also possible. We call this number the fractionation score.
For each homology set, each of its genes was annotated by submitting it to Blast2GO . Then all the annotations from all the genes in this set were considered as annotations for the set as a whole. No account was taken of the multiplicity of "hits" of a single annotation within the set. Of course, for every annotation, each of the higher-level terms of each hit was also counted as an annotation.
Among all the homology sets we constructed, approximately 90% hit at least one GO term, resulting in 10,688 monocot, 6360 rosid and 4638 asterid homology sets for further analysis.
The GO terms are divided at the highest level into "Biological Process", "Molecular Function" and "Cellular Component" and there are a further 67 terms at the next level, which we call "high-level terms". Homology sets with large fractionation scores, i.e., which contain more than one paralog in all or most genomes, tend to have a higher total number of annotations, simply by virtue of having a larger number of genes. This leads to the artifactual observation that almost all functional categories are more favored by homology sets with high fractionation scores. To correct for this bias, we use a normalized proportion of hits for each term for each fractionation score. This is calculated as the number of hits of the term over all homology sets with this fractionation score, divided by the total number of sets with hits for any terms within the appropriate highest-level term. Thus, if "organelle" received 100 hits in all homology sets with fractionation score 3, and if the number of sets hitting any "Cellular Component" term is 300, the normalized "proportion" is 33.3%.
These normalized proportions could then be plotted against fractionation score as in Figure 1. By considering every combination of homology set and functional category as a data point with X-coordinate its fractionation score and its Y-coordinate 1 or 0, depending on whether the homology set was a hit (1) or not (0) for that category, we could then calculate a regression score for the functional category. In Figure 1, the functional categories with significant negative slopes are black and those with significant positive score are red or orange.
Trends for high-level categories
Thus, the three patterns are surprisingly parallel. Recall that the monocot ancestor and its WGD occurred in a completely different part of the flowering plant phylogeny from the eudicots. And though the two groups of core eudicots descend from the same triplication event, the common ancestor of the asterid lamiids evolved for a lengthy period before radiating into the present-day families. Similarly for the common ancestor of peach and cacao, although grape belongs to an early branching rosid order.
The Gene Balance Hypothesis predicts that genes involved in multi-unit protein complexes, or genes involved in cascades in which the downstream genes are involved in multi-unit protein complexes, are more likely to be retained after whole genome duplication. This prediction can in turn predict what GO terms are expected to be fractionation resistant or fractionation prone. The predicted terms are largely the same as the significant in Figure 3, especially "binding", "biological regulation", "response to stimulus" and "developmental process".
Of the significant terms, two stand out as being fractionation prone: "metabolic process" and "catalytic activity". These two terms describe processes or functions that, as traditionally conceived, involve the interactions of single enzymes with substrates, coenzymes and cofactors that are not themselves proteins. This is in contrast with such processes as gene regulation, which explicitly involve the stoichiometry of two or more distinct proteins. On this basis, the Gene Balance Hypothesis would predict that "metabolic process" and "catalytic activity" be fractionation prone. There is much recent enthusiasm, however, for protein complexes in metabolic and other catalytic reactions, for metabolons , protein heteromers and other potential sources of abundance constraints among different proteins, so such predictions might have to be attenuated according to how widespread and how constraining these structures turn out to be (cf [25, 26] as examples of more moderate opinions).
The "mobile" genome
All homology sets involved in the above analysis contain at least one gene from each species To test whether this requirement biases the analysis towards genes with some special functional properties, we also carried out our analysis on homology sets where no ortholog was detected in one of the genomes. The absence of this gene from a synteny block does not necessarily mean that it is absent from the genome; it may have moved to some other location on the same or different chromosome.
Comparison of "stable" genome and the general gene complement
Trends for more specific categories
Finally, we search for tendencies within subcategories of some of the large functional categories. We find that the negative slope within the catalytic activity category reflects in large measure the consistent behaviour, across all three data sets, of isomerase, hydrolase, and some transferase subcategories, while other transferases notably the kinases and those transferring phosphorus-containing groups, as well as several subcategories of oxidoreductase activity are actually fractionation resistant. The details:
Oxidoreductase: No overall tendency, but 5 of 88 terms one level lower are fractionation resistant, significantly so in at least one of the data sets:.
- monooxygenase activity
oxidoreductase activity, acting on diphenols and related substances as donors
oxidoreductase activity, acting on X-H and Y-H to form an X-Y bond
arsenate reductase activity
oxidoreductase activity, acting on the aldehyde or oxo group of donors
Transferase: No overall tendency,
- (prone) transferase activity, transferring one-carbon groups
(prone) methyltransferase activity
(resistant) transferase activity, transferring phosphorus-containing groups
(resistant) kinase activity
(prone) nucleotidyltransferase activity
Hydrolase: Fractionation prone on the general enzyme class level and 2 of 21 lower terms. "Peptidase activity" is fractionation prone (in conflict with ).
(prone) hydrolase activity
(prone) peptidase activity
(prone) hydrolase activity, acting on carbon-nitrogen (but not peptide) bonds
Isomerase: Fractionation prone on the general enzyme class level and 2 of 14 lower terms.
- (prone) isomerase activity
(prone) cis-trans isomerase activity
(prone) intramolecular transferase activity
Lyase: Has conflicting slopes among the three data sets at the general enzyme class level. No fractionation prone or resistance from 13 lower terms.
Ligase: No tendency detected
The apparent conflict between the trends for "organelle" and "organelle part" turns out to be due to be the 95% concentration of the latter in chloroplast terms, where only 30% of the former are so annotated.
Connection between fractionation rates and and paralog retention patterns
Numbers of triples and pairs after fractionation in six rosids.
frequencies of gene family sizes
an original paralogy triplet would survive intact is .
an original triplet would manifest as a pair is , and
an original triplet would be reduced to a single copy is .
Because this model is inspired by concepts of the stable genome, the frequency of single-copy genes, which tend to be part of the mobile genome, are not part of the input data in Table 1, and so the original number of triplicates in the stable moiety, before fractionation, has to be inferred statistically.
We refrained from trying to model the annotation process, so the results in Figure 6 are more contrastive of the two population of genes than the real data. Nevertheless, it is encouraging that the empirical data in Table 1 contain the signal of two classes of genes with different fractionation rates, reminiscent of the patterns in Figures 1 and 3.
Building on our previous results on the functional determinants of fractionation of rosid genomes, we showed that the independent fractionation patterns of rosids, asterids and grasses are highly parallel. The same functional categories of genes are preferentially fractionation prone or resistant.
That the asterids demonstrate the same patterns as the rosids and monocots confirm that there is a general trend for paralogs to be fractionation prone or resistant, as influenced by their functional categories, despite the additional WGD events in the asterid lineages. This may be attenuated by changing fractionation patterns for some gene copies from one WGD to the next, as reported recently .
Additionally, we tested whether these patterns of fractionation were consistent between the stable and mobile portions of the genome. Even though the additional cohort of genes, part of the more "mobile" genome, is less numerous and thus less conducive to statistical significance of the results, we find that it supports most of the tendencies found for the functionally "stabilized" part of the genome. Our analyses support that there is no fundamental difference between the strictly stable and partly mobile parts of the genome. Rather the same functional categories influence the fractionation fate of genes. Nevertheless, we also demonstrate that the homology sets we study are enriched for some categories and depleted for others, in comparison with random samples of genes from the entire extant genomes. This likely has more to do with lineage-specific expansion of gene families, rather than fractionation dynamics.
Perhaps the most important improvement to be envisaged for our method would be to ensure that the homology sets contain only paralogs created by the initial WGD event. Tandem duplicates or other duplicates produced more recently than the WGD may constitute a major proportion of all the duplicates in a genome. Although SynMap normally excludes tandem duplicates, some of them undoubtedly remain and may be reintroduced when we combine pairwise genome results to form homology sets. Indeed, for future work, it would be worthwhile to contrast the paralog loss behaviour from WGD fractionation with that from other sources of duplicates. This approach was pioneered in  for the Arabidopsis genome, with results very similar to our approach, suggesting our comparative approach might bolster these findings.
Our use of general GO terms could conceivably be improved by using a more focused gene ontology database such as Plant Slim developed by the Arabidopsis Information Resource. Indeed, preliminary tests show that our results on high-level terms could be sharpened using this resource, but unfortunately there is relatively little annotation at present using lower level terms, so this avenue is limited for the moment.
It may advance the understanding of functional associations of fractionation to compare the correlations among functional categories enriched for fractionation-prone genes in contrast to the correlations among these categories for fractionation-resistant genes. It is not a methodological flaw that many categories are correlated among themselves -- indeed it opens up the opportunity to compare changes in these correlations for genes with differing fractionation scores.
Research supported in part by grants from the Natural Sciences and Engineering Research Council of Canada (NSERC). DS holds the Canada Research Chair in Mathematical Genomics.
Publication of this article was supported by the Canada Research Chair in Mathematical Genomics.
This article has been published as part of BMC Bioinformatics Volume 14 Supplement 15, 2013: Proceedings from the Eleventh Annual Research in Computational Molecular Biology (RECOMB) Satellite Workshop on Comparative Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/14/S15.
- Comai L: Genetic and epigenetic interactions in allopolyploid plants. Plant Molecular Biology. 2000, 43: 387-399.View ArticlePubMedGoogle Scholar
- Crow KD, Wagner GP: What is the role of genome duplication in the evolution of complexity and diversity?. Molecular Biology and Evolution. 2006, 23: 887-892. 10.1093/molbev/msj083.View ArticlePubMedGoogle Scholar
- Birchler JA, Veitia RA: The Gene Balance Hypothesis: implications for gene regulation, quantitative traits and evolution. New Phytologist. 2010, 186: 54-62. 10.1111/j.1469-8137.2009.03087.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Zheng C, Chen E, Albert VA, Lyons E, Sankoff D: Ancient eudicot hexaploidy meets ancestral eurosid gene order. BMC Genomics. 2013,Google Scholar
- Jung S, Cestaro A, Troggio M, Main D, Zheng P, Cho I, Folta KM, Sosinski B, A A, Celton JM, Arús P, Shulaev V, Verde I, Morgante M, Rokhsar DS, Velasco R, Sargent DJ: Whole genome comparisons of Fragaria, Prunus and Malus reveal different modes of evolution between Rosaceous subfamilies. BMC Genomics. 2012, 13: 129-10.1186/1471-2164-13-129.PubMed CentralView ArticlePubMedGoogle Scholar
- Argout X, Salse J, Aury JM, Guiltinan MJ, Droc G, Gouzy J, Allegre M, Chaparro C, Legavre T, Maximova SN, Abrouk M, Murat F, Fouet O, Poulain J, Ruiz M, Roguet Y, Rodier-Goud M, Barbosa-Neto JF, Sabot F, Kudrna D, Ammiraju JS, Schuster SC, Carlson JE, Sallet E, Schiex T, Dievart A, Kramer M, Gelley L, Shi Z, Bérard A, Viot C, Boccara M, Risterucci A, Guignon V, Sabau X, Axtell MJ, Ma Z, Zhang Y, Brown S, Bourge M, Golser W, Song X, Clement D, Rivallan R, Tahi M, Akaza JM, Pitollat B, Gramacho K, D'Hont A, Brunel D, Infante D, Kebe I, Costet P, Wing R, McCombie WR, Guiderdoni E, Quétier F, Panaud O, Wincker P, Bocs S, Lanaud C: The genome of Theobroma cacao. Nature Genetics. 2011, 43: 101-108. 10.1038/ng.736.View ArticlePubMedGoogle Scholar
- Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E, Jublot D, Poulain J, Bruyère C, Billault A, Segurens B, Gouyvenoux M, Ugarte E, Cattonaro F, Anthouard V, Vico V, Del Fabbro C, Alaux M, Di Gaspero G, Dumas V, Felice N, Paillard S, Juman I, Moroldo M, Scalabrin S, Canaguier A, Le Clainche I, Malacrida G, Durand E, Pesole G, Laucou V, Chatelet P, Merdinoglu D, Delledonne M, Pezzotti M, Lecharny A, Scarpelli C, Artiguenave F, Pè ME, Valle G, Morgante M, Caboche M, Adam-Blondon AF, Weissenbach J, Quétier F, Wincker P, French-Italian Public Consortium for Grapevine Genome Characterization: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449: 463-467. 10.1038/nature06148.View ArticlePubMedGoogle Scholar
- Freeling M: Bias in plant gene content following different sorts of duplication: tandem, whole-genome, segmental, or by transposition. Annual Review of Plant Biology. 2009, 60: 433-453. 10.1146/annurev.arplant.043008.092122.View ArticlePubMedGoogle Scholar
- Lou P, Wu J, Cheng F, Cressman LG, Wang X, McClung CR: Preferential retention of circadian clock genes during diploidization following whole genome triplication in Brassica rapa. The Plant Cell Online. 2012, 24: 2415-2426. 10.1105/tpc.112.099499.View ArticleGoogle Scholar
- Wu X, Qi X: Genes encoding hub and bottleneck enzymes of the Arabidopsis metabolic network preferentially retain homeologs through whole genome duplication. BMC Evolutionary Biology. 2010, 10: 145-10.1186/1471-2148-10-145.PubMed CentralView ArticlePubMedGoogle Scholar
- Duarte J, Wall PK, Edger P, Landherr LL, Ma H, Pires JC, Leebens-Mack J, dePamphilis CW: Identification of shared single copy nuclear genes in Arabidopsis, Populus, Vitis and Oryza and their phylogenetic utility across various taxonomic levels. BMC Evolutionary Biology. 2010, 10: 61-10.1186/1471-2148-10-61.PubMed CentralView ArticlePubMedGoogle Scholar
- Tomato Genome Consortium: The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012, 485: 635-641. 10.1038/nature11119.View ArticleGoogle Scholar
- Ibarra-Laclette E, Lyons E, Hernández-Guzmán G, Pérez-Torres CA, Carretero-Paulet L, Chang TH, Lan T, Welch AJ, Juárez MJ, Simpson J, Fernández-Cortés A, Arteaga-Vázquez M, Góngora-Castillo E, Acevedo- Hernández G, Schuster SC, Himmelbauer H, Minoche AE, Xu S, Lynch M, Oropeza-Aburto A, Cervantes-Pérez SA, de Jesuś Ortega-Estrada M, Cervantes-Luevano JI, Michael TP, Mockler T, Bryant D, Herrera-Estrella A, Albert VA, Herrera-Estrella L: Architecture and evolution of a minute plant genome. Nature. 2013, 498: 94-98. 10.1038/nature12132.View ArticlePubMedGoogle Scholar
- US Department of Energy JGI: Mimulus version 1. 2010, [http://www.phytozome.net/mimulus]Google Scholar
- Kawahara Y, de la Bastide M, Hamilton J, Kanamori H, McCombie WR, Ouyang S, Schwartz D, Tanaka T, Wu J, Zhou S, Childs K, Davidson R, Lin H, Quesada-Ocampo L, Vaillancourt B, Sakai H, Lee SS, Kim J, Numa H, Itoh T, Buell CR, Matsumoto T: Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice. 2013, 6: 4-10.1186/1939-8433-6-4.View ArticlePubMedGoogle Scholar
- Bennetzen JL, Schmutz J, Wang H, Percifield R, Hawkins J, Pontaroli AC, Estep M, Feng L, Vaughn JN, Grimwood J, Jenkins J, Barry K, Lindquist E, Hellsten U, Deshpande S, Wang X, Wu X, Mitros T, Triplett J, Yang X, Ye CY, Mauro-Herrera M, Wang L, Li P, Sharma M, Sharma R, Ronald PC, Panaud O, Kellogg EA, Brutnell TP, Doust AN, Tuskan GA, Rokhsar D, Devos KM: Reference genome sequence of the model plant Setaria. Nature Biotechnology. 2012, 30: 555-561. 10.1038/nbt.2196.View ArticlePubMedGoogle Scholar
- Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A, Schmutz J, Spannagl M, Tang H, Wang X, Wicker T, Bharti AK, Chapman J, Feltus FA, Gowik U, Lyons E, Maher CA, Martis M, Narechania A, Otillar RP, Penning BW, Salamov AA, Wang Y, Zhang L, Carpita NC, Freeling M, Gingle AR, Hash CT, Keller B, Klein P, Kresovich S, McCann MC, Ming R, Peterson DG, ur Rahman M, Ware D, Westhoff P, Mayer KF, Messing J, Rokhsar DS: The Sorghum bicolor genome and the diversification of grasses. Nature. 2009, 457: 551-556. 10.1038/nature07723.View ArticlePubMedGoogle Scholar
- The International Brachypodium Initiative: Genome sequencing and analysis of the model grass Brachypodium distachyon. Nature. 2010, 463: 763-768. 10.1038/nature08747.View ArticleGoogle Scholar
- Abrouk M, Murat F, Pont C, Messing J, Jackson S, Faraut T, Tannier E, Plomion C, Cooke R, Feuillet C, Salse J: Palaeogenomics of plants: synteny-based modelling of extinct ancestors. Trends in Plant Science. 2010, 15: 479-487. 10.1016/j.tplants.2010.06.001.View ArticlePubMedGoogle Scholar
- Lyons E, Freeling M: How to usefully compare homologous plant genes and chromosomes as DNA sequences. The Plant Journal. 2008, 53: 661-673. 10.1111/j.1365-313X.2007.03326.x.View ArticlePubMedGoogle Scholar
- Lyons E, Pedersen B, Kane J, Alam M, Ming R, Tang H, Wang X, Bowers J, Paterson A, Lisch D, Freeling M: Finding and comparing syntenic regions among Arabidopsis and the outgroups papaya, poplar and grape: CoGe with rosids. Plant Physiology. 2008, 148: 1772-1781. 10.1104/pp.108.124867.PubMed CentralView ArticlePubMedGoogle Scholar
- Schnable JC, Freeling M, Lyons E: Genome-wide analysis of syntenic gene deletion in the grasses. Genome Biology and Evolution. 2012, 4: 265-277. 10.1093/gbe/evs009.PubMed CentralView ArticlePubMedGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Sweetlove LJ, Fernie AR: The Spatial Organization of Metabolism Within the Plant Cell. Annual Review of Plant Biology. 2013, 64: 723-746. 10.1146/annurev-arplant-050312-120233.View ArticlePubMedGoogle Scholar
- Lynch M: Evolutionary diversification of the multimeric states of proteins. Proceedings of the National Academy of Sciences USA. 2013, 110: E2821-E2828. 10.1073/pnas.1310980110.View ArticleGoogle Scholar
- Nussinov R, Ma B, Tsai CJ: A broad view of scaffolding suggests that scaffolding proteins can actively control regulation and signaling of multienzyme complexes through allostery. Biochim Biophys Acta. 2013, 1834: 820-829. 10.1016/j.bbapap.2012.12.014.View ArticlePubMedGoogle Scholar
- Schnable JC, Wang X, Pires JC, Freeling M: Escape from preferential retention following repeated whole genome duplications in plants. Frontiers in Plant Science. 2012, 3: 94-PubMed CentralView ArticlePubMedGoogle Scholar
- Cannon S, Mitra A, Baumgarten A, Young N, May G: The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biology. 2004, 4: 10-10.1186/1471-2229-4-10.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.