Gene-balanced duplications, like tetraploidy, provide predictable drive to increase morphological complexity

Gene balance hypothesis: connecting issues of dosage sensitivity across biological disciplines

We summarize, in this review, the evidence that genomic balance influences gene expression, quantitative traits, dosage compensation, aneuploid syndromes, population dynamics of copy number variants and differential evolutionary fate of genes after partial or whole-genome duplication. Gene balance effects are hypothesized to result from stoichiometric differences among members of macromolecular complexes, the interactome, and signaling pathways. The implications of gene balance are discussed.

Deciphering the complex leaf transcriptome of the allotetraploid species Nicotiana tabacum: a phylogenomic perspective

BACKGROUND

Polyploidization is an important mechanism in plant evolution. By analyzing the leaf transcriptomes taken from the allotetraploid Nicotiana tabacum (tobacco) and parental genome donors, N. sylvesteris (S-Genome) and N. tomentosiformis (T-Genome), a phylogenomic approach was taken to map the fate of homeologous gene pairs in this plant.

RESULTS

A comparison between the genes present in the leaf transcriptomes of N. tabacum and modern day representatives of its progenitor species demonstrated that only 33% of assembled transcripts could be distinguished based on their sequences. A large majority of the genes (83.6% of the non parent distinguishable and 87.2% of the phylogenetic topology analyzed clusters) expressed above background level (more than 5 reads) showed similar overall expression levels. Homeologous sequences could be identified for 968 gene clusters, and 90% (6% of all genes) of the set maintained expression of only one of the tobacco homeologs. When both homeologs were expressed, only 15% (0.5% of the total) showed evidence of differential expression, providing limited evidence of subfunctionalization. Comparing the rate of synonymous nucleotide substitution (Ks) and non-synonymous nucleotide substitution (Kn) provided limited evidence for positive selection during the evolution of tobacco since the polyploidization event took place.

CONCLUSIONS

Polyploidization is a powerful mechanism for plant speciation that can occur during one generation; however millions of generations may be necessary for duplicate genes to acquire a new function. Analysis of the tobacco leaf transcriptome reveals that polyploidization, even in a young tetraploid such as tobacco, can lead to complex changes in gene expression. Gene loss and gene silencing, or subfunctionalization may explain why both homeologs are not expressed by the associated genes. With Whole Genome Duplication (WGD) events, polyploid genomes usually maintain a high percentage of gene duplicates. The data provided little evidence of preferential maintenance of gene expression from either the T- or S-genome. Additionally there was little evidence of neofunctionalization in Nicotiana tabacum suggesting it occurs at a low frequency in young polyploidy.

Developmental, cytological and transcriptional analysis of autotetraploid Arabidopsis

An autopolyploid that contains more than two sets of the same chromosomes causes apparent alterations in morphology, development, physiology and gene expression compared to diploid. However, the mechanisms for these changes remain largely unknown. In the present study, cytological observations of mature embryos and growing cotyledons demonstrated that enlarged organ size of an autotetraploid Arabidopsis was caused by cell size and not by cell number. Quantitative real time PCR (qRT-PCR) analysis of 34 core cell cycle genes revealed a subtle but stable increase in the expression of ICK1, ICK2 and ICK5 in autotetraploid seedlings. Autotetraploid Arabidopsis plants were found to be more sensitive to glucose treatment than diploid with decreased number of rosette leaves and suppressed root elongation. Cytological observations demonstrated that both cell proliferation and cell expansion of autotetraploid were dramatically suppressed under glucose treatment. Expression levels of ICK1, ICK5 together with Cyclin D and Cyclin B was increased under glucose treatment in both diploid and autotetraploid plants. These results suggest that ICK1 and ICK5 may be involved in developmental delay and that the suppressed growth under glucose treatment probably resulted from disturbed mitotic and endoreduplication cycle in autotetraploid Arabidopsis.

Polyploidy and the evolution of complex traits

We explore how whole-genome duplications (WGDs) may have given rise to complex innovations in cellular networks, innovations that could not have evolved through sequential single-gene duplications. We focus on two classical WGD events, one in bakers' yeast and the other at the base of vertebrates (i.e., two rounds of whole-genome duplication: 2R-WGD). Two complex adaptations are discussed in detail: aerobic ethanol fermentation in yeast and the rewiring of the vertebrate developmental regulatory network through the 2R-WGD. These two examples, derived from diverged branches on the eukaryotic tree, boldly underline the evolutionary potential of WGD in facilitating major evolutionary transitions. We close by arguing that the evolutionary importance of WGD may require updating certain aspects of modern evolutionary theory, perhaps helping to synthesize a new evolutionary systems biology.

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

Comparative genome biology has unveiled the polyploid origin of all angiosperms and the role of recurrent polyploidization in the amplification of gene families and the structuring of genomes. Which species share certain ancient polyploidy events, and which do not, is ill defined because of the limited number of sequenced genomes and transcriptomes and their uneven phylogenetic distribution. Previously, it has been suggested that most, but probably not all, of the eudicots have shared an ancient hexaploidy event, referred to as the gamma triplication. In this study, detailed phylogenies of subfamilies of MADS-box genes suggest that the gamma triplication has occurred before the divergence of Gunnerales but after the divergence of Buxales and Trochodendrales. Large-scale phylogenetic and K(S)-based approaches on the inflorescence transcriptomes of Gunnera manicata (Gunnerales) and Pachysandra terminalis (Buxales) provide further support for this placement, enabling us to position the gamma triplication in the stem lineage of the core eudicots. This triplication likely initiated the functional diversification of key regulators of reproductive development in the core eudicots, comprising 75% of flowering plants. Although it is possible that the gamma event triggered early core eudicot diversification, our dating estimates suggest that the event occurred early in the stem lineage, well before the rapid speciation of the earliest core eudicot lineages. The evolutionary significance of this paleopolyploidy event may thus rather lie in establishing a species lineage that was resilient to extinction, but with the genomic potential for later diversification. We consider that the traits generated from this potential characterize extant core eudicots both chemically and morphologically.

Arabidopsis thaliana AUCSIA-1 regulates auxin biology and physically interacts with a kinesin-related protein

Aucsia is a green plant gene family encoding 44–54 amino acids long miniproteins. The sequenced genomes of most land plants contain two Aucsia genes. RNA interference of both tomato (Solanum lycopersicum) Aucsia genes (SlAucsia-1 and SlAucsia-2) altered auxin sensitivity, auxin transport and distribution; it caused parthenocarpic development of the fruit and other auxin-related morphological changes. Here we present data showing that the Aucsia-1 gene of Arabidopsis thaliana alters, by itself, root auxin biology and that the AtAUCSIA-1 miniprotein physically interacts with a kinesin-related protein. The AtAucsia-1 gene is ubiquitously expressed, although its expression is higher in roots and inflorescences in comparison to stems and leaves. Two allelic mutants for AtAucsia-1 gene did not display visible root morphological alterations; however both basipetal and acropetal indole-3-acetic acid (IAA) root transport was reduced as compared with wild-type plants. The transcript steady state levels of the auxin efflux transporters ATP BINDING CASSETTE subfamily B (ABCB) ABCB1, ABCB4 and ABCB19 were reduced in ataucsia-1 plants. In ataucsia-1 mutant, lateral root growth showed an altered response to i) exogenous auxin, ii) an inhibitor of polar auxin transport and iii) ethylene. Overexpression of AtAucsia-1 inhibited primary root growth. In vitro and in vivo protein-protein interaction experiments showed that AtAUCSIA-1 interacts with a 185 amino acids long fragment belonging to a 2712 amino acids long protein of unknown function (At4g31570). Bioinformatics analysis indicates that the AtAUCSIA-1 interacting protein (AtAUCSIA-1IP) clusters with a group of CENP-E kinesin-related proteins. Gene ontology predictions for the two proteins are consistent with the hypothesis that the AtAUCSIA-1/AtAUCSIA-1IP complex is involved in the regulation of the cytoskeleton dynamics underlying auxin biology.

Sequence divergence in the 3'-untranslated region has an effect on the subfunctionalization of duplicate genes

Recent studies demonstrated that sequence divergence in both transcriptional regulatory region and coding region contributes to the subfunctionalization of duplicate gene. However, whether sequence divergence in the 3'-untranslated region (3'-UTR) has an impact on the subfunctionalization of duplicate genes remains unclear. Here, we identified two diverging duplicate vsx1 (visual system homeobox-1) loci in goldfish, named vsx1A1 and vsx1A2. Phylogenetic analysis suggests that vsx1A1 and vsx1A2 may arise from a duplication of vsx1 after the separation of goldfish and zebrafish. Sequence comparison revealed that divergence in both transcriptional and translational regulatory regions is higher than divergence in the introns. vsx1A2 expresses during blastula and gastrula stages and in adult retina but silences from segmentation stage to hatching stage, vsx1A1 starts expression from segmentation onward. Comparing to that zebrafish vsx1 expresses in all the developmental stages and in the adult retina, it appears that goldfish vsx1A1 and vsx1A2 are under going to share the functions of ancestral vsx1. The different but overlapping temporal expression patterns of vsx1A1 and vsx1A2 suggest that sequence divergence in the promoter region of duplicate vsx1 is not sufficient for partitioning the functions of ancestral vsx1. By comparing vsx1A1 and vsx1A2 3'-UTR-linked green fluorescent protein gene expression patterns, we demonstrated that the 3'-UTR of vsx1A1 remains but the 3'-UTR of vsx1A2 has lost the capability of mediating bipolar cell specific expression during retina development. These results indicate that sequence divergence in the 3'-UTRs has a clear effect on subfunctionalization of the duplicate genes.

Retention of triplicated phytoene synthase (PSY) genes in Brassica napus L. and its diploid progenitors during the evolution of the Brassiceae

The extent of genome redundancy exhibited by Brassica species provides a model to study the evolutionary fate of multi-copy genes and the effects of polyploidy in economically important crops. Phytoene synthase (PSY) catalyzes the first committed reaction of the carotenoid biosynthetic pathway, which has been shown to be rate-limiting in Brassica napus seeds. In Arabidopsis thaliana, a single PSY gene (AtPSY) regulates phytoene synthesis in all tissues. Considering that diploid Brassica genomes contain three Arabidopsis-like subgenomes, the objectives of the present work were to determine whether PSY gene families exist in B. napus (AACC) and its diploid progenitor species, Brassica rapa (AA) and Brassica oleracea (CC); to establish the level of retention of Brassica PSY genes; to map PSY gene family members in the A and C genomes and to compare Brassica PSY gene expression patterns. A total of 12 PSY homologues were identified, 6 in B. napus (BnaX.PSY.a-f) and 3 in B. rapa (BraA.PSY.a-c) and B. oleracea (BolC.PSY.a-c). Indeed, with six members, B. napus has the largest PSY gene family described to date. Sequence comparison between AtPSY and Brassica PSY genes revealed a highly conserved gene structure and identity percentages above 85% at the coding sequence (CDS) level. Altogether, our data indicate that PSY gene family expansion preceded the speciation of B. rapa and B. oleracea, dating back to the paralogous subgenome triplication event. In these three Brassica species, all PSY homologues are expressed, exhibiting overlapping redundancy and signs of subfunctionalization among photosynthetic and non-photosynthetic tissues. This evidence supports the hypothesis that functional divergence of PSY gene expression facilitates the accumulation of high levels of carotenoids in chromoplast-rich tissues. Thus, functional retention of triplicated Brassica PSY genes could be at least partially explained by the selective advantage provided by increased levels of gene product in floral organs. A better understanding of carotenogenesis in Brassica will aid in the future development of transgenic and conventional cultivars with carotenoid-enriched oil.

Redundancy and rewiring of genetic networks following genome-wide duplication events

Polyploidy or whole-genome duplication is a frequent phenomenon within the plant kingdom and has been associated with the occurrence of evolutionary novelty and increase in biological complexity. Because genome-wide duplication events duplicate whole molecular networks it is of interest to investigate how these networks evolve subsequent to such events. Although genome duplications are generally followed by massive gene loss, at least part of the network is usually retained in duplicate and can rewire to execute novel functions. Alternatively, the network can remain largely redundant and as such confer robustness against mutations. The increasing availability of high-throughput data makes it possible to study evolution following whole genome duplication events at the network level. Here we discuss how the use of 'omics' data in network analysis can provide novel insights on network redundancy and rewiring and conclude with some directions for future research.

Mammalian X chromosome inactivation evolved as a dosage-compensation mechanism for dosage-sensitive genes on the X chromosome

How and why female somatic X-chromosome inactivation (XCI) evolved in mammals remains poorly understood. It has been proposed that XCI is a dosage-compensation mechanism that evolved to equalize expression levels of X-linked genes in females (2X) and males (1X), with a prior twofold increase in expression of X-linked genes in both sexes ("Ohno's hypothesis"). Whereas the parity of X chromosome expression between the sexes has been clearly demonstrated, tests for the doubling of expression levels globally along the X chromosome have returned contradictory results. However, changes in gene dosage during sex-chromosome evolution are not expected to impact on all genes equally, and should have greater consequences for dosage-sensitive genes. We show that, for genes encoding components of large protein complexes (≥ 7 members)--a class of genes that is expected to be dosage-sensitive--expression of X-linked genes is similar to that of autosomal genes within the complex. These data support Ohno's hypothesis that XCI acts as a dosage-compensation mechanism, and allow us to refine Ohno's model of XCI evolution. We also explore the contribution of dosage-sensitive genes to X aneuploidy phenotypes in humans, such as Turner (X0) and Klinefelter (XXY) syndromes. X aneuploidy in humans is common and is known to have mild effects because most of the supernumerary X genes are inactivated and not affected by aneuploidy. Only genes escaping XCI experience dosage changes in X-aneuploidy patients. We combined data on dosage sensitivity and XCI to compute a list of candidate genes for X-aneuploidy syndromes.

Genome-wide analysis of syntenic gene deletion in the grasses

The grasses, Poaceae, are one of the largest and most successful angiosperm families. Like many radiations of flowering plants, the divergence of the major grass lineages was preceded by a whole-genome duplication (WGD), although these events are not rare for flowering plants. By combining identification of syntenic gene blocks with measures of gene pair divergence and different frequencies of ancient gene loss, we have separated the two subgenomes present in modern grasses. Reciprocal loss of duplicated genes or genomic regions has been hypothesized to reproductively isolate populations and, thus, speciation. However, in contrast to previous studies in yeast and teleost fishes, we found very little evidence of reciprocal loss of homeologous genes between the grasses, suggesting that post-WGD gene loss may not be the cause of the grass radiation. The sets of homeologous and orthologous genes and predicted locations of deleted genes identified in this study, as well as links to the CoGe comparative genomics web platform for analyzing pan-grass syntenic regions, are provided along with this paper as a resource for the grass genetics community.

Genome and gene duplications and gene expression divergence: a view from plants

With many plant genomes sequenced, it is now clear that one distinguishing feature of angiosperm (flowering plant) genomes is their high frequency of whole-genome duplication. Single-gene duplication is also widespread in angiosperm genomes. Following various mechanisms of gene duplication, genes are often retained or lost in a biased manner, which has suggested recent models for gene family evolution, such as functional buffering and the gene balance hypothesis in addition to now-classical models, including neofunctionalization and subfunctionalization. Evolutionary consequences of gene duplication, often studied through analyzing gene expression divergence, have enhanced understanding of the biological significance of different mechanisms of gene duplication.

Dynamics and adaptive benefits of protein domain emergence and arrangements during plant genome evolution

Plant genomes are generally very large, mostly paleopolyploid, and have numerous gene duplicates and complex genomic features such as repeats and transposable elements. Many of these features have been hypothesized to enable plants, which cannot easily escape environmental challenges, to rapidly adapt. Another mechanism, which has recently been well described as a major facilitator of rapid adaptation in bacteria, animals, and fungi but not yet for plants, is modular rearrangement of protein-coding genes. Due to the high precision of profile-based methods, rearrangements can be well captured at the protein level by characterizing the emergence, loss, and rearrangements of protein domains, their structural, functional, and evolutionary building blocks. Here, we study the dynamics of domain rearrangements and explore their adaptive benefit in 27 plant and 3 algal genomes. We use a phylogenomic approach by which we can explain the formation of 88% of all arrangements by single-step events, such as fusion, fission, and terminal loss of domains. We find many domains are lost along every lineage, but at least 500 domains are novel, that is, they are unique to green plants and emerged more or less recently. These novel domains duplicate and rearrange more readily within their genomes than ancient domains and are overproportionally involved in stress response and developmental innovations. Novel domains more often affect regulatory proteins and show a higher degree of structural disorder than ancient domains. Whereas a relatively large and well-conserved core set of single-domain proteins exists, long multi-domain arrangements tend to be species-specific. We find that duplicated genes are more often involved in rearrangements. Although fission events typically impact metabolic proteins, fusion events often create new signaling proteins essential for environmental sensing. Taken together, the high volatility of single domains and complex arrangements in plant genomes demonstrate the importance of modularity for environmental adaptability of plants.

Evidence of function for conserved noncoding sequences in Arabidopsis thaliana

• Whole genome duplication events provide a lineage with a large reservoir of genes that can be molded by evolutionary forces into phenotypes that fit alternative environments. A well-studied whole genome duplication, the α-event, occurred in an ancestor of the model plant Arabidopsis thaliana. Retained segments of the α-event have been defined in recent years in the form of duplicate protein coding sequences (α-pairs) and associated conserved noncoding DNA sequences (CNSs). Our aim was to identify any association between CNSs and α-pair co-functionality at the gene expression level. • Here, we tested for correlation between CNS counts and α-pair co-expression and expression intensity across nine expression datasets: aerial tissue, flowers, leaves, roots, rosettes, seedlings, seeds, shoots and whole plants. • We provide evidence for a putative regulatory role of the CNSs. The association of CNSs with α-pair co-expression and expression intensity varied by gene function, subgene position and the presence of transcription factor binding motifs. A range of possible CNS regulatory mechanisms, including intron-mediated enhancement, messenger RNA fold stability and transcriptional regulation, are discussed. • This study provides a framework to understand how CNS motifs are involved in the maintenance of gene expression after a whole genome duplication event.

The Impact of Genome Triplication on Tandem Gene Evolution in Brassica rapa

Whole genome duplication (WGD) and tandem duplication (TD) are both important modes of gene expansion. However, how WGD influences tandemly duplicated genes is not well studied. We used Brassica rapa, which has undergone an additional genome triplication (WGT) and shares a common ancestor with Arabidopsis thaliana, Arabidopsis lyrata, and Thellungiella parvula, to investigate the impact of genome triplication on tandem gene evolution. We identified 2,137, 1,569, 1,751, and 1,135 tandem gene arrays in B. rapa, A. thaliana, A. lyrata, and T. parvula respectively. Among them, 414 conserved tandem arrays are shared by the three species without WGT, which were also considered as existing in the diploid ancestor of B. rapa. Thus, after genome triplication, B. rapa should have 1,242 tandem arrays according to the 414 conserved tandems. Here, we found 400 out of the 414 tandems had at least one syntenic ortholog in the genome of B. rapa. Furthermore, 294 out of the 400 shared syntenic orthologs maintain tandem arrays (more than one gene for each syntenic hit) in B. rapa. For the 294 tandem arrays, we obtained 426 copies of syntenic paralogous tandems in the triplicated genome of B. rapa. In this study, we demonstrated that tandem arrays in B. rapa were dramatically fractionated after WGT when compared either to non-tandem genes in the B. rapa genome or to the tandem arrays in closely related species that have not experienced a recent whole genome polyploidization event.

The odds of duplicate gene persistence after polyploidization

Background

Gene duplication is an important biological phenomenon associated with genomic redundancy, degeneration, specialization, innovation, and speciation. After duplication, both copies continue functioning when natural selection favors duplicated protein function or expression, or when mutations make them functionally distinct before one copy is silenced.

Results

Here we quantify the degree to which genetic parameters related to gene expression, molecular evolution, and gene structure in a diploid frog - Silurana tropicalis - influence the odds of functional persistence of orthologous duplicate genes in a closely related tetraploid species - Xenopus laevis. Using public databases and 454 pyrosequencing, we obtained genetic and expression data from S. tropicalis orthologs of 3,387 X. laevis paralogs and 4,746 X. laevis singletons - the most comprehensive dataset for African clawed frogs yet analyzed. Using logistic regression, we demonstrate that the most important predictors of the odds of duplicate gene persistence in the tetraploid species are the total gene expression level and evenness of expression across tissues and development in the diploid species. Slow protein evolution and information density (fewer exons, shorter introns) in the diploid are also positively correlated with duplicate gene persistence in the tetraploid.

Conclusions

Our findings suggest that a combination of factors contribute to duplicate gene persistence following whole genome duplication, but that the total expression level and evenness of expression across tissues and through development before duplication are most important. We speculate that these parameters are useful predictors of duplicate gene longevity after whole genome duplication in other taxa.

The Medicago genome provides insight into the evolution of rhizobial symbioses

Legumes (Fabaceae or Leguminosae) are unique among cultivated plants for their ability to carry out endosymbiotic nitrogen fixation with rhizobial bacteria, a process that takes place in a specialized structure known as the nodule. Legumes belong to one of the two main groups of eurosids, the Fabidae, which includes most species capable of endosymbiotic nitrogen fixation. Legumes comprise several evolutionary lineages derived from a common ancestor 60 million years ago (Myr ago). Papilionoids are the largest clade, dating nearly to the origin of legumes and containing most cultivated species. Medicago truncatula is a long-established model for the study of legume biology. Here we describe the draft sequence of the M. truncatula euchromatin based on a recently completed BAC assembly supplemented with Illumina shotgun sequence, together capturing ∼94% of all M. truncatula genes. A whole-genome duplication (WGD) approximately 58 Myr ago had a major role in shaping the M. truncatula genome and thereby contributed to the evolution of endosymbiotic nitrogen fixation. Subsequent to the WGD, the M. truncatula genome experienced higher levels of rearrangement than two other sequenced legumes, Glycine max and Lotus japonicus. M. truncatula is a close relative of alfalfa (Medicago sativa), a widely cultivated crop with limited genomics tools and complex autotetraploid genetics. As such, the M. truncatula genome sequence provides significant opportunities to expand alfalfa's genomic toolbox.

Selection for higher gene copy number after different types of plant gene duplications

The evolutionary origins of the multitude of duplicate genes in the plant genomes are still incompletely understood. To gain an appreciation of the potential selective forces acting on these duplicates, we phylogenetically inferred the set of metabolic gene families from 10 flowering plant (angiosperm) genomes. We then compared the metabolic fluxes for these families, predicted using the Arabidopsis thaliana and Sorghum bicolor metabolic networks, with the families' duplication propensities. For duplications produced by both small scale (small-scale duplications) and genome duplication (whole-genome duplications), there is a significant association between the flux and the tendency to duplicate. Following this global analysis, we made a more fine-scale study of the selective constraints observed on plant sodium and phosphate transporters. We find that the different duplication mechanisms give rise to differing selective constraints. However, the exact nature of this pattern varies between the gene families, and we argue that the duplication mechanism alone does not define a duplicated gene's subsequent evolutionary trajectory. Collectively, our results argue for the interplay of history, function, and selection in shaping the duplicate gene evolution in plants.

Phylogenetic perspectives on the origins of nodulation

Recent refinements to the phylogeny of rosid angiosperms support the conclusion that nodulation has evolved several times in the so-called N(2)-fixing clade (NFC), and provide dates for these origins. The hypothesized predisposition that enabled the evolution of nodulation occurred approximately 100 million years ago (MYA), was retained in the various lineages that radiated rapidly shortly thereafter, and was functional in its non-nodulation role for at least an additional 30 million years in each nodulating lineage. Legumes radiated rapidly shortly after their origin approximately 60 MYA, and nodulation most likely evolved several times during this radiation. The major lineages of papilionoid legumes diverged close to the time of origin of nodulation, accounting for the diversity of nodule biology in the group. Nodulation symbioses exemplify the concept of "deep homology," sharing various homologous components across nonhomologous origins of nodulation, largely due to recruitment from existing functions, notably the older arbuscular mycorrhizal symbiosis. Although polyploidy may have played a role in the origin of papilionoid legume nodules, it did not do so in other legumes, nor did the prerosid whole-genome triplication lead directly to the predisposition of nodulation.

Contrasting patterns of evolution following whole genome versus tandem duplication events in Populus

Comparative analysis of multiple angiosperm genomes has implicated gene duplication in the expansion and diversification of many gene families. However, empirical data and theory suggest that whole-genome and small-scale duplication events differ with respect to the types of genes preserved as duplicate pairs. We compared gene duplicates resulting from a recent whole genome duplication to a set of tandemly duplicated genes in the model forest tree Populus trichocarpa. We used a combination of microarray expression analyses of a diverse set of tissues and functional annotation to assess factors related to the preservation of duplicate genes of both types. Whole genome duplicates are 700 bp longer and are expressed in 20% more tissues than tandem duplicates. Furthermore, certain functional categories are over-represented in each class of duplicates. In particular, disease resistance genes and receptor-like kinases commonly occur in tandem but are significantly under-retained following whole genome duplication, while whole genome duplicate pairs are enriched for members of signal transduction cascades and transcription factors. The shape of the distribution of expression divergence for duplicated pairs suggests that nearly half of the whole genome duplicates have diverged in expression by a random degeneration process. The remaining pairs have more conserved gene expression than expected by chance, consistent with a role for selection under the constraints of gene balance. We hypothesize that duplicate gene preservation in Populus is driven by a combination of subfunctionalization of duplicate pairs and purifying selection favoring retention of genes encoding proteins with large numbers of interactions.

DNA methylation rebalances gene dosage after mammalian gene duplications

Although gene duplication plays a major role in organismal evolution, it may also lead to gene dosage imbalance, thereby having an immediate adverse effect on an organism's fitness. Investigating the evolution of the expression patterns of genes that duplicated after the divergence of rodents and primates, we confirm that adaptive evolution has been involved in dosage rebalance after gene duplication. To understand mechanisms underlying this process, we examined 1) microRNA (miRNA)-mediated gene regulation, 2) cis-regulatory sequence modifications, and 3) DNA methylation. Neither miRNA-mediated regulation nor cis-regulatory changes was found to be associated with expression reduction of duplicate genes. By contrast, duplicate genes, especially lowly expressed copies, were heavily methylated in the upstream region. However, for duplicate genes encoding proteins that are members of macromolecular complexes, heavy methylation in the genic region was not consistently observed. This result held after controlling potential confounding factors, such as enrichment in functional categories. Our results suggest that during mammalian evolution, DNA methylation plays a dominant role in dosage rebalance after gene duplication by inhibiting transcription initiation of duplicate genes.

Evolutionary divergence and limits of conserved non-coding sequence detection in plant genomes

The discovery of regulatory motifs embedded in upstream regions of plants is a particularly challenging bioinformatics task. Previous studies have shown that motifs in plants are short compared with those found in vertebrates. Furthermore, plant genomes have undergone several diversification mechanisms such as genome duplication events which impact the evolution of regulatory motifs. In this article, a systematic phylogenomic comparison of upstream regions is conducted to further identify features of the plant regulatory genomes, the component of genomes regulating gene expression, to enable future de novo discoveries. The findings highlight differences in upstream region properties between major plant groups and the effects of divergence times and duplication events. First, clear differences in upstream region evolution can be detected between monocots and dicots, thus suggesting that a separation of these groups should be made when searching for novel regulatory motifs, particularly since universal motifs such as the TATA box are rare. Second, investigating the decay rate of significantly aligned regions suggests that a divergence time of ~100 mya sets a limit for reliable conserved non-coding sequence (CNS) detection. Insights presented here will set a framework to help identify embedded motifs of functional relevance by understanding the limits of bioinformatics detection for CNSs.

Gene duplication and loss in a MADS box gene transcription factor circuit

Although many models have been proposed that could lead to the maintenance of gene duplicates, the ways in which interacting gene duplicates influence each other's evolution and function remain poorly understood. Here, we focus on duplication and loss of the B class MADS box transcription factor genes in the euasterids I and the ramifications of such changes on paralog evolution and their encoded functions. In core eudicots, the B class genes belong to two paralogous lineages whose products form obligate heterodimers. Based on comparative genomic and phylogenetic analyses, we show that five stepwise B class MADS box gene gain or loss events occurred during the radiation of the euasterids I within core eudicots. Gene loss in one sublineage was correlated with a deficit of other sublineage genes. We also show that the gain or loss of B class MADS box gene paralogs were associated with altered protein-protein interactions among the remaining copies. These altered protein interactions were correlated with asymmetric patterns of sequence diversification and selection, suggesting that compensatory changes were driving the evolution of such genes. Furthermore, these B class MADS box gene gain or loss events were associated with the evolutionary divergence of floral morphology in the euasterids I. Together, these observations point to a cooperative strategy by which gene networks evolve, with selection maintaining the overall logic of a network despite changes in individual components.

High time for a roll call: gene duplication and phylogenetic relationships of TCP-like genes in monocots

BACKGROUND AND AIMS

The TCP family is an ancient group of plant developmental transcription factors that regulate cell division in vegetative and reproductive structures and are essential in the establishment of flower zygomorphy. In-depth research on eudicot TCPs has documented their evolutionary and developmental role. This has not happened to the same extent in monocots, although zygomorphy has been critical for the diversification of Orchidaceae and Poaceae, the largest families of this group. Investigating the evolution and function of TCP-like genes in a wider group of monocots requires a detailed phylogenetic analysis of all available sequence information and a system that facilitates comparing genetic and functional information.

METHODS

The phylogenetic relationships of TCP-like genes in monocots were investigated by analysing sequences from the genomes of Zea mays, Brachypodium distachyon, Oryza sativa and Sorghum bicolor, as well as EST data from several other monocot species.

KEY RESULTS

All available monocot TCP-like sequences are associated in 20 major groups with an average identity ≥64 % and most correspond to well-supported clades of the phylogeny. Their sequence motifs and relationships of orthology were documented and it was found that 67 % of the TCP-like genes of Sorghum, Oryza, Zea and Brachypodium are in microsyntenic regions. This analysis suggests that two rounds of whole genome duplication drove the expansion of TCP-like genes in these species.

CONCLUSIONS

A system of classification is proposed where putative or recognized monocot TCP-like genes are assigned to a specific clade of PCF-, CIN- or CYC/tb1-like genes. Specific biases in sequence data of this family that must be tackled when studying its molecular evolution and phylogeny are documented. Finally, the significant retention of duplicated TCP genes from Zea mays is considered in the context of balanced gene drive.

Robustness and evolvability in the B-system of flower development

BACKGROUND

Gene duplication has often been invoked as a key mechanism responsible for evolution of new morphologies. The floral homeotic B-group gene family has undergone a number of gene duplication events, and yet the functions of these genes appear to be largely conserved. However, detailed comparative analysis has indicated that such duplicate genes have considerable cryptic variability in their functions. In the Solanaceae, two duplicate B-group gene lineages have been retained in three subfamilies. Comparisons of orthologous genes across members of the Solanaceae have demonstrated that the combined function of all four B-gene members is to establish petal and stamen identity, but that this function was partitioned differently in each species. These observations emphasize both the robustness and the evolvability of the B-system.

SCOPE

We provide an overview of how the B-function genes can robustly specify petal and stamen identity and at the same time evolve through changes in protein-protein interaction, gene expression patterns, copy number variation or alterations in the downstream genes they control. By using mathematical models we explore regulatory differences between species and how these impose constraints on downstream gene regulation.

CONCLUSIONS

Evolvability of the B-genes can be understood through the multiple ways in which the B-system can be robust. Quantitative approaches should allow for the incorporation of more biological realism in the representations of these regulatory systems and this should contribute to understanding the constraints under which different B-systems can function and evolve. This, in turn, can provide a better understanding of the ways in which B-genes have contributed to flower diversity.

A mystery unveiled

A recent phylogenomic study has provided new evidence for two ancient whole genome duplications in plants, with potential importance for the evolution of seed and flowering plants.

An autotetraploid linkage map of rose (Rosa hybrida) validated using the strawberry (Fragaria vesca) genome sequence

Polyploidy is a pivotal process in plant evolution as it increase gene redundancy and morphological intricacy but due to the complexity of polysomic inheritance we have only few genetic maps of autopolyploid organisms. A robust mapping framework is particularly important in polyploid crop species, rose included (2n = 4x = 28), where the objective is to study multiallelic interactions that control traits of value for plant breeding. From a cross between the garden, peach red and fragrant cultivar Fragrant Cloud (FC) and a cut-rose yellow cultivar Golden Gate (GG), we generated an autotetraploid GGFC mapping population consisting of 132 individuals. For the map we used 128 sequence-based markers, 141 AFLP, 86 SSR and three morphological markers. Seven linkage groups were resolved for FC (Total 632 cM) and GG (616 cM) which were validated by markers that segregated in both parents as well as the diploid integrated consensus map.The release of the Fragaria vesca genome, which also belongs to the Rosoideae, allowed us to place 70 rose sequenced markers on the seven strawberry pseudo-chromosomes. Synteny between Rosa and Fragaria was high with an estimated four major translocations and six inversions required to place the 17 non-collinear markers in the same order. Based on a verified linear order of the rose markers, we could further partition each of the parents into its four homologous groups, thus providing an essential framework to aid the sequencing of an autotetraploid genome.

Transcriptional robustness and protein interactions are associated in yeast

Background

Robustness to insults, both external and internal, is a characteristic feature of life. One level of biological organization for which noise and robustness have been extensively studied is gene expression. Cells have a variety of mechanisms for buffering noise in gene expression, but it is not completely clear what rules govern whether or not a given gene uses such tools to maintain appropriate expression.

Results

Here, we show a general association between the degree to which yeast cells have evolved mechanisms to buffer changes in gene expression and whether they possess protein-protein interactions. We argue that this effect bears an affinity to epistasis, because yeast appears to have evolved regulatory mechanisms such that distant changes in gene copy number for a protein-protein interaction partner gene can alter a gene's expression. This association is not unexpected given recent work linking epistasis and the deleterious effects of changes in gene dosage (i.e., the dosage balance hypothesis). Using gene expression data from artificial aneuploid strains of bakers' yeast, we found that genes coding for proteins that physically interact with other proteins show less expression variation in response to aneuploidy than do other genes. This effect is even more pronounced for genes whose products interact with proteins encoded on aneuploid chromosomes. We further found that genes targeted by transcription factors encoded on aneuploid chromosomes were more likely to change in expression after aneuploidy.

Conclusions

We suggest that these observations can be best understood as resulting from the higher fitness cost of misexpression in epistatic genes and a commensurate greater regulatory control of them.

Two-phase resolution of polyploidy in the Arabidopsis metabolic network gives rise to relative and absolute dosage constraints

The abundance of detected ancient polyploids in extant genomes raises questions regarding evolution after whole-genome duplication (WGD). For instance, what rules govern the preservation or loss of the duplicated genes created by WGD? We explore this question by contrasting two possible preservation forces: selection on relative and absolute gene dosages. Constraints on the relative dosages of central network genes represent an important force for maintaining duplicates (the dosage balance hypothesis). However, preservation may also result from selection on the absolute abundance of certain gene products. The metabolic network of the model plant Arabidopsis thaliana is a powerful system for comparing these hypotheses. We analyzed the surviving WGD-produced duplicate genes in this network, finding evidence that the surviving duplicates from the most recent WGD (WGD-α) are clustered in the network, as predicted by the dosage balance hypothesis. A flux balance analysis suggests an association between the survival of duplicates from a more ancient WGD (WGD-β) and reactions with high metabolic flux. We argue for an interplay of relative and absolute dosage constraints, such that the relative constraints imposed by the recent WGD are still being resolved by evolution, while they have been essentially fully resolved for the ancient event.

Comparative evolution of photosynthetic genes in response to polyploid and nonpolyploid duplication

The likelihood of duplicate gene retention following polyploidy varies by functional properties (e.g. gene ontologies or protein family domains), but little is known about the effects of whole-genome duplication on gene networks related by a common physiological process. Here, we examined the effects of both polyploid and nonpolyploid duplications on genes encoding the major functional groups of photosynthesis (photosystem I, photosystem II, the light-harvesting complex, and the Calvin cycle) in the cultivated soybean (Glycine max), which has experienced two rounds of whole-genome duplication. Photosystem gene families exhibit retention patterns consistent with dosage sensitivity (preferential retention of polyploid duplicates and elimination of nonpolyploid duplicates), whereas Calvin cycle and light-harvesting complex gene families do not. We observed similar patterns in barrel medic (Medicago truncatula), which shared the older genome duplication with soybean but has evolved independently for approximately 50 million years, and in Arabidopsis (Arabidopsis thaliana), which experienced two nested polyploidy events independent from the legume duplications. In both soybean and Arabidopsis, Calvin cycle gene duplicates exhibit a greater capacity for functional differentiation than do duplicates within the photosystems, which likely explains the greater retention of ancient, nonpolyploid duplicates and larger average gene family size for the Calvin cycle relative to the photosystems.

Divergent evolutionary fates of major photosynthetic gene networks following gene and whole genome duplications

Gene and genome duplication are recurring processes in flowering plants, and elucidating the mechanisms by which duplicated genes are lost or deployed is a key component of understanding plant evolution. Using gene ontologies (GO) or protein family (PFAM) domains, distinct patterns of duplicate retention and loss have been identified depending on gene functional properties and duplication mechanism, but little is known about how gene networks encoding interacting proteins (protein complexes or signaling cascades) evolve in response to duplication. We examined patterns of duplicate retention within four major gene networks involved in photosynthesis (the Calvin cycle, photosystem I, photosystem II, and the light harvesting complex) across three species and four whole genome duplications, as well as small-scale duplications, and showed that photosystem gene family evolution is governed largely by dosage sensitivity. ( 1) In contrast, Calvin cycle gene families are not dosage sensitive, but exhibit a greater capacity for functional differentiation. Here we review these findings, highlight how this study, by analyzing defined gene networks, is complementary to global studies using functional annotations such as GO and PFAM, and elaborate on one example of functional differentiation in the Calvin cycle gene family, transketolase.

Transcriptomic shock generates evolutionary novelty in a newly formed, natural allopolyploid plant

New hybrid species might be expected to show patterns of gene expression intermediate to those shown by parental species. "Transcriptomic shock" may also occur, in which gene expression is disrupted; this may be further modified by whole genome duplication (causing allopolyploidy). "Shock" can include instantaneous partitioning of gene expression between parental copies of genes among tissues. These effects have not previously been studied at a population level in a natural allopolyploid plant species. Here, we survey tissue-specific expression of 144 duplicated gene pairs derived from different parental species (homeologs) in two natural populations of 40-generation-old allotetraploid Tragopogon miscellus (Asteraceae) plants. We compare these results with patterns of allelic expression in both in vitro "hybrids" and hand-crossed F(1) hybrids between the parental diploids T. dubius and T. pratensis, and with patterns of homeolog expression in synthetic (S(1)) allotetraploids. Partitioning of expression was frequent in natural allopolyploids, but F(1) hybrids and S(1) allopolyploids showed less partitioning of expression than the natural allopolyploids and the in vitro "hybrids" of diploid parents. Our results suggest that regulation of gene expression is relaxed in a concerted manner upon hybridization, and new patterns of partitioned expression subsequently emerge over the generations following allopolyploidization.

Seventy million years of concerted evolution of a homoeologous chromosome pair, in parallel, in major Poaceae lineages

Whole genome duplication ~70 million years ago provided raw material for Poaceae (grass) diversification. Comparison of rice (Oryza sativa), sorghum (Sorghum bicolor), maize (Zea mays), and Brachypodium distachyon genomes revealed that one paleo-duplicated chromosome pair has experienced very different evolution than all the others. For tens of millions of years, the two chromosomes have experienced illegitimate recombination that has been temporally restricted in a stepwise manner, producing structural stratification in the chromosomes. These strata formed independently in different grass lineages, with their similarities (low sequence divergence between paleo-duplicated genes) preserved in parallel for millions of years since the divergence of these lineages. The pericentromeric region of this homeologous chromosome pair accounts for two-thirds of the gene content differences between the modern chromosomes. Both intriguing and perplexing is a distal chromosomal region with the greatest DNA similarity between surviving duplicated genes but also with the highest concentration of lineage-specific gene pairs found anywhere in these genomes and with a significantly elevated gene evolutionary rate. Intragenomic similarity near this chromosomal terminus may be important in hom(e)ologous chromosome pairing. Chromosome structural stratification, together with enrichment of autoimmune response-related (nucleotide binding site-leucine-rich repeat) genes and accelerated DNA rearrangement and gene loss, confer a striking resemblance of this grass chromosome pair to the sex chromosomes of other taxa.

2R and remodeling of vertebrate signal transduction engine

Background

Whole genome duplication (WGD) is a special case of gene duplication, observed rarely in animals, whereby all genes duplicate simultaneously through polyploidisation. Two rounds of WGD (2R-WGD) occurred at the base of vertebrates, giving rise to an enormous wave of genetic novelty, but a systematic analysis of functional consequences of this event has not yet been performed.

Results

We show that 2R-WGD affected an overwhelming majority (74%) of signalling genes, in particular developmental pathways involving receptor tyrosine kinases, Wnt and transforming growth factor-β ligands, G protein-coupled receptors and the apoptosis pathway. 2R-retained genes, in contrast to tandem duplicates, were enriched in protein interaction domains and multifunctional signalling modules of Ras and mitogen-activated protein kinase cascades. 2R-WGD had a fundamental impact on the cell-cycle machinery, redefined molecular building blocks of the neuronal synapse, and was formative for vertebrate brains. We investigated 2R-associated nodes in the context of the human signalling network, as well as in an inferred ancestral pre-2R (AP2R) network, and found that hubs (particularly involving negative regulation) were preferentially retained, with high connectivity driving retention. Finally, microarrays and proteomics demonstrated a trend for gradual paralog expression divergence independent of the duplication mechanism, but inferred ancestral expression states suggested preferential subfunctionalisation among 2R-ohnologs (2ROs).

Conclusions

The 2R event left an indelible imprint on vertebrate signalling and the cell cycle. We show that 2R-WGD preferentially retained genes are associated with higher organismal complexity (for example, locomotion, nervous system, morphogenesis), while genes associated with basic cellular functions (for example, translation, replication, splicing, recombination; with the notable exception of cell cycle) tended to be excluded. 2R-WGD set the stage for the emergence of key vertebrate functional novelties (such as complex brains, circulatory system, heart, bone, cartilage, musculature and adipose tissue). A full explanation of the impact of 2R on evolution, function and the flow of information in vertebrate signalling networks is likely to have practical consequences for regenerative medicine, stem cell therapies and cancer treatment.

How the vertebrates were made: selective pruning of a double-duplicated genome

Vertebrates are the result of an ancient double duplication of the genome. A new study published in BMC Biology explores the selective retention of genes after this event, finding an extensive enrichment of signaling proteins and transcription factors. Analysis of their expression patterns, interactions and subsequent history reflect the forces that drove their evolution, and with it the evolution of vertebrate complexity.

See research article: http://www.biomedcentral.com/1741-7007/8/146/abstract

Copy number alterations among mammalian enzymes cluster in the metabolic network

Using two high-quality human metabolic networks, we employed comparative genomics techniques to infer metabolic network structures for seven other mammals. We then studied copy number alterations (CNAs) in these networks. Using a graph-theoretic approach, we show that the pattern of CNAs is distinctly different from the random distributions expected under genetic drift. Instead, we find that changes in copy number are most common among transporter genes and that the CNAs differ depending on the mammalian lineage in question. Thus, we find an excess of transporter genes in cattle involved in the milk production, secretion, and regulation. These results suggest a potential role for dosage selection in the evolution of mammalian metabolic networks.

Nonrandom survival of gene conversions among yeast ribosomal proteins duplicated through genome doubling

By comparing the patterns of evolution in the coding and upstream noncoding regions of yeast ribosomal protein (RP) genes duplicated in a genome duplication, we find that although nonsynonymous sites in the coding sequences show strong evidence for the fixation of recent gene conversion events, similar patterns are less evident among the synonymous positions and noncoding regulatory elements. This result suggests a potential explanation for the somewhat puzzling fact that duplicated RP genes are not functionally redundant despite their very high protein sequence identity. An analysis of the patterns of regulatory network evolution after genome duplication also indicates that the duplicated proteins have diverged considerably in expression despite their similar protein sequences.

The evolution of vertebrate tetraspanins: gene loss, retention, and massive positive selection after whole genome duplications

BACKGROUND

The vertebrate tetraspanin family has many features which make it suitable for preserving the imprint of ancient sequence evolution and amenable for phylogenomic analysis. So we believe that an in-depth analysis of the tetraspanin evolution not only provides more complete understanding of tetraspanin biology, but offers new insights into the influence of the two rounds of whole genome duplication (2R-WGD) at the origin of vertebrates.

RESULTS

A detailed phylogeny of vertebrate tetraspanins was constructed by using multiple lines of information, including sequence-based phylogenetics, key structural features, intron configuration and genomic synteny. In particular, a total of 38 modern tetraspanin ortholog lineages in bony vertebrates have been identified and subsequently classified into 17 ancestral lineages existing before 2R-WGD. Based on this phylogeny, we found that the ohnolog retention rate of tetraspanins after 2R-WGD was three times as the average (a rate similar to those of transcription factors and protein kinases). This high rate didn't increase the tetrapanin family size, but changed the family composition, possibly by displacing vertebrate-specific gene lineages with the lineages conserved across deuterostomes. We also found that the period from 2R-WGD to recent time is controlled by gene losses. Meanwhile, positive selection has been detected on 80% of the branches right after 2R-WGDs, which declines significantly on both magnitude and extensity on the following speciation branches. Notably, the loss of mammalian RDS2 is accompanied by strong positive selection on mammalian ROM1, possibly due to gene loss-induced compensatory evolution.

CONCLUSIONS

First, different from transcription factors and kinases, high duplicate retention rate after 2R-WGD didn't increase the tetraspanin family size but just reshaped the family composition. Second, the evolution of tetraspanins right after 2R-WGD had been impacted by a massive wave of gene loss and positive selection on coding sequences. Third, the lingering effect of 2R-WGD on tetraspanin gene loss and positive selection might last for 300-400 million years.

Hidden variability of floral homeotic B genes in Solanaceae provides a molecular basis for the evolution of novel functions

B-class MADS box genes specify petal and stamen identities in several core eudicot species. Members of the Solanaceae possess duplicate copies of these genes, allowing for diversification of function. To examine the changing roles of such duplicate orthologs, we assessed the functions of B-class genes in Nicotiana benthamiana and tomato (Solanum lycopersicum) using virus-induced gene silencing and RNA interference approaches. Loss of function of individual duplicates can have distinct phenotypes, yet complete loss of B-class gene function results in extreme homeotic transformations of petal and stamen identities. We also show that these duplicate gene products have qualitatively different protein-protein interaction capabilities and different regulatory roles. Thus, compensatory changes in B-class MADS box gene duplicate function have occurred in the Solanaceae, in that individual gene roles are distinct, but their combined functions are equivalent. Furthermore, we show that species-specific differences in the stamen regulatory network are associated with differences in the expression of the microRNA miR169. Whereas there is considerable plasticity in individual B-class MADS box transcription factor function, there is overall conservation in the roles of the multimeric MADS box B-class protein complexes, providing robustness in the specification of petal and stamen identities. Such hidden variability in gene function as we observe for individual B-class genes can provide a molecular basis for the evolution of regulatory functions that result in novel morphologies.

Genome-wide phylogenetic comparative analysis of plant transcriptional regulation: a timeline of loss, gain, expansion, and correlation with complexity

Evolutionary retention of duplicated genes encoding transcription-associated proteins (TAPs, comprising transcription factors and other transcriptional regulators) has been hypothesized to be positively correlated with increasing morphological complexity and paleopolyploidizations, especially within the plant kingdom. Here, we present the most comprehensive set of classification rules for TAPs and its application for genome-wide analyses of plants and algae. Using a dated species tree and phylogenetic comparative (PC) analyses, we define the timeline of TAP loss, gain, and expansion among Viridiplantae and find that two major bursts of gain/expansion occurred, coinciding with the water-to-land transition and the radiation of flowering plants. For the first time, we provide PC proof for the long-standing hypothesis that TAPs are major driving forces behind the evolution of morphological complexity, the latter in Plantae being shaped significantly by polyploidization and subsequent biased paleolog retention. Principal component analysis incorporating the number of TAPs per genome provides an alternate and significant proxy for complexity, ideally suited for PC genomics. Our work lays the ground for further interrogation of the shaping of gene regulatory networks underlying the evolution of organism complexity.

Polyploidy did not predate the evolution of nodulation in all legumes

BACKGROUND

Several lines of evidence indicate that polyploidy occurred by around 54 million years ago, early in the history of legume evolution, but it has not been known whether this event was confined to the papilionoid subfamily (Papilionoideae; e.g. beans, medics, lupins) or occurred earlier. Determining the timing of the polyploidy event is important for understanding whether polyploidy might have contributed to rapid diversification and radiation of the legumes near the origin of the family; and whether polyploidy might have provided genetic material that enabled the evolution of a novel organ, the nitrogen-fixing nodule. Although symbioses with nitrogen-fixing partners have evolved in several lineages in the rosid I clade, nodules are widespread only in legume taxa, being nearly universal in the papilionoids and in the mimosoid subfamily (e.g., mimosas, acacias)--which diverged from the papilionoid legumes around 58 million years ago, soon after the origin of the legumes.

METHODOLOGY/PRINCIPAL FINDINGS

Using transcriptome sequence data from Chamaecrista fasciculata, a nodulating member of the mimosoid clade, we tested whether this species underwent polyploidy within the timeframe of legume diversification. Analysis of gene family branching orders and synonymous-site divergence data from C. fasciculata, Glycine max (soybean), Medicago truncatula, and Vitis vinifera (grape; an outgroup to the rosid taxa) establish that the polyploidy event known from soybean and Medicago occurred after the separation of the mimosoid and papilionoid clades, and at or shortly before the Papilionoideae radiation.

CONCLUSIONS

The ancestral legume genome was not fundamentally polyploid. Moreover, because there has not been an independent instance of polyploidy in the Chamaecrista lineage there is no necessary connection between polyploidy and nodulation in legumes. Chamaecrista may serve as a useful model in the legumes that lacks a paleopolyploid history, at least relative to the widely studied papilionoid models.

Quantifying whole transcriptome size, a prerequisite for understanding transcriptome evolution across species: an example from a plant allopolyploid

Evolutionary biologists are increasingly comparing gene expression patterns across species. Due to the way in which expression assays are normalized, such studies provide no direct information about expression per gene copy (dosage responses) or per cell and can give a misleading picture of genes that are differentially expressed. We describe an assay for estimating relative expression per cell. When used in conjunction with transcript profiling data, it is possible to compare the sizes of whole transcriptomes, which in turn makes it possible to compare expression per cell for each gene in the transcript profiling data set. We applied this approach, using quantitative reverse transcriptase-polymerase chain reaction and high throughput RNA sequencing, to a recently formed allopolyploid and showed that its leaf transcriptome was approximately 1.4-fold larger than either progenitor transcriptome (70% of the sum of the progenitor transcriptomes). In contrast, the allopolyploid genome is 94.3% as large as the sum of its progenitor genomes and retains > or =93.5% of the sum of its progenitor gene complements. Thus, "transcriptome downsizing" is greater than genome downsizing. Using this transcriptome size estimate, we inferred dosage responses for several thousand genes and showed that the majority exhibit partial dosage compensation. Homoeologue silencing is nonrandomly distributed across dosage responses, with genes showing extreme responses in either direction significantly more likely to have a silent homoeologue. This experimental approach will add value to transcript profiling experiments involving interspecies and interploidy comparisons by converting expression per transcriptome to expression per genome, eliminating the need for assumptions about transcriptome size.

Life and death among plant lysophosphatidic acid acyltransferases

The tetraploid Brassica napus possesses several seed-expressed microsomal lysophosphatidic acid acyltransferases (LPAAT ) including BAT1.5, which has been retained after genome duplication as a consequence of a subfunctionalisation of the gene encoding the ubiquitously expressed Kennedy pathway enzyme BAT1.13. Next, cDNA BAT1.3, encoding a LPAAT was subsequently isolated from an embryo library. The rapeseed LPAAT encoded by BAT1.3 is orthologous to the Arabidopsis thaliana At1g51260 gene product possibly associated with tapetum development and male fertility. However, BAT1.3 expression is predominant during the mid stages of embryo development in seeds of Brassica napus. Functional characterisation of BAT1.3 provides further support for a hypothesis of gene dosage sensitivity of LPAATs as does an analysis of the chromosomal localisation of LPAAT genes in Arabidopsis thaliana. The pattern of retention or loss of LPAAT genes after polyploidisation or segmental duplication is consistent with a model of balanced gene drive.

Evidence for the involvement of Globosa-like gene duplications and expression divergence in the evolution of floral morphology in the Zingiberales

*The MADS box transcription factor family has long been identified as an important contributor to the control of floral development. It is often hypothesized that the evolution of floral development across angiosperms and within specific lineages may occur as a result of duplication, functional diversification, and changes in regulation of MADS box genes. Here we examine the role of Globosa (GLO)-like genes, members of the B-class MADS box gene lineage, in the evolution of floral development within the monocot order Zingiberales. *We assessed changes in perianth and stamen whorl morphology in a phylogenetic framework. We identified GLO homologs (ZinGLO1-4) from 50 Zingiberales species and investigated the evolution of this gene lineage. Expression of two GLO homologs was assessed in Costus spicatus and Musa basjoo. *Based on the phylogenetic data and expression results, we propose several family-specific losses and gains of GLO homologs that appear to be associated with key morphological changes. The GLO-like gene lineage has diversified concomitant with the evolution of the dimorphic perianth and the staminodial labellum. *Duplications and expression divergence within the GLO-like gene lineage may have played a role in floral diversification in the Zingiberales.

Following tetraploidy in maize, a short deletion mechanism removed genes preferentially from one of the two homologs

Following genome duplication and selfish DNA expansion, maize used a heretofore unknown mechanism to shed redundant genes and functionless DNA with bias toward one of the parental genomes.

Previous work in Arabidopsis showed that after an ancient tetraploidy event, genes were preferentially removed from one of the two homeologs, a process known as fractionation. The mechanism of fractionation is unknown. We sought to determine whether such preferential, or biased, fractionation exists in maize and, if so, whether a specific mechanism could be implicated in this process. We studied the process of fractionation using two recently sequenced grass species: sorghum and maize. The maize lineage has experienced a tetraploidy since its divergence from sorghum approximately 12 million years ago, and fragments of many knocked-out genes retain enough sequence similarity to be easily identifiable. Using sorghum exons as the query sequence, we studied the fate of both orthologous genes in maize following the maize tetraploidy. We show that genes are predominantly lost, not relocated, and that single-gene loss by deletion is the rule. Based on comparisons with orthologous sorghum and rice genes, we also infer that the sequences present before the deletion events were flanked by short direct repeats, a signature of intra-chromosomal recombination. Evidence of this deletion mechanism is found 2.3 times more frequently on one of the maize homeologs, consistent with earlier observations of biased fractionation. The over-fractionated homeolog is also a greater than 3-fold better target for transposon removal, but does not have an observably higher synonymous base substitution rate, nor could we find differentially placed methylation domains. We conclude that fractionation is indeed biased in maize and that intra-chromosomal or possibly a similar illegitimate recombination is the primary mechanism by which fractionation occurs. The mechanism of intra-chromosomal recombination explains the observed bias in both gene and transposon loss in the maize lineage. The existence of fractionation bias demonstrates that the frequency of deletion is modulated. Among the evolutionary benefits of this deletion/fractionation mechanism is bulk DNA removal and the generation of novel combinations of regulatory sequences and coding regions.

All genomes can accumulate dispensable DNA in the form of duplications of individual genes or even partial or whole genome duplications. Genomes also can accumulate selfish DNA elements. Duplication events specifically are often followed by extensive gene loss. The maize genome is particularly extreme, having become tetraploid 10 million years ago and played host to massive transposon amplifications. We compared the genome of sorghum (which is homologous to the pre-tetraploid maize genome) with the two identifiable parental genomes retained in maize. The two maize genomes differ greatly: one of the parental genomes has lost 2.3 times more genes than the other, and the selfish DNA regions between genes were even more frequently lost, suggesting maize can distinguish between the parental genomes present in the original tetraploid. We show that genes are actually lost, not simply relocated. Deletions were rarely longer than a single gene, and occurred between repeated DNA sequences, suggesting mis-recombination as a mechanism of gene removal. We hypothesize an epigenetic mechanism of genome distinction to account for the selective loss. To the extent that the rate of base substitutions tracks time, we neither support nor refute claims of maize allotetraploidy. Finally, we explain why it makes sense that purifying selection in mammals does not operate at all like the gene and genome deletion program we describe here.

The hidden duplication past of the plant pathogen Phytophthora and its consequences for infection

BACKGROUND

Oomycetes of the genus Phytophthora are pathogens that infect a wide range of plant species. For dicot hosts such as tomato, potato and soybean, Phytophthora is even the most important pathogen. Previous analyses of Phytophthora genomes uncovered many genes, large gene families and large genome sizes that can partially be explained by significant repeat expansion patterns.

RESULTS

Analysis of the complete genomes of three different Phytophthora species, using a newly developed approach, unveiled a large number of small duplicated blocks, mainly consisting of two or three consecutive genes. Further analysis of these duplicated genes and comparison with the known gene and genome duplication history of ten other eukaryotes including parasites, algae, plants, fungi, vertebrates and invertebrates, suggests that the ancestor of P. infestans, P. sojae and P. ramorum most likely underwent a whole genome duplication (WGD). Genes that have survived in duplicate are mainly genes that are known to be preferentially retained following WGDs, but also genes important for pathogenicity and infection of the different hosts seem to have been retained in excess. As a result, the WGD might have contributed to the evolutionary and pathogenic success of Phytophthora.

CONCLUSIONS

The fact that we find many small blocks of duplicated genes indicates that the genomes of Phytophthora species have been heavily rearranged following the WGD. Most likely, the high repeat content in these genomes have played an important role in this rearrangement process. As a consequence, the paucity of retained larger duplicated blocks has greatly complicated previous attempts to detect remnants of a large-scale duplication event in Phytophthora. However, as we show here, our newly developed strategy to identify very small duplicated blocks might be a useful approach to uncover ancient polyploidy events, in particular for heavily rearranged genomes.