Search

Transcriptional divergence of duplicated genes after whole genome duplication (WGD) has been described in many plant lineages and is often associated with subgenome dominance, a genome-wide mechanism. However, it is unknown what underlies the transcriptional divergence of duplicated genes in polyploid species that lack subgenome dominance. Soybean is a paleotetraploid with a WGD that occurred 5 to 13 Mya. Approximately 50% of the duplicated genes retained from this WGD exhibit transcriptional divergence. We developed accessible chromatin region (ACR) datasets from leaf, flower, and seed tissues using MNase-hypersensitivity sequencing. We validated enhancer function of several ACRs associated with known genes using CRISPR/Cas9-mediated genome editing. The ACR datasets were used to examine and correlate the transcriptional patterns of 17,111 pairs of duplicated genes in different tissues. We demonstrate that ACR dynamics are correlated with divergence of both expression level and tissue specificity of individual gene pairs. Gain or loss of flanking ACRs and mutation of cis -regulatory elements (CREs) within the ACRs can change the balance of the expression level and/or tissue specificity of the duplicated genes. Analysis of DNA sequences associated with ACRs revealed that the extensive sequence rearrangement after the WGD reshaped the CRE landscape, which appears to play a key role in the transcriptional divergence of duplicated genes in soybean. This may represent a general mechanism for transcriptional divergence of duplicated genes in polyploids that lack subgenome dominance.

The evolutionary mechanisms underlying duplicate gene maintenance and divergence remain highly debated. Epigenetic modifications, such as DNA methylation, may contribute to duplicate gene evolution by facilitating tissue-specific regulation. However, the role of epigenetic divergence on duplicate gene evolution remains little understood. Here we show, using comprehensive data across 10 diverse human tissues, that DNA methylation plays critical roles in several aspects of duplicate gene evolution. We first demonstrate that duplicate genes are initially heavily methylated, before gradually losing DNA methylation as they age. Within each pair, DNA methylation divergence between duplicate partners increases with evolutionary age. Importantly, tissue-specific DNA methylation of duplicates correlates with tissue-specific expression, implicating DNA methylation as a causative factor for functional divergence of duplicate genes. These patterns are apparent in promoters but not in gene bodies, in accord with the complex relationship between gene-body DNA methylation and transcription. Remarkably, many duplicate gene pairs exhibit consistent division of DNA methylation across multiple, divergent tissues: For the majority (73%) of duplicate gene pairs, one partner is always hypermethylated compared with the other. This is indicative of a common underlying determinant of DNA methylation. The division of DNA methylation is also consistent with their chromatin accessibility profiles. Moreover, at least two sequence motifs known to interact with the Sp1 transcription factor mark promoters of more hypomethylated duplicate partners. These results demonstrate critical roles of DNA methylation, as well as complex interaction between genome and epigenome, on duplicate gene evolution.

Gene duplication is a key source of genetic innovation that plays a role in the evolution of phenotypic complexity. Although several evolutionary processes can result in the long-term retention of duplicate genes, their relative contributions in nature are unknown. Here we develop a phylogenetic approach for comparing genome-wide expression profiles of closely related species to quantify the roles of conservation, neofunctionalization, subfunctionalization, and specialization in the preservation of duplicate genes. Application of our method to pairs of young duplicates in Drosophila shows that neofunctionalization, the gain of a novel function in one copy, accounts for the retention of almost two-thirds of duplicate genes. Surprisingly, novel functions nearly always originate in younger (child) copies, whereas older (parent) copies possess functions similar to those of ancestral genes. Further examination of such pairs reveals a strong bias toward RNA-mediated duplication events, implicating asymmetric duplication and positive selection in the evolution of new functions. Moreover, we show that young duplicate genes are expressed primarily in testes and that their expression breadth increases over evolutionary time. This finding supports the "out-of-testes" hypothesis, which posits that testes are a catalyst for the emergence of new genes that ultimately evolve functions in other tissues. Thus, our study highlights the importance of neofunctionalization and positive selection in the retention of young duplicates in Drosophila and illustrates how duplicates become incorporated into novel functional networks over evolutionary time.

We report the genome sequence of melon, an important horticultural crop worldwide. We assembled 375 Mb of the double-haploid line DHL92, representing 83.3% of the estimated melon genome. We predicted 27,427 protein-coding genes, which we analyzed by reconstructing 22,218 phylogenetic trees, allowing mapping of the orthology and paralogy relationships of sequenced plant genomes. We observed the absence of recent whole-genome duplications in the melon lineage since the ancient eudicot triplication, and our data suggest that transposon amplification may in part explain the increased size of the melon genome compared with the close relative cucumber. A low number of nucleotide-binding site-leucine-rich repeat disease resistance genes were annotated, suggesting the existence of specific defense mechanisms in this species. The DHL92 genome was compared with that of its parental lines allowing the quantification of sequence variability in the species. The use of the genome sequence in future investigations will facilitate the understanding of evolution of cucurbits and the improvement of breeding strategies.

Thiol-dependent reductase I (TDR1), an enzyme found in parasitic Leishmania species and Trypanosoma cruzi, is implicated in deglutathionylation and activation of antimonial prodrugs used to treat leishmaniasis. The 2.3 Å resolution structure of TDR1 reveals a unique trimer of subunits each containing two glutathione-S-transferase (GST) domains. The similarities of individual domains and comparisons with GST classes suggest that TDR1 evolved by gene duplication, diversification, and gene fusion; a combination of events previously unknown in the GST protein superfamily and potentially explaining the distinctive enzyme properties of TDR1. The deglutathionylation activity of TDR1 implies that glutathione itself has regulatory intracellular roles in addition to being a precursor for trypanothione, the major low mass thiol present in trypanosomatids. We propose that activation of antiparasite Sb(V)-drugs is a legacy of the deglutathionylation activity of TDR1 and involves processing glutathione adducts with concomitant reduction of the metalloid to active Sb(III) species.

Division of labor among functionally specialized modules occurs at all levels of biological organization in both animals and plants. Well-known examples include the evolution of specialized enzymes after gene duplication, the evolution of specialized cell types, limb diversification in arthropods, and the evolution of specialized colony members in many taxa of marine invertebrates and social insects. Here, we identify conditions favoring the evolution of division of labor by means of a general mathematical model. Our starting point is the assumption that modules contribute to two different biological tasks and that the potential of modules to contribute to these tasks is traded off. Our results are phrased in terms of properties of performance functions that map the phenotype of modules to measures of performance. We show that division of labor is favored by three factors: positional effects that predispose modules for one of the tasks, accelerating performance functions, and synergistic interactions between modules. If modules can be lost or damaged, selection for robustness can counteract selection for functional specialization. To illustrate our theory we apply it to the evolution of specialized enzymes coded by duplicated genes.

Gene duplication plays key roles in organismal evolution. Duplicate genes, if they survive, tend to diverge in regulatory and coding regions. Divergences in coding regions, especially those that can change the function of the gene, can be caused by amino acid-altering substitutions and/or alterations in exon-intron structure. Much has been learned about the mode, tempo, and consequences of nucleotide substitutions, yet relatively little is known about structural divergences. In this study, by analyzing 612 pairs of sibling paralogs from seven representative gene families and 300 pairs of one-to-one orthologs from different species, we investigated the occurrence and relative importance of structural divergences during the evolution of duplicate and nonduplicate genes. We found that structural divergences have been very prevalent in duplicate genes and, in many cases, have led to the generation of functionally distinct paralogs. Comparisons of the genomic sequences of these genes further indicated that the differences in exon-intron structure were actually accomplished by three main types of mechanisms (exon/intron gain/loss, exonization/pseudoexonization, and insertion/deletion), each of which contributed differently to structural divergence. Like nucleotide substitutions, insertion/deletion and exonization/pseudoexonization occurred more or less randomly, with the number of observable mutational events per gene pair being largely proportional to evolutionary time. Notably, however, compared with paralogs with similar evolutionary times, orthologs have accumulated significantly fewer structural changes, whereas the amounts of amino acid replacements accumulated did not show clear differences. This finding suggests that structural divergences have played a more important role during the evolution of duplicate than nonduplicate genes.

Polyploidy has been a common process during the evolution of eukaryotes, especially plants, leading to speciation and the evolution of new gene functions. Gene expression levels and patterns can change, and gene silencing can occur in allopolyploids--phenomena sometimes referred to as "transcriptome shock." Alternative splicing (AS) creates multiple mature mRNAs from a single type of precursor mRNA. Here we examined the evolution of AS patterns after polyploidy, with natural and two resynthesized allotetraploid Brassica napus lines, using RT-PCR and sequencing assays of 82 AS events in duplicated gene pairs (homeologs). Comparing the AS patterns between the two homeologs in natural B. napus revealed that many of the gene pairs show different AS patterns, with a few showing variation that was organ specific or induced by abiotic stress treatments. In the resynthesized allotetraploids, 26-30% of the duplicated genes showed changes in AS compared with the parents, including many cases of AS event loss after polyploidy. Parallel losses of many AS events after allopolyploidy were detected in the two independently resynthesized lines. More changes occurred in parallel between the two lines than changes specific to each line. The PASTICCINO gene showed partitioning of two AS events between the two homeologs in the resynthesized allopolyploids. AS changes after allopolyploidy were much more common than homeolog silencing. Our findings indicate that AS patterns can change rapidly after polyploidy, that many genes are affected, and that AS changes are an important component of the transcriptome shock experienced by new allopolyploids.

Aldehyde dehydrogenases (ALDHs) catabolize toxic aldehydes and process the vitamin A-derived retinaldehyde into retinoic acid (RA), a small diffusible molecule and a pivotal chordate morphogen. In this study, we combine phylogenetic, structural, genomic, and developmental gene expression analyses to examine the evolutionary origins of ALDH substrate preference. Structural modeling reveals that processing of small aldehydes, such as acetaldehyde, by ALDH2, versus large aldehydes, including retinaldehyde, by ALDH1A is associated with small versus large substrate entry channels (SECs), respectively. Moreover, we show that metazoan ALDH1s and ALDH2s are members of a single ALDH1/2 clade and that during evolution, eukaryote ALDH1/2s often switched between large and small SECs after gene duplication, transforming constricted channels into wide opened ones and vice versa. Ancestral sequence reconstructions suggest that during the evolutionary emergence of RA signaling, the ancestral, narrow-channeled metazoan ALDH1/2 gave rise to large ALDH1 channels capable of accommodating bulky aldehydes, such as retinaldehyde, supporting the view that retinoid-dependent signaling arose from ancestral cellular detoxification mechanisms. Our analyses also indicate that, on a more restricted evolutionary scale, ALDH1 duplicates from invertebrate chordates (amphioxus and ascidian tunicates) underwent switches to smaller and narrower SECs. When combined with alterations in gene expression, these switches led to neofunctionalization from ALDH1-like roles in embryonic patterning to systemic, ALDH2-like roles, suggesting functional shifts from signaling to detoxification.

Genetic incompatibility is a barrier contributing to species isolation and is caused by genetic interactions. We made a whole genome survey of two-way interacting loci acting within the gametophyte or zygote using independence tests of marker segregations in an F(2) population from an intersubspecific cross between O. sativa subspecies indica and japonica. We detected only one reproducible interaction, and identified paralogous hybrid incompatibility genes, DOPPELGANGER1 (DPL1) and DOPPELGANGER2 (DPL2), by positional cloning. Independent disruptions of DPL1 and DPL2 occurred in indica and japonica, respectively. DPLs encode highly conserved, plant-specific small proteins (∼10 kDa) and are highly expressed in mature anther. Pollen carrying two defective DPL alleles became nonfunctional and did not germinate, suggesting an essential role for DPLs in pollen germination. Although rice has many duplicated genes resulting from ancient whole genome duplication, the origin of this gene duplication was in recent small-scale gene duplication, occurring after Oryza-Brachypodium differentiation. Comparative analyses suggested the geographic and phylogenetic distribution of these two defective alleles, showing that loss-of-function mutations of DPL1 genes emerged multiple times in indica and its wild ancestor, O. rufipogon, and that the DPL2 gene defect is specific to japonica cultivars.

About 30% of protein-coding genes in the human genome are related through two whole genome duplication (WGD) events. Although WGD is often credited with great evolutionary importance, the processes governing the retention of these genes and their biological significance remain unclear. One increasingly popular hypothesis is that dosage balance constraints are a major determinant of duplicate gene retention. We test this hypothesis and show that WGD-duplicated genes (ohnologs) have rarely experienced subsequent small-scale duplication (SSD) and are also refractory to copy number variation (CNV) in human populations and are thus likely to be sensitive to relative quantities (i.e., they are dosage-balanced). By contrast, genes that have experienced SSD in the vertebrate lineage are more likely to also display CNV. This supports the hypothesis of biased retention of dosage-balanced genes after WGD. We also show that ohnologs have a strong association with human disease. In particular, Down Syndrome (DS) caused by trisomy 21 is widely assumed to be caused by dosage effects, and 75% of previously reported candidate genes for this syndrome are ohnologs that experienced no other copy number changes. We propose the remaining dosage-balanced ohnologs on chromosome 21 as candidate DS genes. These observations clearly show a persistent resistance to dose changes in genes duplicated by WGD. Dosage balance constraints simultaneously explain duplicate gene retention and essentiality after WGD.

Gene and genome duplications create novel genetic material on which evolution can work and have therefore been recognized as a major source of innovation for many eukaryotic lineages. Following duplication, the most likely fate is gene loss; however, a considerable fraction of duplicated genes survive. Not all genes have the same probability of survival, but it is not fully understood what evolutionary forces determine the pattern of gene retention. Here, we use genome sequence data as well as large-scale phosphoproteomics data from the baker's yeast Saccharomyces cerevisiae, which underwent a whole-genome duplication approximately 100 mya, and show that the number of phosphorylation sites on the proteins they encode is a major determinant of gene retention. Protein phosphorylation motifs are short amino acid sequences that are usually embedded within unstructured and rapidly evolving protein regions. Reciprocal loss of those ancestral sites and the gain of new ones are major drivers in the retention of the two surviving duplicates and in their acquisition of distinct functions. This way, small changes in the sequences of unstructured regions in proteins can contribute to the rapid rewiring and adaptation of regulatory networks.

The question of how duplicate genes are retained in a population remains controversial. The duplication-degeneration-complementation model, which involves no positive selection, stipulates a higher retention rate of duplicate genes in a small population than in a large one. This model has been accepted by many evolutionists. However, we found considerably more retentions and fewer losses of duplicate genes in the mouse genome than in the human genome, although the population size of rodents is in general larger than that of primates. Indeed, in nearly every interval of synonymous divergence between duplicate genes, the number of gene retentions in mouse is larger than that in human. Our findings suggest a more important role of positive selection in duplicate retention than duplication-degeneration-complementation. In addition, certain functional categories show a higher tendency of lineage-specific expansion than expected, suggesting lineage-specific selection or functional bias in retained duplicates.

Much recent work has explored molecular and population-genetic constraints on the rate of protein sequence evolution. The best predictor of evolutionary rate is expression level, for reasons that have remained unexplained. Here, we hypothesize that selection to reduce the burden of protein misfolding will favor protein sequences with increased robustness to translational missense errors. Pressure for translational robustness increases with expression level and constrains sequence evolution. Using several sequenced yeast genomes, global expression and protein abundance data, and sets of paralogs traceable to an ancient whole-genome duplication in yeast, we rule out several confounding effects and show that expression level explains roughly half the variation in Saccharomyces cerevisiae protein evolutionary rates. We examine causes for expression's dominant role and find that genome-wide tests favor the translational robustness explanation over existing hypotheses that invoke constraints on function or translational efficiency. Our results suggest that proteins evolve at rates largely unrelated to their functions and can explain why highly expressed proteins evolve slowly across the tree of life.

Homologous proteins occurring through gene duplication may give rise to novel functions through mutations affecting protein sequence or expression. Comparison of such homologues allows insight into how morphological traits evolve. However, it is often unclear which changes are key to determining new functions. To address these ideas, we have studied a system where two homologues have evolved clear and opposite functions in controlling a major developmental switch. In plants, flowering is a major developmental transition that is critical to reproductive success. Arabidopsis phosphatidylethanolamine-binding protein homologues TERMINAL FLOWER 1 (TFL1) and FLOWERING LOCUS T (FT) are key controllers of flowering, determining when and where flowers are made, but as opposing functions: TFL1 is a repressor, FT is an activator. We have uncovered a striking molecular basis for how these homologous proteins have diverged. Although <60% identical, we have shown that swapping a single amino acid is sufficient to convert TFL1 to FT function and vice versa. Therefore, these key residues may have strongly contributed to the selection of these important functions over plant evolution. Further, our results suggest that TFL1 and FT are highly conserved in biochemical function and that they act as repressors or activators of flowering through discrimination of structurally related interactors by a single residue.

Recent analysis of complete eukaryotic genome sequences has revealed that gene duplication has been rampant. Moreover, next to a continuous mode of gene duplication, in many eukaryotic organisms the complete genome has been duplicated in their evolutionary past. Such large-scale gene duplication events have been associated with important evolutionary transitions or major leaps in development and adaptive radiations of species. Here, we present an evolutionary model that simulates the duplication dynamics of genes, considering genome-wide duplication events and a continuous mode of gene duplication. Modeling the evolution of the different functional categories of genes assesses the importance of different duplication events for gene families involved in specific functions or processes. By applying our model to the Arabidopsis genome, for which there is compelling evidence for three whole-genome duplications, we show that gene loss is strikingly different for large-scale and small-scale duplication events and highly biased toward certain functional classes. We provide evidence that some categories of genes were almost exclusively expanded through large-scale gene duplication events. In particular, we show that the three whole-genome duplications in Arabidopsis have been directly responsible for >90% of the increase in transcription factors, signal transducers, and developmental genes in the last 350 million years. Our evolutionary model is widely applicable and can be used to evaluate different assumptions regarding small- or large-scale gene duplication events in eukaryotic genomes.

Mobile element dynamics in seven alleles of the chalcone synthase D locus (CHS-D) of the common morning glory (Ipomoea purpurea) are analyzed in the context of synonymous nucleotide sequence distances for CHS-D exons. By using a nucleotide sequence of CHS-D from the sister species Ipomoea nil (Japanese morning glory [Johzuka-Hisatomi, Y., Hoshino, A., Mori, T., Habu, Y. & Iida, S. (1999) Genes Genet. Syst. 74, 141-147], it is also possible to determine the relative frequency of insertion and loss of elements within the CHS-D locus between these two species. At least four different types of transposable elements exist upstream of the coding region, or within the single intron of the CHS-D locus in I. purpurea. There are three distinct families of miniature inverted-repeat transposable elements (MITES), and some recent transpositions of Activator/Dissociation (Ac/Ds)-like elements (Tip100), of some short interspersed repetitive elements (SINEs), and of an insertion sequence (InsIpCHSD) found in the neighborhood of this locus. The data provide no compelling evidence of the transposition of the mites since the separation of I. nil and I. purpurea roughly 8 million years ago. Finally, it is shown that the number and frequency of mobile elements are highly heterogeneous among different duplicate CHS loci, suggesting that the dynamics observed at CHS-D are locus-specific.