Rapid evolution of expression and regulatory divergences after yeast gene duplication

A Rigorous Framework to Classify the Postduplication Fate of Paralogous Genes

Gene duplication has a central role in evolution; still, little is known on the fates of the duplicated copies, their relative frequency, and on how environmental conditions affect them. Moreover, the lack of rigorous definitions concerning the fate of duplicated genes hinders the development of a global vision of this process. In this paper, we present a new framework aiming at characterizing and formally differentiating the fate of duplicated genes. Our framework has been tested via simulations, where the evolution of populations has been simulated using aevol, an in silico experimental evolution platform. Our results show several patterns that confirm some of the conclusions from previous studies, while also exhibiting new tendencies; this may open up new avenues to better understand the role of duplications as a driver of evolution.

The evolutionary history of Brachyury genes in Hydrozoa involves duplications, divergence, and neofunctionalization

Brachyury, a member of T-box gene family, is widely known for its major role in mesoderm specification in bilaterians. It is also present in non-bilaterian metazoans, such as cnidarians, where it acts as a component of an axial patterning system. In this study, we present a phylogenetic analysis of Brachyury genes within phylum Cnidaria, investigate differential expression and address a functional framework of Brachyury paralogs in hydrozoan Dynamena pumila. Our analysis indicates two duplication events of Brachyury within the cnidarian lineage. The first duplication likely appeared in the medusozoan ancestor, resulting in two copies in medusozoans, while the second duplication arose in the hydrozoan ancestor, resulting in three copies in hydrozoans. Brachyury1 and 2 display a conservative expression pattern marking the oral pole of the body axis in D. pumila. On the contrary, Brachyury3 expression was detected in scattered presumably nerve cells of the D. pumila larva. Pharmacological modulations indicated that Brachyury3 is not under regulation of cWnt signaling in contrast to the other two Brachyury genes. Divergence in expression patterns and regulation suggest neofunctionalization of Brachyury3 in hydrozoans.

Evolutionary conservation of sequence motifs at sites of protein modification

Gene duplications are common in biology and are likely to be an important source of functional diversification and specialization. The yeast Saccharomyces cerevisiae underwent a whole genome duplication event early in evolution, and a substantial number of duplicated genes have been retained. We identified more than 3,500 instances where only one of two paralogous proteins undergoes post-translational modification despite having retained the same amino acid residue in both. We also developed a web-based search algorithm (CoSMoS.c.) that scores conservation of amino acid sequences based on 1011 wild and domesticated yeast isolates and used it to compare differentially-modified pairs of paralogous proteins. We found that the most common modifications - phosphorylation, ubiquitylation and acylation but not N-glycosylation - occur in regions of high sequence conservation. Such conservation is evident even for ubiquitylation and succinylation, where there is no established 'consensus site' for modification. Differences in phosphorylation were not associated with predicted secondary structure or solvent accessibility, but did mirror known differences in kinase-substrate interactions. Thus, differences in post-translational modification likely result from differences in adjoining amino acids and their interactions with modifying enzymes. By integrating data from large scale proteomics and genomics analysis, in a system with such substantial genetic diversity, we obtained a more comprehensive understanding of the functional basis for genetic redundancies that have persisted for 100 million years.

Placing human gene families into their evolutionary context

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

The Heterogeneity in the Landscape of Gene Dominance in Maize is Accompanied by Unique Chromatin Environments

Subgenome dominance after whole-genome duplication (WGD) has been observed in many plant species. However, the degree to which the chromatin environment affects this bias has not been explored. Here, we compared the dominant subgenome (maize1) and the recessive subgenome (maize2) with respect to patterns of sequence substitutions, genes expression, transposable element accumulation, small interfering RNAs, DNA methylation, histone modifications, and accessible chromatin regions (ACRs). Our data show that the degree of bias between subgenomes for all the measured variables does not vary significantly when both of the WGD genes are located in pericentromeric regions. Our data further indicate that the location of maize1 genes in chromosomal arms is pivotal for maize1 to maintain its dominance, but location has a less effect on maize2 homoeologs. In addition to homoeologous genes, we compared ACRs, which often harbor cis-regulatory elements, between the two subgenomes and demonstrate that maize1 ACRs have a higher level of chromatin accessibility, a lower level of sequence substitution, and are enriched in chromosomal arms. Furthermore, we find that a loss of maize1 ACRs near their nearby genes is associated with a reduction in purifying selection and expression of maize1 genes relative to their maize2 homoeologs. Taken together, our data suggest that chromatin environment and cis-regulatory elements are important determinants shaping the divergence and evolution of duplicated genes.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom forming fungi (Agaricomycetes)

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a ‘developmental hourglass,’ act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Multi-Omics Analysis of Multiple Glucose-Sensing Receptor Systems in Yeast

The yeast Saccharomyces cerevisiae has long been used to produce alcohol from glucose and other sugars. While much is known about glucose metabolism, relatively little is known about the receptors and signaling pathways that indicate glucose availability. Here, we compare the two glucose receptor systems in S. cerevisiae. The first is a heterodimer of transporter-like proteins (transceptors), while the second is a seven-transmembrane receptor coupled to a large G protein (Gpa2) that acts in coordination with two small G proteins (Ras1 and Ras2). Through comprehensive measurements of glucose-dependent transcription and metabolism, we demonstrate that the two receptor systems have distinct roles in glucose signaling: the G-protein-coupled receptor directs carbohydrate and energy metabolism, while the transceptors regulate ancillary processes such as ribosome, amino acids, cofactor and vitamin metabolism. The large G-protein transmits the signal from its cognate receptor, while the small G-protein Ras2 (but not Ras1) integrates responses from both receptor pathways. Collectively, our analysis reveals the molecular basis for glucose detection and the earliest events of glucose-dependent signal transduction in yeast.

Retention of duplicated genes in evolution

Gene duplication is a prevalent phenomenon across the tree of life. The processes that lead to the retention of duplicated genes are not well understood. Functional genomics approaches in model organisms, such as yeast, provide useful tools to test the mechanisms underlying retention with functional redundancy and divergence of duplicated genes, including fates associated with neofunctionalization, subfunctionalization, back-up compensation, and dosage amplification. Duplicated genes may also be retained as a consequence of structural and functional entanglement. Advances in human gene editing have enabled the interrogation of duplicated genes in the human genome, providing new tools to evaluate the relative contributions of each of these factors to duplicate gene retention and the evolution of genome structure.

Lineage-Specific Genes and Family Expansions in Dictyostelid Genomes Display Expression Bias and Evolutionary Diversification during Development

Gene duplications generate new genes that can contribute to expression changes and the evolution of new functions. Genomes often consist of gene families that undergo expansions, some of which occur in specific lineages that reflect recent adaptive diversification. In this study, lineage-specific genes and gene family expansions were studied across five dictyostelid species to determine when and how they are expressed during multicellular development. Lineage-specific genes were found to be enriched among genes with biased expression (predominant expression in one developmental stage) in each species and at most developmental time points, suggesting independent functional innovations of new genes throughout the phylogeny. Biased duplicate genes had greater expression divergence than their orthologs and paralogs, consistent with subfunctionalization or neofunctionalization. Lineage-specific expansions in particular had biased genes with both molecular signals of positive selection and high expression, suggesting adaptive genetic and transcriptional diversification following duplication. Our results present insights into the potential contributions of lineage-specific genes and families in generating species-specific phenotypes during multicellular development in dictyostelids.

Horizontally acquired cysteine synthase genes undergo functional divergence in lepidopteran herbivores

Horizontal gene transfer (HGT) plays an important role in evolutionary processes as organisms adapt to their environments, and now cases of gene duplication after HGT in eukaryotes are emerging at an increasing rate. However, the fate and roles of the duplicated genes over time in eukaryotes remain unclear. Here we conducted a comprehensive analysis of the evolution of cysteine synthase (CYS) in lepidopteran insects. Our results indicate that HGT-derived CYS genes are widespread and have undergone duplication following horizontal transfer in many lepidopteran insects. Moreover, lepidopteran CYS proteins not only have β-cyanoalanine synthase activity but also possess cysteine synthase activity that is involved in sulfur amino acid biosynthesis. Duplicated CYS genes show marked divergence in gene expression patterns and enzymatic properties, suggesting that they probably have undergone subfunctionalization and/or neofunctionalization in Lepidoptera. The gene transfer of CYS genes and subsequent duplication appears to have facilitated the adaptation of lepidopteran insects to different diets and promoted their ecological diversification. Overall, this study provides valuable insights into the ecological and evolutionary contributions of CYS in lepidopteran insects.

Innovation, conservation, and repurposing of gene function in root cell type development

Plant species have evolved myriads of solutions, including complex cell type development and regulation, to adapt to dynamic environments. To understand this cellular diversity, we profiled tomato root cell type translatomes. Using xylem differentiation in tomato, examples of functional innovation, repurposing, and conservation of transcription factors are described, relative to the model plant Arabidopsis. Repurposing and innovation of genes are further observed within an exodermis regulatory network and illustrate its function. Comparative translatome analyses of rice, tomato, and Arabidopsis cell populations suggest increased expression conservation of root meristems compared with other homologous populations. In addition, the functions of constitutively expressed genes are more conserved than those of cell type/tissue-enriched genes. These observations suggest that higher order properties of cell type and pan-cell type regulation are evolutionarily conserved between plants and animals.

Diverse Molecular Mechanisms Contribute to Differential Expression of Human Duplicated Genes

Emerging evidence links genes within human-specific segmental duplications (HSDs) to traits and diseases unique to our species. Strikingly, despite being nearly identical by sequence (>98.5%), paralogous HSD genes are differentially expressed across human cell and tissue types, though the underlying mechanisms have not been examined. We compared cross-tissue mRNA levels of 75 HSD genes from 30 families between humans and chimpanzees and found expression patterns consistent with relaxed selection on or neofunctionalization of derived paralogs. In general, ancestral paralogs exhibited greatest expression conservation with chimpanzee orthologs, though exceptions suggest certain derived paralogs may retain or supplant ancestral functions. Concordantly, analysis of long-read isoform sequencing data sets from diverse human tissues and cell lines found that about half of derived paralogs exhibited globally lower expression. To understand mechanisms underlying these differences, we leveraged data from human lymphoblastoid cell lines (LCLs) and found no relationship between paralogous expression divergence and post-transcriptional regulation, sequence divergence, or copy-number variation. Considering cis-regulation, we reanalyzed ENCODE data and recovered hundreds of previously unidentified candidate CREs in HSDs. We also generated large-insert ChIP-sequencing data for active chromatin features in an LCL to better distinguish paralogous regions. Some duplicated CREs were sufficient to drive differential reporter activity, suggesting they may contribute to divergent cis-regulation of paralogous genes. This work provides evidence that cis-regulatory divergence contributes to novel expression patterns of recent gene duplicates in humans.

Genome-Wide Identification and Analysis of Chitinase GH18 Gene Family in Mycogone perniciosa

Mycogone perniciosa causes wet bubble disease in Agaricus bisporus and various Agaricomycetes species. In a previous work, we identified 41 GH18 chitinase genes and other pathogenicity-related genes in the genome of M. perniciosa Hp10. Chitinases are enzymes that degrade chitin, and they have diverse functions in nutrition, morphogenesis, and pathogenesis. However, these important genes in M. perniciosa have not been fully characterized, and their functions remain unclear. Here, we performed a genome-wide analysis of M. perniciosa GH18 genes and analyzed the transcriptome profiles and GH18 expression patterns in M. perniciosa during the time course of infection in A. bisporus. Phylogenetic analysis of the 41 GH18 genes with those of 15 other species showed that the genes were clustered into three groups and eight subgroups based on their conserved domains. The GH18 genes clustered in the same group shared different gene structures but had the same protein motifs. All GH18 genes were localized in different organelles, were unevenly distributed on 11 contigs, and had orthologs in the other 13 species. Twelve duplication events were identified, and these had undergone both positive and purifying selection. The transcriptome analyses revealed that numerous genes, including transporters, cell wall degrading enzymes (CWDEs), cytochrome P450, pathogenicity-related genes, secondary metabolites, and transcription factors, were significantly upregulated at different stages of M. perniciosa Hp10 infection of A. bisporus. Twenty-three out of the 41 GH18 genes were differentially expressed. The expression patterns of the 23 GH18 genes were different and were significantly expressed from 3 days post-inoculation of M. perniciosa Hp10 in A. bisporus. Five differentially expressed GH18 genes were selected for RT-PCR and gene cloning to verify RNA-seq data accuracy. The results showed that those genes were successively expressed in different infection stages, consistent with the previous sequencing results. Our study provides a comprehensive analysis of pathogenicity-related and GH18 chitinase genes’ influence on M. perniciosa mycoparasitism of A. bisporus. Our findings may serve as a basis for further studies of M. perniciosa mycoparasitism, and the results have potential value for improving resistance in A. bisporus and developing efficient disease-management strategies to mitigate wet bubble disease.

The Effects of Gene Duplication Modes on the Evolution of Regulatory Divergence in Wild and Cultivated Soybean

Regulatory changes include divergence in both cis-elements and trans-factors, which play roles in organismal evolution. Whole genome duplications (WGD) followed by diploidization are a recurrent feature in the evolutionary history of angiosperms. Prior studies have shown that duplicated genes have different evolutionary fates due to variable selection constraints and results in genomic compositions with hallmarks of paleopolyploidy. The recent sequential WGDs and post-WGD evolution in the common ancestor of cultivated soybean (Glycine max) and wild soybean (Glycine soja), together with other models of gene duplication, have resulted in a highly duplicated genome. In this study, we investigated the transcriptional changes in G. soja and G. max. We identified a sizable proportion of interspecific differentially expressed genes (DEGs) and found parental expression level dominance of G. max in their F1 hybrids. By classifying genes into different regulatory divergence types, we found the trans-regulatory changes played a predominant role in transcriptional divergence between wild and cultivated soybean. The same gene ontology (GO) and protein family (Pfam) terms were found to be over-represented in DEGs and genes of cis-only between JY47 and GS, suggesting the substantial contribution of cis-regulatory divergences to the evolution of wild and cultivated soybeans. By further dissecting genes into five different duplication modes, we found genes in different duplication modes tend to accumulate different types of regulatory differences. A relatively higher proportion of cis-only regulatory divergences was detected in singleton, dispersed, proximal, and tandem duplicates than WGD duplicates and genome-wide level, which is in line with the prediction of gene balance hypothesis for the differential fates of duplicated genes post-WGD. The numbers of cis-only and trans-only regulated genes were similar for singletons, whereas there were more genes of trans-only than cis-only in the rest duplication types, especially in WGD in which there were two times more trans-only genes than that in cis-only type. Tandem duplicates showed the highest proportion of trans-only genes probably due to some special features of this class. In summary, our results demonstrate that genes in different duplication modes have different fates in transcriptional evolution underpinned by cis- or trans-regulatory divergences in soybean and likely in other paleopolyploid higher organisms.

Evolution of Homeologous Gene Expression in Polyploid Wheat

Polyploidization has played a prominent role in the evolutionary history of plants. Two recent and sequential allopolyploidization events have resulted in the formation of wheat species with different ploidies, and which provide a model to study the effects of polyploidization on the evolution of gene expression. In this study, we identified differentially expressed genes (DEGs) between four BBAA tetraploid wheats of three different ploidy backgrounds. DEGs were found to be unevenly distributed among functional categories and duplication modes. We observed more DEGs in the extracted tetraploid wheat (ETW) than in natural tetraploid wheats (TD and TTR13) as compared to a synthetic tetraploid (AT2). Furthermore, DEGs showed higher Ka/Ks ratios than those that did not show expression changes (non-DEGs) between genotypes, indicating DEGs and non-DEGs experienced different selection pressures. For A-B homeolog pairs with DEGs, most of them had only one differentially expressed copy, however, when both copies of a homeolog pair were DEGs, the A and B copies were more likely to be regulated to the same direction. Our results suggest that both cis- and inter-subgenome trans-regulatory changes are important drivers in the evolution of homeologous gene expression in polyploid wheat, with ploidy playing a significant role in the process.

Functional Constitutional Dynamic Networks Revealing Evolutionary Reproduction/Variation/Selection Principles

Within the broad research efforts to engineer chemical pathways to yield high-throughput evolutionary synthesis of genes and their screening for dictated functionalities, we introduce the evolution of nucleic-acid-based constitutional dynamic networks (CDNs) that follow reproduction/variation/selection principles. These fundamental principles are demonstrated by assembling a library of nucleic-acid strands and hairpins as functional modules for evolving networks. Primary T1-initiated selection of components from the library assembles a parent CDN X, where the evolved constituents exhibit catalytic properties to cleave the hairpins in the library. Cleavage of the hairpins yields fragments, which reproduces T1 to replicate CDN X, whereas the other fragments T2 and T3 select other components to evolve two other CDNs, Y and Z (variation). By applying appropriate counter triggers, we demonstrate the guided selection of networks from the evolved CDNs. By integrating additional hairpin substrates into the system, CDN-dictated emergent catalytic transformations are accomplished. The study provides pathways to construct evolutionary dynamic networks revealing enhanced gated and cascaded functions.

Out of the testis, into the ovary: biased outcomes of gene duplication and deletion in Drosophila

Gene turnover is a key source of adaptive variation. Yet most evolutionary studies have focused on gene duplication, dismissing gene deletion as a mechanism that simply eradicates redundancy. Here, I use genome-scale sequence and multi-tissue expression data from Drosophila melanogaster and Drosophila pseudoobscura to simultaneously assess the evolutionary outcomes of gene duplication and deletion in Drosophila. I find that gene duplication is more frequent than gene deletion in both species, indicating that it may play a more important role in Drosophila evolution. However, examination of several genic properties reveals that genes likely possess distinct functions after duplication that diverge further before deletion, suggesting that loss of redundancy cannot explain a majority of gene deletion events in Drosophila. Moreover, in addition to providing support for the well-known "out of the testis" origin of young duplicate genes, analyses of gene expression profiles uncover a preferential bias against deletion of old ovary-expressed genes. Therefore, I propose a novel "into the ovary" hypothesis for gene deletion in Drosophila, in which gene deletion may promote adaptation by salvaging genes that contribute to the evolution of female reproductive phenotypes. Under this combined "out of the testis, into the ovary" evolutionary model, gene duplication and deletion work in concert to generate and maintain a balanced repertoire of genes that promote sex-specific adaptation in Drosophila.

Exploring whole-genome duplicate gene retention with complex genetic interaction analysis

Whole-genome duplication has played a central role in the genome evolution of many organisms, including the human genome. Most duplicated genes are eliminated, and factors that influence the retention of persisting duplicates remain poorly understood. We describe a systematic complex genetic interaction analysis with yeast paralogs derived from the whole-genome duplication event. Mapping of digenic interactions for a deletion mutant of each paralog, and of trigenic interactions for the double mutant, provides insight into their roles and a quantitative measure of their functional redundancy. Trigenic interaction analysis distinguishes two classes of paralogs: a more functionally divergent subset and another that retained more functional overlap. Gene feature analysis and modeling suggest that evolutionary trajectories of duplicated genes are dictated by combined functional and structural entanglement factors.

Duplication and diversification of trehalase confers evolutionary advantages on lepidopteran insects

Gene duplication provides a major source of new genes for evolutionary novelty and ecological adaptation. However, the maintenance of duplicated genes and their relevance to adaptive evolution has long been debated. Insect trehalase (Treh) plays key roles in energy metabolism, growth, and stress recovery. Here, we show that the duplication of Treh in Lepidoptera (butterflies and moths) is linked with their adaptation to various environmental stresses. Generally, two Treh genes are present in insects: Treh1 and Treh2. We report three distinct forms of Treh in lepidopteran insects, where Treh1 was duplicated into two gene clusters (Treh1a and Treh1b). These gene clusters differ in gene expression patterns, enzymatic properties, and subcellular localizations, suggesting that the enzymes probably underwent sub- and/or neofunctionalization in the lepidopteran insects. Interestingly, selective pressure analysis provided significant evidence of positive selection on duplicate Treh1b gene in lepidopteran insect lineages. Most positively selected sites were located in the alpha-helical region, and several sites were close to the trehalose binding and catalytic sites. Subcellular adaptation of duplicate Treh1b driven by positive selection appears to have occurred as a result of selected changes in specific sequences, allowing for rapid reprogramming of duplicated Treh during evolution. Our results suggest that gene duplication of Treh and subsequent functional diversification could increase the survival rate of lepidopteran insects through various regulations of intracellular trehalose levels, facilitating their adaptation to diverse habitats. This study provides evidence regarding the mechanism by which gene family expansion can contribute to species adaptation through gene duplication and subsequent functional diversification.

Duplication of toll-like receptor 22 in teleost fishes

The TLRs of teleost fishes have distinct features and are highly diverse, but the duplication characteristics and expression patterns of the tlr22 gene remain unclear. Here, we identified paralogous tlr22 genes in 13 teleost fishes by screening available fish genomic resources and using molecular cloning. We then conducted comprehensive bioinformatics analyses and investigated spatiotemporal differences in the expression patterns of the tlr22 genes in G. eckloni. The results indicated that more than three paralogous tlr22 genes were possessed by some teleost fishes. Of these, tlr22c is specific to some subfamilies of the Cyprinidae (e.g., Barbinae, Cyprininae, Schizothoracinae, and Leuciscinae). Phylogenetic and syntenic analyses showed that the paralogous tlr22 genes originated from two single-gene duplication events. Molecular clock calculations dated the two gene duplication events at 49.5 and 39.3 MYA, which is before the common carp-specific genome duplication event and well after the fish-specific genome duplication. Gene duplication of tlr22 was followed by gene loss or pseudogene events in certain lineages. Spatiotemporal expression differences between the three duplicated tlr22 genes from G. eckloni suggested that these genes diverged functionally after gene duplication.

Evolution of Nucleic-Acid-Based Constitutional Dynamic Networks Revealing Adaptive and Emergent Functions

The evolution of networks is a fundamental unresolved issue in developing the area of systems chemistry. We introduce a versatile rewiring mechanism that leads to the emergence of nucleic-acid-based constitutional dynamic networks (CDNs). A two-component constituent AA' functionalized with a Mg²⁺ -ion-dependent DNAzyme activator unit forms a complex with an intact hairpin H_BB' composed of B and B' sequences. Cleavage of H_BB' leads to the two-component constituent BB', and its rewiring with AA' yields CDN X composed of the equilibrated constituents AA', AB', BA', and BB'. In analogy, subjecting AA' to an intact hairpin H_CC' leads to the formation of CDN Y consisting of AA', AC', CA', and CC'. Subjecting AA' to the mixture of H_BB' and H_CC' evolves the [3×3] CDN Z, composed of nine constituents, thus demonstrating hierarchical adaptive properties. Furthermore, the DNAzyme units associated with the constituents are applied to tailor emerging catalytic functions from the different CDNs.

Genome-wide identification, molecular evolution, and expression analysis of auxin response factor (ARF) gene family in Brachypodium distachyon L

BACKGROUND

The auxin response factor (ARF) gene family is involved in plant development and hormone regulation. Although the ARF gene family has been studied in some plant species, its structural features, molecular evolution, and expression profiling in Brachypodium distachyon L. are still not clear.

RESULTS

Genome-wide analysis identified 19 ARF genes in B. distachyon. A phylogenetic tree constructed with 182 ARF genes from seven plant species revealed three different clades, and the ARF genes from within a clade exhibited structural conservation, although certain divergences occurred in different clades. The branch-site model identified some sites where positive selection may have occurred, and functional divergence analysis found more Type II divergence sites than Type I. In particular, both positive selection and functional divergence may have occurred in 241H, 243G, 244 L, 310 T, 340G and 355 T. Subcellular localization prediction and experimental verification indicated that BdARF proteins were present in the nucleus. Transcript expression analysis revealed that BdARFs were mainly expressed in the leaf and root tips, stems, and developing seeds. Some BdARF genes exhibited significantly upregulated expression under various abiotic stressors. Particularly, BdARF4 and BdARF8 were significantly upregulated in response to abiotic stress factors such as salicylic acid and heavy metals.

CONCLUSION

The ARF gene family in B. distachyon was highly conserved. Several important amino acid sites were identified where positive selection and functional divergence occurred, and they may play important roles in functional differentiation. BdARF genes had clear tissue and organ expression preference and were involved in abiotic stress response, suggesting their roles in plant growth and stress resistance.

Genome-wide Estimation of Evolutionary Distance and Phylogenetic Analysis of Homologous Genes

Homologous genes, including paralogs and orthologs, are genes that share sequence homologies within or between different species. Homologous genes originate from a common origin through speciation, genetic duplication or horizontal gene transfer. Estimation of the sequence divergence of homologous genes help us to understand divergence time, which makes it possible to understand the evolutionary patterns of speciation, gene duplication and gene transfer events. This protocol will provide a detailed bioinformatics pipeline on how to identify the homologous genes, compare their sequence divergence and phylogenetic relationships, focusing on homologous genes that show syntenic relationships using soybean (Glycine max) and common bean (Phaseolus vulgaris) as example species.

Comparative Transcriptomics Highlights New Features of the Iron Starvation Response in the Human Pathogen Candida glabrata

In this work, we used comparative transcriptomics to identify regulatory outliers (ROs) in the human pathogen Candida glabrata. ROs are genes that have very different expression patterns compared to their orthologs in other species. From comparative transcriptome analyses of the response of eight yeast species to toxic doses of selenite, a pleiotropic stress inducer, we identified 38 ROs in C. glabrata. Using transcriptome analyses of C. glabrata response to five different stresses, we pointed out five ROs which were more particularly responsive to iron starvation, a process which is very important for C. glabrata virulence. Global chromatin Immunoprecipitation and gene profiling analyses showed that four of these genes are actually new targets of the iron starvation responsive Aft2 transcription factor in C. glabrata. Two of them (HBS1 and DOM34b) are required for C. glabrata optimal growth in iron limited conditions. In S. cerevisiae, the orthologs of these two genes are involved in ribosome rescue by the NO GO decay (NGD) pathway. Hence, our results suggest a specific contribution of NGD co-factors to the C. glabrata adaptation to iron starvation.

Patterns and Consequences of Subgenome Differentiation Provide Insights into the Nature of Paleopolyploidy in Plants

Polyploidy is an important feature of plant genomes, but the nature of many polyploidization events remains to be elucidated. Here, we demonstrate that the evolutionary fates of the subgenomes in maize (Zea mays) and soybean (Glycine max) have followed different trajectories. One subgenome has been subject to relaxed selection, lower levels of gene expression, higher rates of transposable element accumulation, more small interfering RNAs and DNA methylation around genes, and higher rates of gene loss in maize, whereas none of these features were observed in soybean. Nevertheless, individual gene pairs exhibit differentiation with respect to these features in both species. In addition, we observed a higher number of chromosomal rearrangements and higher frequency of retention of duplicated genes in soybean than in maize. Furthermore, soybean "singletons" were found to be more frequently tandemly duplicated than "duplicates" in soybean, which may, to some extent, counteract the genome imbalance caused by gene loss. We propose that unlike in maize, in which two subgenomes were distinct prior to the allotetraploidization event and thus experienced global differences in selective constraints, in soybean, the two subgenomes were far less distinct prior to polyploidization, such that individual gene pairs, rather than subgenomes, experienced stochastic differences over longer periods of time, resulting in retention of the majority of duplicates.

Intrinsic protein disorder reduces small-scale gene duplicability

Whereas the rate of gene duplication is relatively high, only certain duplications survive the filter of natural selection and can contribute to genome evolution. However, the reasons why certain genes can be retained after duplication whereas others cannot remain largely unknown. Many proteins contain intrinsically disordered regions (IDRs), whose structures fluctuate between alternative conformational states. Due to their high flexibility, IDRs often enable protein–protein interactions and are the target of post-translational modifications. Intrinsically disordered proteins (IDPs) have characteristics that might either stimulate or hamper the retention of their encoding genes after duplication. On the one hand, IDRs may enable functional diversification, thus promoting duplicate retention. On the other hand, increased IDP availability is expected to result in deleterious unspecific interactions. Here, we interrogate the proteomes of human, Drosophila melanogaster, Caenorhabditis elegans, Saccharomyces cerevisiae, Arabidopsis thaliana and Escherichia coli, in order to ascertain the impact of protein intrinsic disorder on gene duplicability. We show that, in general, proteins encoded by duplicated genes tend to be less disordered than those encoded by singletons. The only exception is proteins encoded by ohnologs, which tend to be more disordered than those encoded by singletons or genes resulting from small-scale duplications. Our results indicate that duplication of genes encoding IDPs outside the context of whole-genome duplication (WGD) is often deleterious, but that IDRs facilitate retention of duplicates in the context of WGD. We discuss the potential evolutionary implications of our results.

Expansion of TALE homeobox genes and the evolution of spiralian development

Spiralians, including molluscs, annelids and platyhelminths, share a unique development process that includes the typical geometry of early cleavage and early segregation of cell fate in blastomeres along the animal-vegetal axis. However, the molecular mechanisms underlying this early cell fate segregation are largely unknown. Here, we report spiralian-specific expansion of the three-amino-acid loop extension (TALE) class of homeobox genes. During early development, some of these TALE genes are expressed in staggered domains along the animal-vegetal axis in the limpet Nipponacmea fuscoviridis and the polychaete Spirobranchus kraussii. Inhibition or overexpression of these genes alters the developmental fate of blastomeres, as predicted by the gene expression patterns. These results suggest that the expansion of novel TALE genes plays a critical role in the establishment of a novel cell fate segregation mechanism in spiralians.

Diversification of Transcriptional Regulation Determines Subfunctionalization of Paralogous Branched Chain Aminotransferases in the Yeast Saccharomyces cerevisiae

Saccharomyces cerevisiae harbors BAT1 and BAT2 paralogous genes that encode branched chain aminotransferases and have opposed expression profiles and physiological roles . Accordingly, in primary nitrogen sources such as glutamine, BAT1 expression is induced, supporting Bat1-dependent valine-isoleucine-leucine (VIL) biosynthesis, while BAT2 expression is repressed. Conversely, in the presence of VIL as the sole nitrogen source, BAT1 expression is hindered while that of BAT2 is activated, resulting in Bat2-dependent VIL catabolism. The presented results confirm that BAT1 expression is determined by transcriptional activation through the action of the Leu3-α-isopropylmalate (α-IPM) active isoform, and uncovers the existence of a novel α-IPM biosynthetic pathway operating in a put3Δ mutant grown on VIL, through Bat2-Leu2-Leu1 consecutive action. The classic α-IPM biosynthetic route operates in glutamine through the action of the leucine-sensitive α-IPM synthases. The presented results also show that BAT2 repression in glutamine can be alleviated in a ure2Δ mutant or through Gcn4-dependent transcriptional activation. Thus, when S. cerevisiae is grown on glutamine, VIL biosynthesis is predominant and is preferentially achieved through BAT1; while on VIL as the sole nitrogen source, catabolism prevails and is mainly afforded by BAT2.

The ancestral levels of transcription and the evolution of sexual phenotypes in filamentous fungi

Changes in gene expression have been hypothesized to play an important role in the evolution of divergent morphologies. To test this hypothesis in a model system, we examined differences in fruiting body morphology of five filamentous fungi in the Sordariomycetes, culturing them in a common garden environment and profiling genome-wide gene expression at five developmental stages. We reconstructed ancestral gene expression phenotypes, identifying genes with the largest evolved increases in gene expression across development. Conducting knockouts and performing phenotypic analysis in two divergent species typically demonstrated altered fruiting body development in the species that had evolved increased expression. Our evolutionary approach to finding relevant genes proved far more efficient than other gene deletion studies targeting whole genomes or gene families. Combining gene expression measurements with knockout phenotypes facilitated the refinement of Bayesian networks of the genes underlying fruiting body development, regulation of which is one of the least understood processes of multicellular development.

Brownian model of transcriptome evolution and phylogenetic network visualization between tissues

While phylogenetic analysis of transcriptomes of the same tissue is usually congruent with the species tree, the controversy emerges when multiple tissues are included, that is, whether species from the same tissue are clustered together, or different tissues from the same species are clustered together. Recent studies have suggested that phylogenetic network approach may shed some lights on our understanding of multi-tissue transcriptome evolution; yet the underlying evolutionary mechanism remains unclear. In this paper we develop a Brownian-based model of transcriptome evolution under the phylogenetic network that can statistically distinguish between the patterns of species-clustering and tissue-clustering. Our model can be used as a null hypothesis (neutral transcriptome evolution) for testing any correlation in tissue evolution, can be applied to cancer transcriptome evolution to study whether two tumors of an individual appeared independently or via metastasis, and can be useful to detect convergent evolution at the transcriptional level.

How evolution learns to generalise: Using the principles of learning theory to understand the evolution of developmental organisation

One of the most intriguing questions in evolution is how organisms exhibit suitable phenotypic variation to rapidly adapt in novel selective environments. Such variability is crucial for evolvability, but poorly understood. In particular, how can natural selection favour developmental organisations that facilitate adaptive evolution in previously unseen environments? Such a capacity suggests foresight that is incompatible with the short-sighted concept of natural selection. A potential resolution is provided by the idea that evolution may discover and exploit information not only about the particular phenotypes selected in the past, but their underlying structural regularities: new phenotypes, with the same underlying regularities, but novel particulars, may then be useful in new environments. If true, we still need to understand the conditions in which natural selection will discover such deep regularities rather than exploiting ‘quick fixes’ (i.e., fixes that provide adaptive phenotypes in the short term, but limit future evolvability). Here we argue that the ability of evolution to discover such regularities is formally analogous to learning principles, familiar in humans and machines, that enable generalisation from past experience. Conversely, natural selection that fails to enhance evolvability is directly analogous to the learning problem of over-fitting and the subsequent failure to generalise. We support the conclusion that evolving systems and learning systems are different instantiations of the same algorithmic principles by showing that existing results from the learning domain can be transferred to the evolution domain. Specifically, we show that conditions that alleviate over-fitting in learning systems successfully predict which biological conditions (e.g., environmental variation, regularity, noise or a pressure for developmental simplicity) enhance evolvability. This equivalence provides access to a well-developed theoretical framework from learning theory that enables a characterisation of the general conditions for the evolution of evolvability.

A striking feature of evolving organisms is their ability to acquire novel characteristics that help them adapt in new environments. The origin and the conditions of such ability remain elusive and is a long-standing question in evolutionary biology. Recent theory suggests that organisms can evolve designs that help them generate novel features that are more likely to be beneficial. Specifically, this is possible when the environments that organisms are exposed to share common regularities. However, the organisms develop robust designs that tend to produce what had been selected in the past and might be inflexible for future environments. The resolution comes from a recent theory introduced by Watson and Szathmáry that suggests a deep analogy between learning and evolution. Accordingly, here we utilise learning theory to explain the conditions that lead to more evolvable designs. We successfully demonstrate this by equating evolvability to the way humans and machines generalise to previously-unseen situations. Specifically, we show that the same conditions that enhance generalisation in learning systems have biological analogues and help us understand why environmental noise and the reproductive and maintenance costs of gene-regulatory connections can lead to more evolvable designs.

Parallel reorganization of protein function in the spindle checkpoint pathway through evolutionary paths in the fitness landscape that appear neutral in laboratory experiments

Regulatory networks often increase in complexity during evolution through gene duplication and divergence of component proteins. Two models that explain this increase in complexity are: 1) adaptive changes after gene duplication, such as resolution of adaptive conflicts, and 2) non-adaptive processes such as duplication, degeneration and complementation. Both of these models predict complementary changes in the retained duplicates, but they can be distinguished by direct fitness measurements in organisms with short generation times. Previously, it has been observed that repeated duplication of an essential protein in the spindle checkpoint pathway has occurred multiple times over the eukaryotic tree of life, leading to convergent protein domain organization in its duplicates. Here, we replace the paralog pair in S. cerevisiae with a single-copy protein from a species that did not undergo gene duplication. Surprisingly, using quantitative fitness measurements in laboratory conditions stressful for the spindle-checkpoint pathway, we find no evidence that reorganization of protein function after gene duplication is beneficial. We then reconstruct several evolutionary intermediates from the inferred ancestral network to the extant one, and find that, at the resolution of our assay, there exist stepwise mutational paths from the single protein to the divergent pair of extant proteins with no apparent fitness defects. Parallel evolution has been taken as strong evidence for natural selection, but our results suggest that even in these cases, reorganization of protein function after gene duplication may be explained by neutral processes.

Exploring xylose metabolism in Spathaspora species: XYL1.2 from Spathaspora passalidarum as the key for efficient anaerobic xylose fermentation in metabolic engineered Saccharomyces cerevisiae

BACKGROUND

The production of ethanol and other fuels and chemicals from lignocellulosic materials is dependent of efficient xylose conversion. Xylose fermentation capacity in yeasts is usually linked to xylose reductase (XR) accepting NADH as cofactor. The XR from Scheffersomyces stipitis, which is able to use NADH as cofactor but still prefers NADPH, has been used to generate recombinant xylose-fermenting Saccharomyces cerevisiae. Novel xylose-fermenting yeasts species, as those from the Spathaspora clade, have been described and are potential sources of novel genes to improve xylose fermentation in S. cerevisiae.

RESULTS

Xylose fermentation by six strains from different Spathaspora species isolated in Brazil, plus the Sp. passalidarum type strain (CBS 10155(T)), was characterized under two oxygen-limited conditions. The best xylose-fermenting strains belong to the Sp. passalidarum species, and their highest ethanol titers, yields, and productivities were correlated to higher XR activity with NADH than with NADPH. Among the different Spathaspora species, Sp. passalidarum appears to be the sole harboring two XYL1 genes: XYL1.1, similar to the XYL1 found in other Spathaspora and yeast species and XYL1.2, with relatively higher expression level. XYL1.1p and XYL1.2p from Sp. passalidarum were expressed in S. cerevisiae TMB 3044 and XYL1.1p was confirmed to be strictly NADPH-dependent, while XYL1.2p to use both NADPH and NADH, with higher activity with the later. Recombinant S. cerevisiae strains expressing XYL1.1p did not show anaerobic growth in xylose medium. Under anaerobic xylose fermentation, S. cerevisiae TMB 3504, which expresses XYL1.2p from Sp. passalidarum, revealed significant higher ethanol yield and productivity than S. cerevisiae TMB 3422, which harbors XYL1p N272D from Sc. stipitis in the same isogenic background (0.40 vs 0.34 g gCDW (-1) and 0.33 vs 0.18 g gCDW (-1) h(-1), respectively).

CONCLUSION

This work explored a new clade of xylose-fermenting yeasts (Spathaspora species) towards the engineering of S. cerevisiae for improved xylose fermentation. The new S. cerevisiae TMB 3504 displays higher XR activity with NADH than with NADPH, with consequent improved ethanol yield and productivity and low xylitol production. This meaningful advance in anaerobic xylose fermentation by recombinant S. cerevisiae (using the XR/XDH pathway) paves the way for the development of novel industrial pentose-fermenting strains.

The divergence and positive selection of the plant-specific BURP-containing protein family

BURP domain-containing proteins belong to a plant-specific protein family and have diverse roles in plant development and stress responses. However, our understanding about the genetic divergence patterns and evolutionary rates of these proteins remain inadequate. In this study, 15 plant genomes were explored to elucidate the genetic origins, divergence, and functions of these proteins. One hundred and twenty-five BURP protein-encoding genes were identified from four main plant lineages, including 13 higher plant species. The absence of BURP family genes in unicellular and multicellular algae suggests that this family (1) appeared when plants shifted from relatively stable aquatic environments to land, where conditions are more variable and stressful, and (2) is critical in the adaptation of plants to adverse environments. Promoter analysis revealed that several responsive elements to plant hormones and external environment stresses are concentrated in the promoter region of BURP protein-encoding genes. This finding confirms that these genes influence plant stress responses. Several segmentally and tandem-duplicated gene pairs were identified from eight plant species. Thus, in general, BURP domain-containing genes have been subject to strong positive selection, even though these genes have conformed to different expansion models in different species. Our study also detected certain critical amino acid sites that may have contributed to functional divergence among groups or subgroups. Unexpectedly, all of the critical amino acid residues of functional divergence and positive selection were exclusively located in the C-terminal region of the BURP domain. In conclusion, our results contribute novel insights into the genetic divergence patterns and evolutionary rates of BURP proteins.

Drosophila duplicate genes evolve new functions on the fly

Gene duplication is thought to play a key role in phenotypic innovation. While several processes have been hypothesized to drive the retention and functional evolution of duplicate genes, their genomic contributions have never been determined. We recently developed the first genome-wide method to classify these processes by comparing distances between expression profiles of duplicate genes and their ancestral single-copy orthologs. Application of our approach to spatial gene expression profiles in two Drosophila species revealed that a majority of young duplicate genes possess new functions, and that new functions are acquired rapidly-often within a few million years. Surprisingly, new functions tend to arise in younger copies of duplicate gene pairs. Moreover, we found that young duplicates are often specifically expressed in testes, whereas old duplicates are broadly expressed across several tissues, providing strong support for the hypothetical "out-of-testes" origin of new genes. In this Extra View, I discuss our findings in the context of theoretical predictions about gene duplication, with a particular emphasis on the importance of natural selection in the evolution of novel phenotypes.

A Minimal Set of Glycolytic Genes Reveals Strong Redundancies in Saccharomyces cerevisiae Central Metabolism

As a result of ancestral whole-genome and small-scale duplication events, the genomes of Saccharomyces cerevisiae and many eukaryotes still contain a substantial fraction of duplicated genes. In all investigated organisms, metabolic pathways, and more particularly glycolysis, are specifically enriched for functionally redundant paralogs. In ancestors of the Saccharomyces lineage, the duplication of glycolytic genes is purported to have played an important role leading to S. cerevisiae's current lifestyle favoring fermentative metabolism even in the presence of oxygen and characterized by a high glycolytic capacity. In modern S. cerevisiae strains, the 12 glycolytic reactions leading to the biochemical conversion from glucose to ethanol are encoded by 27 paralogs. In order to experimentally explore the physiological role of this genetic redundancy, a yeast strain with a minimal set of 14 paralogs was constructed (the "minimal glycolysis" [MG] strain). Remarkably, a combination of a quantitative systems approach and semiquantitative analysis in a wide array of growth environments revealed the absence of a phenotypic response to the cumulative deletion of 13 glycolytic paralogs. This observation indicates that duplication of glycolytic genes is not a prerequisite for achieving the high glycolytic fluxes and fermentative capacities that are characteristic of S. cerevisiae and essential for many of its industrial applications and argues against gene dosage effects as a means of fixing minor glycolytic paralogs in the yeast genome. The MG strain was carefully designed and constructed to provide a robust prototrophic platform for quantitative studies and has been made available to the scientific community.

Rapid divergence and diversification of mammalian duplicate gene functions

BACKGROUND

Gene duplication provides raw material for the evolution of functional innovation. We recently developed a phylogenetic method that classifies evolutionary processes driving the retention of duplicate genes by quantifying divergence between their spatial gene expression profiles and that of their single-copy orthologous gene in a closely related sister species.

RESULTS

Here, we apply our classification method to pairs of duplicate genes in eight mammalian genomes, using data from 11 tissues to construct spatial gene expression profiles. We find that young mammalian duplicates are often functionally conserved, and that expression divergence rapidly increases over evolutionary time. Moreover, expression divergence results in increased tissue specificity, with an overrepresentation of expression in male kidney, underrepresentation of expression in female liver, and strong underrepresentation of expression in testis. Thus, duplicate genes acquire a diversity of new tissue-specific functions outside of the testis, possibly contributing to the origin of a multitude of complex phenotypes during mammalian evolution.

CONCLUSIONS

Our findings reveal that mammalian duplicate genes are initially functionally conserved, and then undergo rapid functional divergence over evolutionary time, acquiring diverse tissue-specific biological roles. These observations are in stark contrast to the much faster expression divergence and acquisition of broad housekeeping roles we previously observed in Drosophila duplicate genes. Due to the smaller effective population sizes of mammals relative to Drosophila, these analyses implicate natural selection in the functional evolution of duplicate genes.

Understanding tissue expression evolution: from expression phylogeny to phylogenetic network

Our understanding of tissue expression evolution in multi-cellular model organisms has been considerably advanced with the help of high-throughput technologies from EST, microarray to RNA-seq. Yet, many controversies remained unsolved, ranging from the evolutionary patterns of tissue expressions to expression phylogenetic analysis. Moreover, despite numerous reports published, it is desirable to have a general framework for study of tissue expression evolution. In this article, we first provide an up-to-date and concise review for the study of tissue expression evolution in multi-cellular organisms. While the expression phylogeny of the same tissues sampled from closely or intermediately related species largely reflects the species phylogeny, we demonstrate that phylogenetic network approach may shed some lights for our understanding of the developmental similarity and evolutionary relatedness during the multi-tissue evolution.

Dosage compensation can buffer copy-number variation in wild yeast

Aneuploidy is linked to myriad diseases but also facilitates organismal evolution. It remains unclear how cells overcome the deleterious effects of aneuploidy until new phenotypes evolve. Although laboratory strains are extremely sensitive to aneuploidy, we show here that aneuploidy is common in wild yeast isolates, which show lower-than-expected expression at many amplified genes. We generated diploid strain panels in which cells carried two, three, or four copies of the affected chromosomes, to show that gene-dosage compensation functions at >30% of amplified genes. Genes subject to dosage compensation are under higher expression constraint in wild populations—but they show elevated rates of gene amplification, suggesting that copy-number variation is buffered at these genes. We find that aneuploidy provides a clear ecological advantage to oak strain YPS1009, by amplifying a causal gene that escapes dosage compensation. Our work presents a model in which dosage compensation buffers gene amplification through aneuploidy to provide a natural, but likely transient, route to rapid phenotypic evolution.

DOI:http://dx.doi.org/10.7554/eLife.05462.001

Evolution is driven by changes to the genes and other genetic information found in the DNA of an organism. These changes might, for example, alter the physical characteristics of the organism, or change how efficiently crucial tasks are carried out inside cells. Whatever the change, if it makes it easier for the organism to survive and reproduce, it is more likely to be passed on to future generations.

DNA is organized inside cells in structures called chromosomes. Most of the cells in animals, plants, and fungi contain two copies of each chromosome. However, sometimes mistakes happen during cell division and extra copies of a chromosome—and hence the genes contained within it—may end up in a cell. These extra copies of genes might help to speed up the rate at which a species evolves, as the ‘spare’ copies are free to adapt to new roles. However, having extra copies of genes can also often be harmful, and in humans can cause genetic disorders such as Down syndrome.

In the laboratory, chromosomes are commonly studied in a species of yeast called Saccharomyces cerevisiae. This species consists of several groups—or strains—that are genetically distinct from each other. Over the years, breeding the yeast for experiments has created laboratory strains that have lost some of the characteristics seen in wild strains. Earlier studies suggested that these cells fail to grow properly if they contain extra copies of chromosomes. Now, Hose et al. have studied nearly 50 wild strains of Saccharomyces cerevisiae. In these, extra copies of chromosomes are commonplace, and seemingly have no detrimental effect on growth. Instead, Hose et al. found that cells with too many copies of a gene use many of those genes less often than would be expected. This process is known as ‘dosage compensation’. This dosage compensation has not been observed in laboratory strains, in part because the extra gene copies make them sickly and hard to study.

Together, the results provide examples of how dosage compensation could help new traits to evolve in a species by reducing the negative effects of duplicated genes. This knowledge may have broad application, from suggesting methods to alleviate human disorders to implicating new ways to engineer useful traits in yeast and other microbes.

DOI:http://dx.doi.org/10.7554/eLife.05462.002

How do regulatory networks evolve and expand throughout evolution?

Throughout evolution, regulatory networks need to expand and adapt to accommodate novel genes and gene functions. However, the molecular details explaining how gene networks evolve remain largely unknown. Recent studies demonstrate that changes in transcription factors contribute to the evolution of regulatory networks. In particular, duplication of transcription factors followed by specific mutations in their DNA-binding or interaction domains propels the divergence and emergence of new networks. The innate promiscuity and modularity of regulatory networks contributes to their evolvability: duplicated promiscuous regulators and their target promoters can acquire mutations that lead to gradual increases in specificity, allowing neofunctionalization or subfunctionalization.

Gene age predicts the strength of purifying selection acting on gene expression variation in humans

Gene expression levels can be subject to selection. We hypothesized that the age of gene origin is associated with expression constraints, given that it affects the level of gene integration into the functional cellular environment. By studying the genetic variation affecting gene expression levels (cis expression quantitative trait loci [cis-eQTLs]) and protein levels (cis protein QTLs [cis-pQTLs]), we determined that young, primate-specific genes are enriched in cis-eQTLs and cis-pQTLs. Compared to cis-eQTLs of old genes originating before the zebrafish divergence, cis-eQTLs of young genes have a higher effect size, are located closer to the transcription start site, are more significant, and tend to influence genes in multiple tissues and populations. These results suggest that the expression constraint of each gene increases throughout its lifespan. We also detected a positive correlation between expression constraints (approximated by cis-eQTL properties) and coding constraints (approximated by Ka/Ks) and observed that this correlation might be driven by gene age. To uncover factors associated with the increase in gene-age-related expression constraints, we demonstrated that gene connectivity, gene involvement in complex regulatory networks, gene haploinsufficiency, and the strength of posttranscriptional regulation increase with gene age. We also observed an increase in heritability of gene expression levels with age, implying a reduction of the environmental component. In summary, we show that gene age shapes key gene properties during evolution and is therefore an important component of genome function.

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.

Duplication of a promiscuous transcription factor drives the emergence of a new regulatory network

The emergence of new genes throughout evolution requires rewiring and extension of regulatory networks. However, the molecular details of how the transcriptional regulation of new gene copies evolves remain largely unexplored. Here we show how duplication of a transcription factor gene allowed the emergence of two independent regulatory circuits. Interestingly, the ancestral transcription factor was promiscuous and could bind different motifs in its target promoters. After duplication, one paralogue evolved increased binding specificity so that it only binds one type of motif, whereas the other copy evolved a decreased activity so that it only activates promoters that contain multiple binding sites. Interestingly, only a few mutations in both the DNA-binding domains and in the promoter binding sites were required to gradually disentangle the two networks. These results reveal how duplication of a promiscuous transcription factor followed by concerted cis and trans mutations allows expansion of a regulatory network.

Expression evolution facilitated the convergent neofunctionalization of a sodium channel gene

Ion channels have played a substantial role in the evolution of novel traits across all of the domains of life. A fascinating example of a novel adaptation is the convergent evolution of electric organs in the Mormyroid and Gymnotiform electric fishes. The regulated currents that flow through ion channels directly generate the electrical signals which have evolved in these fish. Here, we investigated how the expression evolution of two sodium channel paralogs (Scn4aa and Scn4ab) influenced their convergent molecular evolution following the teleost-specific whole-genome duplication. We developed a reliable assay to accurately measure the expression stoichiometry of these genes and used this technique to analyze relative expression of the duplicate genes in a phylogenetic context. We found that before a major shift in expression from skeletal muscle and neofunctionalization in the muscle-derived electric organ, Scn4aa was first downregulated in the ancestors of both electric lineages. This indicates that underlying the convergent evolution of this gene, there was a greater propensity toward neofunctionalization due to its decreased expression relative to its paralog Scn4ab. We investigated another derived muscle tissue, the sonic organ of Porichthys notatus, and show that, as in the electric fishes, Scn4aa again shows a radical shift in expression away from the ancestral muscle cells into the evolutionarily novel muscle-derived tissue. This study presents evidence that expression downregulation facilitates neofunctionalization after gene duplication, a pattern that may often set the stage for novel trait evolution after gene duplication.