Gamma paleohexaploidy in the stem lineage of core eudicots: significance for MADS-box gene and species diversification

Factors Influencing Gene Family Size Variation Among Related Species in a Plant Family, Solanaceae

Gene duplication and loss contribute to gene content differences as well as phenotypic divergence across species. However, the extent to which gene content varies among closely related plant species and the factors responsible for such variation remain unclear. Here, using the Solanaceae family as a model and Pfam domain families as a proxy for gene families, we investigated variation in gene family sizes across species and the likely factors contributing to the variation. We found that genes in highly variable families have high turnover rates and tend to be involved in processes that have diverged between Solanaceae species, whereas genes in low-variability families tend to have housekeeping roles. In addition, genes in high- and low-variability gene families tend to be duplicated by tandem and whole genome duplication, respectively. This finding together with the observation that genes duplicated by different mechanisms experience different selection pressures suggest that duplication mechanism impacts gene family turnover. We explored using pseudogene number as a proxy for gene loss but discovered that a substantial number of pseudogenes are actually products of pseudogene duplication, contrary to the expectation that most plant pseudogenes are remnants of once-functional duplicates. Our findings reveal complex relationships between variation in gene family size, gene functions, duplication mechanism, and evolutionary rate. The patterns of lineage-specific gene family expansion within the Solanaceae provide the foundation for a better understanding of the genetic basis underlying phenotypic diversity in this economically important family.

Two Likely Auto-Tetraploidization Events Shaped Kiwifruit Genome and Contributed to Establishment of the Actinidiaceae Family

The genome of kiwifruit (Actinidia chinensis) was sequenced previously, the first in the Actinidiaceae family. It was shown to have been affected by polyploidization events, the nature of which has been elusive. Here, we performed a reanalysis of the genome and found clear evidence of 2 tetraploidization events, with one occurring ∼50–57 million years ago (Mya) and the other ∼18–20 Mya. Two subgenomes produced by each event have been under balanced fractionation. Moreover, genes were revealed to express in a balanced way between duplicated copies of chromosomes. Besides, lowered evolutionary rates of kiwifruit genes were observed. These findings could be explained by the likely auto-tetraploidization nature of the polyploidization events. Besides, we found that polyploidy contributed to the expansion of key functional genes, e.g., vitamin C biosynthesis genes. The present work also provided an important comparative genomics resource in the Actinidiaceae and related families.

•

Two independent paleo-tetraploidization events may have occurred in Actinidiaceae

•

The tetraploidization events are likely autotetraploid in nature

•

These events contribute to the expansion of key trait genes

•

Hierarchical deconvolution allowed analysis of the kiwifruit genome interweaving homology

TM8 represses developmental timing in Nicotiana benthamiana and has functionally diversified in angiosperms

BACKGROUND

MADS-box genes are key regulators of plant reproductive development and members of most lineages of this gene family have been extensively studied. However, the function and diversification of the ancient TM8 lineage remains elusive to date. The available data suggest a possible function in flower development in tomato and fast evolution through numerous gene loss events in flowering plants.

RESULTS

We show the broad conservation of TM8 within angiosperms and find that in contrast to other MADS-box gene lineages, no gene duplicates have been retained after major whole genome duplication events. Through knock-down of NbTM8 by virus induced gene silencing in Nicotiana benthamiana, we show that NbTM8 represses miR172 together with another MADS-box gene, SHORT VEGETATIVE PHASE (NbSVP). In the closely related species Petunia hybrida, PhTM8 is not expressed under the conditions we investigated and consistent with this, a knock-out mutant did not show a phenotype. Finally, we generated transgenic tomato plants in which TM8 was silenced or ectopically expressed, but these plants did not display a clear phenotype. Therefore, no clear function could be confirmed for Solanum lycopersium.

CONCLUSIONS

While the presence of TM8 is generally conserved, it remains difficult to propose a general function in angiosperms. Based on all the available data to date, supplemented with our own results, TM8 function seems to have diversified quickly throughout angiosperms and acts as repressor of miR172 in Nicotiana benthamiana, together with NbSVP.

Genomic Changes Associated with the Evolutionary Transitions of Nostoc to a Plant Symbiont

Cyanobacteria belonging to the genus Nostoc comprise free-living strains and also facultative plant symbionts. Symbiotic strains can enter into symbiosis with taxonomically diverse range of host plants. Little is known about genomic changes associated with evolutionary transition of Nostoc from free-living to plant symbiont. Here, we compared the genomes derived from 11 symbiotic Nostoc strains isolated from different host plants and infer phylogenetic relationships between strains. Phylogenetic reconstructions of 89 Nostocales showed that symbiotic Nostoc strains with a broad host range, entering epiphytic and intracellular or extracellular endophytic interactions, form a monophyletic clade indicating a common evolutionary history. A polyphyletic origin was found for Nostoc strains which enter only extracellular symbioses, and inference of transfer events implied that this trait was likely acquired several times in the evolution of the Nostocales. Symbiotic Nostoc strains showed enriched functions in transport and metabolism of organic sulfur, chemotaxis and motility, as well as the uptake of phosphate, branched-chain amino acids, and ammonium. The genomes of the intracellular clade differ from that of other Nostoc strains, with a gain/enrichment of genes encoding proteins to generate l-methionine from sulfite and pathways for the degradation of the plant metabolites vanillin and vanillate, and of the macromolecule xylan present in plant cell walls. These compounds could function as C-sources for members of the intracellular clade. Molecular clock analysis indicated that the intracellular clade emerged ca. 600 Ma, suggesting that intracellular Nostoc symbioses predate the origin of land plants and the emergence of their extant hosts.

Macroevolutionary Patterns of Flowering Plant Speciation and Extinction

Species diversity is remarkably unevenly distributed among flowering plant lineages. Despite a growing toolbox of research methods, the reasons underlying this patchy pattern have continued to perplex plant biologists for the past two decades. In this review, we examine the present understanding of transitions in flowering plant evolution that have been proposed to influence speciation and extinction. In particular, ploidy changes, transitions between tropical and nontropical biomes, and shifts in floral form have received attention and have offered some surprises in terms of which factors influence speciation and extinction rates. Mating systems and dispersal characteristics once predominated as determining factors, yet recent evidence suggests that these changes are not as influential as previously thought or are important only when paired with range shifts. Although range extent is an important correlate of speciation, it also influences extinction and brings an applied focus to diversification research. Recent studies that find that past diversification can predict present-day extinction risk open an exciting avenue for future research to help guide conservation prioritization.

Phylotranscriptomic analysis and genome evolution of the Cypripedioideae (Orchidaceae)

PREMISE OF THE STUDY

The slipper orchids (Cypripedioideae) are a morphologically distinct subfamily of Orchidaceae. They also have some of the largest genomes in the orchids, which may be due to polyploidy or some other mechanism of genome evolution. We generated 10 transcriptomes and incorporated existing RNA-seq data to infer a multilocus nuclear phylogeny of the Cypripedioideae and to determine whether a whole-genome duplication event (WGD) correlated with the large genome size of this subfamily. Knowing more about timing of ancient polyploidy events can help us understand the evolution of one of the most species-rich plant families.

METHODS

Transcriptome data were used to identify low-copy orthologous genes to infer a phylogeny of Orchidaceae and to identify paralogs to place any WGD events on the species tree.

KEY RESULTS

Our transcriptome phylogeny confirmed relationships published in previous studies that used fewer markers but incorporated more taxa. We did not find a WGD event at the base of the slipper orchids; however, we did identify one on the Orchidaceae stem lineage. We also confirmed the presence of a previously identified WGD event deeper in the monocot phylogeny.

CONCLUSIONS

Although WGD has played a role in the evolution of Orchidaceae, polyploidy does not appear to be responsible for the large genome size of slipper orchids. The conserved set of 775 largely single-copy nuclear genes identified in this study should prove useful in future studies of orchid evolution.

Phylogenomics reveals an extensive history of genome duplication in diatoms (Bacillariophyta)

PREMISE OF THE STUDY

Diatoms are one of the most species-rich lineages of microbial eukaryotes. Similarities in clade age, species richness, and primary productivity motivate comparisons to angiosperms, whose genomes have been inordinately shaped by whole-genome duplication (WGD). WGDs have been linked to speciation, increased rates of lineage diversification, and identified as a principal driver of angiosperm evolution. We synthesized a large but scattered body of evidence that suggests polyploidy may be common in diatoms as well.

METHODS

We used gene counts, gene trees, and distributions of synonymous divergence to carry out a phylogenomic analysis of WGD across a diverse set of 37 diatom species.

KEY RESULTS

Several methods identified WGDs of varying age across diatoms. Determining the occurrence, exact number, and placement of events was greatly impacted by uncertainty in gene trees. WGDs inferred from synonymous divergence of paralogs varied depending on how redundancy in transcriptomes was assessed, gene families were assembled, and synonymous distances (Ks) were calculated. Our results highlighted a need for systematic evaluation of key methodological aspects of Ks-based approaches to WGD inference. Gene tree reconciliations supported allopolyploidy as the predominant mode of polyploid formation, with strong evidence for ancient allopolyploid events in the thalassiosiroid and pennate diatom clades.

CONCLUSIONS

Our results suggest that WGD has played a major role in the evolution of diatom genomes. We outline challenges in reconstructing paleopolyploid events in diatoms that, together with these results, offer a framework for understanding the impact of genome duplication in a group that likely harbors substantial genomic diversity.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

Evolution of the process underlying floral zygomorphy development in pentapetalous angiosperms

PREMISE OF THE STUDY

Observations of floral ontogeny indicated that floral organ initiation in pentapetalous flowers most commonly results in a median-abaxial (MAB) petal during early development, a median-adaxial (MAD) petal being less common. Such different patterns of floral organ initiation might be linked with different morphologies of floral zygomorphy that have evolved in Asteridae. Here, we provide the first study of zygomorphy in pentapetalous angiosperms placed in a phylogenetic framework, the goal being to find if the different patterns of floral organ initiation are connected with particular patterns of zygomorphy.

METHODS

We analyzed patterns of floral organ initiation and displays of zygomorphy, extracted from floral diagrams representing 405 taxa in 330 genera, covering 83% of orders (30 out of 36) and 37% of families (116 out of 313) in core eudicots in the context of a phylogeny using ancestral state reconstructions.

KEY RESULTS

The MAB petal initiation is the ancestral state of the pattern of floral organ initiation in pentapetalous angiosperms. Taxa with MAD petal initiation represent ∼30 independent origins from the ancestral MAB initiation. There are distinct developmental processes that give rise to zygomorphy in different lineages of pentapetalous angiosperms, closely related lineages being likely to share similar developmental processes.

CONCLUSIONS

We have demonstrated that development indeed constrains the processes that give rise to floral zygomorphy, while phylogenetic distance allows relaxation of these constraints, which provides novel insights on the role that development plays in the evolution of floral zygomorphy.

Evolution of floral diversity: genomics, genes and gamma

A salient feature of flowering plant diversification is the emergence of a novel suite of floral features coinciding with the origin of the most species-rich lineage, Pentapetalae. Advances in phylogenetics, developmental genetics and genomics, including new analyses presented here, are helping to reconstruct the specific evolutionary steps involved in the evolution of this clade. The enormous floral diversity among Pentapetalae appears to be built on a highly conserved ground plan of five-parted (pentamerous) flowers with whorled phyllotaxis. By contrast, lability in the number and arrangement of component parts of the flower characterize the early-diverging eudicot lineages subtending Pentapetalae. The diversification of Pentapetalae also coincides closely with ancient hexaploidy, referred to as the gamma whole-genome triplication, for which the phylogenetic timing, mechanistic details and molecular evolutionary consequences are as yet not fully resolved. Transcription factors regulating floral development often persist in duplicate or triplicate in gamma-derived genomes, and both individual genes and whole transcriptional programmes exhibit a shift from broadly overlapping to tightly defined expression domains in Pentapetalae flowers. Investigations of these changes associated with the origin of Pentapetalae can lead to a more comprehensive understanding of what is arguably one of the most important evolutionary diversification events within terrestrial plants.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Revisiting ancestral polyploidy in plants

Whole-genome duplications (WGDs) or polyploidy events have been studied extensively in plants. In a now widely cited paper, Jiao et al. presented evidence for two ancient, ancestral plant WGDs predating the origin of flowering and seed plants, respectively. This finding was based primarily on a bimodal age distribution of gene duplication events obtained from molecular dating of almost 800 phylogenetic gene trees. We reanalyzed the phylogenomic data of Jiao et al. and found that the strong bimodality of the age distribution may be the result of technical and methodological issues and may hence not be a "true" signal of two WGD events. By using a state-of-the-art molecular dating algorithm, we demonstrate that the reported bimodal age distribution is not robust and should be interpreted with caution. Thus, there exists little evidence for two ancient WGDs in plants from phylogenomic dating.

The evolutionary significance of polyploidy

Polyploidy, or the duplication of entire genomes, has been observed in prokaryotic and eukaryotic organisms, and in somatic and germ cells. The consequences of polyploidization are complex and variable, and they differ greatly between systems (clonal or non-clonal) and species, but the process has often been considered to be an evolutionary 'dead end'. Here, we review the accumulating evidence that correlates polyploidization with environmental change or stress, and that has led to an increased recognition of its short-term adaptive potential. In addition, we discuss how, once polyploidy has been established, the unique retention profile of duplicated genes following whole-genome duplication might explain key longer-term evolutionary transitions and a general increase in biological complexity.

Resolution of deep eudicot phylogeny and their temporal diversification using nuclear genes from transcriptomic and genomic datasets

Explosive diversification is widespread in eukaryotes, making it difficult to resolve phylogenetic relationships. Eudicots contain c. 75% of extant flowering plants, are important for human livelihood and terrestrial ecosystems, and have probably experienced explosive diversifications. The eudicot phylogenetic relationships, especially among those of the Pentapetalae, remain unresolved. Here, we present a highly supported eudicot phylogeny and diversification rate shifts using 31 newly generated transcriptomes and 88 other datasets covering 70% of eudicot orders. A highly supported eudicot phylogeny divided Pentapetalae into two groups: one with rosids, Saxifragales, Vitales and Santalales; the other containing asterids, Caryophyllales and Dilleniaceae, with uncertainty for Berberidopsidales. Molecular clock analysis estimated that crown eudicots originated c. 146 Ma, considerably earlier than earliest tricolpate pollen fossils and most other molecular clock estimates, and Pentapetalae sequentially diverged into eight major lineages within c. 15 Myr. Two identified increases of diversification rate are located in the stems leading to Pentapetalae and asterids, and lagged behind the gamma hexaploidization. The nuclear genes from newly generated transcriptomes revealed a well-resolved eudicot phylogeny, sequential separation of major core eudicot lineages and temporal mode of diversifications, providing new insights into the evolutionary trend of morphologies and contributions to the diversification of eudicots.

Protein interaction evolution from promiscuity to specificity with reduced flexibility in an increasingly complex network

A key question regarding protein evolution is how proteins adapt to the dynamic environment in which they function and how in turn their evolution shapes the protein interaction network. We used extant and resurrected ancestral plant MADS-domain transcription factors to understand how SEPALLATA3, a protein with hub and glue properties, evolved and takes part in network organization. Although the density of dimeric interactions was saturated in the network, many new interactions became mediated by SEPALLATA3 after a whole genome triplication event. By swapping SEPALLATA3 and its ancestors between dimeric networks of different ages, we found that the protein lost the capacity of promiscuous interaction and acquired specificity in evolution. This was accompanied with constraints on conformations through proline residue accumulation, which made the protein less flexible. SHORT VEGETATIVE PHASE on the other hand (non-hub) was able to gain protein-protein interactions due to a C-terminal domain insertion, allowing for a larger interaction interface. These findings illustrate that protein interaction evolution occurs at the level of conformational dynamics, when the binding mechanism concerns an induced fit or conformational selection. Proteins can evolve towards increased specificity with reduced flexibility when the complexity of the protein interaction network requires specificity.

The evolution of plant microRNAs: insights from a basal eudicot sacred lotus

microRNAs (miRNAs) are important noncoding small RNAs that regulate mRNAs in eukaryotes. However, under which circumstances different miRNAs/miRNA families exhibit different evolutionary trajectories in plants remains unclear. In this study, we sequenced the small RNAs and degradome from a basal eudicot, sacred lotus (Nelumbo nucifera or lotus), to identify miRNAs and their targets. Combining with public miRNAs, we predicted 57 pre-eudicot miRNA families from different evolutionary stages. We found that miRNA families featuring older age, higher copy and target number tend to show lower propensity for miRNA family loss (PGL) and stronger signature of purifying selection during divergence of temperate and tropical lotus. Further analyses of lotus genome revealed that there is an association between loss of miRNA families in descendent plants and in duplicated genomes. Gene dosage balance is crucial in maintaining those preferentially retained MIRNA duplicates by imposing stronger purifying selection. However, these factors and selection influencing miRNA family evolution are not applicable to the putative MIRNA-likes. Additionally, the MIRNAs participating in lotus pollen-pistil interaction, a conserved process in angiosperms, also have a strong signature of purifying selection. Functionally, sequence divergence in MIRNAs escalates expression divergence of their target genes between temperate and tropical lotus during rhizome and leaf growth. Overall, our study unravels several important factors and selection that determine the miRNA family distribution in plants and duplicated genomes, and provides evidence for functional impact of MIRNA sequence evolution.

Most Compositae (Asteraceae) are descendants of a paleohexaploid and all share a paleotetraploid ancestor with the Calyceraceae

PREMISE OF THE STUDY

Like many other flowering plants, members of the Compositae (Asteraceae) have a polyploid ancestry. Previous analyses found evidence for an ancient duplication or possibly triplication in the early evolutionary history of the family. We sought to better place this paleopolyploidy in the phylogeny and assess its nature.

METHODS

We sequenced new transcriptomes for Barnadesia, the lineage sister to all other Compositae, and four representatives of closely related families. Using a recently developed algorithm, MAPS, we analyzed nuclear gene family phylogenies for evidence of paleopolyploidy.

KEY RESULTS

We found that the previously recognized Compositae paleopolyploidy is also in the ancestry of the Calyceraceae. Our phylogenomic analyses uncovered evidence for a successive second round of genome duplication among all sampled Compositae except Barnadesia.

CONCLUSIONS

Our analyses of new samples with new tools provide a revised view of paleopolyploidy in the Compositae. Together with results from a high density Lactuca linkage map, our results suggest that the Compositae and Calyceraceae have a common paleotetraploid ancestor and that most Compositae are descendants of a paleohexaploid. Although paleohexaploids have been previously identified, this is the first example where the paleotetraploid and paleohexaploid lineages have survived over tens of millions of years. The complex polyploidy in the ancestry of the Compositae and Calyceraceae represents a unique opportunity to study the long-term evolutionary fates and consequences of different ploidal levels.

Epigenetic regulation of subgenome dominance following whole genome triplication in Brassica rapa

Subgenome dominance is an important phenomenon observed in allopolyploids after whole genome duplication, in which one subgenome retains more genes as well as contributes more to the higher expressing gene copy of paralogous genes. To dissect the mechanism of subgenome dominance, we systematically investigated the relationships of gene expression, transposable element (TE) distribution and small RNA targeting, relating to the multicopy paralogous genes generated from whole genome triplication in Brassica rapa. The subgenome dominance was found to be regulated by a relatively stable factor established previously, then inherited by and shared among B. rapa varieties. In addition, we found a biased distribution of TEs between flanking regions of paralogous genes. Furthermore, the 24-nt small RNAs target TEs and are negatively correlated to the dominant expression of individual paralogous gene pairs. The biased distribution of TEs among subgenomes and the targeting of 24-nt small RNAs together produce the dominant expression phenomenon at a subgenome scale. Based on these findings, we propose a bucket hypothesis to illustrate subgenome dominance and hybrid vigor. Our findings and hypothesis are valuable for the evolutionary study of polyploids, and may shed light on studies of hybrid vigor, which is common to most species.

Single nucleotide polymorphism discovery via genotyping by sequencing to assess population genetic structure and recurrent polyploidization in Andropogon gerardii

PREMISE OF THE STUDY

Autopolyploidy, genome duplication within a single lineage, can result in multiple cytotypes within a species. Geographic distributions of cytotypes may reflect the evolutionary history of autopolyploid formation and subsequent population dynamics including stochastic (drift) and deterministic (differential selection among cytotypes) processes. Here, we used a population genomic approach to investigate whether autopolyploidy occurred once or multiple times in Andropogon gerardii, a widespread, North American grass with two predominant cytotypes.

METHODS

Genotyping by sequencing was used to identify single nucleotide polymorphisms (SNPs) in individuals collected from across the geographic range of A. gerardii. Two independent approaches to SNP calling were used: the reference-free UNEAK pipeline and a reference-guided approach based on the sequenced Sorghum bicolor genome. SNPs generated using these pipelines were analyzed independently with genetic distance and clustering.

KEY RESULTS

Analyses of the two SNP data sets showed very similar patterns of population-level clustering of A. gerardii individuals: a cluster of A. gerardii individuals from the southern Plains, a northern Plains cluster, and a western cluster. Groupings of individuals corresponded to geographic localities regardless of cytotype: 6x and 9x individuals from the same geographic area clustered together.

CONCLUSIONS

SNPs generated using reference-guided and reference-free pipelines in A. gerardii yielded unique subsets of genomic data. Both data sets suggest that the 9x cytotype in A. gerardii likely evolved multiple times from 6x progenitors across the range of the species. Genomic approaches like GBS and diverse bioinformatics pipelines used here facilitate evolutionary analyses of complex systems with multiple ploidy levels.

Large-Scale Gene Relocations following an Ancient Genome Triplication Associated with the Diversification of Core Eudicots

Different modes of gene duplication including whole-genome duplication (WGD), and tandem, proximal and dispersed duplications are widespread in angiosperm genomes. Small-scale, stochastic gene relocations and transposed gene duplications are widely accepted to be the primary mechanisms for the creation of dispersed duplicates. However, here we show that most surviving ancient dispersed duplicates in core eudicots originated from large-scale gene relocations within a narrow window of time following a genome triplication (γ) event that occurred in the stem lineage of core eudicots. We name these surviving ancient dispersed duplicates as relocated γ duplicates. In Arabidopsis thaliana, relocated γ, WGD and single-gene duplicates have distinct features with regard to gene functions, essentiality, and protein interactions. Relative to γ duplicates, relocated γ duplicates have higher non-synonymous substitution rates, but comparable levels of expression and regulation divergence. Thus, relocated γ duplicates should be distinguished from WGD and single-gene duplicates for evolutionary investigations. Our results suggest large-scale gene relocations following the γ event were associated with the diversification of core eudicots.

Of dups and dinos: evolution at the K/Pg boundary

Fifteen years into sequencing entire plant genomes, more than 30 paleopolyploidy events could be mapped on the tree of flowering plants (and many more when also transcriptome data sets are considered). While some genome duplications are very old and have occurred early in the evolution of dicots and monocots, or even before, others are more recent and seem to have occurred independently in many different plant lineages. Strikingly, a majority of these duplications date somewhere between 55 and 75 million years ago (mya), and thus likely correlate with the K/Pg boundary. If true, this would suggest that plants that had their genome duplicated at that time, had an increased chance to survive the most recent mass extinction event, at 66mya, which wiped out a majority of plant and animal life, including all non-avian dinosaurs. Here, we review several processes, both neutral and adaptive, that might explain the establishment of polyploid plants, following the K/Pg mass extinction.

Flower development: the evolutionary history and functions of the AGL6 subfamily MADS-box genes

AGL6 is an ancient subfamily of MADS-box genes found in both gymnosperms and angiosperms. Its functions remained elusive despite the fact that the MADS-box genes and the ABC model have been studied for >20 years. Nevertheless, recent discoveries in petunia, rice, and maize support its involvement in the 'E' function of floral development, very similar to the closely related AGL2 (SEPALLATA) subfamily which has been well characterized. The known functions of AGL6 span from ancient conserved roles to new functions acquired in specific plant families. The AGL6 genes are involved in floral meristem regulation, in floral organs, and ovule (integument) and seed development, and have possible roles in both male and female germline and gametophyte development. In grasses, they are also important for the development of the first whorl of the flower, whereas in Arabidopsis they may play additional roles before floral meristem formation. This review covers these recent insights and some other aspects that are not yet fully elucidated, which deserve more studies in the future.

Polyploidy and genome evolution in plants

Plant genomes vary in size and complexity, fueled in part by processes of whole-genome duplication (WGD; polyploidy) and subsequent genome evolution. Despite repeated episodes of WGD throughout the evolutionary history of angiosperms in particular, the genomes are not uniformly large, and even plants with very small genomes carry the signatures of ancient duplication events. The processes governing the evolution of plant genomes following these ancient events are largely unknown. Here, we consider mechanisms of diploidization, evidence of genome reorganization in recently formed polyploid species, and macroevolutionary patterns of WGD in plant genomes and propose that the ongoing genomic changes observed in recent polyploids may illustrate the diploidization processes that result in ancient signatures of WGD over geological timescales.

Large-Scale Evolutionary Analysis of Genes and Supergene Clusters from Terpenoid Modular Pathways Provides Insights into Metabolic Diversification in Flowering Plants

An important component of plant evolution is the plethora of pathways producing more than 200,000 biochemically diverse specialized metabolites with pharmacological, nutritional and ecological significance. To unravel dynamics underlying metabolic diversification, it is critical to determine lineage-specific gene family expansion in a phylogenomics framework. However, robust functional annotation is often only available for core enzymes catalyzing committed reaction steps within few model systems. In a genome informatics approach, we extracted information from early-draft gene-space assemblies and non-redundant transcriptomes to identify protein families involved in isoprenoid biosynthesis. Isoprenoids comprise terpenoids with various roles in plant-environment interaction, such as pollinator attraction or pathogen defense. Combining lines of evidence provided by synteny, sequence homology and Hidden-Markov-Modelling, we screened 17 genomes including 12 major crops and found evidence for 1,904 proteins associated with terpenoid biosynthesis. Our terpenoid genes set contains evidence for 840 core terpene-synthases and 338 triterpene-specific synthases. We further identified 190 prenyltransferases, 39 isopentenyl-diphosphate isomerases as well as 278 and 219 proteins involved in mevalonate and methylerithrol pathways, respectively. Assessing the impact of gene and genome duplication to lineage-specific terpenoid pathway expansion, we illustrated key events underlying terpenoid metabolic diversification within 250 million years of flowering plant radiation. By quantifying Angiosperm-wide versatility and phylogenetic relationships of pleiotropic gene families in terpenoid modular pathways, our analysis offers significant insight into evolutionary dynamics underlying diversification of plant secondary metabolism. Furthermore, our data provide a blueprint for future efforts to identify and more rapidly clone terpenoid biosynthetic genes from any plant species.

Adaptation in flower form: a comparative evodevo approach

Evolutionary developmental biology (evodevo) attempts to explain how the process of organismal development evolves, utilizing a comparative approach to investigate changes in developmental pathways and processes that occur during the evolution of a given lineage. Evolutionary genetics uses a population approach to understand how organismal changes in form or function are linked to underlying genetics, focusing on changes in gene and genotype frequencies within populations and the fixation of genotypic variation into traits that define species or evoke speciation events. Microevolutionary processes, including mutation, genetic drift, natural selection and gene flow, can provide the foundation for macroevolutionary patterns observed as morphological evolution and adaptation. The temporal element linking microevolutionary processes to macroevolutionary patterns is development: an organism's genotype is converted to phenotype by ontogenetic processes. Because selection acts upon the phenotype, the connection between evolutionary genetics and developmental evolution becomes essential to understanding adaptive evolution in organismal form and function. Here, we discuss how developmental genetic studies focused on key developmental processes could be linked within a comparative framework to study the developmental genetics of adaptive evolution, providing examples from research on two key processes of plant evodevo - floral symmetry and organ fusion - and their role in the adaptation of floral form.

Polyploidy-associated genome modifications during land plant evolution

The occurrence of polyploidy in land plant evolution has led to an acceleration of genome modifications relative to other crown eukaryotes and is correlated with key innovations in plant evolution. Extensive genome resources provide for relating genomic changes to the origins of novel morphological and physiological features of plants. Ancestral gene contents for key nodes of the plant family tree are inferred. Pervasive polyploidy in angiosperms appears likely to be the major factor generating novel angiosperm genes and expanding some gene families. However, most gene families lose most duplicated copies in a quasi-neutral process, and a few families are actively selected for single-copy status. One of the great challenges of evolutionary genomics is to link genome modifications to speciation, diversification and the morphological and/or physiological innovations that collectively compose biodiversity. Rapid accumulation of genomic data and its ongoing investigation may greatly improve the resolution at which evolutionary approaches can contribute to the identification of specific genes responsible for particular innovations. The resulting, more 'particulate' understanding of plant evolution, may elevate to a new level fundamental knowledge of botanical diversity, including economically important traits in the crop plants that sustain humanity.

Tangled up in two: a burst of genome duplications at the end of the Cretaceous and the consequences for plant evolution

Genome sequencing has demonstrated that besides frequent small-scale duplications, large-scale duplication events such as whole genome duplications (WGDs) are found on many branches of the evolutionary tree of life. Especially in the plant lineage, there is evidence for recurrent WGDs, and the ancestor of all angiosperms was in fact most likely a polyploid species. The number of WGDs found in sequenced plant genomes allows us to investigate questions about the roles of WGDs that were hitherto impossible to address. An intriguing observation is that many plant WGDs seem associated with periods of increased environmental stress and/or fluctuations, a trend that is evident for both present-day polyploids and palaeopolyploids formed around the Cretaceous–Palaeogene (K–Pg) extinction at 66 Ma. Here, we revisit the WGDs in plants that mark the K–Pg boundary, and discuss some specific examples of biological innovations and/or diversifications that may be linked to these WGDs. We review evidence for the processes that could have contributed to increased polyploid establishment at the K–Pg boundary, and discuss the implications on subsequent plant evolution in the Cenozoic.

A complex interplay of tandem- and whole-genome duplication drives expansion of the L-type lectin receptor kinase gene family in the brassicaceae

The comparative analysis of plant gene families in a phylogenetic framework has greatly accelerated due to advances in next generation sequencing. In this study, we provide an evolutionary analysis of the L-type lectin receptor kinase and L-type lectin domain proteins (L-type LecRKs and LLPs) that are considered as components in plant immunity, in the plant family Brassicaceae and related outgroups. We combine several lines of evidence provided by sequence homology, HMM-driven protein domain annotation, phylogenetic analysis, and gene synteny for large-scale identification of L-type LecRK and LLP genes within nine core-eudicot genomes. We show that both polyploidy and local duplication events (tandem duplication and gene transposition duplication) have played a major role in L-type LecRK and LLP gene family expansion in the Brassicaceae. We also find significant differences in rates of molecular evolution based on the mode of duplication. Additionally, we show that LLPs share a common evolutionary origin with L-type LecRKs and provide a consistent gene family nomenclature. Finally, we demonstrate that the largest and most diverse L-type LecRK clades are lineage-specific. Our evolutionary analyses of these plant immune components provide a framework to support future plant resistance breeding.

A novel approach for multi-domain and multi-gene family identification provides insights into evolutionary dynamics of disease resistance genes in core eudicot plants

Background

Recent advances in DNA sequencing techniques resulted in more than forty sequenced plant genomes representing a diverse set of taxa of agricultural, energy, medicinal and ecological importance. However, gene family curation is often only inferred from DNA sequence homology and lacks insights into evolutionary processes contributing to gene family dynamics. In a comparative genomics framework, we integrated multiple lines of evidence provided by gene synteny, sequence homology and protein-based Hidden Markov Modelling to extract homologous super-clusters composed of multi-domain resistance (R)-proteins of the NB-LRR type (for NUCLEOTIDE BINDING/LEUCINE-RICH REPEATS), that are involved in plant innate immunity.

Results

To assess the diversity of R-proteins within and between species, we screened twelve eudicot plant genomes including six major crops and found a total of 2,363 NB-LRR genes. Our curated R-proteins set shows a 50% average for tandem duplicates and a 22% fraction of gene copies retained from ancient polyploidy events (ohnologs). We provide evidence for strong positive selection and show significant differences in molecular evolution rates (Ka/Ks-ratio) among tandem- (mean = 1.59), ohnolog (mean = 1.36) and singleton (mean = 1.22) R-gene duplicates. To foster the process of gene-edited plant breeding, we report species-specific presence/absence of all 140 NB-LRR genes present in the model plant Arabidopsis and describe four distinct clusters of NB-LRR “gatekeeper” loci sharing syntenic orthologs across all analyzed genomes.

Conclusion

By curating a near-complete set of multi-domain R-protein clusters in an eudicot-wide scale, our analysis offers significant insight into evolutionary dynamics underlying diversification of the plant innate immune system. Furthermore, our methods provide a blueprint for future efforts to identify and more rapidly clone functional NB-LRR genes from any plant species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-966) contains supplementary material, which is available to authorized users.

Modeling the evolution of molecular systems from a mechanistic perspective

Systems biology-inspired genotype-phenotype mapping models are increasingly being used to study the evolutionary properties of molecular biological systems, in particular the general emergent properties of evolving systems, such as modularity, robustness, and evolvability. However, the level of abstraction at which many of these models operate might not be sufficient to capture all relevant intricacies of biological evolution in sufficient detail. Here, we argue that in particular gene and genome duplications, both evolutionary mechanisms of potentially major importance for the evolution of molecular systems and of special relevance to plant evolution, are not adequately accounted for in most GPM modeling frameworks, and that more fine-grained mechanistic models may significantly advance understanding of how gen(om)e duplication impacts molecular systems evolution.

Functional and evolutionary analysis of the AP1/SEP/AGL6 superclade of MADS-box genes in the basal eudicot Epimedium sagittatum

BACKGROUND AND AIMS

MADS-box transcriptional regulators play important roles during plant development. Based on phylogenetic reconstruction, the AP1/SEP/AGL6 superclade of floral MADS-box genes underwent one or two duplication events in the common ancestor of the core eudicots. However, the functional evolution of the AP1/SEP/AGL6 superclade in basal eudicots remains uncharacterized. Epimedium sagittatum is a basal eudicot species valued for its medicinal properties and showing unique floral morphology. In this study, structural and functional variation of FUL-like (AP1 subfamily), SEP-like and AGL6-like genes in this species was investigated to further our understanding of flower evolution in angiosperms. Detailed investigations into the microsynteny and evolutionary history of the floral A and E class MADS-box genes in eudicots were undertaken and used to trace their genomic rearrangements.

METHODS

One AP1-like gene, two SEP-like genes and one AGL6-like gene were cloned from E. sagittatum. Their expression patterns were examined using quantitative RT-PCR in different vegetative and reproductive organs at two developmental stages. Yeast two-hybrid assays were carried out among AP1/SEP/AGL6 superclade, AP3/PI and AGAMOUS subfamily members for elucidation of dimerization patterns. In addition, possible formation of a ternary complex involving B class proteins with the A class protein EsFUL-like, the E class SEP-like protein EsAGL2-1 or the AGL6-class protein EsAGL6 were detected using yeast three-hybrid assays. Transgenic Arabidopsis or tobacco plants expressing EsFUL-like, EsAGL2-1 and EsAGL6-like under the cauliflower mosaic virus (CaMV) 35S promoter were generated and analysed. Genomic studies of AP1 syntenic regions in arabidopsis, columbine, strawberry, papaya, peach, grapevine and tomato were conducted for microsyntenic analyses.

KEY RESULTS

Sequence and phylogenetic analyses showed that EsFUL-like is a member of the AP1 (A class) subfamily, EsAGL2-1 and EsAGL2-2 belong to the SEP-like (E class) subfamily, and EsAGL6-like belongs to the AGL6 (AGL6 class) subfamily. Quantitative RT-PCR analyses revealed that the transcripts of the four genes are absent, or minimal, in vegetative tissues and are most highly expressed in floral organs. Yeast two-hybrid results revealed that of the eight MADS-box proteins tested, only EsAGL6-like, EsAGL2-1 and EsAGL2 were able to form strong homo- and heterodimers, with EsAGL6-like and EsAGL2-1 showing similar interaction patterns. Yeast three-hybrid analysis revealed that EsFUL1-like, EsAGL6-like and EsAGL2-1 (representing the three major lineages of the Epimedium AGL/SEP/ALG6 superclade) could act as bridging proteins in ternary complexes with both EsAP3-2 (B class) and EsPI (B class), which do not heterodimerize themselves. Syntenic analyses of sequenced basal eudicots, rosids and asterids showed that most AP1-like and SEP-like genes have been tightly associated as neighbours since the origin of basal eudicots. Ectopic expression of EsFUL-like in arabidopsis caused early flowering through endogenous high-level expression of AP1 and formation of secondary flowers between the first and second whorls. Tobacco plants with ectopic expression of EsAGL2-1 showed shortened pistils and styles, as well as axillary and extra petals in the initial flower.

CONCLUSIONS

This study provides a description of EsFUL-like, EsAGL2-1, EsAGL2-2 and EsAGL6-like function divergence and conservation in comparison with a selection of model core eudicots. The study also highlights how organization in genomic segments containing A and E class genes in sequenced model species has resulted in similar topologies of AP1 and SEP-like gene trees.

TRAPID: an efficient online tool for the functional and comparative analysis of de novo RNA-Seq transcriptomes

Transcriptome analysis through next-generation sequencing technologies allows the generation of detailed gene catalogs for non-model species, at the cost of new challenges with regards to computational requirements and bioinformatics expertise. Here, we present TRAPID, an online tool for the fast and efficient processing of assembled RNA-Seq transcriptome data, developed to mitigate these challenges. TRAPID offers high-throughput open reading frame detection, frameshift correction and includes a functional, comparative and phylogenetic toolbox, making use of 175 reference proteomes. Benchmarking and comparison against state-of-the-art transcript analysis tools reveals the efficiency and unique features of the TRAPID system. TRAPID is freely available at http://bioinformatics.psb.ugent.be/webtools/trapid/.

The sacred lotus genome provides insights into the evolution of flowering plants

Sacred lotus (Nelumbo nucifera) is an ornamental plant that is also used for food and medicine. This basal eudicot species is especially important from an evolutionary perspective, as it occupies a critical phylogenetic position in flowering plants. Here we report the draft genome of a wild strain of sacred lotus. The assembled genome is 792 Mb, which is approximately 85-90% of genome size estimates. We annotated 392 Mb of repeat sequences and 36,385 protein-coding genes within the genome. Using these sequence data, we constructed a phylogenetic tree and confirmed the basal location of sacred lotus within eudicots. Importantly, we found evidence for a relatively recent whole-genome duplication event; any indication of the ancient paleo-hexaploid event was, however, absent. Genomic analysis revealed evidence of positive selection within 28 embryo-defective genes and one annexin gene that may be related to the long-term viability of sacred lotus seed. We also identified a significant expansion of starch synthase genes, which probably elevated starch levels within the rhizome of sacred lotus. Sequencing this strain of sacred lotus thus provided important insights into the evolution of flowering plant and revealed genetic mechanisms that influence seed dormancy and starch synthesis.

Assessing duplication and loss of APETALA1/FRUITFULL homologs in Ranunculales

Gene duplication and loss provide raw material for evolutionary change within organismal lineages as functional diversification of gene copies provide a mechanism for phenotypic variation. Here we focus on the APETALA1/FRUITFULL MADS-box gene lineage evolution. AP1/FUL genes are angiosperm-specific and have undergone several duplications. By far the most significant one is the core-eudicot duplication resulting in the euAP1 and euFUL clades. Functional characterization of several euAP1 and euFUL genes has shown that both function in proper floral meristem identity, and axillary meristem repression. Independently, euAP1 genes function in floral meristem and sepal identity, whereas euFUL genes control phase transition, cauline leaf growth, compound leaf morphogenesis and fruit development. Significant functional variation has been detected in the function of pre-duplication basal-eudicot FUL-like genes, but the underlying mechanisms for change have not been identified. FUL-like genes in the Papaveraceae encode all functions reported for euAP1 and euFUL genes, whereas FUL-like genes in Aquilegia (Ranunculaceae) function in inflorescence development and leaf complexity, but not in flower or fruit development. Here we isolated FUL-like genes across the Ranunculales and used phylogenetic approaches to analyze their evolutionary history. We identified an early duplication resulting in the RanFL1 and RanFL2 clades. RanFL1 genes were present in all the families sampled and are mostly under strong negative selection in the MADS, I and K domains. RanFL2 genes were only identified from Eupteleaceae, Papaveraceae s.l., Menispermaceae and Ranunculaceae and show relaxed purifying selection at the I and K domains. We discuss how asymmetric sequence diversification, new motifs, differences in codon substitutions and likely protein-protein interactions resulting from this Ranunculiid-specific duplication can help explain the functional differences among basal-eudicot FUL-like genes.

Combining phylogenetic and syntenic analyses for understanding the evolution of TCP ECE genes in eudicots

TCP ECE genes encode transcription factors which have received much attention for their repeated recruitment in the control of floral symmetry in core eudicots, and more recently in monocots. Major duplications of TCP ECE genes have been described in core eudicots, but the evolutionary history of this gene family is unknown in basal eudicots. Reconstructing the phylogeny of ECE genes in basal eudicots will help set a framework for understanding the functional evolution of these genes. TCP ECE genes were sequenced in all major lineages of basal eudicots and Gunnera which belongs to the sister clade to all other core eudicots. We show that in these lineages they have a complex evolutionary history with repeated duplications. We estimate the timing of the two major duplications already identified in the core eudicots within a timeframe before the divergence of Gunnera and after the divergence of Proteales. We also use a synteny-based approach to examine the extent to which the expansion of TCP ECE genes in diverse eudicot lineages may be due to genome-wide duplications. The three major core-eudicot specific clades share a number of collinear genes, and their common evolutionary history may have originated at the γ event. Genomic comparisons in Arabidopsis thaliana and Solanumlycopersicum highlight their separate polyploid origin, with syntenic fragments with and without TCP ECE genes showing differential gene loss and genomic rearrangements. Comparison between recently available genomes from two basal eudicots Aquilegiacoerulea and Nelumbonucifera suggests that the two TCP ECE paralogs in these species are also derived from large-scale duplications. TCP ECE loci from basal eudicots share many features with the three main core eudicot loci, and allow us to infer the makeup of the ancestral eudicot locus.

The Tarenaya hassleriana genome provides insight into reproductive trait and genome evolution of crucifers

The Brassicaceae, including Arabidopsis thaliana and Brassica crops, is unmatched among plants in its wealth of genomic and functional molecular data and has long served as a model for understanding gene, genome, and trait evolution. However, genome information from a phylogenetic outgroup that is essential for inferring directionality of evolutionary change has been lacking. We therefore sequenced the genome of the spider flower (Tarenaya hassleriana) from the Brassicaceae sister family, the Cleomaceae. By comparative analysis of the two lineages, we show that genome evolution following ancient polyploidy and gene duplication events affect reproductively important traits. We found an ancient genome triplication in Tarenaya (Th-α) that is independent of the Brassicaceae-specific duplication (At-α) and nested Brassica (Br-α) triplication. To showcase the potential of sister lineage genome analysis, we investigated the state of floral developmental genes and show Brassica retains twice as many floral MADS (for minichromosome maintenance1, AGAMOUS, DEFICIENS and serum response factor) genes as Tarenaya that likely contribute to morphological diversity in Brassica. We also performed synteny analysis of gene families that confer self-incompatibility in Brassicaceae and found that the critical serine receptor kinase receptor gene is derived from a lineage-specific tandem duplication. The T. hassleriana genome will facilitate future research toward elucidating the evolutionary history of Brassicaceae genomes.

Timing and consequences of recurrent polyploidy in meadow-rues (thalictrum, ranunculaceae)

The discovery of ancient whole-genome duplications in eukaryotic lineages has renewed the interest in polyploidy and its effects on the diversification of organisms. Polyploidy has large-scale effects on both genotype and phenotype and has been linked to the evolution of genome size, dioecy, and changes in ecological interactions, such as pollinator visitation. Here, we take a molecular systematics approach to examine the evolution of polyploidy in the plant genus Thalictrum (Ranunculaceae) and test its correlation to changes in genome size, sexual system, and pollination mode. Thalictrum is an ideal study system due to its extensive ploidy range and floral diversity. Phylogenetic analyses were used for character reconstructions, correlation tests, and dating estimates. Our results suggest that polyploidization occurred frequently and recently in the evolution of Thalictrum, mostly within the last 10.6-5.8 My, coinciding with the diversification of particular clades. In spite of an overall trend of genomic downsizing accompanying polyploidy in angiosperms and proportional increases observed at finer scales, our genome size estimates for Thalictrum show no correlation with chromosome number. Instead, we observe genomic expansion in diploids and genomic contraction in polyploids with increased age. Additionally, polyploidy is not correlated with dioecy in Thalictrum; therefore, other factors must have influenced the evolution of separate sexes in this group. A novel finding from our study is the association of polyploidy with shifts to wind pollination, in particular, during a time period of global cooling and mountain uplift in the Americas.

Genome of the long-living sacred lotus (Nelumbo nucifera Gaertn.)

Background

Sacred lotus is a basal eudicot with agricultural, medicinal, cultural and religious importance. It was domesticated in Asia about 7,000 years ago, and cultivated for its rhizomes and seeds as a food crop. It is particularly noted for its 1,300-year seed longevity and exceptional water repellency, known as the lotus effect. The latter property is due to the nanoscopic closely packed protuberances of its self-cleaning leaf surface, which have been adapted for the manufacture of a self-cleaning industrial paint, Lotusan.

Results

The genome of the China Antique variety of the sacred lotus was sequenced with Illumina and 454 technologies, at respective depths of 101× and 5.2×. The final assembly has a contig N50 of 38.8 kbp and a scaffold N50 of 3.4 Mbp, and covers 86.5% of the estimated 929 Mbp total genome size. The genome notably lacks the paleo-triplication observed in other eudicots, but reveals a lineage-specific duplication. The genome has evidence of slow evolution, with a 30% slower nucleotide mutation rate than observed in grape. Comparisons of the available sequenced genomes suggest a minimum gene set for vascular plants of 4,223 genes. Strikingly, the sacred lotus has 16 COG2132 multi-copper oxidase family proteins with root-specific expression; these are involved in root meristem phosphate starvation, reflecting adaptation to limited nutrient availability in an aquatic environment.

Conclusions

The slow nucleotide substitution rate makes the sacred lotus a better resource than the current standard, grape, for reconstructing the pan-eudicot genome, and should therefore accelerate comparative analysis between eudicots and monocots.

Genome-wide expression analysis of soybean MADS genes showing potential function in the seed development

The MADS family is an ancient and best-studied transcription factor and plays fundamental roles in almost every developmental process in plants. In the plant evolutionary history, the whole genome duplication (WGD) events are important not only to the plant species evolution, but to expansion of members of the gene families. Soybean as a model legume crop has experience three rounds of WGD events. Members of some MIKC(C) subfamilies, such as SOC, AGL6, SQUA, SVP, AGL17 and DEF/GLO, were expanded after soybean three rounds of WGD events. And some MIKC(C) subfamilies, MIKC* and type I MADS families had experienced faster birth-and-death evolution and their traces before the Glycine WGD event were not found. Transposed duplication played important roles in tandem arrangements among the members of different subfamilies. According to the expression profiles of type I and MIKC paralog pair genes, the fates of MIKC paralog gene pairs were subfunctionalization, and the fates of type I MADS paralog gene pairs were nonfunctionalization. 137 out of 163 MADS genes were close to 186 loci within 2 Mb genomic regions associated with seed-relative QTLs, among which 115 genes expressed during the seed development. Although MIKC(C) genes kept the important and conserved functions of the flower development, most MIKC(C) genes showed potentially essential roles in the seed development as well as the type I MADS.

Convergent gene loss following gene and genome duplications creates single-copy families in flowering plants

The importance of gene gain through duplication has long been appreciated. In contrast, the importance of gene loss has only recently attracted attention. Indeed, studies in organisms ranging from plants to worms and humans suggest that duplication of some genes might be better tolerated than that of others. Here we have undertaken a large-scale study to investigate the existence of duplication-resistant genes in the sequenced genomes of 20 flowering plants. We demonstrate that there is a large set of genes that is convergently restored to single-copy status following multiple genome-wide and smaller scale duplication events. We rule out the possibility that such a pattern could be explained by random gene loss only and therefore propose that there is selection pressure to preserve such genes as singletons. This is further substantiated by the observation that angiosperm single-copy genes do not comprise a random fraction of the genome, but instead are often involved in essential housekeeping functions that are highly conserved across all eukaryotes. Furthermore, single-copy genes are generally expressed more highly and in more tissues than non-single-copy genes, and they exhibit higher sequence conservation. Finally, we propose different hypotheses to explain their resistance against duplication.