A burst of protein sequence evolution and a prolonged period of asymmetric evolution follow gene duplication in yeast

Plant genome evolution in the genus Eucalyptus is driven by structural rearrangements that promote sequence divergence

Genomes have a highly organized architecture (nonrandom organization of functional and nonfunctional genetic elements within chromosomes) that is essential for many biological functions, particularly gene expression and reproduction. Despite the need to conserve genome architecture, a high level of structural variation has been observed within species. As species separate and diverge, genome architecture also diverges, becoming increasingly poorly conserved as divergence time increases. However, within plant genomes, the processes of genome architecture divergence are not well described. Here we use long-read sequencing and de novo assembly of 33 phylogenetically diverse, wild and naturally evolving Eucalyptus species, covering 1-50 million years of diverging genome evolution to measure genome architectural conservation and describe architectural divergence. The investigation of these genomes revealed that following lineage divergence, genome architecture is highly fragmented by rearrangements. As genomes continue to diverge, the accumulation of mutations and the subsequent divergence beyond recognition of rearrangements become the primary driver of genome divergence. The loss of syntenic regions also contribute to genome divergence but at a slower pace than that of rearrangements. We hypothesize that duplications and translocations are potentially the greatest contributors to Eucalyptus genome divergence.

Frustration can Limit the Adaptation of Promiscuous Enzymes Through Gene Duplication and Specialisation

Virtually all enzymes catalyse more than one reaction, a phenomenon known as enzyme promiscuity. It is unclear whether promiscuous enzymes are more often generalists that catalyse multiple reactions at similar rates or specialists that catalyse one reaction much more efficiently than other reactions. In addition, the factors that shape whether an enzyme evolves to be a generalist or a specialist are poorly understood. To address these questions, we follow a three-pronged approach. First, we examine the distribution of promiscuity in empirical enzymes reported in the BRENDA database. We find that the promiscuity distribution of empirical enzymes is bimodal. In other words, a large fraction of promiscuous enzymes are either generalists or specialists, with few intermediates. Second, we demonstrate that enzyme biophysics is not sufficient to explain this bimodal distribution. Third, we devise a constraint-based model of promiscuous enzymes undergoing duplication and facing selection pressures favouring subfunctionalization. The model posits the existence of constraints between the catalytic efficiencies of an enzyme for different reactions and is inspired by empirical case studies. The promiscuity distribution predicted by our constraint-based model is consistent with the empirical bimodal distribution. Our results suggest that subfunctionalization is possible and beneficial only in certain enzymes. Furthermore, the model predicts that conflicting constraints and selection pressures can cause promiscuous enzymes to enter a 'frustrated' state, in which competing interactions limit the specialisation of enzymes. We find that frustration can be both a driver and an inhibitor of enzyme evolution by duplication and subfunctionalization. In addition, our model predicts that frustration becomes more likely as enzymes catalyse more reactions, implying that natural selection may prefer catalytically simple enzymes. In sum, our results suggest that frustration may play an important role in enzyme evolution.

Parallel Nonfunctionalization of CK1δ/ε Kinase Ohnologs Following a Whole-Genome Duplication Event

Whole-genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in Saccharomyces cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post-WGD. In S. cerevisiae, HRR25 encodes the casein kinase 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post-WGD and non-WGD species, and have nonconserved cellular localization, consistent with their ongoing loss of function. The analysis in Naumovozyma castellii shows that the noncomplementing ohnolog is expressed at a lower level and has become nonessential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

A high-quality reference genome for the fission yeast Schizosaccharomyces osmophilus

Fission yeasts are an ancient group of fungal species that diverged from each other from tens to hundreds of million years ago. Among them is the preeminent model organism Schizosaccharomyces pombe, which has significantly contributed to our understandings of molecular mechanisms underlying fundamental cellular processes. The availability of the genomes of S. pombe and 3 other fission yeast species S. japonicus, S. octosporus, and S. cryophilus has enabled cross-species comparisons that provide insights into the evolution of genes, pathways, and genomes. Here, we performed genome sequencing on the type strain of the recently identified fission yeast species S. osmophilus and obtained a complete mitochondrial genome and a nuclear genome assembly with gaps only at rRNA gene arrays. A total of 5,098 protein-coding nuclear genes were annotated and orthologs for more than 95% of them were identified. Genome-based phylogenetic analysis showed that S. osmophilus is most closely related to S. octosporus and these 2 species diverged around 16 million years ago. To demonstrate the utility of this S. osmophilus reference genome, we conducted cross-species comparative analyses of centromeres, telomeres, transposons, the mating-type region, Cbp1 family proteins, and mitochondrial genomes. These analyses revealed conservation of repeat arrangements and sequence motifs in centromere cores, identified telomeric sequences composed of 2 types of repeats, delineated relationships among Tf1/sushi group retrotransposons, characterized the evolutionary origins and trajectories of Cbp1 family domesticated transposases, and discovered signs of interspecific transfer of 2 types of mitochondrial selfish elements.

Deletion of QDR genes in a bioethanol-producing yeast strain reduces propagation of contaminating lactic acid bacteria

Bacterial contaminations in yeast fermentation tanks are a recurring problem for the bioethanol production industry. Lactic acid bacteria (LAB), particularly of the genus Lactobacillus, are the most common contaminants. Their proliferation can reduce fermentation efficiency or even impose premature shutdown for cleaning. We have previously reported that laboratory yeast strains naturally excrete amino acids via transporters of the Drug: H⁺ Antiporter-1 (DHA1) family. This excretion allows yeast to cross-feed LAB, which are most often unable to grow without an external amino acid supply. Whether industrial yeast strains used in bioethanol production likewise promote LAB proliferation through cross-feeding has not been investigated. In this study, we first show that the yeast strain Ethanol Red used in ethanol production supports growth of Lactobacillus fermentum in an amino-acid-free synthetic medium. This effect was markedly reduced upon homozygous deletion of the QDR3 gene encoding a DHA1-family amino acid exporter. We further show that cultivation of Ethanol Red in a nonsterile sugarcane-molasses-based medium is associated with an increase in lactic acid due to LAB growth. When Ethanol Red lacked the QDR1, QDR2, and QDR3 genes, this lactic acid production was not observed and ethanol production was not significantly reduced. Our results indicate that Ethanol Red cultivated in synthetic or molasses medium sustains LAB proliferation in a manner that depends on its ability to excrete amino acids via Qdr transporters. They further suggest that using mutant industrial yeast derivatives lacking DHA1-family amino acid exporters may be a way to reduce the risk of bacterial contaminations during fermentation.

Comparative gene retention analysis in barley, wild emmer, and bread wheat pangenome lines reveals factors affecting gene retention following gene duplication

BACKGROUND

Gene duplication is a prevalent phenomenon and a major driving force underlying genome evolution. The process leading to the fixation of gene duplicates following duplication is critical to understand how genome evolves but remains fragmentally understood. Most previous studies on gene retention are based on gene duplicate analyses in single reference genome. No population-based comparative gene retention analysis has been performed to date.

RESULTS

Taking advantage of recently published genomic data in Triticeae, we dissected a divergent homogentisate phytyltransferase (HPT2) lineage caught in the middle stage of gene fixation following duplication. The presence/absence of HPT2 in barley (diploid), wild emmer (tetraploid), and bread wheat (hexaploid) pangenome lines appears to be associated with gene dosage constraint and environmental adaption. Based on these observations, we adopted a phylogeny-based orthology inference approach and performed comparative gene retention analyses across barley, wild emmer, and bread wheat. This led to the identification of 326 HPT2-pattern-like genes at whole genome scale, representing a pool of gene duplicates in the middle stage of gene fixation. Majority of these HPT2-pattern-like genes were identified as small-scale duplicates, such as dispersed, tandem, and proximal duplications. Natural selection analyses showed that HPT2-pattern-like genes have experienced relaxed selection pressure, which is generally accompanied with partial positive selection and transcriptional divergence. Functional enrichment analyses showed that HPT2-pattern-like genes are over-represented with molecular-binding and defense response functions, supporting the potential role of environmental adaption during gene retention. We also observed that gene duplicates from larger gene family are more likely to be lost, implying a gene dosage constraint effect. Further comparative gene retention analysis in barley and bread wheat pangenome lines revealed combined effects of species-specific selection and gene dosage constraint.

CONCLUSIONS

Comparative gene retention analyses at the population level support gene dosage constraint, environmental adaption, and species-specific selection as three factors that may affect gene retention following gene duplication. Our findings shed light on the evolutionary process leading to the retention of newly formed gene duplicates and will greatly improve our understanding on genome evolution via duplication.

Evolution of TOP1 and TOP1MT Topoisomerases in Chordata

Type IB topoisomerases relax the torsional stress associated with DNA metabolism in the nucleus and mitochondria and constitute important molecular targets of anticancer drugs. Vertebrates stand out among eukaryotes by having two Type IB topoisomerases acting specifically in the nucleus (TOP1) and mitochondria (TOP1MT). Despite their major importance, the origin and evolution of these paralogues remain unknown. Here, we examine the molecular evolutionary processes acting on both TOP1 and TOP1MT in Chordata, taking advantage of the increasing number of available genome sequences. We found that both TOP1 and TOP1MT evolved under strong purifying selection, as expected considering their essential biological functions. Critical active sites, including those associated with resistance to anticancer agents, were found particularly conserved. However, TOP1MT presented a higher rate of molecular evolution than TOP1, possibly related with its specialized activity on the mitochondrial genome and a less critical role in cells. We could place the duplication event that originated the TOP1 and TOP1MT paralogues early in the radiation of vertebrates, most likely associated with the first round of vertebrate tetraploidization (1R). Moreover, our data suggest that cyclostomes present a specialized mitochondrial Type IB topoisomerase. Interestingly, we identified two missense mutations replacing amino acids in the Linker region of TOP1MT in Neanderthals, which appears as a rare event when comparing the genome of both species. In conclusion, TOP1 and TOP1MT differ in their rates of evolution, and their evolutionary histories allowed us to better understand the evolution of chordates.

How gene duplication diversifies the landscape of protein oligomeric state and function

Oligomeric proteins are central to cellular life and the duplication and divergence of their genes is a key driver of evolutionary innovations. The duplication of a gene coding for an oligomeric protein has numerous possible outcomes, which motivates questions on the relationship between structural and functional divergence. How do protein oligomeric states diversify after gene duplication? In the simple case of duplication of a homo-oligomeric protein gene, what properties can influence the fate of descendant paralogs toward forming independent homomers or maintaining their interaction as a complex? Furthermore, how are functional innovations associated with the diversification of oligomeric states? Here, we review recent literature and present specific examples in an attempt to illustrate and answer these questions.

A Population-Genetic Lens into the Process of Gene Loss Following Whole-Genome Duplication

Whole-genome duplications (WGDs) have occurred in many eukaryotic lineages. However, the underlying evolutionary forces and molecular mechanisms responsible for the long-term retention of gene duplicates created by WGDs are not well understood. We employ a population-genomic approach to understand the selective forces acting on paralogs and investigate ongoing duplicate-gene loss in multiple species of Paramecium that share an ancient WGD. We show that mutations that abolish protein function are more likely to be segregating in retained WGD paralogs than in single-copy genes, most likely because of ongoing nonfunctionalization post-WGD. This relaxation of purifying selection occurs in only one WGD paralog, accompanied by the gradual fixation of nonsynonymous mutations and reduction in levels of expression, and occurs over a long period of evolutionary time, "marking" one locus for future loss. Concordantly, the fitness effects of new nonsynonymous mutations and frameshift-causing indels are significantly more deleterious in the highly expressed copy compared with their paralogs with lower expression. Our results provide a novel mechanistic model of gene duplicate loss following WGDs, wherein selection acts on the sum of functional activity of both duplicate genes, allowing the two to wander in expression and functional space, until one duplicate locus eventually degenerates enough in functional efficiency or expression that its contribution to total activity is too insignificant to be retained by purifying selection. Retention of duplicates by such mechanisms predicts long times to duplicate-gene loss, which should not be falsely attributed to retention due to gain/change in function.

Evolutionary History of TOPIIA Topoisomerases in Animals

TOPIIA topoisomerases are required for the regulation of DNA topology by DNA cleavage and re-ligation and are important targets of antibiotic and anticancer agents. Humans possess two TOPIIA paralogue genes (TOP2A and TOP2B) with high sequence and structural similarity but distinct cellular functions. Despite their functional and clinical relevance, the evolutionary history of TOPIIA is still poorly understood. Here we show that TOPIIA is highly conserved in Metazoa. We also found that TOPIIA paralogues from jawed and jawless vertebrates had different origins related with tetraploidization events. After duplication, TOP2B evolved under a stronger purifying selection than TOP2A, perhaps promoted by the more specialized role of TOP2B in postmitotic cells. We also detected genetic signatures of positive selection in the highly variable C-terminal domain (CTD), possibly associated with adaptation to cellular interactions. By comparing TOPIIA from modern and archaic humans, we found two amino acid substitutions in the TOP2A CTD, suggesting that TOP2A may have contributed to the evolution of present-day humans, as proposed for other cell cycle-related genes. Finally, we identified six residues conferring resistance to chemotherapy differing between TOP2A and TOP2B. These six residues could be targets for the development of TOP2A-specific inhibitors that would avoid the side effects caused by inhibiting TOP2B. Altogether, our findings clarify the origin, diversification and selection pressures governing the evolution of animal TOPIIA.

Special Care Is Needed in Applying Phylogenetic Comparative Methods to Gene Trees with Speciation and Duplication Nodes

How gene function evolves is a central question of evolutionary biology. It can be investigated by comparing functional genomics results between species and between genes. Most comparative studies of functional genomics have used pairwise comparisons. Yet it has been shown that this can provide biased results, as genes, like species, are phylogenetically related. Phylogenetic comparative methods should be used to correct for this, but they depend on strong assumptions, including unbiased tree estimates relative to the hypothesis being tested. Such methods have recently been used to test the “ortholog conjecture,” the hypothesis that functional evolution is faster in paralogs than in orthologs. Although pairwise comparisons of tissue specificity (τ) provided support for the ortholog conjecture, phylogenetic independent contrasts did not. Our reanalysis on the same gene trees identified problems with the time calibration of duplication nodes. We find that the gene trees used suffer from important biases, due to the inclusion of trees with no duplication nodes, to the relative age of speciations and duplications, to systematic differences in branch lengths, and to non-Brownian motion of tissue specificity on many trees. We find that incorrect implementation of phylogenetic method in empirical gene trees with duplications can be problematic. Controlling for biases allows successful use of phylogenetic methods to study the evolution of gene function and provides some support for the ortholog conjecture using three different phylogenetic approaches.

Selection following Gene Duplication Shapes Recent Genome Evolution in the Pea Aphid Acyrthosiphon pisum

Ecology of insects is as wide as their diversity, which reflects their high capacity of adaptation in most of the environments of our planet. Aphids, with over 4,000 species, have developed a series of adaptations including a high phenotypic plasticity and the ability to feed on the phloem sap of plants, which is enriched in sugars derived from photosynthesis. Recent analyses of aphid genomes have indicated a high level of shared ancestral gene duplications that might represent a basis for genetic innovation and broad adaptations. In addition, there are a large number of recent, species-specific gene duplications whose role in adaptation remains poorly understood. Here, we tested whether duplicates specific to the pea aphid Acyrthosiphon pisum are related to genomic innovation by combining comparative genomics, transcriptomics, and chromatin accessibility analyses. Consistent with large levels of neofunctionalization, we found that most of the recent pairs of gene duplicates evolved asymmetrically, showing divergent patterns of positive selection and gene expression. Genes under selection involved a plethora of biological functions, suggesting that neofunctionalization and tissue specificity, among other evolutionary mechanisms, have orchestrated the evolution of recent paralogs in the pea aphid and may have facilitated host-symbiont cooperation. Our comprehensive phylogenomics analysis allowed us to tackle the history of duplicated genes to pave the road toward understanding the role of gene duplication in ecological adaptation.

Genes derived from ancient polyploidy have higher genetic diversity and are associated with domestication in Brassica rapa

Many crops are polyploid or have a polyploid ancestry. Recent phylogenetic analyses have found that polyploidy often preceded the domestication of crop plants. One explanation for this observation is that increased genetic diversity following polyploidy may have been important during the strong artificial selection that occurs during domestication. In order to test the connection between domestication and polyploidy, we identified and examined candidate genes associated with the domestication of the diverse crop varieties of Brassica rapa. Like all 'diploid' flowering plants, B. rapa has a diploidized paleopolyploid genome and experienced many rounds of whole genome duplication (WGD). We analyzed transcriptome data of more than 100 cultivated B. rapa accessions. Using a combination of approaches, we identified > 3000 candidate genes associated with the domestication of four major B. rapa crop varieties. Consistent with our expectation, we found that the candidate genes were significantly enriched with genes derived from the Brassiceae mesohexaploidy. We also observed that paleologs were significantly more diverse than non-paleologs. Our analyses find evidence for that genetic diversity derived from ancient polyploidy played a key role in the domestication of B. rapa and provide support for its importance in the success of modern agriculture.

Genetic, transcriptional, and regulatory landscape of monolignol biosynthesis pathway in Miscanthus × giganteus

BACKGROUND

Miscanthus × giganteus is widely recognized as a promising lignocellulosic biomass crop due to its advantages of high biomass production, low environmental impacts, and the potential to be cultivated on marginal land. However, the high costs of bioethanol production still limit the current commercialization of lignocellulosic bioethanol. The lignin in the cell wall and its by-products released in the pretreatment step is the main component inhibiting the enzymatic reactions in the saccharification and fermentation processes. Hence, genetic modification of the genes involved in lignin biosynthesis could be a feasible strategy to overcome this barrier by manipulating the lignin content and composition of M. × giganteus. For this purpose, the essential knowledge of these genes and understanding the underlying regulatory mechanisms in M. × giganteus is required.

RESULTS

In this study, MgPAL1, MgPAL5, Mg4CL1, Mg4CL3, MgHCT1, MgHCT2, MgC3'H1, MgCCoAOMT1, MgCCoAOMT3, MgCCR1, MgCCR2, MgF5H, MgCOMT, and MgCAD were identified as the major monolignol biosynthetic genes in M. × giganteus based on genetic and transcriptional evidence. Among them, 12 genes were cloned and sequenced. By combining transcription factor binding site prediction and expression correlation analysis, MYB46, MYB61, MYB63, WRKY24, WRKY35, WRKY12, ERF021, ERF058, and ERF017 were inferred to regulate the expression of these genes directly. On the basis of these results, an integrated model was summarized to depict the monolignol biosynthesis pathway and the underlying regulatory mechanism in M. × giganteus.

CONCLUSIONS

This study provides a list of potential gene targets for genetic improvement of lignocellulosic biomass quality of M. × giganteus, and reveals the genetic, transcriptional, and regulatory landscape of the monolignol biosynthesis pathway in M. × giganteus.

Patterns of gene evolution following duplications and speciations in vertebrates

BACKGROUND

Eukaryotic genes typically form independent evolutionary lineages through either speciation or gene duplication events. Generally, gene copies resulting from speciation events (orthologs) are expected to maintain similarity over time with regard to sequence, structure and function. After a duplication event, however, resulting gene copies (paralogs) may experience a broader set of possible fates, including partial (subfunctionalization) or complete loss of function, as well as gain of new function (neofunctionalization). This assumption, known as the Ortholog Conjecture, is prevalent throughout molecular biology and notably plays an important role in many functional annotation methods. Unfortunately, studies that explicitly compare evolutionary processes between speciation and duplication events are rare and conflicting.

METHODS

To provide an empirical assessment of ortholog/paralog evolution, we estimated ratios of nonsynonymous to synonymous substitutions (ω = dN/dS) for 251,044 lineages in 6,244 gene trees across 77 vertebrate taxa.

RESULTS

Overall, we found ω to be more similar between lineages descended from speciation events (p < 0.001) than lineages descended from duplication events, providing strong support for the Ortholog Conjecture. The asymmetry in ω following duplication events appears to be largely driven by an increase along one of the paralogous lineages, while the other remains similar to the parent. This trend is commonly associated with neofunctionalization, suggesting that gene duplication is a significant mechanism for generating novel gene functions.

A Single, Shared Triploidy in Three Species of Parasitic Nematodes

The root-knot nematodes of the genus Meloidogyne are important and damaging parasites capable of infecting most flowering plants. Within this genus, several species of the Meloidogyne incognita group show evidence of paleopolyploidy in their genomes. We used our software tool POInT, the Polyploidy Orthology Inference Tool, to phylogenetically model the gene losses that followed that polyploidy. These models, and simulations based on them, show that three of these species (M. incognita, M. arenaria and M. javanica) descend from a single common hybridization event that yielded triplicated genomes with three distinguishable subgenomes. While one of the three subgenomes shows elevated gene loss rates relative to the other two, this subgenome does not show elevated sequence divergence. In all three species, ancestral loci where two of the three gene copies have been lost are less likely to have orthologs in Caenorhabditis elegans that are lethal when knocked down than are ancestral loci with surviving duplicate copies.

Contribution of Functional Divergence Through Copy Number Variations to the Inter-Species and Intra-Species Diversity in Specialized Metabolites

There is considerable diversity in the specialized metabolites within a single plant species (intra-species) and among different plant species (inter-species). The functional divergence associated with gene duplications largely contributes to the inter-species diversity in the specialized metabolites, whereas the intra-species diversity is due to gene dosage changes via gene duplications [i.e., copy number variants (CNVs)] at the intra-species level of evolution. This is because CNVs are thought to undergo associated with less functional divergence at the intra-species level of evolution. However, functional divergence caused by CNVs may induce specialized metabolite diversity at the intra-species and inter-species levels of evolution. We herein discuss the functional divergence of CNVs in metabolic quantitative trait genes (mQTGs). We focused on 5,654 previously identified mQTGs in 270 Arabidopsis thaliana accessions. The ratio of nonsynonymous to synonymous variations tends to be higher for mQTGs with CNVs than for mQTGs without CNVs within A. thaliana accessions, suggesting that CNVs are responsible for the functional divergence among mQTGs at the intra-species level of evolution. To evaluate the contribution of CNVs to inter-species diversity, we calculated the ratio of nonsynonymous to synonymous substitutions in the Arabidopsis lineage. The ratio tends to be higher for the mQTGs with CNVs than for the mQTGs without CNVs. Additionally, we determined that mQTGs with CNVs are subject to positive selection in the Arabidopsis lineage. Our data suggest that CNVs are closely related to functional divergence contributing to adaptations via the production of diverse specialized metabolites at the intra-species and inter-species levels of evolution.

Rapid and pervasive development- and tissue-specific homeolog expression partitioning in newly formed inter-subspecific rice segmental allotetraploids

Background

In diverse plant taxa, whole-genome duplication (WGD) events are major sources of phenotypic novelty. Studies of gene expression in synthetic polyploids have shown immediate expression and functional partitioning of duplicated genes among different tissues. Many studies of the tissue-specific homeolog expression partitioning have focused on allopolyploids that have very different parental genomes, while few studies have focused on autopolyploids or allopolyploids that have similar parental genomes.

Results

In this study, we used a set of reciprocal F1 hybrids and synthetic tetraploids constructed from subspecies (japonica and indica) of Asian rice (Oryza sativa L.) as a model to gain insights into the expression partitioning of homeologs among tissues in a developmental context. We assayed the tissue-specific silencing (TSS) of the parental homeologs of 30 key genes in the hybrids and tetraploids relative to the in vitro “hybrids” (parental mixes) using Sequenom MassARRAY. We found that the parental mix and synthetic tetraploids had higher frequencies of homeolog TSS than the F1, revealing an instantaneous role of WGD on homeolog expression partitioning.

Conclusions

Our observations contradicted those of previous studies in which newly formed allopolyploids had a low TSS frequency, similar to that of F1 hybrids, suggesting that the impact of WGD on homeolog expression requires a longer time to manifest. In addition, we found that the TSS frequency in the tetraploids varied at different growth stages and that roots had a much higher frequency of TSS than leaves, which indicated that developmental and metabolic traits may influence the expression states of duplicated genes in newly formed plant polyploids.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5150-7) contains supplementary material, which is available to authorized users.

Patterns of Population Variation in Two Paleopolyploid Eudicot Lineages Suggest That Dosage-Based Selection on Homeologs Is Long-Lived

Genes that are inherently subject to strong selective constraints tend to be overretained in duplicate after polyploidy. They also continue to experience similar, but somewhat relaxed, constraints after that polyploidy event. We sought to assess for how long the influence of polyploidy is felt on these genes’ selective pressures. We analyzed two nested polyploidy events in Brassicaceae: the At-α genome duplication that is the most recent polyploidy in the model plant Arabidopsis thaliana and a more recent hexaploidy shared by the genus Brassica and its relatives. By comparing the strength and direction of the natural selection acting at the population and at the species level, we find evidence for continued intensified purifying selection acting on retained duplicates from both polyploidies even down to the present. The constraint observed in preferentially retained genes is not a result of the polyploidy event: the orthologs of such genes experience even stronger constraint in nonpolyploid outgroup genomes. In both the Arabidopsis and Brassica lineages, we further find evidence for segregating mildly deleterious variants, confirming that the population-level data uncover patterns not visible with between-species comparisons. Using the A. thaliana metabolic network, we also explored whether network position was correlated with the measured selective constraint. At both the population and species level, nodes/genes tended to show similar constraints to their neighbors. Our results paint a picture of the long-lived effects of polyploidy on plant genomes, suggesting that even yesterday’s polyploids still have distinct evolutionary trajectories.

Expression properties exhibit correlated patterns with the fate of duplicated genes, their divergence, and transcriptional plasticity in Saccharomycotina

Gene duplication is an important source of novelties and genome complexity. What genes are preserved as duplicated through long evolutionary times can shape the evolution of innovations. Identifying factors that influence gene duplicability is therefore an important aim in evolutionary biology. Here, we show that in the yeast Saccharomyces cerevisiae the levels of gene expression correlate with gene duplicability, its divergence, and transcriptional plasticity. Genes that were highly expressed before duplication are more likely to be preserved as duplicates for longer evolutionary times and wider phylogenetic ranges than genes that were lowly expressed. Duplicates with higher expression levels exhibit greater divergence between their gene copies. Duplicates that exhibit higher expression divergence are those enriched for TATA-containing promoters. These duplicates also show transcriptional plasticity, which seems to be involved in the origin of adaptations to environmental stresses in yeast. While the expression properties of genes strongly affect their duplicability, divergence and transcriptional plasticity are enhanced after gene duplication. We conclude that highly expressed genes are more likely to be preserved as duplicates due to their promoter architectures, their greater tolerance to expression noise, and their ability to reduce the noise-plasticity conflict.

Position Matters: Network Centrality Considerably Impacts Rates of Protein Evolution in the Human Protein-Protein Interaction Network

The proteins of any organism evolve at disparate rates. A long list of factors affecting rates of protein evolution have been identified. However, the relative importance of each factor in determining rates of protein evolution remains unresolved. The prevailing view is that evolutionary rates are dominantly determined by gene expression, and that other factors such as network centrality have only a marginal effect, if any. However, this view is largely based on analyses in yeasts, and accurately measuring the importance of the determinants of rates of protein evolution is complicated by the fact that the different factors are often correlated with each other, and by the relatively poor quality of available functional genomics data sets. Here, we use correlation, partial correlation and principal component regression analyses to measure the contributions of several factors to the variability of the rates of evolution of human proteins. For this purpose, we analyzed the entire human protein–protein interaction data set and the human signal transduction network—a network data set of exceptionally high quality, obtained by manual curation, which is expected to be virtually free from false positives. In contrast with the prevailing view, we observe that network centrality (measured as the number of physical and nonphysical interactions, betweenness, and closeness) has a considerable impact on rates of protein evolution. Surprisingly, the impact of centrality on rates of protein evolution seems to be comparable, or even superior according to some analyses, to that of gene expression. Our observations seem to be independent of potentially confounding factors and from the limitations (biases and errors) of interactomic data sets.

Increased rates of protein evolution and asymmetric deceleration after the whole-genome duplication in yeasts

Background

Whole-genome duplication (WGD) events have shaped the genomes of eukaryotic organisms. Relaxed selection after duplication along with inherent functional constraints are thought to determine the fate of the paralogs and, ultimately, the evolution of gene function. Here, we investigated the rate of protein evolution (as measured by dN/dS ratios) before and after the WGD in the hemiascomycete yeasts, and the way in which changes in such rates relate to molecular and biological function.

Results

For most groups of orthologous genes (81%) we observed a change in the rates of evolution after genome duplication. Genes with atypically-low dN/dS ratio before the WGD were prone to increase their rates of evolution after duplication. Importantly, the paralogs were often different in their rates of evolution after the WGD (50% cases), however, this was more consistent with an asymmetric deceleration in the protein-evolution rates, rather than an asymmetric increase of the initial rates. Functional-category analysis showed that regulatory proteins such as protein kinases and transcription factors were enriched in genes that increase their rates of evolution after the WGD. While changes in the rate of protein-sequence evolution were associated to protein abundance, content of disordered regions, and contribution to fitness, these features were an attribute of specific functional classes.

Conclusions

Our results indicate that strong purifying selection in ancestral pre-duplication sequences is a strong predictor of increased rates after the duplication in yeasts and that asymmetry in evolution rate is established during the deceleration phase. In addition, changes in the rates at which paralogous sequences evolve before and after WGD are different for specific protein functions; increased rates of protein evolution after duplication occur preferentially in specific protein functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0895-1) contains supplementary material, which is available to authorized users.

Evolution of Gene Duplication in Plants

Ancient duplication events and a high rate of retention of extant pairs of duplicate genes have contributed to an abundance of duplicate genes in plant genomes. These duplicates have contributed to the evolution of novel functions, such as the production of floral structures, induction of disease resistance, and adaptation to stress. Additionally, recent whole-genome duplications that have occurred in the lineages of several domesticated crop species, including wheat (Triticum aestivum), cotton (Gossypium hirsutum), and soybean (Glycine max), have contributed to important agronomic traits, such as grain quality, fruit shape, and flowering time. Therefore, understanding the mechanisms and impacts of gene duplication will be important to future studies of plants in general and of agronomically important crops in particular. In this review, we survey the current knowledge about gene duplication, including gene duplication mechanisms, the potential fates of duplicate genes, models explaining duplicate gene retention, the properties that distinguish duplicate from singleton genes, and the evolutionary impact of gene duplication.

Comparing the Statistical Fate of Paralogous and Orthologous Sequences

For several decades, sequence alignment has been a widely used tool in bioinformatics. For instance, finding homologous sequences with a known function in large databases is used to get insight into the function of nonannotated genomic regions. Very efficient tools like BLAST have been developed to identify and rank possible homologous sequences. To estimate the significance of the homology, the ranking of alignment scores takes a background model for random sequences into account. Using this model we can estimate the probability to find two exactly matching subsequences by chance in two unrelated sequences. For two homologous sequences, the corresponding probability is much higher, which allows us to identify them. Here we focus on the distribution of lengths of exact sequence matches between protein-coding regions of pairs of evolutionarily distant genomes. We show that this distribution exhibits a power-law tail with an exponent [Formula: see text] Developing a simple model of sequence evolution by substitutions and segmental duplications, we show analytically and computationally that paralogous and orthologous gene pairs contribute differently to this distribution. Our model explains the differences observed in the comparison of coding and noncoding parts of genomes, thus providing a better understanding of statistical properties of genomic sequences and their evolution.

The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons

To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences identified in human genome-wide association studies. Transcriptomic analyses showed that the sums of expression domains and expression levels for duplicated teleost genes often approximate the patterns and levels of expression for gar genes, consistent with subfunctionalization. The gar genome provides a resource for understanding evolution after genome duplication, the origin of vertebrate genomes and the function of human regulatory sequences.

Binding interface change and cryptic variation in the evolution of protein-protein interactions

Background

Physical interactions between proteins are essential for almost all biological functions and systems. To understand the evolution of function it is therefore important to understand the evolution of molecular interactions. Of key importance is the evolution of binding specificity, the set of interactions made by a protein, since change in specificity can lead to “rewiring” of interaction networks. Unfortunately, the interfaces through which proteins interact are complex, typically containing many amino-acid residues that collectively must contribute to binding specificity as well as binding affinity, structural integrity of the interface and solubility in the unbound state.

Results

In order to study the relationship between interface composition and binding specificity, we make use of paralogous pairs of yeast proteins. Immediately after duplication these paralogues will have identical sequences and protein products that make an identical set of interactions. As the sequences diverge, we can correlate amino-acid change in the interface with any change in the specificity of binding. We show that change in interface regions correlates only weakly with change in specificity, and many variants in interfaces are functionally equivalent. We show that many of the residue replacements within interfaces are silent with respect to their contribution to binding specificity.

Conclusions

We conclude that such functionally-equivalent change has the potential to contribute to evolutionary plasticity in interfaces by creating cryptic variation, which in turn may provide the raw material for functional innovation and coevolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0608-1) contains supplementary material, which is available to authorized users.

Compensatory Drift and the Evolutionary Dynamics of Dosage-Sensitive Duplicate Genes

Dosage-balance selection preserves functionally redundant duplicates (paralogs) at the optimum for their combined expression. Here we present a model of the dynamics of duplicate genes coevolving under dosage-balance selection. We call this the compensatory drift model. Results show that even when strong dosage-balance selection constrains total expression to the optimum, expression of each duplicate can diverge by drift from its original level. The rate of divergence slows as the strength of stabilizing selection, the size of the mutation effect, and/or the size of the population increases. We show that dosage-balance selection impedes neofunctionalization early after duplication but can later facilitate it. We fit this model to data from sodium channel duplicates in 10 families of teleost fish; these include two convergent lineages of electric fish in which one of the duplicates neofunctionalized. Using the model, we estimated the strength of dosage-balance selection for these genes. The results indicate that functionally redundant paralogs still may undergo radical functional changes after a prolonged period of compensatory drift.

MglC, a Paralog of Myxococcus xanthus GTPase-Activating Protein MglB, Plays a Divergent Role in Motility Regulation

In order to optimize interactions with their environment and one another, bacteria regulate their motility. In the case of the rod-shaped cells of Myxococcus xanthus, regulated motility is essential for social behaviors. M. xanthus moves over surfaces using type IV pilus-dependent motility and gliding motility. These two motility systems are coordinated by a protein module that controls cell polarity and consists of three polarly localized proteins, the small G protein MglA, the cognate MglA GTPase-activating protein MglB, and the response regulator RomR. Cellular reversals are induced by the Frz chemosensory system, and the output response regulator of this system, FrzZ, interfaces with the MglA/MglB/RomR module to invert cell polarity. Using a computational approach, we identify a paralog of MglB, MXAN_5770 (MglC). Genetic epistasis experiments demonstrate that MglC functions in the same pathway as MglA, MglB, RomR, and FrzZ and is important for regulating cellular reversals. Like MglB, MglC localizes to the cell poles asymmetrically and with a large cluster at the lagging pole. Correct polar localization of MglC depends on RomR and MglB. Consistently, MglC interacts directly with MglB and the C-terminal output domain of RomR, and we identified a surface of MglC that is necessary for the interaction with MglB and for MglC function. Together, our findings identify an additional member of the M. xanthus polarity module involved in regulating motility and demonstrate how gene duplication followed by functional divergence can add a layer of control to the complex cellular processes of motility and motility regulation.

IMPORTANCE Gene duplication and the subsequent divergence of the duplicated genes are important evolutionary mechanisms for increasing both biological complexity and regulation of biological processes. The bacterium Myxococcus xanthus is a soil bacterium with an unusually large genome that carries out several social processes, including predation of other bacterial species and formation of multicellular, spore-filled fruiting bodies. One feature of the large M. xanthus genome is that it contains many gene duplications. Here, we compare the products of one example of gene duplication and divergence, in which a paralog of the cognate MglA GTPase-activating protein MglB has acquired a different and opposing role in the regulation of cellular polarity and motility, processes critical to the bacterium's social behaviors.

Drosophila innate immunity: regional and functional specialization of prophenoloxidases

Background

The diversification of immune systems during evolution involves the expansion of particular gene families in given phyla. A better understanding of the metazoan immune system requires an analysis of the logic underlying such immune gene amplification. This analysis is now within reach due to the ease with which we can generate multiple mutations in an organism. In this paper, we analyze the contribution of the three Drosophila prophenoloxidases (PPOs) to host defense by generating single, double and triple mutants. PPOs are enzymes that catalyze the production of melanin at the site of infection and around parasites. They are the rate-limiting enzymes that contribute to the melanization reaction, a major immune mechanism of arthropods. The number of PPO-encoding genes is variable among insects, ranging from one in the bee to ten in the mosquito.

Results

By analyzing mutations alone and in combination, we ascribe a specific function to each of the three PPOs of Drosophila. Our study confirms that two PPOs produced by crystal cells, PPO1 and PPO2, contribute to the bulk of melanization in the hemolymph, upon septic or clean injury. In contrast, PPO3, a PPO restricted to the D. melanogaster group, is expressed in lamellocytes and contributes to melanization during the encapsulation process. Interestingly, another overlapping set of PPOs, PPO2 and PPO3, achieve melanization of the capsule upon parasitoid wasp infection.

Conclusions

The use of single or combined mutations allowed us to show that each PPO mutant has a specific phenotype, and that knocking out two of three genes is required to abolish fully a particular function. Thus, Drosophila PPOs have partially overlapping functions to optimize melanization in at least two conditions: following injury or during encapsulation. Since PPO3 is restricted to the D. melanogaster group, this suggests that production of PPO by lamellocytes emerged as a recent defense mechanism against parasitoid wasps. We conclude that differences in spatial localization, immediate or late availability, and mode of activation underlie the functional diversification of the three Drosophila PPOs, with each of them having non-redundant but overlapping functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0193-6) contains supplementary material, which is available to authorized users.

Key Role of Amino Acid Repeat Expansions in the Functional Diversification of Duplicated Transcription Factors

The high regulatory complexity of vertebrates has been related to two rounds of whole genome duplication (2R-WGD) that occurred before the divergence of the major vertebrate groups. Following these events, many developmental transcription factors (TFs) were retained in multiple copies and subsequently specialized in diverse functions, whereas others reverted to their singleton state. TFs are known to be generally rich in amino acid repeats or low-complexity regions (LCRs), such as polyalanine or polyglutamine runs, which can evolve rapidly and potentially influence the transcriptional activity of the protein. Here we test the hypothesis that LCRs have played a major role in the diversification of TF gene duplicates. We find that nearly half of the TF gene families originated during the 2R-WGD contains LCRs. The number of gene duplicates with LCRs is 155 out of 550 analyzed (28%), about twice as many as the number of single copy genes with LCRs (15 out of 115, 13%). In addition, duplicated TFs preferentially accumulate certain LCR types, the most prominent of which are alanine repeats. We experimentally test the role of alanine-rich LCRs in two different TF gene families, PHOX2A/PHOX2B and LHX2/LHX9. In both cases, the presence of the alanine-rich LCR in one of the copies (PHOX2B and LHX2) significantly increases the capacity of the TF to activate transcription. Taken together, the results provide strong evidence that LCRs are important driving forces of evolutionary change in duplicated genes.

Differential gene retention as an evolutionary mechanism to generate biodiversity and adaptation in yeasts

The evolutionary history of the characters underlying the adaptation of microorganisms to food and biotechnological uses is poorly understood. We undertook comparative genomics to investigate evolutionary relationships of the dairy yeast Geotrichum candidum within Saccharomycotina. Surprisingly, a remarkable proportion of genes showed discordant phylogenies, clustering with the filamentous fungus subphylum (Pezizomycotina), rather than the yeast subphylum (Saccharomycotina), of the Ascomycota. These genes appear not to be the result of Horizontal Gene Transfer (HGT), but to have been specifically retained by G. candidum after the filamentous fungi-yeasts split concomitant with the yeasts' genome contraction. We refer to these genes as SRAGs (Specifically Retained Ancestral Genes), having been lost by all or nearly all other yeasts, and thus contributing to the phenotypic specificity of lineages. SRAG functions include lipases consistent with a role in cheese making and novel endoglucanases associated with degradation of plant material. Similar gene retention was observed in three other distantly related yeasts representative of this ecologically diverse subphylum. The phenomenon thus appears to be widespread in the Saccharomycotina and argues that, alongside neo-functionalization following gene duplication and HGT, specific gene retention must be recognized as an important mechanism for generation of biodiversity and adaptation in yeasts.

Molecular Polymorphism and Divergence of Duplicated Genes in Tetraploid African Clawed Frogs (Xenopus)

Genome duplication creates redundancy in proteins and their interaction networks, and subsequent smaller-scale gene duplication can further amplify genetic redundancy. Mutations then lead to the loss, maintenance or functional divergence of duplicated genes. Genome duplication occurred many times in African clawed frogs (genus Xenopus), and almost all extant species in this group evolved from a polyploid ancestor. To better understand the nature of selective constraints in a polyploid genome, we examined molecular polymorphism and divergence of duplicates and single-copy genes in 2 tetraploid African clawed frog species, Xenopus laevis and X. victorianus. We found that molecular polymorphism in the coding regions of putative duplicated genes was higher than in singletons, but not significantly so. Our findings also suggest that transcriptome evolution in polyploids is influenced by variation in the genome-wide mutation rate, and do not reject the hypothesis that gene dosage balance is also important.

Evolution of new functions de novo and from preexisting genes

How the enormous structural and functional diversity of new genes and proteins was generated (estimated to be 10(10)-10(12) different proteins in all organisms on earth [Choi I-G, Kim S-H. 2006. Evolution of protein structural classes and protein sequence families. Proc Natl Acad Sci 103: 14056-14061] is a central biological question that has a long and rich history. Extensive work during the last 80 years have shown that new genes that play important roles in lineage-specific phenotypes and adaptation can originate through a multitude of different mechanisms, including duplication, lateral gene transfer, gene fusion/fission, and de novo origination. In this review, we focus on two main processes as generators of new functions: evolution of new genes by duplication and divergence of pre-existing genes and de novo gene origination in which a whole protein-coding gene evolves from a noncoding sequence.

Lipids containing medium-chain fatty acids are specific to post-whole genome duplication Saccharomycotina yeasts

Background

Yeasts belonging to the subphylum Saccharomycotina have been used for centuries in food processing and, more recently, biotechnology. Over the past few decades, these yeasts have also been studied in the interest of their potential to produce oil to replace fossil resources. Developing yeasts for massive oil production requires increasing yield and modifying the profiles of the fatty acids contained in the oil to satisfy specific technical requirements. For example, derivatives of medium-chain fatty acids (MCFAs, containing 6–14 carbons) are used for the production of biodiesels, cleaning products, lubricants and cosmetics. Few studies are available in the literature on the production of MCFAs in yeasts.

Results

We analyzed the MCFA content in Saccharomyces cerevisiae grown in various conditions. The results revealed that MCFAs preferentially accumulated when cells were grown on synthetic media with a high C/N ratio at low temperature (23 °C). Upon screening deletion mutant strains for genes encoding lipid droplet-associated proteins, we found two genes, LOA1 and TGL3, involved in MCFA homeostasis. A phylogenetic analysis on 16 Saccharomycotina species showed that fatty acid profiles differed drastically among yeasts. Interestingly, MCFAs are only present in post-whole genome duplication yeast species.

Conclusions

In this study, we produced original data on fatty acid diversity in yeasts. We demonstrated that yeasts are amenable to genetic and metabolic engineering to increase their MCFA production. Furthermore, we revealed that yeast lipid biodiversity has not been fully explored, but that yeasts likely harbor as-yet-undiscovered strains or enzymes that can contribute to the production of high-value fatty acids for green chemistry.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0369-2) contains supplementary material, which is available to authorized users.

The origins of mutational robustness

Biological systems are resistant to genetic changes; a property known as mutational robustness, the origin of which remains an open question. In recent years, researchers have explored emergent properties of biological systems and mechanisms of genetic redundancy to reveal how mutational robustness emerges and persists. Several mechanisms have been proposed to explain the origin of mutational robustness, including molecular chaperones and gene duplication. The latter has received much attention, but its role in robustness remains controversial. Here, I examine recent findings linking genetic redundancy through gene duplication and mutational robustness. Experimental evolution and genome resequencing have made it possible to test the role of gene duplication in tolerating mutations at both the coding and regulatory levels. This evidence as well as previous findings on regulatory reprogramming of duplicates support the role of gene duplication in the origin of robustness.

Multispecies Analysis of Expression Pattern Diversification in the Recently Expanded Insect Ly6 Gene Family

Gene families often consist of members with diverse expression domains reflecting their functions in a wide variety of tissues. However, how the expression of individual members, and thus their tissue-specific functions, diversified during the course of gene family expansion is not well understood. In this study, we approached this question through the analysis of the duplication history and transcriptional evolution of a rapidly expanding subfamily of insect Ly6 genes. We analyzed different insect genomes and identified seven Ly6 genes that have originated from a single ancestor through sequential duplication within the higher Diptera. We then determined how the original embryonic expression pattern of the founding gene diversified by characterizing its tissue-specific expression in the beetle Tribolium castaneum, the butterfly Bicyclus anynana, and the mosquito Anopheles stephensi and those of its duplicates in three higher dipteran species, representing various stages of the duplication history (Megaselia abdita, Ceratitis capitata, and Drosophila melanogaster). Our results revealed that frequent neofunctionalization episodes contributed to the increased expression breadth of this subfamily and that these events occurred after duplication and speciation events at comparable frequencies. In addition, at each duplication node, we consistently found asymmetric expression divergence. One paralog inherited most of the tissue-specificities of the founder gene, whereas the other paralog evolved drastically reduced expression domains. Our approach attests to the power of combining a well-established duplication history with a comprehensive coverage of representative species in acquiring unequivocal information about the dynamics of gene expression evolution in gene families.

Detecting functional divergence after gene duplication through evolutionary changes in posttranslational regulatory sequences

Gene duplication is an important evolutionary mechanism that can result in functional divergence in paralogs due to neo-functionalization or sub-functionalization. Consistent with functional divergence after gene duplication, recent studies have shown accelerated evolution in retained paralogs. However, little is known in general about the impact of this accelerated evolution on the molecular functions of retained paralogs. For example, do new functions typically involve changes in enzymatic activities, or changes in protein regulation? Here we study the evolution of posttranslational regulation by examining the evolution of important regulatory sequences (short linear motifs) in retained duplicates created by the whole-genome duplication in budding yeast. To do so, we identified short linear motifs whose evolutionary constraint has relaxed after gene duplication with a likelihood-ratio test that can account for heterogeneity in the evolutionary process by using a non-central chi-squared null distribution. We find that short linear motifs are more likely to show changes in evolutionary constraints in retained duplicates compared to single-copy genes. We examine changes in constraints on known regulatory sequences and show that for the Rck1/Rck2, Fkh1/Fkh2, Ace2/Swi5 paralogs, they are associated with previously characterized differences in posttranslational regulation. Finally, we experimentally confirm our prediction that for the Ace2/Swi5 paralogs, Cbk1 regulated localization was lost along the lineage leading to SWI5 after gene duplication. Our analysis suggests that changes in posttranslational regulation mediated by short regulatory motifs systematically contribute to functional divergence after gene duplication.

How a protein is controlled is intimately linked to its function. Therefore, evolution can drive the functional divergence of proteins by tweaking their regulation, even if enzymatic capacities are preserved. Changes in posttranslational regulation (protein phosphorylation, degradation, subcellular localization, etc.) could therefore represent key mechanisms in functional divergence and lead to different phenotypic outcomes. Since disordered protein regions contain sites of protein modification and interaction (known as short linear motifs) and evolve rapidly relative to domains encoding enzymatic functions, these regions are good candidates to harbour sequence changes that underlie changes in function. In this study, we develop a statistical framework to identify changes in rate of evolution specific to protein regulatory sequences and identify hundreds of short linear motifs in disordered regions that are likely to have diverged after the whole-genome duplication in budding yeast. We show that these divergent motifs are much more frequent in paralogs than in single-copy proteins, and that they are more frequent in duplicate pairs that have functionally diverged. Our analysis suggests that changes in short linear motifs in disordered protein regions could be important molecular mechanisms of functional divergence after gene duplication.

Comparative genomics as a time machine: how relative gene dosage and metabolic requirements shaped the time-dependent resolution of yeast polyploidy

Using a phylogenetic model of evolution after genome duplication (i.e., polyploidy) and 12 yeast genomes with a shared genome duplication, I show that the loss of duplicate genes after that duplication occurred in three phases. First, losses that occurred immediately after the event were biased toward genes functioning in DNA repair and organellar functions. Then, the main group of duplicate losses appear to have been shaped by a requirement to maintain balance in protein levels: There is a strong statistical association between the number of protein interactions a gene's product is involved in and its propensity to have remained in duplicate. Moreover, when duplicated genes with interactions were lost, it was more common than expected for both members of an interaction pair to have been lost on the same branch of the phylogeny. Finally, in the third phase of the resolution process, overretention of duplicated enzymes carrying high flux and of duplicated genes involved in transcriptional regulation became dominant. I speculate that initial retention of such genes by a requirement to maintain gene dosage set the stage for the later functional changes that then maintained these duplicates for long periods.

Analysis of 41 plant genomes supports a wave of successful genome duplications in association with the Cretaceous-Paleogene boundary

Ancient whole-genome duplications (WGDs), also referred to as paleopolyploidizations, have been reported in most evolutionary lineages. Their attributed role remains a major topic of discussion, ranging from an evolutionary dead end to a road toward evolutionary success, with evidence supporting both fates. Previously, based on dating WGDs in a limited number of plant species, we found a clustering of angiosperm paleopolyploidizations around the Cretaceous-Paleogene (K-Pg) extinction event about 66 million years ago. Here we revisit this finding, which has proven controversial, by combining genome sequence information for many more plant lineages and using more sophisticated analyses. We include 38 full genome sequences and three transcriptome assemblies in a Bayesian evolutionary analysis framework that incorporates uncorrelated relaxed clock methods and fossil uncertainty. In accordance with earlier findings, we demonstrate a strongly nonrandom pattern of genome duplications over time with many WGDs clustering around the K-Pg boundary. We interpret these results in the context of recent studies on invasive polyploid plant species, and suggest that polyploid establishment is promoted during times of environmental stress. We argue that considering the evolutionary potential of polyploids in light of the environmental and ecological conditions present around the time of polyploidization could mitigate the stark contrast in the proposed evolutionary fates of polyploids.

Long-term asymmetrical acceleration of protein evolution after gene duplication

Rapid divergence of gene copies after duplication is thought to determine the fate of the copies and evolution of novel protein functions. However, data on how long the gene copies continue to experience an elevated rate of evolution remain scarce. Standard theory of gene duplications based on some level of genetic redundancy of gene copies predicts that the period of accelerated evolution must end relatively quickly. Using a maximum-likelihood approach we estimate preduplication, initial postduplication, and recent postduplication rates of evolution that occurred in the mammalian lineage. We find that both gene copies experience a similar in magnitude acceleration in their rate of evolution. The copy located in the original genomic position typically returns to the preduplication rates of evolution in a short period of time. The burst of faster evolution of the copy that is located in a new genomic position typically lasts longer. Furthermore, the fast-evolving copies on average continue to evolve faster than the preduplication rates far longer than predicted by standard theory of gene duplications. We hypothesize that the prolonged elevated rates of evolution are determined by functional properties that were acquired during, or soon after, the gene duplication event.

Dosage, duplication, and diploidization: clarifying the interplay of multiple models for duplicate gene evolution over time

Requirements to maintain dosage balance shape many genome-scale patterns in organisms, including the resolution of whole genome duplications (WGD), as well as the varied effects of aneuploidy, segmental duplications, tandem duplications, gene copy number variations (CNV), and epigenetic marks. Like neofunctionalization and subfunctionalization, the impact of absolute and relative dosage varies over time. These variations are of particular importance in understanding the role of dosage in the evolution of polyploid organisms. Numerous investigations have found the consequences of polyploidy remain distinct from small-scale duplications (SSD). This observation is significant as all flowering plants have experienced at least two ancient polyploid events, and many angiosperm lineages have undergone additional rounds of polyploidy. Intriguingly, recent studies indicate a link between how epigenetic marks in recent allopolyploids may induce immediate changes in gene expression and the longer-term patterns of biased fractionation and chromosomal evolution. We argue that dosage effects represent one aspect of an emerging pluralistic framework, a framework that will use biophysics, genomic technologies, and systems-level models of cells to broaden our view of how genomes evolve.

Insights into molecular evolution from yeast genomics

Enabled by comparative genomics, yeasts have increasingly developed into a powerful model system for molecular evolution. Here we survey several areas in which yeast studies have made important contributions, including regulatory evolution, gene duplication and divergence, evolution of gene order and evolution of complexity. In each area we highlight key studies and findings based on techniques ranging from statistical analysis of large datasets to direct laboratory measurements of fitness. Future work will combine traditional evolutionary genetics analysis and experimental evolution with tools from systems biology to yield mechanistic insight into complex phenotypes.

Consequences of Whole-Genome Triplication as Revealed by Comparative Genomic Analyses of the Wild Radish Raphanus raphanistrum and Three Other Brassicaceae Species

Polyploidization events are frequent among flowering plants, and the duplicate genes produced via such events contribute significantly to plant evolution. We sequenced the genome of wild radish (Raphanus raphanistrum), a Brassicaceae species that experienced a whole-genome triplication event prior to diverging from Brassica rapa. Despite substantial gene gains in these two species compared with Arabidopsis thaliana and Arabidopsis lyrata, ∼70% of the orthologous groups experienced gene losses in R. raphanistrum and B. rapa, with most of the losses occurring prior to their divergence. The retained duplicates show substantial divergence in sequence and expression. Based on comparison of A. thaliana and R. raphanistrum ortholog floral expression levels, retained radish duplicates diverged primarily via maintenance of ancestral expression level in one copy and reduction of expression level in others. In addition, retained duplicates differed significantly from genes that reverted to singleton state in function, sequence composition, expression patterns, network connectivity, and rates of evolution. Using these properties, we established a statistical learning model for predicting whether a duplicate would be retained postpolyploidization. Overall, our study provides new insights into the processes of plant duplicate loss, retention, and functional divergence and highlights the need for further understanding factors controlling duplicate gene fate.

Extensive rewiring and complex evolutionary dynamics in a C. elegans multiparameter transcription factor network

Gene duplication results in two identical paralogs that diverge through mutation, leading to loss or gain of interactions with other biomolecules. Here, we comprehensively characterize such network rewiring for C. elegans transcription factors (TFs) within and across four newly delineated molecular networks. Remarkably, we find that even highly similar TFs often have different interaction degrees and partners. In addition, we find that most TF families have a member that is highly connected in multiple networks. Further, different TF families have opposing correlations between network connectivity and phylogenetic age, suggesting that they are subject to different evolutionary pressures. Finally, TFs that have similar partners in one network generally do not in another, indicating a lack of pressure to retain cross-network similarity. Our multiparameter analyses provide unique insights into the evolutionary dynamics that shaped TF networks.

Accelerated evolution after gene duplication: a time-dependent process affecting just one copy

Gene duplication is widely regarded as a major mechanism modeling genome evolution and function. However, the mechanisms that drive the evolution of the two, initially redundant, gene copies are still ill defined. Many gene duplicates experience evolutionary rate acceleration, but the relative contribution of positive selection and random drift to the retention and subsequent evolution of gene duplicates, and for how long the molecular clock may be distorted by these processes, remains unclear. Focusing on rodent genes that duplicated before and after the mouse and rat split, we find significantly increased sequence divergence after duplication in only one of the copies, which in nearly all cases corresponds to the novel daughter copy, independent of the mechanism of duplication. We observe that the evolutionary rate of the accelerated copy, measured as the ratio of nonsynonymous to synonymous substitutions, is on average 5-fold higher in the period spanning 4-12 My after the duplication than it was before the duplication. This increase can be explained, at least in part, by the action of positive selection according to the results of the maximum likelihood-based branch-site test. Subsequently, the rate decelerates until purifying selection completely returns to preduplication levels. Reversion to the original rates has already been accomplished 40.5 My after the duplication event, corresponding to a genetic distance of about 0.28 synonymous substitutions per site. Differences in tissue gene expression patterns parallel those of substitution rates, reinforcing the role of neofunctionalization in explaining the evolution of young gene duplicates.