- Open Access
The Aquilegia genome reveals a hybrid origin of core eudicots
Genome Biology volume 20, Article number: 256 (2019)
Whole-genome duplications (WGDs) have dominated the evolutionary history of plants. One consequence of WGD is a dramatic restructuring of the genome as it undergoes diploidization, a process under which deletions and rearrangements of various sizes scramble the genetic material, leading to a repacking of the genome and eventual return to diploidy. Here, we investigate the history of WGD in the columbine genus Aquilegia, a basal eudicot, and use it to illuminate the origins of the core eudicots.
Within-genome synteny confirms that columbines are ancient tetraploids, and comparison with the grape genome reveals that this tetraploidy appears to be shared with the core eudicots. Thus, the ancient gamma hexaploidy found in all core eudicots must have involved a two-step process: first, tetraploidy in the ancestry of all eudicots, then hexaploidy in the ancestry of core eudicots. Furthermore, the precise pattern of synteny sharing suggests that the latter involved allopolyploidization and that core eudicots thus have a hybrid origin.
Novel analyses of synteny sharing together with the well-preserved structure of the columbine genome reveal that the gamma hexaploidy at the root of core eudicots is likely a result of hybridization between a tetraploid and a diploid species.
Whole-genome duplication (WGD) is common in the evolutionary history of plants (reviewed in [1, 2]). All flowering plants are descended from a polyploid ancestor, which in turn shows evidence of an even older WGD shared by all seed plants . These repeated cycles of polyploidy dramatically restructure plant genomes. Presumably driven by the “diploidization” process, whereby genomes are returned to an effectively diploid state, chromosomes are scrambled via fusions and fissions, lose both repetitive and genic sequences, or are lost entirely [4,5,6,7,8,9,10,11]. Intriguingly, gene loss after WGD is non-random: there is a bias not only against the retention of certain genes [12, 13] but also against the retention of one of the WGD-derived paralog chromosomes [6, 9, 14,15,16].
We investigated the history of WGDs in the columbine genus Aquilegia for two reasons. The first is related to its phylogenetic position: columbines have been referred to as “basal” eudicots because they appear to be an outgroup to most other “core” eudicot lineages [17, 18]. This matters because our understanding of eudicot karyotype evolution is limited to the heavily sampled core eudicots. Using the recently published Aquilegia coerulea genome , we were able to address key questions about the history of polyploidization in all eudicots. Second, we traced the origins of the columbine chromosomes with a particular focus on the strange chromosome 4, which differs from the rest of the genome in many ways. In particular, it harbors more genetic polymorphism and transposable elements, has lower gene density and reduced gene expression, and appears to migrate more, including between species. It also carries the rDNA clusters, and there is reason to believe that knowing the history of the chromosome could help explain its aberrant behavior .
2.1.Within-genome synteny confirms columbine paleotetraploidy
Ancient WGDs have commonly been inferred from the distribution of divergences between gene duplicates. The simultaneous generation of gene duplicates via WGD is expected to produce a peak in the age distribution relative to the background age distribution of single-gene duplicates [20,21,22]. Such a spike of ancient gene birth was the first evidence of paleotetraploidy in columbines  and was later supported by gene count-based modeling .
Given an assembled genome, a more direct method to infer ancient polyploidy is to look for regions with conserved gene order [25, 26]. Such conservation (a.k.a. synteny) decreases over time due to gene loss and rearrangements, but will still be detectable unless the extent of change is too extensive. We detected a total of 121 synteny block pairs harboring at least five paralogous gene pairs within the columbine genome. The distribution of these blocks across the seven columbine chromosomes indicates pairwise synteny between large genomic regions (Fig. 1). This 1:1 relationship suggests a single round of WGD in columbines and is further supported by similar levels of divergence between synteny pairs (Additional file 1: Figure S1 and S2).
2.2.Columbines share ancient tetraploidy with core eudicots
All sequenced core eudicots appear to share a triplicate genome structure due to paleohexaploidy postdating the separation of monocots and eudicots ([9, 28,29,30,31,32] and Supplementary Note 5 in ). The tetraploidy in columbines, a basal eudicot, might be independent of this ancient “gamma” hexaploidy (Scenarios 1 and 2 in Fig. 2) or might be a remnant of a WGD at the base of all eudicots, which formed the first step of the gamma hexaploidy in core eudicots (Scenario 3 in Fig. 2).
We used the grape (Vitis vinifera) genome as a representative of the core eudicots to distinguish between the three scenarios in Fig. 2. Grape has experienced a relatively small number of chromosomal rearrangements post-gamma and thus strongly resembles the ancestral pre-hexaploid genome . Given the ploidy level of columbine under each scenario, we can predict the ratio of haploid chromosome sets in grape to that in a columbine. If tetraploidy in columbines is lineage-specific and superimposed on the gamma hexaploidy (Scenario 1), we would expect to find a 3:6 ratio of grape and columbine synteny blocks. Instead, we observe a 3:2 relationship (Fig. 3 and Additional file 1: Figure S3) as expected under Scenarios 2 or 3. A similar 3:2 pattern is found in comparisons between grape and sacred lotus . This strongly suggests that basal eudicots do not share the triplicate genome structure of core eudicots, ruling out Scenario 1.
Distinguishing between the two remaining scenarios is more difficult. We began by comparing the divergence at synonymous sites (Ks) between columbine paralogs, grape paralogs, and columbine-grape homologs. In agreement with the analysis of Jiao et al. , the Ks distribution for grape paralogs appears to two major peaks, as expected under the two-step model for gamma hexaploidy (Fig. 4). However, columbine paralogs and columbine-grape homologs each show a single peak of divergence—and the peaks overlap each other and the “older” divergence peak of grape paralogs. This suggests that columbine tetraploidy is derived from the tetraploidy that eventually led to gamma hexaploidy in core eudicots (Scenario 3).
Next we considered the gene order similarity between pairs of columbine and grape chromosomes. If columbine and grape have descended from a common tetraploid ancestor (Scenario 3), they should have inherited all diploidization-driven changes to gene order between the paralogous chromosomes of the tetraploid ancestor. As a result, we expect to see alternative paralogous gene orders to be uniquely shared between two pairs of columbine and grape chromosomes (Fig. 5)— whereas if tetraploidization occurred twice (Scenario 2), no such pattern is expected.
We used two different approaches to detect this pattern. First, we clustered homologous segments based on gene order similarity (“Materials and methods”). The pairwise comparisons show that each member of columbine paralogs matches a different grape chromosome (Additional file 1: Figure S6–S8). Reshuffling genes on grape chromosomes further indicates that this pattern of clustering is highly unlikely to be produced by chance (p = 0–0.05). Second, we attempted to corroborate the clustering based on gene order similarity by clustering homologous regions based on similarity in protein sequence (“Materials and methods”). Because of the deep history of shared tetraploidy, only a small fraction of all the informative gene trees (0.016–0.044) show the “expected” pairings (Additional file 1: Figure S6–S8), and it is thus not possible to infer history from individual trees. Although it is possible that more sophisticated tree-building methods would perform better, the order of the homologous genes that do show the expected pairwise clustering (based on sequence divergence) again recaptures the clustering pattern inferred from synteny alone (compare Additional file 1: Figure S7 and S9). Thus, the clustering pattern inferred from synteny is mirrored in clustering based on sequence divergence.
An eudicot-wide WGD is further supported by the observation that a chromosomal fusion, presumably experienced by the common tetraploid ancestor, is still detectable in the genomes of columbine and grape despite their separation of around 125 million years . The first hint comes from the composition of the chromosomes: columbine chromosome 5 and grape chromosome 7 both appear to be fusions of the same ancestral chromosomes (Fig. 6). If these chromosomes were created by a single fusion event in the common tetraploid ancestor of eudicots, they should match each other with respect to gene order on each of the two homologous portions (“orange” and “green” portions in Fig. 7). This is what we see: columbine chromosome 5 and grape chromosome 7 cluster together with respect to their gene order on the “orange” portion (Additional file 1: Figure S7). For the “green” portion, columbine chromosome 5 matches grape chromosome 4 (Additional file 1: Figure S8), which used to be fused to grape chromosome 7 . The shared ancestral fusion can also be seen using the cacao (Theobroma cacao) genome . The first chromosome of cacao does not only show a similar pattern of chromosomal ancestry [38, 39] but also shares the gene order exclusively with the grape chromosomes 4 and 7 on the corresponding homologous portions (Additional file 1: Figure S10). In summary, the columbine fusion clusters with that of grape, which, in turn, clusters with that of cacao, strongly favoring a common origin of the fusion between “orange” and “green” ancestral chromosomes (Fig. 6).
2.3.The core eudicots have a hybrid origin
Our inference of shared tetraploidy between basal and core eudicots makes use of the signals presumably generated by diploidization (Figs. 5 and 7). However, hybridization of unreduced gametes from two divergent diploid genomes, “allotetraploidy,” would also lead to gene order-based clustering between two different pairs of grape and columbine chromosomes (Additional file 1: Figure S11–S12). In this case, the alternative paralogous gene orders in the tetraploid ancestor reflect the gene orders on the progenitor chromosomes. Thus, the clustering pattern does not depend on whether the eudicot tetraploid genome evolved via “auto-” or “allopolyploidy.” The same is not true for the second part of the process leading to hexaploidy. In this case, autopolyploidy would lead to the duplication of one of the existing chromosomes, whereas allohexaploidy would lead to one of the three paralogous grape chromosomes being an “outlier” with respect to the two grape-columbine pairing (Additional file 1: Figure S11–S12)— which is what we see in our data (Additional file 1: Figure S6–S8).
If our interpretation is correct and all core eudicots have a hybrid origin, the pattern of gene order-based clustering should be conserved. That is, we should be able to identify the hexaploidy-derived “outlier” chromosomes in other core eudicot genomes as well. To check this expectation, we again used the cacao genome, one of the most conserved genomes after grape [9, 39]. Pairwise alignment between the homologous regions of columbine and cacao confirms our expectation: each member of columbine paralogs pairs up with a single cacao chromosome, leaving one of the cacao paralogs as an outlier (Additional file 1: Figure S13–S14). Furthermore, as shown in Fig. 8 (see also Additional file 1: Figure S10), the cacao regions putatively derived from tetraploidy and hexaploidy, respectively, show a very clear one-to-one match to those in grape (detected in the grape-columbine comparison). As expected, the putatively orthologous pairs of cacao and grape regions show similar levels of synteny conservation with their paralogous counterparts, with the “outlier” regions being the most divergent . Thus, although not an independent source of evidence, the cacao genome supports a hybrid origin and highlights the key role of the columbine genome in unraveling the history of the eudicot genome.
2.4.Current columbine chromosomes have mostly been generated via fusions
It is widely accepted that genome shuffling post-WGD has shaped the present-day karyotypes of all plant genomes . Nevertheless, the extent of genome shuffling as a part of the “re-diploidization” process seems to vary widely: only 3 chromosomal rearrangements post-gamma are enough to explain the current structure of the grape genome (Fig. 6) while almost 150 chromosomal rearrangements were necessary for the sunflower genome to reach its current karyotype after several rounds of WGD . To check where columbine falls in this spectrum, we identified chromosomal rearrangements likely to have happened after the tetraploidy shared by all eudicots: if the pre-WGD ancestral eudicot karyotype had a haploid number of 7 chromosomes , only seven columbine-specific fusions and a single fission are enough to explain the reduction in columbine chromosome number from n = 13 to n = 7 after the ancestral fusion event (Fig. 6). These rearrangements involve all the chromosomes in columbine apart from chromosomes 4 and 6, the former of which paradoxically shows the greatest erosion of synteny with grape chromosomes (Fig. 6 and Additional file 1: Figure S3). Given all the evidence suggesting a “decaying” nature of columbine chromosome 4 , we repeated the analysis of grape-columbine synteny detection with relaxed parameter settings. We did this by decreasing the minimum number of aligned gene pairs within a block (from 5 to 3) and increasing the maximum genic distance between matches (from 20 to 30). This allowed us to extend the synteny blocks towards more proximal regions (Additional file 1: Figure S15). Further zooming into the synteny relationship between grape chromosomes that are homologous to columbine chromosome 4 confirmed that there is no evidence of a fusion event (Fig. 9).
The lack of a fusion event on columbine chromosome 6 might explain the fact that it is the smallest chromosome of columbine (Fig. 6). However, chromosome 4 is comparable in size to the remaining chromosomes, all of which are products of ancient fusion events. The observations that chromosome 4 has a higher proportion of genes in tandem duplicates (0.37 versus genome-wide mean of 0.22) and a greater extent of intra-chromosomal synteny (indicative of segmental duplications) (Additional file 1: Figure S17) suggest that chromosome 4 has reached a comparable size partly due to numerous tandem and segmental duplications and partly due to an expansion of repetitive DNA . These results reinforce the idea that chromosome 4 has followed a distinct evolutionary path from the rest of the genome .
Fusion-dominated genome shuffling  is not the only facet of diploidization . Following WGDs, gene duplicates get lost and this happens in a non-random manner. Genes involved in connected molecular functions like kinases, transcription factors, and ribosomal proteins are retained in pairs [41,42,43,44,45] potentially due to dosage-related constraints : losing or duplicating some, but not all of these dosage-sensitive genes, might upset the stoichiometric relationship between their protein products [47,48,49]. Consistent with this dosage balance hypothesis, columbine genes potentially retained post WGD (1302 genes across 76 syntenic regions; Additional file 2: Table S1) are enriched for the GO categories “structural constituent of ribosome,” “transcription factor activity,” “translation” (p < 0.001), and “protein tyrosine kinase activity” (p < 0.01). Tandemly duplicated genes (n = 6972), on the other hand, are depleted for the GO categories “structural constituent of ribosome” and “translation” (p = 10−17), reflecting the role of dosage-related purifying selection.
All flowering plants are descended from a polyploid ancestor, and with only a few exceptions (e.g., Amborella from basal angiosperms ), all of them experienced at least one further round of WGD. Within-genome synteny (Fig. 1) shows that the columbine is an example of the latter, confirming the conclusions from other studies [23, 24, 51]. Here we show that this columbine tetraploidy is a remnant of a WGD at the base of all eudicots and is thus far more ancient than previously thought [23, 51]. Furthermore, we use this observation to argue that the hexaploidy shared by the core eudicots must have involved allopolyploidy, i.e., presumably hybridization between the ancestral tetraploid and a diploid species.
A eudicot-wide WGD has been suggested by several studies [26, 35, 36, 52]. Our synteny-based approach solidifies these findings by demonstrating that the columbine and grape genomes have inherited the genome structure of a common tetraploid ancestor. That we can trace such an ancient tetraploidy is due to two facts. First, the genome structure of columbine is well-preserved and free from recent WGDs. Second, genomes provide much greater information than genes alone. Indeed, a recent study on another basal eudicot , the opium poppy, highlights these two facts. Having experienced a recent WGD (~ 8 million years ago), the genome of the poppy is dominated by syntenic gene pairs of low divergence (Fig. 1C in ), although it also carries highly diverged paralogs whose Ks values nicely overlap with our estimates for columbine and grape, consistent with an eudicot-wide WGD (compare Fig. S13D in  to Fig. 4). However, the strength of the signal from recent polyploidization largely obscures the much weaker signal of ancient polyploidization (interpreted as segmental duplications by Guo et al. ). In fact, although overlooked by the authors, the intergenomic synteny between columbine and poppy provides a clear signature of an eudicot-wide WGD (Fig. S9D in ). Differing from columbines with only one additional genome duplication, the poppy genome aligns to the columbine genome in a 4:2 manner, with 4 paralogous regions of poppy syntenic to 2 paralogous regions of columbine derived from the ancient shared tetraploidy.
Our approach also helps us shed light on the nature of the gamma hexaploidy found in all core eudicots ([9, 28,29,30,31,32], and Supplementary Note 5 in ). WGDs have often been discussed as if they were “events,” ignoring the process by which they originated. We show here that core eudicot hexaploidy is the result of two processes: an ancient tetraploidization shared by all eudicots, followed by allopolyploidization leading to the core eudicots. In other words, all core eudicots have a hybrid origin. An allohexaploid origin has indeed been previously suggested by Murat et al. , who identified the three subgenomes of grape using differential patterns of gene loss on “dominant” versus “sensitive” subgenomes. Their classification assumes that the most recently added set of paralogous chromosomes will be “dominant,” because they have spent a shorter amount of time in the polyploid genome and thus experienced fewer gene losses. Contrary to this, our results suggest that the most recently added grape chromosomes (chromosomes 3, 8, 9, and 14) largely correspond to the “sensitive” grape chromosomes identified by Murat et al. . Instead, we argue that the extensive gene loss in the most recently added subgenome reflects its divergence from the other two subgenomes at the time of hexaploid formation, perhaps similar to the situation in the allotetraploid Arabidopsis suecica, which is a hybrid between the more ancestral-like (n = 8) genome of A. arenosa and the heavily reduced (n = 5) genome of A. thaliana . Another example is hexaploid wheat, which is a hybrid between tetraploid emmer wheat and wild diploid grass, Aegilops tauschii ( and references therein]).
Our findings reveal the hybrid structure of core eudicot genomes and will hopefully help us understand what hybridization has meant for core eudicots—a group which comprises more than 70% of all living flowering plants . What are the hybridization-coupled changes that have led to the current patterns of gene expression, methylation, or transposable element density/distribution? All these questions call for additional genomes from basal eudicots which—as this study illustrates—have great values as outgroup to the core eudicots. More data will also allow the development of sophisticated analysis methods based on explicit models of the evolution of gene order, which our results suggest is a very powerful source of information about the past.
5.Materials and methods
We performed all genes (CDS)-against-all genes (CDS) BLAST for the latest version of Aquilegia coerulea reference genome (v3.1) using the SynMap tool  in the online CoGe portal . We also looked at the synteny within Vitis vinifera (v12) and between A. coerulea and V. vinifera using default parameter combinations in DAGChainer. We filtered the raw output files for both within-grape and columbine-to-grape synteny. For the former, we only kept the blocks that are syntenic between the polyploidy-derived paralogous chromosomes of grape as identified by Jaillon et al.  (Additional file 4: Table S3). For the latter, we required that a given columbine chromosome is overall syntenic to all the three paralogous chromosomes of grape (Additional file 3: Table S2). So, for a given pair of columbine and grape chromosomes, we only kept the blocks if the columbine chromosome also matches to the other members of paralogous grape chromosomes.
The raw output files can be regenerated at the CoGe portal  using the id numbers provided below for each species (Availability of data and materials) and changing the default parameter combination in DAGChainer (D:A = 20:5) when needed. D and A specify the maximum genic distance between two matches and the minimum number of aligned gene pairs, respectively, to form a collinear syntenic block.
5.2.Estimating the divergence between synteny block pairs
We used Ks (the number of synonymous substitutions per synonymous site) values provided for each homolog gene pair by the CoGe portal . We estimated the median Ks of homologous genes in a synteny block after filtering gene pairs with Ks > 10 due to a saturation effect . Both values are provided in Supplementary Data 1–3 for within columbine, columbine-to-grape, and within-grape synteny, respectively.
5.3.Quantifying gene order similarity
We “reconstructed” a given set of columbine and grape chromosomes at their homologous regions (color-coded in Fig. 6). We seeded this reconstruction by focusing on at least three consecutive genes aligning between a pair of columbine and grape chromosomes (D:A = 0:3). We particularly chose three genes since it is the most stringent value we could use to detect homologous synteny blocks; we detected almost nothing when we required 4 consecutive genes (D:A = 0:4). This stringent criteria aim to minimize the effect of gene movement on the homology between columbine and grape chromosomes. Once we had the list of genes, we then looked for their paralogous counterparts on the remaining columbine and grape chromosomes using intragenomic gene-to-gene blast (D:A = 0:1). Having chromosomes represented by syntenic gene sets and reminiscent of these sets (Additional file 1: Figure S5 and Additional file 5: Table S4), we assigned a unique word to each synteny block and the genes forming the block to be able to use the text alignment provided by the R package align_local . We then quantified the gene (“word”) similarity as such: for an initial N number of words on a columbine chromosome (N = window size), we did a pairwise alignment between these N words and all the words a grape chromosome (match = 4, gap = − 1). We repeated the same analysis with the inverted order of N words and picked the maximum alignment score. We repeated these steps by sliding the window by one word and keeping the N constant to get a distribution of scores as in Additional file 1: Figure S6–S8. We used different N values ranging from 4 to 15. Note that we excluded columbine chromosomes 3, 4, and 7 from this analysis since they all have a complex history of lineage-specific chromosomal reshuffling events (Figs. 1, 6, Additional file 1: Figure S4 and S17).
We applied the same stringent criteria (D:A = 0:3) to detect the homologous regions between grape and cacao (Theobroma cacao, v1). The same criteria led to very few homologous regions between columbine and cacao. So, we relaxed the parameters for the synteny detection between these two genomes (D:A = 0:2) and quantified the gene order similarity with greater window sizes (N = 20, 30, 35, 40, and 50). Note that we focused on the triplicated regions distributed across 3 different cacao chromosomes (Fig. 8, Additional file 1: Figure S13–S14), which are rather unaffected by lineage-specific shuffling .
5.4.Statistical testing of gene order similarity
Given the gene order similarity between the two different pairs of columbine and grape chromosomes harboring homologous regions, we performed permutation tests to estimate the probability of observing such a clustering just by chance. To do so, we first combined all the grape genes and sampled the same number of genes (“words”) as we observe to reconstruct each of the paralogous grape chromosome. We repeated the quantification step as above to get a permuted distribution of alignment score between a pair of columbine and grape chromosomes. We used the Wilcoxon rank sum test (W-statistic) to quantify the shift in the distribution of alignment scores between one of the members of columbine paralogous chromosomes and its best grape hit when combined with the alignment scores between the same columbine chromosome and other grape chromosomes. We repeated the same analysis for the other member of columbine paralogous chromosomes as well. Having these observed W-statistics, we counted the number of cases (out of 100) where the permuted distributions generate W-statistics as high as or higher than the observed ones. We ran permutation tests for the columbine-cacao pairing as well (Additional file 1: Figure S13–S14).
5.5.Building gene trees
We built upgma trees for the homologous genes (Additional file 6: Table S5) distributed across a given set of columbine and grape chromosomes (color-coded in Fig. 6). We first detected homologous genes aligning between a pair of columbine and grape chromosomes (D:A = 0:1). We then searched for their paralogous counterparts using intragenomic blasts (D:A = 0:1). For protein alignment, we required at least five homologous genes, each from a single chromosome in the given set and ran ClustalW2 (v2.1) with the options -TREE -KIMURA -CLUSTERING = UPGMA -OUTPUTTREE = dist . Of all the trees generated by ClustalW2 (informative trees), we only focused on the ones that support the synteny-based pairings (Additional file 1: Figure S6–S8), which are detected by the subtrees function in R package ape [60, 61]. Once we had the sets of homologous genes from this subset of trees, we assigned a unique word to the each set and quantified the gene order similarity between pairs of columbine and grape chromosomes as mentioned above.
For protein sequences, we used the annotations provided by JGI and Ensembl for columbine  and grape , respectively. Note that CoGe  outputs grape genes with the “PAC” tag while they are tagged with “VIT” in the Ensembl database. To match these different ids, we used two intermediary files. The first one is a gff file provided by CoGe (available at https://genomevolution.org/coge/GenomeInfo.pl?gid=19990). The second one is a conversion file provided by the Grape Genome Database  that lists the correspondence between different gene ids (can be downloaded from http://genomes.cribi.unipd.it/DATA/). These two files contain the common tag “GSVIVT” which bridges the “PAC” and “VIT” tags.
5.6.GO enrichment analysis
We used gene annotations provided by JGI  to test the null hypothesis that the property for a gene to be retained post-WGD and to belong to a given GO category are independent. We created a 2 × 2 contingency table (Table 1) and applied Fisher’s exact test for each GO category independently. We repeated the same analysis for tandem gene duplicates as identified by SynMap [29, 56]; this time, testing the null hypotheses that the property for a gene to be tandemly duplicated and to belong to a given GO category are independent. We excluded genes on scaffolds and reported enriched/depleted categories if they remain significant (p < 0.05) after multiple test correction (fdr).
Availability of data and materials
The columbine, grape, and cacao genomes are available at the CoGE portal for the synteny analyses with the id numbers 28620, 19990, and 25287, respectively .
Wendel JF. The wondrous cycles of polyploidy in plants. Am J Bot. 2015;102:1753–6.
Van de Peer Y, Mizrachi E, Marchal K. The evolutionary significance of polyploidy. Nat Rev Genet. 2017;18:411–24.
Jiao Y, Wickett NJ, Ayyampalayam S, Chanderbali AS, Landherr L, Ralph PE, et al. Ancestral polyploidy in seed plants and angiosperms. Nature. 2011;473:97–100.
Leitch IJ, Bennett MD. Genome downsizing in polyploid plants. Biol J Linn Soc Lond. Oxford University Press. 2004;82:651–63.
Comai L. The advantages and disadvantages of being polyploid. Nat Rev Genet. 2005;6:836–46.
Thomas BC, Pedersen B, Freeling M. Following tetraploidy in an Arabidopsis ancestor, genes were removed preferentially from one homeolog leaving clusters enriched in dose-sensitive genes. Genome Res. 2006;16:934–46.
Otto SP. The evolutionary consequences of polyploidy. Cell. 2007;131:452–62.
Renny-Byfield S, Chester M, Kovařík A, Le Comber SC, Grandbastien M-A, Deloger M, et al. Next generation sequencing reveals genome downsizing in allotetraploid Nicotiana tabacum, predominantly through the elimination of paternally derived repetitive DNAs. Mol Biol Evol. 2011;28:2843–54.
Murat F, Zhang R, Guizard S, Gavranović H, Flores R, Steinbach D, et al. Karyotype and gene order evolution from reconstructed extinct ancestors highlight contrasts in genome plasticity of modern rosid crops. Genome Biol Evol. 2015;7:735–49.
Murat F, Armero A, Pont C, Klopp C, Salse J. Reconstructing the genome of the most recent common ancestor of flowering plants. Nat Genet. 2017;49:490–6.
Badouin H, Gouzy J, Grassa CJ, Murat F, Staton SE, Cottret L, et al. The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution. Nature. 2017;546:148–52.
Blanc G, Wolfe KH. Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution. Plant Cell. 2004;16:1679–91.
De Smet R, Adams KL, Vandepoele K, Van Montagu MCE, Maere S, Van de Peer Y. Convergent gene loss following gene and genome duplications creates single-copy families in flowering plants. Proc Natl Acad Sci U S A. 2013;110:2898–903.
Woodhouse MR, Schnable JC, Pedersen BS, Lyons E, Lisch D, Subramaniam S, et al. Following Tetraploidy in Maize, a Short Deletion Mechanism Removed Genes Preferentially from One of the Two Homeologs. PLoS Biol. Public Library of Science; 2010;8:e1000409.
Schnable JC, Springer NM, Freeling M. Differentiation of the maize subgenomes by genome dominance and both ancient and ongoing gene loss. Proc Natl Acad Sci U S A. 2011;108:4069–74.
Tang H, Woodhouse MR, Cheng F, Schnable JC, Pedersen BS, Conant G, et al. Altered patterns of fractionation and exon deletions in Brassica rapa support a two-step model of paleohexaploidy. Genetics. 2012;190:1563–74.
Hodges SA, Kramer EM. Columbines. Curr Biol. 2007;17:R992–4.
Soltis DE, Smith SA, Cellinese N, Wurdack KJ, Tank DC, Brockington SF, et al. Angiosperm phylogeny: 17 genes, 640 taxa. Am J Bot. 2011;98:704–30.
Filiault DL, Ballerini ES, Mandáková T, Aköz G, Derieg NJ, Schmutz J, et al. The Aquilegia genome provides insight into adaptive radiation and reveals an extraordinarily polymorphic chromosome with a unique history. Elife. 2018;7. https://doi.org/10.7554/eLife.36426.
Lynch M, Conery JS. The evolutionary fate and consequences of duplicate genes. Science. 2000;290:1151–5.
Blanc G, Wolfe KH. Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes. Plant Cell. 2004;16:1667–78.
Cui L, Wall PK, Leebens-Mack JH, Lindsay BG, Soltis DE, Doyle JJ, et al. Widespread genome duplications throughout the history of flowering plants. Genome Res. 2006;16:738–49.
Vanneste K, Baele G, Maere S, Van de Peer Y. Analysis of 41 plant genomes supports a wave of successful genome duplications in association with the Cretaceous–Paleogene boundary. Genome Res. 2014;24:1334–47.
Tiley GP, Ané C, Burleigh JG. Evaluating and Characterizing Ancient Whole-Genome Duplications in Plants with Gene Count Data. Genome Biol Evol. 2016;8:1023–37.
Doyle JJ, Egan AN. Dating the origins of polyploidy events. New Phytol. 2010;186:73–85.
Jiao Y, Paterson AH. Polyploidy-associated genome modifications during land plant evolution. Philos Trans R Soc Lond B Biol Sci. 2014;369. https://doi.org/10.1098/rstb.2013.0355.
Melters DP, Bradnam KR, Young HA, Telis N, May MR, Ruby JG, et al. Comparative analysis of tandem repeats from hundreds of species reveals unique insights into centromere evolution. Genome Biol. 2013;14:R10.
Jaillon O, Aury J-M, Noel B, Policriti A, Clepet C, Casagrande A, et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007;449:463–7.
Lyons E, Pedersen B, Kane J, Freeling M. The Value of Nonmodel Genomes and an Example Using SynMap Within CoGe to Dissect the Hexaploidy that Predates the Rosids. Trop Plant Biol. 2008;1:181–90.
Potato Genome Sequencing Consortium, Xu X, Pan S, Cheng S, Zhang B, Mu D, et al. Genome sequence and analysis of the tuber crop potato. Nature. 2011;475:189–95.
Truco MJ, Ashrafi H, Kozik A, van Leeuwen H, Bowers J, Wo SRC, et al. An Ultra-High-Density, Transcript-Based, Genetic Map of Lettuce. G3. 2013;3:617–31.
Denoeud F, Carretero-Paulet L, Dereeper A, Droc G, Guyot R, Pietrella M, et al. The coffee genome provides insight into the convergent evolution of caffeine biosynthesis. Science. 2014;345:1181–4.
Bombarely A, Moser M, Amrad A, Bao M, Bapaume L, Barry CS, et al. Insight into the evolution of the Solanaceae from the parental genomes of Petunia hybrida. Nat Plants. 2016;2:16074.
Salse J. Ancestors of modern plant crops. Curr Opin Plant Biol. 2016;30:134–42.
Ming R, VanBuren R, Liu Y, Yang M, Han Y, Li L-T, et al. Genome of the long-living sacred lotus (Nelumbo nucifera Gaertn.). Genome Biol. 2013;14:R41.
Jiao Y, Leebens-Mack J, Ayyampalayam S, Bowers JE, McKain MR, McNeal J, et al. A genome triplication associated with early diversification of the core eudicots. Genome Biol. 2012;13:R3.
Zeng L, Zhang Q, Sun R, Kong H, Zhang N, Ma H. Resolution of deep angiosperm phylogeny using conserved nuclear genes and estimates of early divergence times. Nat Commun. 2014;5:4956.
Zheng C, Chen E, Albert VA, Lyons E, Sankoff D. Ancient eudicot hexaploidy meets ancestral eurosid gene order. BMC Genomics. 2013;14(Suppl 7):S3.
Argout X, Salse J, Aury J-M, Guiltinan MJ, Droc G, Gouzy J, et al. The genome of Theobroma cacao. Nat Genet. 2011;43:101–8.
Soltis DE, Visger CJ, Marchant DB, Soltis PS. Polyploidy: Pitfalls and paths to a paradigm. Am J Bot. 2016;103:1146–66.
Seoighe C, Wolfe KH. Yeast genome evolution in the post-genome era. Curr Opin Microbiol. 1999;2:548–54.
Tian C-G, Xiong Y-Q, Liu T-Y, Sun S-H, Chen L-B, Chen M-S. Evidence for an ancient whole-genome duplication event in rice and other cereals. Yi Chuan Xue Bao. 2005;32:519–27.
Aury J-M, Jaillon O, Duret L, Noel B, Jubin C, Porcel BM, et al. Global trends of whole-genome duplications revealed by the ciliate Paramecium tetraurelia. Nature. 2006;444:171–8.
Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, et al. The B73 maize genome: complexity, diversity, and dynamics. Science. 2009;326:1112–5.
Rodgers-Melnick E, Mane SP, Dharmawardhana P, Slavov GT, Crasta OR, Strauss SH, et al. Contrasting patterns of evolution following whole genome versus tandem duplication events in Populus. Genome Res. 2012;22:95–105.
Birchler JA, Bhadra U, Bhadra MP, Auger DL. Dosage-dependent gene regulation in multicellular eukaryotes: implications for dosage compensation, aneuploid syndromes, and quantitative traits. Dev Biol. 2001;234:275–88.
Papp B, Pál C, Hurst LD. Dosage sensitivity and the evolution of gene families in yeast. Nature. 2003;424:194–7.
Freeling M, Thomas BC. Gene-balanced duplications, like tetraploidy, provide predictable drive to increase morphological complexity. Genome Res. 2006;16:805–14.
Birchler JA, Veitia RA. Gene balance hypothesis: connecting issues of dosage sensitivity across biological disciplines. Proc Natl Acad Sci U S A. 2012;109:14746–53.
Amborella Genome Project. The Amborella genome and the evolution of flowering plants. Science. 2013;342:1241089.
Guo L, Winzer T, Yang X, Li Y, Ning Z, He Z, et al. The opium poppy genome and morphinan production. Science. 2018. https://doi.org/10.1126/science.aat4096.
Malacarne G, Perazzolli M, Cestaro A, Sterck L, Fontana P, Van de Peer Y, et al. Deconstruction of the (paleo)polyploid grapevine genome based on the analysis of transposition events involving NBS resistance genes. PLoS One. 2012;7:e29762.
O’Kane SL, Schaal BA, Al-Shehbaz IA. The Origins of Arabidopsis suecica (Brassicaceae) as Indicated by Nuclear rDNA Sequences. Syst Bot. 1996;21:559–66.
Matsuoka Y. Evolution of polyploid triticum wheats under cultivation: the role of domestication, natural hybridization and allopolyploid speciation in their diversification. Plant Cell Physiol. 2011;52:750–64.
Friis EM, Pedersen KR, Crane PR. The emergence of core eudicots: new floral evidence from the earliest Late Cretaceous. Proc Biol Sci. 2016;283. https://doi.org/10.1098/rspb.2016.1325.
Lyons E, Freeling M. How to usefully compare homologous plant genes and chromosomes as DNA sequences. Plant J. 2008;53:661–73.
Blanc G, Hokamp K, Wolfe KH. A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome. Genome Res. 2003;13:137–44.
Smith DA, Cordell R, Dillon EM. Infectious texts: Modeling text reuse in nineteenth-century newspapers. 2013 IEEE International Conference on Big Data; 2013. p. 86–94.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013. Available from: https://www.R-project.org/
Paradis E, Claude J, Strimmer K. APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics. 2004;20:289–90.
Kersey PJ, Allen JE, Allot A, Barba M, Boddu S, Bolt BJ, et al. Ensembl Genomes 2018: an integrated omics infrastructure for non-vertebrate species. Nucleic Acids Res. 2018;46:D802–8.
Vitulo N, Forcato C, Carpinelli EC, Telatin A, Campagna D, D’Angelo M, et al. A deep survey of alternative splicing in grape reveals changes in the splicing machinery related to tissue, stress condition and genotype. BMC Plant Biol. 2014;14:99.
We thank Robin Burns and Claus Vogl for their comments on the manuscript; Daniel Gómez Sánchez and Benjamin Jaegle for the fruitful discussions.
Peer review information
Andrew Cosgrove and Barbara Cheifet were the primary editors on this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.
G.A. was supported by the Vienna Graduate School of Population Genetics (Austrian Science Fund, FWF: DK W1225-B20).
Ethics approval and consent to participate
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Aköz, G., Nordborg, M. The Aquilegia genome reveals a hybrid origin of core eudicots. Genome Biol 20, 256 (2019). https://doi.org/10.1186/s13059-019-1888-8