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

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.

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.

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.

Systematic analysis of the Hippo pathway organization and oncogenic alteration in evolution

The Hippo pathway is a central regulator of organ size and a key tumor suppressor via coordinating cell proliferation and death. Initially discovered in Drosophila, the Hippo pathway has been implicated as an evolutionarily conserved pathway in mammals; however, how this pathway was evolved to be functional from its origin is still largely unknown. In this study, we traced the Hippo pathway in premetazoan species, characterized the intrinsic functions of its ancestor components, and unveiled the evolutionary history of this key signaling pathway from its unicellular origin. In addition, we elucidated the paralogous gene history for the mammalian Hippo pathway components and characterized their cancer-derived somatic mutations from an evolutionary perspective. Taken together, our findings not only traced the conserved function of the Hippo pathway to its unicellular ancestor components, but also provided novel evolutionary insights into the Hippo pathway organization and oncogenic alteration.

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.

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.

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.

Insights into Biological Complexity from Simple Foundations

We discuss an overtly "simple approach" to complex biological systems borrowing selectively from theoretical physics. The approach is framed by three maxims, and we show examples of its success in two different applications: investigating cellular robustness at the level of gene regulatory networks and quantifying rare events of DNA replication errors.

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.