2R and remodeling of vertebrate signal transduction engine

Gene Function is a Driver of Activin Signaling Pathway Evolution Following Whole-Genome Duplication in Rainbow Trout (Oncorhynchus mykiss)

The genomes of plant and animal species are influenced by ancestral whole-genome duplication (WGD) events, which have profound impacts on the regulation and function of gene networks. To gain insight into the consequences of WGD events, we characterized the sequence conservation and expression patterns of ohnologs in the highly duplicated activin receptor signaling pathway in rainbow trout (RBT). The RBT activin receptor signaling pathway is defined by tissue-specific expression of inhibitors and ligands and broad expression of receptors and Co-Smad signaling molecules. Signaling pathway ligands exhibited shared expression, while inhibitors and Smad signaling molecules primarily express a single dominant ohnolog. Our findings suggest that gene function influences ohnolog evolution following duplication of the activin signaling pathway in RBT.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

Paralog dispensability shapes homozygous deletion patterns in tumor genomes

Genomic instability is a hallmark of cancer, resulting in tumor genomes having large numbers of genetic aberrations, including homozygous deletions of protein coding genes. That tumor cells remain viable in the presence of such gene loss suggests high robustness to genetic perturbation. In model organisms and cancer cell lines, paralogs have been shown to contribute substantially to genetic robustness-they are generally more dispensable for growth than singletons. Here, by analyzing copy number profiles of > 10,000 tumors, we test the hypothesis that the increased dispensability of paralogs shapes tumor genome evolution. We find that genes with paralogs are more likely to be homozygously deleted and that this cannot be explained by other factors known to influence copy number variation. Furthermore, features that influence paralog dispensability in cancer cell lines correlate with paralog deletion frequency in tumors. Finally, paralogs that are broadly essential in cancer cell lines are less frequently deleted in tumors than non-essential paralogs. Overall, our results suggest that homozygous deletions of paralogs are more frequently observed in tumor genomes because paralogs are more dispensable.

Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes

The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.

Chromosomal Instability in Genome Evolution: From Cancer to Macroevolution

The integrity of the genome is crucial for the survival of all living organisms. However, genomes need to adapt to survive certain pressures, and for this purpose use several mechanisms to diversify. Chromosomal instability (CIN) is one of the main mechanisms leading to the creation of genomic heterogeneity by altering the number of chromosomes and changing their structures. In this review, we will discuss the different chromosomal patterns and changes observed in speciation, in evolutional biology as well as during tumor progression. By nature, the human genome shows an induction of diversity during gametogenesis but as well during tumorigenesis that can conclude in drastic changes such as the whole genome doubling to more discrete changes as the complex chromosomal rearrangement chromothripsis. More importantly, changes observed during speciation are strikingly similar to the genomic evolution observed during tumor progression and resistance to therapy. The different origins of CIN will be treated as the importance of double-strand breaks (DSBs) or the consequences of micronuclei. We will also explain the mechanisms behind the controlled DSBs, and recombination of homologous chromosomes observed during meiosis, to explain how errors lead to similar patterns observed during tumorigenesis. Then, we will also list several diseases associated with CIN, resulting in fertility issues, miscarriage, rare genetic diseases, and cancer. Understanding better chromosomal instability as a whole is primordial for the understanding of mechanisms leading to tumor progression.

Evidence for Multilevel Chemopreventive Activities of Natural Phenols from Functional Genomic Studies of Curcumin, Resveratrol, Genistein, Quercetin, and Luteolin

Herein, I present an updated and contextualized literature review of functional genomic studies of natural phenols in the context of cancer. I suggest multilevel chemopreventive and anticancer mechanisms of action, which are shared by multiple dietary natural phenols. Specifically, I cite evidence that curcumin and resveratrol have multilevel anti-cancer effects through: (1) inducing either p53-dependent or p53-independent apoptosis in cancer cell lines, (2) acting as potent regulators of expression of oncogenic and anti-oncogenic microRNAs, and (3) inducing complex epigenetic changes that can switch off oncogenes/switch on anti-oncogenes. There is no simple reductionist explanation for anti-cancer effects of curcumin and resveratrol. More generally, multilevel models of chemoprevention are suggested for related natural phenols and flavonoids such as genistein, quercetin, or luteolin.

Differential roles and regulation of the protein kinases PAK4, PAK5 and PAK6 in melanoma cells

The protein kinases PAK4, PAK5 and PAK6 comprise a family of ohnologues. In multiple cancers including melanomas PAK5 most frequently carries non-synonymous mutations; PAK6 and PAK4 have fewer; and PAK4 is often amplified. To help interpret these genomic data, initially we compared the cellular regulation of the sister kinases and their roles in melanoma cells. In common with many ohnologue protein kinases, PAK4, PAK5 and PAK6 each have two 14-3-3-binding phosphosites of which phosphoSer99 is conserved. PAK4 localises to the leading edge of cells in response to phorbol ester-stimulated binding of 14-3-3 to phosphoSer99 and phosphoSer181, which are phosphorylated by two different PKCs or PKDs. These phosphorylations of PAK4 are essential for its phorbol ester-stimulated phosphorylation of downstream substrates. In contrast, 14-3-3 interacts with PAK5 in response to phorbol ester-stimulated phosphorylation of Ser99 and epidermal growth factor-stimulated phosphorylation of Ser288; whereas PAK6 docks onto 14-3-3 and is prevented from localising to cell–cell junctions when Ser133 is phosphorylated in response to cAMP-elevating agents via PKA and insulin-like growth factor 1 via PKB/Akt. Silencing of PAK4 impairs viability, migration and invasive behaviour of melanoma cells carrying BRAF^V600E or NRAS^Q61K mutations. These defects are rescued by ectopic expression of PAK4, more so by a 14-3-3-binding deficient PAK4, and barely by PAK5 or PAK6. Together these genomic, biochemical and cellular data suggest that the oncogenic properties of PAK4 are regulated by PKC–PKD signalling in melanoma, while PAK5 and PAK6 are dispensable in this cancer.

Comprehensive expression analysis for the core cell cycle regulators in the chicken embryo reveals novel tissue-specific synexpression groups and similarities and differences with expression in mouse, frog and zebrafish

The core cell cycle machinery is conserved from yeast to humans, and hence it is assumed that all vertebrates share the same set of players. Yet during vertebrate evolution, the genome was duplicated twice, followed by a further genome duplication in teleost fish. Thereafter, distinct genes were retained in different vertebrate lineages; some individual gene duplications also occurred. To which extent these diversifying tendencies were compensated by retaining the same expression patterns across homologous genes is not known. This study for the first time undertook a comprehensive expression analysis for the core cell cycle regulators in the chicken, focusing in on early neurula and pharyngula stages of development, with the latter representing the vertebrate phylotypic stage. We also compared our data with published data for the mouse, Xenopus and zebrafish, the other established vertebrate models. Our work shows that, while many genes are expressed widely, some are upregulated or specifically expressed in defined tissues of the chicken embryo, forming novel synexpression groups with markers for distinct developmental pathways. Moreover, we found that in the neural tube and in the somite, mRNAs of some of the genes investigated accumulate in a specific subcellular localisation, pointing at a novel link between the site of mRNA translation, cell cycle control and interkinetic nuclear movements. Finally, we show that expression patterns of orthologous genes may differ in the four vertebrate models. Thus, for any study investigating cell proliferation, cell differentiation, tissue regeneration, stem cell behaviour and cancer/cancer therapy, it has to be carefully examined which of the observed effects are due to the specific model organism used, and which can be generalised.

MicroRNAs as Indicators into the Causes and Consequences of Whole-Genome Duplication Events

Whole-genome duplications (WGDs) have long been considered the causal mechanism underlying dramatic increases to morphological complexity due to the neo-functionalization of paralogs generated during these events. Nonetheless, an alternative hypothesis suggests that behind the retention of most paralogs is not neo-functionalization, but instead the degree of the inter-connectivity of the intended gene product, as well as the mode of the WGD itself. Here, we explore both the causes and consequences of WGD by examining the distribution, expression, and molecular evolution of microRNAs (miRNAs) in both gnathostome vertebrates as well as chelicerate arthropods. We find that although the number of miRNA paralogs tracks the number of WGDs experienced within the lineage, few of these paralogs experienced changes to the seed sequence, and thus are functionally equivalent relative to their mRNA targets. Nonetheless, in gnathostomes, although the retention of paralogs following the 1R autotetraploidization event is similar across the two subgenomes, the paralogs generated by the gnathostome 2R allotetraploidization event are retained in higher numbers on one subgenome relative to the second, with the miRNAs found on the preferred subgenome showing both higher expression of mature miRNA transcripts and slower molecular evolution of the precursor miRNA sequences. Importantly, WGDs do not result in the creation of miRNA novelty, nor do WGDs correlate to increases in complexity. Instead, it is the number of miRNA seed sequences in the genome itself that not only better correlate to instances in complexification, but also mechanistically explain why complexity increases when new miRNA families are established.

The impact of whole genome duplications on the human gene regulatory networks

This work studies the effects of the two rounds of Whole Genome Duplication (WGD) at the origin of the vertebrate lineage on the architecture of the human gene regulatory networks. We integrate information on transcriptional regulation, miRNA regulation, and protein-protein interactions to comparatively analyse the role of WGD and Small Scale Duplications (SSD) in the structural properties of the resulting multilayer network. We show that complex network motifs, such as combinations of feed-forward loops and bifan arrays, deriving from WGD events are specifically enriched in the network. Pairs of WGD-derived proteins display a strong tendency to interact both with each other and with common partners and WGD-derived transcription factors play a prominent role in the retention of a strong regulatory redundancy. Combinatorial regulation and synergy between different regulatory layers are in general enhanced by duplication events, but the two types of duplications contribute in different ways. Overall, our findings suggest that the two WGD events played a substantial role in increasing the multi-layer complexity of the vertebrate regulatory network by enhancing its combinatorial organization, with potential consequences on its overall robustness and ability to perform high-level functions like signal integration and noise control. Lastly, we discuss in detail the RAR/RXR pathway as an illustrative example of the evolutionary impact of WGD duplications in human.

Gene duplication is one of the main mechanisms driving genome evolution. The duplication of a genomic segment can be the result of a local event, involving only a small portion of the genome, or of a dramatic duplication of the whole genome, which is however only rarely retained. All vertebrates descend from two rounds of Whole-Genome Duplication (WGD) that occurred approximately 500 Mya. We show that these events influenced in unique ways the evolution of different human gene regulatory networks, with sizeable effects on their current structure. We find that WGDs statistically increased the presence of specific classes of simple genetic circuits, considered to be fundamental building blocks of more sophisticated circuitry and commonly associated to complex functions. Our findings support the hypothesis that these rare, large-scale events have played a substantial role in the emergence of complex traits in vertebrates.

Models of the Gene Must Inform Data-Mining Strategies in Genomics

The gene is a fundamental concept of genetics, which emerged with the Mendelian paradigm of heredity at the beginning of the 20th century. However, the concept has since diversified. Somewhat different narratives and models of the gene developed in several sub-disciplines of genetics, that is in classical genetics, population genetics, molecular genetics, genomics, and, recently, also, in systems genetics. Here, I ask how the diversity of the concept impacts data-integration and data-mining strategies for bioinformatics, genomics, statistical genetics, and data science. I also consider theoretical background of the concept of the gene in the ideas of empiricism and experimentalism, as well as reductionist and anti-reductionist narratives on the concept. Finally, a few strategies of analysis from published examples of data-mining projects are discussed. Moreover, the examples are re-interpreted in the light of the theoretical material. I argue that the choice of an optimal level of abstraction for the gene is vital for a successful genome analysis.

The 14-3-3 protein is an essential component of cyclic AMP signaling for regulation of chemotaxis and development in Dictyostelium

The evolutionarily-conserved 14-3-3 proteins regulate many cellular processes through binding to various phosphorylated targets in eukaryotes. It first appears in Dictyostelium, however its role in this organism is poorly understood. Here we show that down-regulation of the 14-3-3 impairs chemotaxis and causes multiple-tip formation in Dictyostelium. Mechanistically, the 14-3-3 is a critical component of cyclic AMP (cAMP) signaling and binds to nearly a hundred of proteins in Dictyostelium, including a number of evolutionarily-conserved proteins. 14-3-3 - interaction with its targets is up-regulated in response to developmental cues/regulators including starvation, osmotic stress and cAMP. cAMP stimulates 14-3-3 - binding to phospho-Ser431 on a guanine nucleotide exchange factor Gef-Q. Interestingly, overexpression of Gef-Q^Ser431Ala mutant but not wild-type Gef-Q protein causes a multiple-tip phenotype in Dictyostelium, which partially resembles phenotypes of the 14-3-3 - deficient mutant. Collectively, these data demonstrate that the 14-3-3 plays an important role in Dictyostelium and may help to deepen our understanding of the evolution of 14-3-3 - interactomes in eukaryotes.

A Contemporary Message from Mendel's Logical Empiricism

The gene is one of the most fundamental concepts in life sciences, having been developed in the mold of the Mendelian paradigm of heredity, which shaped genetics across 150 years. How could Mendel possibly be so prophetic in the middle of 19^th century, using only the small garden of the monastery as his experimental breeding field? I believe that we are indebted to Mendel's mastery of the scientific method, which was far ahead of his time. Although his experimental technology was literally garden-variety, Mendel's excellence in the method of science, algebra, and logical analysis helped him in designing the right experiment and in interpreting the results insightfully. This may be valuable to recall in today's technology-focused culture, where the center of interest tends to be on the generation and description of high-throughput datasets from specialized genomics screens. As Mendel's story suggests, progress in 21^st century genetics may also depend on the development of robust concepts and generalizations.

OHNOLOGS v2: a comprehensive resource for the genes retained from whole genome duplication in vertebrates

All vertebrates including human have evolved from an ancestor that underwent two rounds of whole genome duplication (2R-WGD). In addition, teleost fish underwent an additional third round of genome duplication (3R-WGD). The genes retained from these genome duplications, so-called ohnologs, have been instrumental in the evolution of vertebrate complexity, development and susceptibility to genetic diseases. However, the identification of vertebrate ohnologs has been challenging, due to lineage specific genome rearrangements since 2R- and 3R-WGD. We previously identified vertebrate ohnologs using a novel synteny comparison across multiple genomes. Here, we refine and apply this approach on 27 vertebrate genomes to identify ohnologs from both 2R- and 3R-WGD, while taking into account the phylogenetically biased sampling of available species. We assemble vertebrate ohnolog pairs and families in an expanded OHNOLOGS v2 database. We find that teleost fish have retained more 2R-WGD ohnologs than mammals and sauropsids, and that these 2R-ohnologs have retained significantly more ohnologs from the subsequent 3R-WGD than genes without 2R-ohnologs. Interestingly, species with fewer extant genes, such as sauropsids, have retained similar or higher proportions of ohnologs. OHNOLOGS v2 should allow deeper evolutionary genomic analysis of the impact of WGD on vertebrates and can be freely accessed at http://ohnologs.curie.fr.

Magic roundabout is an endothelial-specific ohnolog of ROBO1 which neo-functionalized to an essential new role in angiogenesis

BACKGROUND

Magic roundabout (ROBO4) is an unusual endothelial-specific paralog of the family of neuronally-expressed axon guidance receptors called roundabouts. Endothelial cells (ECs), whose uninterrupted sheet delimits the lumen of all vertebrate blood vessels and which are absent from invertebrate species, are a vertebrate-specific evolutionary novelty.

RESULTS

Herein, the evolutionary mechanism of the duplication, retention and divergence of ROBO4 was investigated for the first time. Phylogenetic analyses carried out suggested that ROBO4 is a fast-evolving paralog of ROBO1 formed at the base of vertebrates. The ancestral expression pattern was neuronal. ROBO4 dramatically shifted its expression and became exceptionally specific to ECs. The data-mining of FANTOM5 and ENCODE reveals that ROBO4's endothelial expression arises from a single transcription start site (TSS), conserved in mouse, controlled by a proximal promoter with a complex architecture suggestive of regulatory neo-functionalization. (An analysis of promoter probabilities suggested the architecture was not due to a chance arrangement of TFBSes). Further evidence for the neo-functionalization of ROBO4 comes from the analysis of its protein interactions, the rates of protein evolution, and of positively selected sites.

CONCLUSIONS

The neo-functionalization model explains why ROBO4 protein acquired new context-specific biological functions in the control of angiogenesis. This endothelial-specific roundabout receptor is an illustrative example of the emergence of an essential vertebrate molecular novelty and an endothelial-specific signaling sub-network through 2R-WGD. The emergence of novel cell types, such as ECs, might be a neglected evolutionary force contributing to the high rate of retention of duplicates post-2R-WGD. Crucially, expression neo-functionalization to evolutionarily novel sites of expression conceptually extends the classical model of neo-functionalization.

Pax3/7 duplicated and diverged independently in amphioxus, the basal chordate lineage

The Pax3/7 transcription factor family is integral to developmental gene networks contributing to important innovations in vertebrate evolution, including the neural crest. The basal chordate lineage of amphioxus is ideally placed to understand the dynamics of the gene regulatory network evolution that produced these novelties. We report here the discovery that the cephalochordate lineage possesses two Pax3/7 genes, Pax3/7a and Pax3/7b. The tandem duplication is ancestral to all extant amphioxus, occurring in both Asymmetron and Branchiostoma, but originated after the split from the lineage leading to vertebrates. The two paralogues are differentially expressed during embryonic development, particularly in neural and somitic tissues, suggesting distinct regulation. Our results have implications for the study of amphioxus regeneration, neural plate and crest evolution, and differential tandem paralogue evolution.

The effects of repeated whole genome duplication events on the evolution of cytokinin signaling pathway

BACKGROUND

It is thought that after whole-genome duplications (WGDs), a large fraction of the duplicated gene copies is lost over time while few duplicates are retained. Which factors promote survival or death of a duplicate remains unclear and the underlying mechanisms are poorly understood. According to the model of gene dosage balance, genes encoding interacting proteins are predicted to be preferentially co-retained after WGDs. Among these are genes encoding proteins involved in complexes or in signal transduction.

RESULTS

We have investigated the way that repeated WGDs during land plant evolution have affected cytokinin signaling to study patterns of gene duplicability and co-retention in this important signal transduction pathway. Through the integration of phylogenetic analyses with comparisons of genome collinearity, we have found that signal input mediated by cytokinin receptors proved to be highly conserved over long evolutionary time-scales, with receptors showing predominantly gene loss after repeated WGDs. However, the downstream elements, e,g. response regulators, were mainly retained after WGDs and thereby formed gene families in most plant lineages.

CONCLUSIONS

Gene dosage balance between the interacting components indicated by co-retention after WGDs seems to play a minor role in the evolution of cytokinin signaling pathway. Overall, core genes of cytokinin signaling show a highly heterogeneous pattern of gene retention after WGD, reflecting complex relationships between the various factors that shape the long-term fate of a duplicated gene.

Teleost Fish-Specific Preferential Retention of Pigmentation Gene-Containing Families After Whole Genome Duplications in Vertebrates

Vertebrate pigmentation is a highly diverse trait mainly determined by neural crest cell derivatives. It has been suggested that two rounds (1R/2R) of whole-genome duplications (WGDs) at the basis of vertebrates allowed changes in gene regulation associated with neural crest evolution. Subsequently, the teleost fish lineage experienced other WGDs, including the teleost-specific Ts3R before teleost radiation and the more recent Ss4R at the basis of salmonids. As the teleost lineage harbors the highest number of pigment cell types and pigmentation diversity in vertebrates, WGDs might have contributed to the evolution and diversification of the pigmentation gene repertoire in teleosts. We have compared the impact of the basal vertebrate 1R/2R duplications with that of the teleost-specific Ts3R and salmonid-specific Ss4R WGDs on 181 gene families containing genes involved in pigmentation. We show that pigmentation genes (PGs) have been globally more frequently retained as duplicates than other genes after Ts3R and Ss4R but not after the early 1R/2R. This is also true for non-pigmentary paralogs of PGs, suggesting that the function in pigmentation is not the sole key driver of gene retention after WGDs. On the long-term, specific categories of PGs have been repeatedly preferentially retained after ancient 1R/2R and Ts3R WGDs, possibly linked to the molecular nature of their proteins (e.g., DNA binding transcriptional regulators) and their central position in protein-protein interaction networks. Taken together, our results support a major role of WGDs in the diversification of the pigmentation gene repertoire in the teleost lineage, with a possible link with the diversity of pigment cell lineages observed in these animals compared to other vertebrates.

Identification and Evolution of TGF-β Signaling Pathway Members in Twenty-Four Animal Species and Expression in Tilapia

Transforming growth factor β (TGF-β) signaling controls diverse cellular processes during embryogenesis as well as in mature tissues of multicellular animals. Here we carried out a comprehensive analysis of TGF-β pathway members in 24 representative animal species. The appearance of the TGF-β pathway was intrinsically linked to the emergence of metazoan. The total number of TGF-β ligands, receptors, and smads changed slightly in all invertebrates and jawless vertebrates analyzed. In contrast, expansion of the pathway members, especially ligands, was observed in jawed vertebrates most likely due to the second round of whole genome duplication (2R) and additional rounds in teleosts. Duplications of TGFB2, TGFBR2, ACVR1, SMAD4 and SMAD6, which were resulted from 2R, were first isolated. Type II receptors may be originated from the ACVR2-like ancestor. Interestingly, AMHR2 was not identified in Chimaeriformes and Cypriniformes even though they had the ligand AMH. Based on transcriptome data, TGF-β ligands exhibited a tissue-specific expression especially in the heart and gonads. However, most receptors and smads were expressed in multiple tissues indicating they were shared by different ligands. Spatial and temporal expression profiles of 8 genes in gonads of different developmental stages provided a fundamental clue for understanding their important roles in sex determination and reproduction. Taken together, our findings provided a global insight into the phylogeny and expression patterns of the TGF-β pathway genes, and hence contribute to the greater understanding of their biological roles in the organism especially in teleosts.

Selective Constraints on Coding Sequences of Nervous System Genes Are a Major Determinant of Duplicate Gene Retention in Vertebrates

The evolutionary history of vertebrates is marked by three ancient whole-genome duplications: two successive rounds in the ancestor of vertebrates, and a third one specific to teleost fishes. Biased loss of most duplicates enriched the genome for specific genes, such as slow evolving genes, but this selective retention process is not well understood. To understand what drives the long-term preservation of duplicate genes, we characterized duplicated genes in terms of their expression patterns. We used a new method of expression enrichment analysis, TopAnat, applied to in situ hybridization data from thousands of genes from zebrafish and mouse. We showed that the presence of expression in the nervous system is a good predictor of a higher rate of retention of duplicate genes after whole-genome duplication. Further analyses suggest that purifying selection against the toxic effects of misfolded or misinteracting proteins, which is particularly strong in nonrenewing neural tissues, likely constrains the evolution of coding sequences of nervous system genes, leading indirectly to the preservation of duplicate genes after whole-genome duplication. Whole-genome duplications thus greatly contributed to the expansion of the toolkit of genes available for the evolution of profound novelties of the nervous system at the base of the vertebrate radiation.

Keeping the home fires burning: AMP-activated protein kinase

Living cells obtain energy either by oxidizing reduced compounds of organic or mineral origin or by absorbing light. Whichever energy source is used, some of the energy released is conserved by converting adenosine diphosphate (ADP) to adenosine triphosphate (ATP), which are analogous to the chemicals in a rechargeable battery. The energy released by the conversion of ATP back to ADP is used to drive most energy-requiring processes, including cell growth, cell division, communication and movement. It is clearly essential to life that the production and consumption of ATP are always maintained in balance, and the AMP-activated protein kinase (AMPK) is one of the key cellular regulatory systems that ensures this. In eukaryotic cells (cells with nuclei and other internal membrane-bound structures, including human cells), most ATP is produced in mitochondria, which are thought to have been derived by the engulfment of oxidative bacteria by a host cell not previously able to use molecular oxygen. AMPK is activated by increasing AMP or ADP (AMP being generated from ADP whenever ADP rises) coupled with falling ATP. Relatives of AMPK are found in essentially all eukaryotes, and it may have evolved to allow the host cell to monitor the output of the newly acquired mitochondria and step their ATP production up or down according to the demand. Structural studies have illuminated how AMPK achieves the task of detecting small changes in AMP and ADP, despite the presence of much higher concentrations of ATP. Recently, it has been shown that AMPK can also sense the availability of glucose, the primary carbon source for most eukaryotic cells, via a mechanism independent of changes in AMP or ADP. Once activated by energy imbalance or glucose lack, AMPK modifies many target proteins by transferring phosphate groups to them from ATP. By this means, numerous ATP-producing processes are switched on (including the production of new mitochondria) and ATP-consuming processes are switched off, thus restoring energy homeostasis. Drugs that modulate AMPK have great potential in the treatment of metabolic disorders such as obesity and Type 2 diabetes, and even cancer. Indeed, some existing drugs such as metformin and aspirin, which were derived from traditional herbal remedies, appear to work, in part, by activating AMPK.

Recent advances in understanding the roles of whole genome duplications in evolution

Ancient whole-genome duplications (WGDs)- paleopolyploidy events-are key to solving Darwin's 'abominable mystery' of how flowering plants evolved and radiated into a rich variety of species. The vertebrates also emerged from their invertebrate ancestors via two WGDs, and genomes of diverse gymnosperm trees, unicellular eukaryotes, invertebrates, fishes, amphibians and even a rodent carry evidence of lineage-specific WGDs. Modern polyploidy is common in eukaryotes, and it can be induced, enabling mechanisms and short-term cost-benefit assessments of polyploidy to be studied experimentally. However, the ancient WGDs can be reconstructed only by comparative genomics: these studies are difficult because the DNA duplicates have been through tens or hundreds of millions of years of gene losses, mutations, and chromosomal rearrangements that culminate in resolution of the polyploid genomes back into diploid ones (rediploidisation). Intriguing asymmetries in patterns of post-WGD gene loss and retention between duplicated sets of chromosomes have been discovered recently, and elaborations of signal transduction systems are lasting legacies from several WGDs. The data imply that simpler signalling pathways in the pre-WGD ancestors were converted via WGDs into multi-stranded parallelised networks. Genetic and biochemical studies in plants, yeasts and vertebrates suggest a paradigm in which different combinations of sister paralogues in the post-WGD regulatory networks are co-regulated under different conditions. In principle, such networks can respond to a wide array of environmental, sensory and hormonal stimuli and integrate them to generate phenotypic variety in cell types and behaviours. Patterns are also being discerned in how the post-WGD signalling networks are reconfigured in human cancers and neurological conditions. It is fascinating to unpick how ancient genomic events impact on complexity, variety and disease in modern life.

Evolution of the vertebrate insulin receptor substrate (Irs) gene family

BACKGROUND

Insulin receptor substrate (Irs) proteins are essential for insulin signaling as they allow downstream effectors to dock with, and be activated by, the insulin receptor. A family of four Irs proteins have been identified in mice, however the gene for one of these, IRS3, has been pseudogenized in humans. While it is known that the Irs gene family originated in vertebrates, it is not known when it originated and which members are most closely related to each other. A better understanding of the evolution of Irs genes and proteins should provide insight into the regulation of metabolism by insulin.

RESULTS

Multiple genes for Irs proteins were identified in a wide variety of vertebrate species. Phylogenetic and genomic neighborhood analyses indicate that this gene family originated very early in vertebrae evolution. Most Irs genes were duplicated and retained in fish after the fish-specific genome duplication. Irs genes have been lost of various lineages, including Irs3 in primates and birds and Irs1 in most fish. Irs3 and Irs4 experienced an episode of more rapid protein sequence evolution on the ancestral mammalian lineage. Comparisons of the conservation of the proteins sequences among Irs paralogs show that domains involved in binding to the plasma membrane and insulin receptors are most strongly conserved, while divergence has occurred in sequences involved in interacting with downstream effector proteins.

CONCLUSIONS

The Irs gene family originated very early in vertebrate evolution, likely through genome duplications, and in parallel with duplications of other components of the insulin signaling pathway, including insulin and the insulin receptor. While the N-terminal sequences of these proteins are conserved among the paralogs, changes in the C-terminal sequences likely allowed changes in biological function.

Lineage-specific rediploidization is a mechanism to explain time-lags between genome duplication and evolutionary diversification

BACKGROUND

The functional divergence of duplicate genes (ohnologues) retained from whole genome duplication (WGD) is thought to promote evolutionary diversification. However, species radiation and phenotypic diversification are often temporally separated from WGD. Salmonid fish, whose ancestor underwent WGD by autotetraploidization ~95 million years ago, fit such a 'time-lag' model of post-WGD radiation, which occurred alongside a major delay in the rediploidization process. Here we propose a model, 'lineage-specific ohnologue resolution' (LORe), to address the consequences of delayed rediploidization. Under LORe, speciation precedes rediploidization, allowing independent ohnologue divergence in sister lineages sharing an ancestral WGD event.

RESULTS

Using cross-species sequence capture, phylogenomics and genome-wide analyses of ohnologue expression divergence, we demonstrate the major impact of LORe on salmonid evolution. One-quarter of each salmonid genome, harbouring at least 4550 ohnologues, has evolved under LORe, with rediploidization and functional divergence occurring on multiple independent occasions >50 million years post-WGD. We demonstrate the existence and regulatory divergence of many LORe ohnologues with functions in lineage-specific physiological adaptations that potentially facilitated salmonid species radiation. We show that LORe ohnologues are enriched for different functions than 'older' ohnologues that began diverging in the salmonid ancestor.

CONCLUSIONS

LORe has unappreciated significance as a nested component of post-WGD divergence that impacts the functional properties of genes, whilst providing ohnologues available solely for lineage-specific adaptation. Under LORe, which is predicted following many WGD events, the functional outcomes of WGD need not appear 'explosively', but can arise gradually over tens of millions of years, promoting lineage-specific diversification regimes under prevailing ecological pressures.

Evolution of complexity in the zebrafish synapse proteome

The proteome of human brain synapses is highly complex and is mutated in over 130 diseases. This complexity arose from two whole-genome duplications early in the vertebrate lineage. Zebrafish are used in modelling human diseases; however, its synapse proteome is uncharacterized, and whether the teleost-specific genome duplication (TSGD) influenced complexity is unknown. We report the characterization of the proteomes and ultrastructure of central synapses in zebrafish and analyse the importance of the TSGD. While the TSGD increases overall synapse proteome complexity, the postsynaptic density (PSD) proteome of zebrafish has lower complexity than mammals. A highly conserved set of ∼1,000 proteins is shared across vertebrates. PSD ultrastructural features are also conserved. Lineage-specific proteome differences indicate that vertebrate species evolved distinct synapse types and functions. The data sets are a resource for a wide range of studies and have important implications for the use of zebrafish in modelling human synaptic diseases.

Robust Yet Fragile: Expression Noise, Protein Misfolding, and Gene Dosage in the Evolution of Genomes

The complex manner in which organisms respond to changes in their gene dosage has long fascinated geneticists. Oddly, although the existence of dominance implies that dosage reductions often have mild phenotypes, extra copies of whole chromosomes (aneuploidy) are generally strongly deleterious. Even more paradoxically, an extra copy of the genome is better tolerated than is aneuploidy. We review the resolution of this paradox, highlighting the roles of biochemistry, protein aggregation, and disruption of cellular microstructure in that explanation. Returning to life's curious combination of robustness and sensitivity to dosage changes, we argue that understanding how biological robustness evolved makes these observations less inexplicable. We propose that noise in gene expression and evolutionary strategies for its suppression play a role in generating dosage phenotypes. Finally, we outline an unappreciated mechanism for the preservation of duplicate genes, namely preservation to limit expression noise, arguing that it is particularly relevant in polyploid organisms.

AMP-activated protein kinase: a cellular energy sensor that comes in 12 flavours

The AMP-activated protein kinase (AMPK) is a sensor of cellular energy status that is expressed in essentially all eukaryotic cells, suggesting that it arose during early eukaryotic evolution. It occurs universally as heterotrimeric complexes containing catalytic α subunits and regulatory β and γ subunits. Although Drosophila melanogaster contains single genes encoding each subunit, in mammals, each subunit exists as multiple isoforms encoded by distinct genes, giving rise to up to 12 heterotrimeric combinations. The multiple isoforms of each subunit are 2R-ohnologues generated by the two rounds of whole genome duplication that occurred at the evolutionary origin of the vertebrates. Although the differential roles of these isoform combinations remain only partly understood, there are indications that they may have different subcellular locations, different inputs and outputs, and different functions. The multiple isoforms are of particular interest with respect to the roles of AMPK in cancer because the genes encoding some isoforms, such as PRKAA1 and PRKAB2 (encoding α1 and β2), are quite frequently amplified in tumour cells, whereas the genes encoding others, such as PRKAA2 (encoding α2), tend to be mutated, which, in some but not all cases, may result in a loss of function. Thus, although AMPK acts downstream of the tumour suppressor liver kinase B1, and some of its isoform combinations may act as tumour suppressors that restrain the growth and proliferation of tumour cells, other isoform combinations may paradoxically act as oncogenes, perhaps by aiding the survival of tumour cells undergoing environmental stresses such as hypoxia or nutrient deprivation.

Spatial Colocalization of Human Ohnolog Pairs Acts to Maintain Dosage-Balance

Ohnologs -paralogous gene pairs generated by whole genome duplication- are enriched for dosage sensitive genes, that is, genes that have a phenotype due to copy number changes. Dosage sensitive genes frequently occur in the same metabolic pathway and in physically interacting proteins. Accumulating evidence reveals that functionally related genes tend to co-localize in the three-dimensional (3D) arrangement of chromosomes. We query whether the spatial distribution of ohnologs has implications for their dosage balance. We analyzed the colocalization frequency of ohnologs based on chromatin interaction datasets of seven human cell lines and found that ohnolog pairs exhibit higher spatial proximity in 3D nuclear organization than other paralog pairs and than randomly chosen ohnologs in the genome. We also found that colocalized ohnologs are more resistant to copy number variations and more likely to be disease-associated genes, which indicates a stronger dosage balance in ohnologs with high spatial proximity. This phenomenon is further supported by the stronger similarity of gene co-expression and of gene ontology terms of colocalized ohnologs. In addition, for a large fraction of ohnologs, the spatial colocalization is conserved in mouse cells, suggestive of functional constraint on their 3D positioning in the nucleus.

Whole Genome Duplications Shaped the Receptor Tyrosine Kinase Repertoire of Jawed Vertebrates

The receptor tyrosine kinase (RTK) gene family, involved primarily in cell growth and differentiation, comprises proteins with a common enzymatic tyrosine kinase intracellular domain adjacent to a transmembrane region. The amino-terminal portion of RTKs is extracellular and made of different domains, the combination of which characterizes each of the 20 RTK subfamilies among mammals. We analyzed a total of 7,376 RTK sequences among 143 vertebrate species to provide here the first comprehensive census of the jawed vertebrate repertoire. We ascertained the 58 genes previously described in the human and mouse genomes and established their phylogenetic relationships. We also identified five additional RTKs amounting to a total of 63 genes in jawed vertebrates. We found that the vertebrate RTK gene family has been shaped by the two successive rounds of whole genome duplications (WGD) called 1R and 2R (1R/2R) that occurred at the base of the vertebrates. In addition, the Vegfr and Ephrin receptor subfamilies were expanded by single gene duplications. In teleost fish, 23 additional RTK genes have been retained after another expansion through the fish-specific third round (3R) of WGD. Several lineage-specific gene losses were observed. For instance, birds have lost three RTKs, and different genes are missing in several fish sublineages. The RTK gene family presents an unusual high gene retention rate from the vertebrate WGDs (58.75% after 1R/2R, 64.4% after 3R), resulting in an expansion that might be correlated with the evolution of complexity of vertebrate cellular communication and intracellular signaling.

Evolution and functions of Oct4 homologs in non-mammalian vertebrates

PouV class transcription factor Oct4/Pou5f1 is a central regulator of indefinite pluripotency in mammalian embryonic stem cells (ESCs) but also participates in cell lineage specification in mouse embryos and in differentiating cell cultures. The molecular basis for this versatility, which is shared between Oct4 and its non-mammalian homologs Pou5f1 and Pou5f3, is not yet completely understood. Here, I review the current understanding of the evolution of PouV class transcription factors and discuss equivalent and diverse roles of Oct4 homologs in pluripotency, differentiation, and cell behavior in different vertebrate embryos. This article is part of a Special Issue entitled: The Oct Transcription Factor Family, edited by Dr. Dean Tantin.

Early origin and adaptive evolution of the GW182 protein family, the key component of RNA silencing in animals

The GW182 proteins are a key component of the miRNA-dependent post-transcriptional silencing pathway in animals. They function as scaffold proteins to mediate the interaction of Argonaute (AGO)-containing complexes with cytoplasmic poly(A)-binding proteins (PABP) and PAN2-PAN3 and CCR4-NOT deadenylases. The AGO-GW182 complexes mediate silencing of the target mRNA through induction of translational repression and/or mRNA degradation. Although the GW182 proteins are a subject of extensive experimental research in the recent years, very little is known about their origin and evolution. Here, based on complex functional annotation and phylogenetic analyses, we reveal 448 members of the GW182 protein family from the earliest animals to humans. Our results indicate that a single-copy GW182/TNRC6C progenitor gene arose with the emergence of multicellularity and it multiplied in the last common ancestor of vertebrates in 2 rounds of whole genome duplication (WGD) resulting in 3 genes. Before the divergence of vertebrates, both the AGO- and CCR4-NOT-binding regions of GW182s showed significant acceleration in the accumulation of amino acid changes, suggesting functional adaptation toward higher specificity to the molecules of the silencing complex. We conclude that the silencing ability of the GW182 proteins improves with higher position in the taxonomic classification and increasing complexity of the organism. The first reconstruction of the molecular journey of GW182 proteins from the ancestral metazoan protein to the current mammalian configuration provides new insight into development of the miRNA-dependent post-transcriptional silencing pathway in animals.

Gene Duplicability of Core Genes Is Highly Consistent across All Angiosperms

Gene duplication is an important mechanism for adding to genomic novelty. Hence, which genes undergo duplication and are preserved following duplication is an important question. It has been observed that gene duplicability, or the ability of genes to be retained following duplication, is a nonrandom process, with certain genes being more amenable to survive duplication events than others. Primarily, gene essentiality and the type of duplication (small-scale versus large-scale) have been shown in different species to influence the (long-term) survival of novel genes. However, an overarching view of "gene duplicability" is lacking, mainly due to the fact that previous studies usually focused on individual species and did not account for the influence of genomic context and the time of duplication. Here, we present a large-scale study in which we investigated duplicate retention for 9178 gene families shared between 37 flowering plant species, referred to as angiosperm core gene families. For most gene families, we observe a strikingly consistent pattern of gene duplicability across species, with gene families being either primarily single-copy or multicopy in all species. An intermediate class contains gene families that are often retained in duplicate for periods extending to tens of millions of years after whole-genome duplication, but ultimately appear to be largely restored to singleton status, suggesting that these genes may be dosage balance sensitive. The distinction between single-copy and multicopy gene families is reflected in their functional annotation, with single-copy genes being mainly involved in the maintenance of genome stability and organelle function and multicopy genes in signaling, transport, and metabolism. The intermediate class was overrepresented in regulatory genes, further suggesting that these represent putative dosage-balance-sensitive genes.

The Constrained Maximal Expression Level Owing to Haploidy Shapes Gene Content on the Mammalian X Chromosome

X chromosomes are unusual in many regards, not least of which is their nonrandom gene content. The causes of this bias are commonly discussed in the context of sexual antagonism and the avoidance of activity in the male germline. Here, we examine the notion that, at least in some taxa, functionally biased gene content may more profoundly be shaped by limits imposed on gene expression owing to haploid expression of the X chromosome. Notably, if the X, as in primates, is transcribed at rates comparable to the ancestral rate (per promoter) prior to the X chromosome formation, then the X is not a tolerable environment for genes with very high maximal net levels of expression, owing to transcriptional traffic jams. We test this hypothesis using The Encyclopedia of DNA Elements (ENCODE) and data from the Functional Annotation of the Mammalian Genome (FANTOM5) project. As predicted, the maximal expression of human X-linked genes is much lower than that of genes on autosomes: on average, maximal expression is three times lower on the X chromosome than on autosomes. Similarly, autosome-to-X retroposition events are associated with lower maximal expression of retrogenes on the X than seen for X-to-autosome retrogenes on autosomes. Also as expected, X-linked genes have a lesser degree of increase in gene expression than autosomal ones (compared to the human/Chimpanzee common ancestor) if highly expressed, but not if lowly expressed. The traffic jam model also explains the known lower breadth of expression for genes on the X (and the Z of birds), as genes with broad expression are, on average, those with high maximal expression. As then further predicted, highly expressed tissue-specific genes are also rare on the X and broadly expressed genes on the X tend to be lowly expressed, both indicating that the trend is shaped by the maximal expression level not the breadth of expression per se. Importantly, a limit to the maximal expression level explains biased tissue of expression profiles of X-linked genes. Tissues whose tissue-specific genes are very highly expressed (e.g., secretory tissues, tissues abundant in structural proteins) are also tissues in which gene expression is relatively rare on the X chromosome. These trends cannot be fully accounted for in terms of alternative models of biased expression. In conclusion, the notion that it is hard for genes on the Therian X to be highly expressed, owing to transcriptional traffic jams, provides a simple yet robustly supported rationale of many peculiar features of X’s gene content, gene expression, and evolution.

Laurence Hurst, Lukasz Huminiecki, and the FANTOM5 consortium propose a new explanation for the peculiar expression properties of genes on the human X chromosome, based on the premise that very high expression levels cannot be achieved on a haploid-expressed chromosome.

Genes located on the human X chromosome are not a random mix of genes: they tend to be expressed in relatively few tissues or are specific for a particular set of tissues, e.g., brain regions. Prior attempts to explain this skewed gene content have hypothesized that the X chromosome might be peculiar because it has to balance mutations that are advantageous to one sex but deleterious to the other, or because it has to shut down during the process of sperm manufacture in males. Here we suggest and test a third possible explanation: that genes on the X chromosome are limited in their transcription levels and thus tend to be genes that are lowly or specifically expressed. We consider the suggestion that since these genes can only be expressed from one chromosome, as males only have one X, the ability to express a gene at very high rates is limited owing to potential transcriptional traffic jams. As predicted, we find that human X-located genes have maximal expression rates far below that of genes residing on autosomes. When we look at genes that have moved onto or off the X chromosome during recent evolution, we find the maximal expression is higher when not on the X chromosome. We also find that X-located genes that are relatively highly expressed are not able to increase their expression level further. Our model explains both the enrichment for tissue specificity and the paucity of certain tissues with X-located genes. Genes underrepresented on the X are either expressed in many tissues—such genes tend to have high maximal expression—or are from tissues that require a lot of transcription (e.g., fast secreting tissues like the liver). Just as many of the findings cannot be explained by the two earlier models, neither can the traffic jam model explain all the peculiar features of the genes found on the X chromosome. Indeed, we find evidence of a reproduction-related bias in X-located genes, even after allowing for the traffic jam problem.

Dissecting signalling by individual Akt/PKB isoforms, three steps at once

The serine/threonine kinase Akt/PKB (protein kinase B) is key for mammalian cell growth, survival, metabolism and oncogenic transformation. The diverse level and tissue expression of its three isoforms, Akt1/PKBα, Akt2/PKBβ and Akt3/PKBγ, make it daunting to identify isoform-specific actions in vivo and even in isolated tissues/cells. To date, isoform-specific knockout and knockdown have been the best strategies to dissect their individual overall functions. In a recent article in the Biochemical Journal, Kajno et al. reported a new strategy to study isoform selectivity in cell lines. Individual Akt/PKB isoforms in 3T3-L1 pre-adipocytes are first silenced via shRNA and stable cellular clones lacking one or the other isoform are selected. The stably silenced isoform is then replaced by a mutant engineered to be refractory to inhibition by MK-2206 (Akt1(W80A) or Akt2(W80A)). Akt1(W80A) or Akt2(W80A) are functional and effectively recruited to the plasma membrane in response to insulin. The system affords the opportunity to acutely control the activity of the endogenous non-silenced isoform through timely addition of MK-2206. Using this approach, it is confirmed that Akt1/PKBα is the preferred isoform sustaining adipocyte differentiation, but both Akt1/PKBα and Akt2/PKBβ can indistinctly support insulin-dependent FoxO1 (forkhead box O1) nuclear exclusion. Surprisingly, either isoform can also support insulin-dependent glucose transporter (GLUT) 4 translocation to the membrane, in contrast with the preferential role of Akt2/PKBβ assessed by knockdown studies. The new strategy should allow analysis of the plurality of Akt/PKB functions in other cells and in response to other stimuli. It should also be amenable to high-throughput studies to speed up advances in signal transmission by this pivotal kinase.

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

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

Identification of Ohnolog Genes Originating from Whole Genome Duplication in Early Vertebrates, Based on Synteny Comparison across Multiple Genomes

Whole genome duplications (WGD) have now been firmly established in all major eukaryotic kingdoms. In particular, all vertebrates descend from two rounds of WGDs, that occurred in their jawless ancestor some 500 MY ago. Paralogs retained from WGD, also coined 'ohnologs' after Susumu Ohno, have been shown to be typically associated with development, signaling and gene regulation. Ohnologs, which amount to about 20 to 35% of genes in the human genome, have also been shown to be prone to dominant deleterious mutations and frequently implicated in cancer and genetic diseases. Hence, identifying ohnologs is central to better understand the evolution of vertebrates and their susceptibility to genetic diseases. Early computational analyses to identify vertebrate ohnologs relied on content-based synteny comparisons between the human genome and a single invertebrate outgroup genome or within the human genome itself. These approaches are thus limited by lineage specific rearrangements in individual genomes. We report, in this study, the identification of vertebrate ohnologs based on the quantitative assessment and integration of synteny conservation between six amniote vertebrates and six invertebrate outgroups. Such a synteny comparison across multiple genomes is shown to enhance the statistical power of ohnolog identification in vertebrates compared to earlier approaches, by overcoming lineage specific genome rearrangements. Ohnolog gene families can be browsed and downloaded for three statistical confidence levels or recompiled for specific, user-defined, significance criteria at http://ohnologs.curie.fr/. In the light of the importance of WGD on the genetic makeup of vertebrates, our analysis provides a useful resource for researchers interested in gaining further insights on vertebrate evolution and genetic diseases.

14-3-3-Pred: improved methods to predict 14-3-3-binding phosphopeptides

MOTIVATION

The 14-3-3 family of phosphoprotein-binding proteins regulates many cellular processes by docking onto pairs of phosphorylated Ser and Thr residues in a constellation of intracellular targets. Therefore, there is a pressing need to develop new prediction methods that use an updated set of 14-3-3-binding motifs for the identification of new 14-3-3 targets and to prioritize the downstream analysis of >2000 potential interactors identified in high-throughput experiments.

RESULTS

Here, a comprehensive set of 14-3-3-binding targets from the literature was used to develop 14-3-3-binding phosphosite predictors. Position-specific scoring matrix, support vector machines (SVM) and artificial neural network (ANN) classification methods were trained to discriminate experimentally determined 14-3-3-binding motifs from non-binding phosphopeptides. ANN, position-specific scoring matrix and SVM methods showed best performance for a motif window spanning from -6 to +4 around the binding phosphosite, achieving Matthews correlation coefficient of up to 0.60. Blind prediction showed that all three methods outperform two popular 14-3-3-binding site predictors, Scansite and ELM. The new methods were used for prediction of 14-3-3-binding phosphosites in the human proteome. Experimental analysis of high-scoring predictions in the FAM122A and FAM122B proteins confirms the predictions and suggests the new 14-3-3-predictors will be generally useful.

AVAILABILITY AND IMPLEMENTATION

A standalone prediction web server is available at http://www.compbio.dundee.ac.uk/1433pred. Human candidate 14-3-3-binding phosphosites were integrated in ANIA: ANnotation and Integrated Analysis of the 14-3-3 interactome database.

AnimalTFDB 2.0: a resource for expression, prediction and functional study of animal transcription factors

Transcription factors (TFs) are key regulators for gene expression. Here we updated the animal TF database AnimalTFDB to version 2.0 (http://bioinfo.life.hust.edu.cn/AnimalTFDB/). Using the improved prediction pipeline, we identified 72 336 TF genes, 21 053 transcription co-factor genes and 6502 chromatin remodeling factor genes from 65 species covering main animal lineages. Besides the abundant annotations (basic information, gene model, protein functional domain, gene ontology, pathway, protein interaction, ortholog and paralog, etc.) in the previous version, we made several new features and functions in the updated version. These new features are: (i) gene expression from RNA-Seq for nine model species, (ii) gene phenotype information, (iii) multiple sequence alignment of TF DNA-binding domains, and the weblogo and phylogenetic tree based on the alignment, (iv) a TF prediction server to identify new TFs from input sequences and (v) a BLAST server to search against TFs in AnimalTFDB. A new nice web interface was designed for AnimalTFDB 2.0 allowing users to browse and search all data in the database. We aim to maintain the AnimalTFDB as a solid resource for TF identification and studies of transcription regulation and comparative genomics.

A simple metric of promoter architecture robustly predicts expression breadth of human genes suggesting that most transcription factors are positive regulators

BACKGROUND

Conventional wisdom holds that, owing to the dominance of features such as chromatin level control, the expression of a gene cannot be readily predicted from knowledge of promoter architecture. This is reflected, for example, in a weak or absent correlation between promoter divergence and expression divergence between paralogs. However, an inability to predict may reflect an inability to accurately measure or employment of the wrong parameters. Here we address this issue through integration of two exceptional resources: ENCODE data on transcription factor binding and the FANTOM5 high-resolution expression atlas.

RESULTS

Consistent with the notion that in eukaryotes most transcription factors are activating, the number of transcription factors binding a promoter is a strong predictor of expression breadth. In addition, evolutionarily young duplicates have fewer transcription factor binders and narrower expression. Nonetheless, we find several binders and cooperative sets that are disproportionately associated with broad expression, indicating that models more complex than simple correlations should hold more predictive power. Indeed, a machine learning approach improves fit to the data compared with a simple correlation. Machine learning could at best moderately predict tissue of expression of tissue specific genes.

CONCLUSIONS

We find robust evidence that some expression parameters and paralog expression divergence are strongly predictable with knowledge of transcription factor binding repertoire. While some cooperative complexes can be identified, consistent with the notion that most eukaryotic transcription factors are activating, a simple predictor, the number of binding transcription factors found on a promoter, is a robust predictor of expression breadth.

Identification of 2R-ohnologue gene families displaying the same mutation-load skew in multiple cancers

The complexity of signalling pathways was boosted at the origin of the vertebrates, when two rounds of whole genome duplication (2R-WGD) occurred. Those genes and proteins that have survived from the 2R-WGD-termed 2R-ohnologues-belong to families of two to four members, and are enriched in signalling components relevant to cancer. Here, we find that while only approximately 30% of human transcript-coding genes are 2R-ohnologues, they carry 42-60% of the gene mutations in 30 different cancer types. Across a subset of cancer datasets, including melanoma, breast, lung adenocarcinoma, liver and medulloblastoma, we identified 673 2R-ohnologue families in which one gene carries mutations at multiple positions, while sister genes in the same family are relatively mutation free. Strikingly, in 315 of the 322 2R-ohnologue families displaying such a skew in multiple cancers, the same gene carries the heaviest mutation load in each cancer, and usually the second-ranked gene is also the same in each cancer. Our findings inspire the hypothesis that in certain cancers, heterogeneous combinations of genetic changes impair parts of the 2R-WGD signalling networks and force information flow through a limited set of oncogenic pathways in which specific non-mutated 2R-ohnologues serve as effectors. The non-mutated 2R-ohnologues are therefore potential therapeutic targets. These include proteins linked to growth factor signalling, neurotransmission and ion channels.

A synergism between adaptive effects and evolvability drives whole genome duplication to fixation

Whole genome duplication has shaped eukaryotic evolutionary history and has been associated with drastic environmental change and species radiation. While the most common fate of WGD duplicates is a return to single copy, retained duplicates have been found enriched for highly interacting genes. This pattern has been explained by a neutral process of subfunctionalization and more recently, dosage balance selection. However, much about the relationship between environmental change, WGD and adaptation remains unknown. Here, we study the duplicate retention pattern postWGD, by letting virtual cells adapt to environmental changes. The virtual cells have structured genomes that encode a regulatory network and simple metabolism. Populations are under selection for homeostasis and evolve by point mutations, small indels and WGD. After populations had initially adapted fully to fluctuating resource conditions re-adaptation to a broad range of novel environments was studied by tracking mutations in the line of descent. WGD was established in a minority (≈30%) of lineages, yet, these were significantly more successful at re-adaptation. Unexpectedly, WGD lineages conserved more seemingly redundant genes, yet had higher per gene mutation rates. While WGD duplicates of all functional classes were significantly over-retained compared to a model of neutral losses, duplicate retention was clearly biased towards highly connected TFs. Importantly, no subfunctionalization occurred in conserved pairs, strongly suggesting that dosage balance shaped retention. Meanwhile, singles diverged significantly. WGD, therefore, is a powerful mechanism to cope with environmental change, allowing conservation of a core machinery, while adapting the peripheral network to accommodate change.

The evolution of eukaryotes is characterized by drastic changes in their genome content. Genome expansions have often occurred by duplication of the entire genome. It is generally not know whether organisms gain any adaptive advantage from these mutations. However, they appear to become fixed in response to environmental change. Many interesting whole genome duplications happened long ago in eukaryotic evolutionary history during periods of turbulent genome and species evolution. Genomic data analysis alone cannot resolve the evolutionary mechanisms and consequences of whole genome duplication. Here, we modeled evolution with whole genome duplications in a Virtual Cell model. Simulating populations that undergo a range of different environmental changes we found that next to often increasing fitness directly, whole genome duplications made lineages more evolvable and hence more able to adapt to harsh new environments. Although most duplicates are deleted in subsequent evolution, genes with many interaction partners were retained preferentially, increasing regulatory complexity. Interestingly however, we found that innovation happened most likely in the more loosely connected and less essential genes.

On the retention of gene duplicates prone to dominant deleterious mutations

Recent studies have shown that gene families from different functional categories have been preferentially expanded either by small scale duplication (SSD) or by whole-genome duplication (WGD). In particular, gene families prone to dominant deleterious mutations and implicated in cancers and other genetic diseases in human have been greatly expanded through two rounds of WGD dating back from early vertebrates. Here, we strengthen this intriguing observation, showing that human oncogenes involved in different primary tumors have retained many WGD duplicates compared to other human genes. In order to rationalize this evolutionary outcome, we propose a consistent population genetics model to analyze the retention of SSD and WGD duplicates taking into account their propensity to acquire dominant deleterious mutations. We solve a deterministic haploid model including initial duplicated loci, their retention through sub-functionalization or their neutral loss-of-function or deleterious gain-of-function at one locus. Extensions to diploid genotypes are presented and population size effects are analyzed using stochastic simulations. The only difference between the SSD and WGD scenarios is the initial number of individuals with duplicated loci. While SSD duplicates need to spread through the entire population from a single individual to reach fixation, WGD duplicates are de facto fixed in the small initial post-WGD population arising through the ploidy incompatibility between post-WGD individuals and the rest of the pre-WGD population. WGD duplicates prone to dominant deleterious mutations are then shown to be indirectly selected through purifying selection in post-WGD species, whereas SSD duplicates typically require positive selection. These results highlight the long-term evolution mechanisms behind the surprising accumulation of WGD duplicates prone to dominant deleterious mutations and are shown to be consistent with cancer genome data on the prevalence of human oncogenes with WGD duplicates.

ANIA: ANnotation and Integrated Analysis of the 14-3-3 interactome

The dimeric 14-3-3 proteins dock onto pairs of phosphorylated Ser and Thr residues on hundreds of proteins, and thereby regulate many events in mammalian cells. To facilitate global analyses of these interactions, we developed a web resource named ANIA: ANnotation and Integrated Analysis of the 14-3-3 interactome, which integrates multiple data sets on 14-3-3-binding phosphoproteins. ANIA also pinpoints candidate 14-3-3-binding phosphosites using predictor algorithms, assisted by our recent discovery that the human 14-3-3-interactome is highly enriched in 2R-ohnologues. 2R-ohnologues are proteins in families of two to four, generated by two rounds of whole genome duplication at the origin of the vertebrate animals. ANIA identifies candidate ‘lynchpins’, which are 14-3-3-binding phosphosites that are conserved across members of a given 2R-ohnologue protein family. Other features of ANIA include a link to the catalogue of somatic mutations in cancer database to find cancer polymorphisms that map to 14-3-3-binding phosphosites, which would be expected to interfere with 14-3-3 interactions. We used ANIA to map known and candidate 14-3-3-binding enzymes within the 2R-ohnologue complement of the human kinome. Our projections indicate that 14-3-3s dock onto many more human kinases than has been realized. Guided by ANIA, PAK4, 6 and 7 (p21-activated kinases 4, 6 and 7) were experimentally validated as a 2R-ohnologue family of 14-3-3-binding phosphoproteins. PAK4 binding to 14-3-3 is stimulated by phorbol ester, and involves the ‘lynchpin’ site phosphoSer99 and a major contribution from Ser181. In contrast, PAK6 and PAK7 display strong phorbol ester-independent binding to 14-3-3, with Ser113 critical for the interaction with PAK6. These data point to differential 14-3-3 regulation of PAKs in control of cell morphology.

Database URL: https://ania-1433.lifesci.dundee.ac.uk/prediction/webserver/index.py

Age distribution patterns of human gene families: divergent for Gene Ontology categories and concordant between different subcellular localizations

The age distribution of gene duplication events within the human genome exhibits two waves of duplications along with an ancient component. However, because of functional constraint differences, genes in different functional categories might show dissimilar retention patterns after duplication. It is known that genes in some functional categories are highly duplicated in the early stage of vertebrate evolution. However, the correlations of the age distribution pattern of gene duplication between the different functional categories are still unknown. To investigate this issue, we developed a robust pipeline to date the gene duplication events in the human genome. We successfully estimated about three-quarters of the duplication events within the human genome, along with the age distribution pattern in each Gene Ontology (GO) slim category. We found that some GO slim categories show different distribution patterns when compared to the whole genome. Further hierarchical clustering of the GO slim functional categories enabled grouping into two main clusters. We found that human genes located in the duplicated copy number variant regions, whose duplicate genes have not been fixed in the human population, were mainly enriched in the groups with a high proportion of recently duplicated genes. Moreover, we used a phylogenetic tree-based method to date the age of duplications in three signaling-related gene superfamilies: transcription factors, protein kinases and G-protein coupled receptors. These superfamilies were expressed in different subcellular localizations. They showed a similar age distribution as the signaling-related GO slim categories. We also compared the differences between the age distributions of gene duplications in multiple subcellular localizations. We found that the distribution patterns of the major subcellular localizations were similar to that of the whole genome. This study revealed the whole picture of the evolution patterns of gene functional categories in the human genome.

The evolutionary analysis reveals domain fusion of proteins with Frizzled-like CRD domain

Frizzleds (FZDs) are transmembrane receptors in the Wnt signaling pathway and they play pivotal roles in developments. The Frizzled-like extracellular Cysteine-rich domain (Fz-CRD) has been identified in FZDs and other proteins. The origin and evolution of these proteins with Fz-CRD is the main interest of this study. We found that the Fz-CRD exists in FZD, SFRP, RTK, MFRP, CPZ, CORIN, COL18A1 and other proteins. Our systematic analysis revealed that the Fz-CRD domain might have originated in protists and then fused with the Frizzled-like seven-transmembrane domain (7TM) to form the FZD receptors, which duplicated and diversified into about 11 members in Vertebrates. The SFRPs and RTKs with the Fz-CRD were found in sponge and expanded in Vertebrates. Other proteins with Fz-CRD may have emerged during Vertebrate evolution through domain fusion. Moreover, we found a glycosylation site and several conserved motifs in FZDs, which may be related to Wnt interaction. Based on these results, we proposed a model showing that the domain fusion and expansion of Fz-CRD genes occurred in Metazoa and Vertebrates. Our study may help to pave the way for further research on the conservation and diversification of Wnt signaling functions during evolution.

Evolution of signal multiplexing by 14-3-3-binding 2R-ohnologue protein families in the vertebrates

14-3-3 proteins regulate cellular responses to stimuli by docking onto pairs of phosphorylated residues on target proteins. The present study shows that the human 14-3-3-binding phosphoproteome is highly enriched in 2R-ohnologues, which are proteins in families of two to four members that were generated by two rounds of whole genome duplication at the origin of the vertebrates. We identify 2R-ohnologue families whose members share a ‘lynchpin’, defined as a 14-3-3-binding phosphosite that is conserved across members of a given family, and aligns with a Ser/Thr residue in pro-orthologues from the invertebrate chordates. For example, the human receptor expression enhancing protein (REEP) 1–4 family has the commonest type of lynchpin motif in current datasets, with a phosphorylatable serine in the –2 position relative to the 14-3-3-binding phosphosite. In contrast, the second 14-3-3-binding sites of REEPs 1–4 differ and are phosphorylated by different kinases, and hence the REEPs display different affinities for 14-3-3 dimers. We suggest a conceptual model for intracellular regulation involving protein families whose evolution into signal multiplexing systems was facilitated by 14-3-3 dimer binding to lynchpins, which gave freedom for other regulatory sites to evolve. While increased signalling complexity was needed for vertebrate life, these systems also generate vulnerability to genetic disorders.

Watching the grin fade: tracing the effects of polyploidy on different evolutionary time scales

Polyploidy, or whole-genome duplication (WGD), is a recurrent mutation both in cell lineages and over evolutionary time. By globally changing the relationship between gene copy number and other cellular entities, it can induce dramatic changes at the cellular and phenotypic level. Perhaps surprisingly, then, the insights that these events can bring to understanding other cellular features are not as well appreciated as they could be. In this review, we draw on examples of polyploidy from animals, plants and yeast to explore how investigations of polyploid cells have improved our understanding of the cell cycle, biological network complexity, metabolic phenotypes and tumor biology. We argue that the study of polyploidy across organisms, cell types, and time scales serves not only as a window into basic cell biology, but also as a basis for a predictive biology with applications ranging from crop improvement to treating cancer.