Insights into bear evolution from a Pleistocene polar bear genome.

The polar bear (Ursus maritimus) has become a symbol of the threat to biodiversity from climate change. Understanding polar bear evolutionary history may provide insights into apex carnivore responses and prospects during periods of extreme environmental perturbations. In recent years, genomic studies have examined bear speciation and population history, including evidence for ancient admixture between polar bears and brown bears (Ursus arctos). Here, we extend our earlier studies of a 130,000- to 115,000-y-old polar bear from the Svalbard Archipelago using a 10× coverage genome sequence and 10 new genomes of polar and brown bears from contemporary zones of overlap in northern Alaska. We demonstrate a dramatic decline in effective population size for this ancient polar bear's lineage, followed by a modest increase just before its demise. A slightly higher genetic diversity in the ancient polar bear suggests a severe genetic erosion over a prolonged bottleneck in modern polar bears. Statistical fitting of data to alternative admixture graph scenarios favors at least one ancient introgression event from brown bears into the ancestor of polar bears, possibly dating back over 150,000 y. Gene flow was likely bidirectional, but allelic transfer from brown into polar bear is the strongest detected signal, which contrasts with other published work. These findings may have implications for our understanding of climate change impacts: Polar bears, a specialist Arctic lineage, may not only have undergone severe genetic bottlenecks but also been the recipient of generalist, boreal genetic variants from brown bears during critical phases of Northern Hemisphere glacial oscillations.

The digestive systems of carnivorous plants.

To survive in the nutrient-poor habitats, carnivorous plants capture small organisms comprising complex substances not suitable for immediate reuse. The traps of carnivorous plants, which are analogous to the digestive systems of animals, are equipped with mechanisms for the breakdown and absorption of nutrients. Such capabilities have been acquired convergently over the past tens of millions of years in multiple angiosperm lineages by modifying plant-specific organs including leaves. The epidermis of carnivorous trap leaves bears groups of specialized cells called glands, which acquire substances from their prey via digestion and absorption. The digestive glands of carnivorous plants secrete mucilage, pitcher fluids, acids, and proteins, including digestive enzymes. The same (or morphologically distinct) glands then absorb the released compounds via various membrane transport proteins or endocytosis. Thus, these glands function in a manner similar to animal cells that are physiologically important in the digestive system, such as the parietal cells of the stomach and intestinal epithelial cells. Yet, carnivorous plants are equipped with strategies that deal with or incorporate plant-specific features, such as cell walls, epidermal cuticles, and phytohormones. In this review, we provide a systematic perspective on the digestive and absorptive capacity of convergently evolved carnivorous plants, with an emphasis on the forms and functions of glands.

Chiral Central Charge from a Single Bulk Wave Function.

A (2+1)-dimensional gapped quantum many-body system can have a topologically protected energy current at its edge. The magnitude of this current is determined entirely by the temperature and the chiral central charge, a quantity associated with the effective field theory of the edge. We derive a formula for the chiral central charge that, akin to the topological entanglement entropy, is completely determined by the many-body ground state wave function in the bulk. According to our formula, nonzero chiral central charge gives rise to a topological obstruction that prevents the ground state wave function from being real valued in any local product basis.

The avocado genome informs deep angiosperm phylogeny, highlights introgressive hybridization, and reveals pathogen-influenced gene space adaptation.

The avocado, Persea americana, is a fruit crop of immense importance to Mexican agriculture with an increasing demand worldwide. Avocado lies in the anciently diverged magnoliid clade of angiosperms, which has a controversial phylogenetic position relative to eudicots and monocots. We sequenced the nuclear genomes of the Mexican avocado race, P. americana var. drymifolia, and the most commercially popular hybrid cultivar, Hass, and anchored the latter to chromosomes using a genetic map. Resequencing of Guatemalan and West Indian varieties revealed that ∼39% of the Hass genome represents Guatemalan source regions introgressed into a Mexican race background. Some introgressed blocks are extremely large, consistent with the recent origin of the cultivar. The avocado lineage experienced 2 lineage-specific polyploidy events during its evolutionary history. Although gene-tree/species-tree phylogenomic results are inconclusive, syntenic ortholog distances to other species place avocado as sister to the enormous monocot and eudicot lineages combined. Duplicate genes descending from polyploidy augmented the transcription factor diversity of avocado, while tandem duplicates enhanced the secondary metabolism of the species. Phenylpropanoid biosynthesis, known to be elicited by Colletotrichum (anthracnose) pathogen infection in avocado, is one enriched function among tandems. Furthermore, transcriptome data show that tandem duplicates are significantly up- and down-regulated in response to anthracnose infection, whereas polyploid duplicates are not, supporting the general view that collections of tandem duplicates contribute evolutionarily recent "tuning knobs" in the genome adaptive landscapes of given species.

TCP and MADS-Box Transcription Factor Networks Regulate Heteromorphic Flower Type Identity in Gerbera hybrida.

The large sunflower family, Asteraceae, is characterized by compressed, flower-like inflorescences that may bear phenotypically distinct flower types. The CYCLOIDEA (CYC)/TEOSINTE BRANCHED1-like transcription factors (TFs) belonging to the TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR (TCP) protein family are known to regulate bilateral symmetry in single flowers. In Asteraceae, they function at the inflorescence level, and were recruited to define differential flower type identities. Here, we identified upstream regulators of GhCYC3, a gene that specifies ray flower identity at the flower head margin in the model plant Gerbera hybrida We discovered a previously unidentified expression domain and functional role for the paralogous CINCINNATA-like TCP proteins. They function upstream of GhCYC3 and affect the developmental delay of marginal ray primordia during their early ontogeny. At the level of single flowers, the Asteraceae CYC genes show a unique function in regulating the elongation of showy ventral ligules that play a major role in pollinator attraction. We discovered that during ligule development, the E class MADS-box TF GRCD5 activates GhCYC3 expression. We propose that the C class MADS-box TF GAGA1 contributes to stamen development upstream of GhCYC3 Our data demonstrate how interactions among and between the conserved floral regulators, TCP and MADS-box TFs, contribute to the evolution of the elaborate inflorescence architecture of Asteraceae.

Tandem gene duplications drive divergent evolution of caffeine and crocin biosynthetic pathways in plants.

BACKGROUND: Plants have evolved a panoply of specialized metabolites that increase their environmental fitness. Two examples are caffeine, a purine psychotropic alkaloid, and crocins, a group of glycosylated apocarotenoid pigments. Both classes of compounds are found in a handful of distantly related plant genera (Coffea, Camellia, Paullinia, and Ilex for caffeine; Crocus, Buddleja, and Gardenia for crocins) wherein they presumably evolved through convergent evolution. The closely related Coffea and Gardenia genera belong to the Rubiaceae family and synthesize, respectively, caffeine and crocins in their fruits. RESULTS: Here, we report a chromosomal-level genome assembly of Gardenia jasminoides, a crocin-producing species, obtained using Oxford Nanopore sequencing and Hi-C technology. Through genomic and functional assays, we completely deciphered for the first time in any plant the dedicated pathway of crocin biosynthesis. Through comparative analyses with Coffea canephora and other eudicot genomes, we show that Coffea caffeine synthases and the first dedicated gene in the Gardenia crocin pathway, GjCCD4a, evolved through recent tandem gene duplications in the two different genera, respectively. In contrast, genes encoding later steps of the Gardenia crocin pathway, ALDH and UGT, evolved through more ancient gene duplications and were presumably recruited into the crocin biosynthetic pathway only after the evolution of the GjCCD4a gene. CONCLUSIONS: This study shows duplication-based divergent evolution within the coffee family (Rubiaceae) of two characteristic secondary metabolic pathways, caffeine and crocin biosynthesis, from a common ancestor that possessed neither complete pathway. These findings provide significant insights on the role of tandem duplications in the evolution of plant specialized metabolism.

Symmetry Breaking and Error Correction in Open Quantum Systems.

Symmetry-breaking transitions are a well-understood phenomenon of closed quantum systems in quantum optics, condensed matter, and high energy physics. However, symmetry breaking in open systems is less thoroughly understood, in part due to the richer steady-state and symmetry structure that such systems possess. For the prototypical open system-a Lindbladian-a unitary symmetry can be imposed in a "weak" or a "strong" way. We characterize the possible Z_{n} symmetry-breaking transitions for both cases. In the case of Z_{2}, a weak-symmetry-broken phase guarantees at most a classical bit steady-state structure, while a strong-symmetry-broken phase admits a partially protected steady-state qubit. Viewing photonic cat qubits through the lens of strong-symmetry breaking, we show how to dynamically recover the logical information after any gap-preserving strong-symmetric error; such recovery becomes perfect exponentially quickly in the number of photons. Our study forges a connection between driven-dissipative phase transitions and error correction.

Long-read sequencing uncovers the adaptive topography of a carnivorous plant genome.

Utricularia gibba, the humped bladderwort, is a carnivorous plant that retains a tiny nuclear genome despite at least two rounds of whole genome duplication (WGD) since common ancestry with grapevine and other species. We used a third-generation genome assembly with several complete chromosomes to reconstruct the two most recent lineage-specific ancestral genomes that led to the modern U. gibba genome structure. Patterns of subgenome dominance in the most recent WGD, both architectural and transcriptional, are suggestive of allopolyploidization, which may have generated genomic novelty and led to instantaneous speciation. Syntenic duplicates retained in polyploid blocks are enriched for transcription factor functions, whereas gene copies derived from ongoing tandem duplication events are enriched in metabolic functions potentially important for a carnivorous plant. Among these are tandem arrays of cysteine protease genes with trap-specific expression that evolved within a protein family known to be useful in the digestion of animal prey. Further enriched functions among tandem duplicates (also with trap-enhanced expression) include peptide transport (intercellular movement of broken-down prey proteins), ATPase activities (bladder-trap acidification and transmembrane nutrient transport), hydrolase and chitinase activities (breakdown of prey polysaccharides), and cell-wall dynamic components possibly associated with active bladder movements. Whereas independently polyploid Arabidopsis syntenic gene duplicates are similarly enriched for transcriptional regulatory activities, Arabidopsis tandems are distinct from those of U. gibba, while still metabolic and likely reflecting unique adaptations of that species. Taken together, these findings highlight the special importance of tandem duplications in the adaptive landscapes of a carnivorous plant genome.

Architecture and evolution of a minute plant genome.

It has been argued that the evolution of plant genome size is principally unidirectional and increasing owing to the varied action of whole-genome duplications (WGDs) and mobile element proliferation. However, extreme genome size reductions have been reported in the angiosperm family tree. Here we report the sequence of the 82-megabase genome of the carnivorous bladderwort plant Utricularia gibba. Despite its tiny size, the U. gibba genome accommodates a typical number of genes for a plant, with the main difference from other plant genomes arising from a drastic reduction in non-genic DNA. Unexpectedly, we identified at least three rounds of WGD in U. gibba since common ancestry with tomato (Solanum) and grape (Vitis). The compressed architecture of the U. gibba genome indicates that a small fraction of intergenic DNA, with few or no active retrotransposons, is sufficient to regulate and integrate all the processes required for the development and reproduction of a complex organism.

Dissecting functions of SEPALLATA-like MADS box genes in patterning of the pseudanthial inflorescence of Gerbera hybrida.

The pseudanthial inflorescences of the sunflower family, Asteraceae, mimic a solitary flower but are composed of multiple flowers. Our studies in Gerbera hybrida indicate functional diversification for SEPALLATA (SEP)-like MADS box genes that often function redundantly in other core eudicots. We conducted phylogenetic and expression analysis for eight SEP-like GERBERA REGULATOR OF CAPITULUM DEVELOPMENT (GRCD) genes, including previously unstudied gene family members. Transgenic gerbera plants were used to infer gene functions. Adding to the previously identified stamen and carpel functions for GRCD1 and GRCD2, two partially redundant genes, GRCD4 and GRCD5, were found to be indispensable for petal development. Stepwise conversion of floral organs into leaves in the most severe RNA interference lines suggest redundant and additive GRCD activities in organ identity regulation. We show conserved and redundant functions for several GRCD genes in regulation of flower meristem maintenance, while functional diversification for three SEP1/2/4 clade genes in regulation of inflorescence meristem patterning was observed. GRCD genes show both specialized and pleiotropic functions contributing to organ differentiation and flower meristem fate, and uniquely, to patterning of the inflorescence meristem. Altogether, we provide an example of how plant reproductive evolution has used conserved genetic modules for regulating the elaborate inflorescence architecture in Asteraceae.

Evolutionary Co-Option of Floral Meristem Identity Genes for Patterning of the Flower-Like Asteraceae Inflorescence.

The evolutionary success of Asteraceae, the largest family of flowering plants, has been attributed to the unique inflorescence architecture of the family, which superficially resembles an individual flower. Here, we show that Asteraceae inflorescences (flower heads, or capitula) resemble solitary flowers not only morphologically but also at the molecular level. By conducting functional analyses for orthologs of the flower meristem identity genes LEAFY (LFY) and UNUSUAL FLORAL ORGANS (UFO) in Gerbera hybrida, we show that GhUFO is the master regulator of flower meristem identity, while GhLFY has evolved a novel, homeotic function during the evolution of head-like inflorescences. Resembling LFY expression in a single flower meristem, uniform expression of GhLFY in the inflorescence meristem defines the capitulum as a determinate structure that can assume floral fate upon ectopic GhUFO expression. We also show that GhLFY uniquely regulates the ontogeny of outer, expanded ray flowers but not inner, compact disc flowers, indicating that the distinction of different flower types in Asteraceae is connected with their independent evolutionary origins from separate branching systems.

Genotyping-by-sequencing provides the first well-resolved phylogeny for coffee (Coffea) and insights into the evolution of caffeine content in its species: GBS coffee phylogeny and the evolution of caffeine content.

A comprehensive and meaningful phylogenetic hypothesis for the commercially important coffee genus (Coffea) has long been a key objective for coffee researchers. For molecular studies, progress has been limited by low levels of sequence divergence, leading to insufficient topological resolution and statistical support in phylogenetic trees, particularly for the major lineages and for the numerous species occurring in Madagascar. We report here the first almost fully resolved, broadly sampled phylogenetic hypothesis for coffee, the result of combining genotyping-by-sequencing (GBS) technology with a newly developed, lab-based workflow to integrate short read next-generation sequencing for low numbers of additional samples. Biogeographic patterns indicate either Africa or Asia (or possibly the Arabian Peninsula) as the most likely ancestral locality for the origin of the coffee genus, with independent radiations across Africa, Asia, and the Western Indian Ocean Islands (including Madagascar and Mauritius). The evolution of caffeine, an important trait for commerce and society, was evaluated in light of our phylogeny. High and consistent caffeine content is found only in species from the equatorial, fully humid environments of West and Central Africa, possibly as an adaptive response to increased levels of pest predation. Moderate caffeine production, however, evolved at least one additional time recently (between 2 and 4Mya) in a Madagascan lineage, which suggests that either the biosynthetic pathway was already in place during the early evolutionary history of coffee, or that caffeine synthesis within the genus is subject to convergent evolution, as is also the case for caffeine synthesis in coffee versus tea and chocolate.

Cat Codes with Optimal Decoherence Suppression for a Lossy Bosonic Channel.

We investigate cat codes that can correct multiple excitation losses and identify two types of logical errors: bit-flip errors due to excessive excitation loss and dephasing errors due to quantum backaction from the environment. We show that selected choices of logical subspace and coherent amplitude significantly reduce dephasing errors. The trade-off between the two major errors enables optimized performance of cat codes in terms of minimized decoherence. With high coupling efficiency, we show that one-way quantum repeaters with cat codes feature a boosted secure communication rate per mode when compared to conventional encoding schemes, showcasing the promising potential of quantum information processing with continuous variable quantum codes.

Solvable Family of Driven-Dissipative Many-Body Systems.

Exactly solvable models have played an important role in establishing the sophisticated modern understanding of equilibrium many-body physics. Conversely, the relative scarcity of solutions for nonequilibrium models greatly limits our understanding of systems away from thermal equilibrium. We study a family of nonequilibrium models, some of which can be viewed as dissipative analogues of the transverse-field Ising model, in that an effectively classical Hamiltonian is frustrated by dissipative processes that drive the system toward states that do not commute with the Hamiltonian. Surprisingly, a broad and experimentally relevant subset of these models can be solved efficiently. We leverage these solutions to compute the effects of decoherence on a canonical trapped-ion-based quantum computation architecture, and to prove a no-go theorem on steady-state phase transitions in a many-body model that can be realized naturally with Rydberg atoms or trapped ions.

Evolution of RLSB, a nuclear-encoded S1 domain RNA binding protein associated with post-transcriptional regulation of plastid-encoded rbcL mRNA in vascular plants.

BACKGROUND: RLSB, an S-1 domain RNA binding protein of Arabidopsis, selectively binds rbcL mRNA and co-localizes with Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) within chloroplasts of C3 and C4 plants. Previous studies using both Arabidopsis (C3) and maize (C4) suggest RLSB homologs are post-transcriptional regulators of plastid-encoded rbcL mRNA. While RLSB accumulates in all Arabidopsis leaf chlorenchyma cells, in C4 leaves RLSB-like proteins accumulate only within Rubisco-containing bundle sheath chloroplasts of Kranz-type species, and only within central compartment chloroplasts in the single cell C4 plant Bienertia. Our recent evidence implicates this mRNA binding protein as a primary determinant of rbcL expression, cellular localization/compartmentalization, and photosynthetic function in all multicellular green plants. This study addresses the hypothesis that RLSB is a highly conserved Rubisco regulatory factor that occurs in the chloroplasts all higher plants. RESULTS: Phylogenetic analysis has identified RLSB orthologs and paralogs in all major plant groups, from ancient liverworts to recent angiosperms. RLSB homologs were also identified in algae of the division Charophyta, a lineage closely related to land plants. RLSB-like sequences were not identified in any other algae, suggesting that it may be specific to the evolutionary line leading to land plants. The RLSB family occurs in single copy across most angiosperms, although a few species with two copies were identified, seemingly randomly distributed throughout the various taxa, although perhaps correlating in some cases with known ancient whole genome duplications. Monocots of the order Poales (Poaceae and Cyperaceae) were found to contain two copies, designated here as RLSB-a and RLSB-b, with only RLSB-a implicated in the regulation of rbcL across the maize developmental gradient. Analysis of microsynteny in angiosperms revealed high levels of conservation across eudicot species and for both paralogs in grasses, highlighting the possible importance of maintaining this gene and its surrounding genomic regions. CONCLUSIONS: Findings presented here indicate that the RLSB family originated as a unique gene in land plant evolution, perhaps in the common ancestor of charophytes and higher plants. Purifying selection has maintained this as a highly conserved single- or two-copy gene across most extant species, with several conserved gene duplications. Together with previous findings, this study suggests that RLSB has been sustained as an important regulatory protein throughout the course of land plant evolution. While only RLSB-a has been directly implicated in rbcL regulation in maize, RLSB-b could have an overlapping function in the co-regulation of rbcL, or may have diverged as a regulator of one or more other plastid-encoded mRNAs. This analysis confirms that RLSB is an important and unique photosynthetic regulatory protein that has been continuously expressed in land plants as they emerged and diversified from their ancient common ancestor.

High Gene Family Turnover Rates and Gene Space Adaptation in the Compact Genome of the Carnivorous Plant Utricularia gibba.

Utricularia gibba is an aquatic carnivorous plant with highly specialized morphology, featuring fibrous floating networks of branches and leaf-like organs, no recognizable roots, and bladder traps that capture and digest prey. We recently described the compressed genome of U. gibba as sufficient to control the development and reproduction of a complex organism. We hypothesized intense deletion pressure as a mechanism whereby most noncoding DNA was deleted, despite evidence for three independent whole-genome duplications (WGDs). Here, we explore the impact of intense genome fractionation in the evolutionary dynamics of U. gibba's functional gene space. We analyze U. gibba gene family turnover by modeling gene gain/death rates under a maximum-likelihood statistical framework. In accord with our deletion pressure hypothesis, we show that the U. gibba gene death rate is significantly higher than those of four other eudicot species. Interestingly, the gene gain rate is also significantly higher, likely reflecting the occurrence of multiple WGDs and possibly also small-scale genome duplications. Gene ontology enrichment analyses of U. gibba-specific two-gene orthogroups, multigene orthogroups, and singletons highlight functions that may represent adaptations in an aquatic carnivorous plant. We further discuss two homeodomain transcription factor gene families (WOX and HDG/HDZIP-IV) showing conspicuous differential expansions and contractions in U. gibba. Our results 1) reconcile the compactness of the U. gibba genome with its accommodation of a typical number of genes for a plant genome, and 2) highlight the role of high gene family turnover in the evolutionary diversification of U. gibba's functional gene space and adaptations to its unique lifestyle and highly specialized body plan.

The Fragaria vesca homolog of suppressor of overexpression of constans1 represses flowering and promotes vegetative growth.

In the annual long-day plant Arabidopsis thaliana, suppressor of overexpression of constans1 (SOC1) integrates endogenous and environmental signals to promote flowering. We analyzed the function and regulation of the SOC1 homolog (Fragaria vesca [Fv] SOC1) in the perennial short-day plant woodland strawberry (Fragaria vesca). We found that Fv SOC1 overexpression represses flower initiation under inductive short days, whereas its silencing causes continuous flowering in both short days and noninductive long days, similar to mutants in the floral repressor Fv terminal flower1 (Fv TFL1). Molecular analysis of these transgenic lines revealed that Fv SOC1 activates Fv TFL1 in the shoot apex, leading to the repression of flowering in strawberry. In parallel, Fv SOC1 regulates the differentiation of axillary buds to runners or axillary leaf rosettes, probably through the activation of gibberellin biosynthetic genes. We also demonstrated that Fv SOC1 is regulated by photoperiod and Fv flowering locus T1, suggesting that it plays a central role in the photoperiodic control of both generative and vegetative growth in strawberry. In conclusion, we propose that Fv SOC1 is a signaling hub that regulates yearly cycles of vegetative and generative development through separate genetic pathways.

Holonomic Quantum Control with Continuous Variable Systems.

Universal computation of a quantum system consisting of superpositions of well-separated coherent states of multiple harmonic oscillators can be achieved by three families of adiabatic holonomic gates. The first gate consists of moving a coherent state around a closed path in phase space, resulting in a relative Berry phase between that state and the other states. The second gate consists of "colliding" two coherent states of the same oscillator, resulting in coherent population transfer between them. The third gate is an effective controlled-phase gate on coherent states of two different oscillators. Such gates should be realizable via reservoir engineering of systems that support tunable nonlinearities, such as trapped ions and circuit QED.

Functional diversification of duplicated CYC2 clade genes in regulation of inflorescence development in Gerbera hybrida (Asteraceae).

The complex inflorescences (capitula) of Asteraceae consist of different types of flowers. In Gerbera hybrida (gerbera), the peripheral ray flowers are bilaterally symmetrical and lack functional stamens while the central disc flowers are more radially symmetrical and hermaphroditic. Proteins of the CYC2 subclade of the CYC/TB1-like TCP domain transcription factors have been recruited several times independently for parallel evolution of bilaterally symmetrical flowers in various angiosperm plant lineages, and have also been shown to regulate flower-type identity in Asteraceae. The CYC2 subclade genes in gerbera show largely overlapping gene expression patterns. At the level of single flowers, their expression domain in petals shows a spatial shift from the dorsal pattern known so far in species with bilaterally symmetrical flowers, suggesting that this change in expression may have evolved after the origin of Asteraceae. Functional analysis indicates that GhCYC2, GhCYC3 and GhCYC4 mediate positional information at the proximal-distal axis of the inflorescence, leading to differentiation of ray flowers, but that they also regulate ray flower petal growth by affecting cell proliferation until the final size and shape of the petals is reached. Moreover, our data show functional diversification for the GhCYC5 gene. Ectopic activation of GhCYC5 increases flower density in the inflorescence, suggesting that GhCYC5 may promote the flower initiation rate during expansion of the capitulum. Our data thus indicate that modification of the ancestral network of TCP factors has, through gene duplications, led to the establishment of new expression domains and to functional diversification.

Syntenic block overlap multiplicities with a panel of reference genomes provide a signature of ancient polyploidization events.

BACKGROUND: Following whole genome duplication (WGD), there is a compact distribution of gene similarities within the genome reflecting duplicate pairs of all the genes in the genome. With time, the distribution broadens and loses volume due to variable decay of duplicate gene similarity and to the process of duplicate gene loss. If there are two WGD, the older one becomes so reduced and broad that it merges with the tail of the distributions resulting from more recent events, and it becomes difficult to distinguish them. The goal of this paper is to advance statistical methods of identifying, or at least counting, the WGD events in the lineage of a given genome. METHODS: For a set of 15 angiosperm genomes, we analyze all 15 × 14 = 210 ordered pairs of target genome versus reference genome, using SynMap to find syntenic blocks. We consider all sets of B ≥ 2 syntenic blocks in the target genome that overlap in the reference genome as evidence of WGD activity in the target, whether it be one event or several. We hypothesize that in fitting an exponential function to the tail of the empirical distribution f (B) of block multiplicities, the size of the exponent will reflect the amount of WGD in the history of the target genome. RESULTS: By amalgamating the results from all reference genomes, a range of values of SynMap parameters, and alternative cutoff points for the tail, we find a clear pattern whereby multiple-WGD core eudicots have the smallest (negative) exponents, followed by core eudicots with only the single "γ" triplication in their history, followed by a non-core eudicot with a single WGD, followed by the monocots, with a basal angiosperm, the WGD-free Amborella having the largest exponent. CONCLUSION: The hypothesis that the exponent of the fit to the tail of the multiplicity distribution is a signature of the amount of WGD is verified, but there is also a clear complicating factor in the monocot clade, where a history of multiple WGD is not reflected in a small exponent.