Teleost fish with specific genome duplication as unique models of vertebrate evolution

Whole-genome duplication (WGD) is believed to be one of the major evolutionary events that shaped the genome organization of vertebrates. Here, we review recent research on vertebrate genome evolution, specifically on WGD and its consequences for gene and genome evolution in teleost fishes. Recent genome analyses confirmed that all vertebrates experienced two rounds of WGD early in their evolution, and that teleosts experienced a subsequent additional third-round (3R)-WGD. The 3R-WGD was estimated to have occurred 320–400 million years ago in a teleost ancestor, but after its divergence from a common ancestor with living non-teleost actinopterygians (Bichir, Sturgeon, Bowfin, and Gar) based on the analyses of teleost-specific duplicate genes. This 3R-WGD was confirmed by synteny analysis and ancestral karyotype inference using the genome sequences of Tetraodon and medaka. Most of the tetrapods, on the other hand, have not experienced an additional WGD; however, they have experienced repeated chromosomal rearrangements throughout the whole genome. Therefore, different types of chromosomal events have characterized the genomes of teleosts and tetrapods, respectively. The 3R-WGD is useful to investigate the consequences of WGD because it is an evolutionarily recent WGD and thus teleost genomes retain many more WGD-derived duplicates and “traces” of their evolution. In addition, the remarkable morphological, physiological, and ecological diversity of teleosts may facilitate understanding of macrophenotypic evolution on the basis of genetic/genomic information. We highlight the teleosts with 3R-WGD as unique models for future studies on ecology and evolution taking advantage of emerging genomics technologies and systems biology environments.     

GPCR genes are preferentially retained after whole genome duplication

One of the most interesting questions in biology is whether certain pathways have been favored during evolution, and if so, what properties could cause such a preference. Due to the lack of experimental evidence, whether select gene families have been preferentially retained over time after duplication in metazoan organisms remains unclear. Here, by syntenic mapping of nonchemosensory G protein-coupled receptor genes (nGPCRs which represent half the receptome for transmembrane signaling) in the vertebrate genomes, we found that, as opposed to the 8-15% retention rate for whole genome duplication (WGD)-derived gene duplicates in the entire genome of pufferfish, greater than 27.8% of WGD-derived nGPCRs which interact with a nonpeptide ligand were retained after WGD in pufferfish Tetraodon nigroviridis. In addition, we show that concurrent duplication of cognate ligand genes by WGD could impose selection of nGPCRs that interact with a polypeptide ligand. Against less than 2.25% probability for parallel retention of a pair of WGD-derived ligands and a pair of cognate receptor duplicates, we found a more than 8.9% retention of WGD-derived ligand-nGPCR pairs--threefold greater than one would surmise. These results demonstrate that gene retention is not uniform after WGD in vertebrates, and suggest a Darwinian selection of GPCR-mediated intercellular communication in metazoan organisms.     

Functional analysis of gene duplications in Saccharomyces cerevisiae

Gene duplication can occur on two scales: whole-genome duplications (WGD) and smaller-scale duplications (SSD) involving individual genes or genomic segments. Duplication may result in functionally redundant genes or diverge in function through neofunctionalization or subfunctionalization. The effect of duplication scale on functional evolution has not yet been explored, probably due to the lack of global knowledge of protein function and different times of duplication events. To address this question, we used integrated Bayesian analysis of diverse functional genomic data to accurately evaluate the extent of functional similarity and divergence between paralogs on a global scale. We found that paralogs resulting from the whole-genome duplication are more likely to share interaction partners and biological functions than smaller-scale duplicates, independent of sequence similarity. In addition, WGD paralogs show lower frequency of essential genes and higher synthetic lethality rate, but instead diverge more in expression pattern and upstream regulatory region. Thus, our analysis demonstrates that WGD paralogs generally have similar compensatory functions but diverging expression patterns, suggesting a potential of distinct evolutionary scenarios for paralogs that arose through different duplication mechanisms. Furthermore, by identifying these functional disparities between the two types of duplicates, we reconcile previous disputes on the relationship between sequence divergence and expression divergence or essentiality.     

Functional divergence of Populus MYB158 and MYB189 gene pair created by whole genome duplication

Whole genome duplication (WGD) providesnew genetic material for genome evolution. After a WGD event, some duplicates are lost, while other duplicates still persist and evolve diverse functions. A particular challenge is to understand how this diversity arises. This study identified two WGD-derived duplicates, MYB158 and MYB189, from Populus tomentosa. Populus MYB158 and MYB189 had expression divergence. Populus tomentosa overexpressing MYB158 or MYB189 had similar phenotypes: creep growth, decreased width of xylem and secondary cell wall thickness. Compared to wild-type, neither myb158 mutant nor myb158 myb189 double mutant showed obvious phenotypic variation in P. tomentosa. Although MYB158 and MYB189 proteins could repress the same structural genes involved in lignin, cellulose, and xylan biosynthesis, the two proteins had their own specific regulatory targets. Populus MYB158 could act as the upstream regulator of secondary cell wall NAC master switch and directly represses the expression of the SND1-B2 gene. Taken together, Populus MYB158 and MYB189 have retained similar functions in negatively regulating secondary cell wall biosynthesis, but have evolved partially distinct functions in direct regulation of NAC master switch, with MYB158 playing a more crucial role. Our findings provide new insights into the evolutionary and functional divergence of WGD-derived duplicate genes.     

Impact of gene gains,losses and duplication modes on the origin and diversification of vertebrates

The study of the evolutionary origin of vertebrates has been linked to the study of genome duplications since Susumo Ohno suggested that the successful diversification of vertebrate innovations was facilitated by two rounds of whole-genome duplication (2R-WGD) in the stem vertebrate. Since then, studies on the functional evolution of many genes duplicated in the vertebrate lineage have provided the grounds to support experimentally this link. This article reviews cases of gene duplications derived either from the 2R-WGD or from local gene duplication events in vertebrates, analyzing their impact on the evolution of developmental innovations. We analyze how gene regulatory networks can be rewired by the activity of transposable elements after genome duplications, discuss how different mechanisms of duplication might affect the fate of duplicated genes, and how the loss of gene duplicates might influence the fate of surviving paralogs. We also discuss the evolutionary relationships between gene duplication and alternative splicing, in particular in the vertebrate lineage. Finally, we discuss the role that the 2R-WGD might have played in the evolution of vertebrate developmental gene networks, paying special attention to those related to vertebrate key features such as neural crest cells, placodes, and the complex tripartite brain. In this context, we argue that current evidences points that the 2R-WGD may not be linked to the origin of vertebrate innovations, but to their subsequent diversification in a broad variety of complex structures and functions that facilitated the successful transition from peaceful filter-feeding non-vertebrate ancestors to voracious vertebrate predators.     

Whole-genome duplications spurred the functional diversification of the globin gene superfamily in vertebrates

It has been hypothesized that two successive rounds of whole-genome duplication (WGD) in the stem lineage of vertebrates provided genetic raw materials for the evolutionary innovation of many vertebrate-specific features. However, it has seldom been possible to trace such innovations to specific functional differences between paralogous gene products that derive from a WGD event. Here, we report genomic evidence for a direct link between WGD and key physiological innovations in the vertebrate oxygen transport system. Specifically, we demonstrate that key globin proteins that evolved specialized functions in different aspects of oxidative metabolism (hemoglobin, myoglobin, and cytoglobin) represent paralogous products of two WGD events in the vertebrate common ancestor. Analysis of conserved macrosynteny between the genomes of vertebrates and amphioxus (subphylum Cephalochordata) revealed that homologous chromosomal segments defined by myoglobin + globin-E, cytoglobin, and the α-globin gene cluster each descend from the same linkage group in the reconstructed proto-karyotype of the chordate common ancestor. The physiological division of labor between the oxygen transport function of hemoglobin and the oxygen storage function of myoglobin played a pivotal role in the evolution of aerobic energy metabolism, supporting the hypothesis that WGDs helped fuel key innovations in vertebrate evolution.     

Pattern of divergence of amino acid sequences encoded by paralogous genes in human and pufferfish

We used phylogenetic analyses of protein families containing two or more pairs of orthologues in the genomes of human and pufferfish (Takifugu rubripes) to test the hypothesis that these sequences show a strong signal of polyploidization events hypothesized to have occurred early in vertebrate history. In order to test for evidence of two distinct rounds of polyploidization (the 2R hypothesis), we compared the pattern of amino acid sequence divergence of proteins encoded by genes duplicated just prior to the most recent common ancestor of human and pufferfish with that of proteins encoded genes duplicated earlier. These sequence divergences were statistically indistinguishable, contrary to the prediction of the 2R hypothesis. The variance of amino acid sequence divergences between paralogues was significantly greater than expected from that of orthologues in the same families. Estimation of gene duplication times assuming a molecular clock provided earlier estimates than expected, suggesting that it may not be appropriate to time the duplication of paralogues using rate estimates derived from orthologous comparisons. Overall, the results indicate that amino acid sequences do not provide a strong signal supporting the hypothesis that gene duplications early in vertebrate history occurred by polyploidization. On the other hand, the data are easily explained under an alternative model that gene duplications occurred at different times in different vertebrate gene families.     

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.     

Fourfold paralogy regions on human HOX-bearing chromosomes: Role of ancient segmental duplications in the evolution of vertebrate genome

BackgroundSusumu Ohno’s idea that modern vertebrates are degenerate polyploids (concept referred as 2R hypothesis) has been the subject of intense debate for past four decades. It was proposed that intra-genomic synteny regions (paralogons) in human genome are remains of ancient polyploidization events that occurred early in the vertebrate history. The quadruplicated paralogon centered on human HOX clusters is taken as evidence that human HOX-bearing chromosomes were structured by two rounds of whole genome duplication (WGD) events.ResultsEvolutionary history of human HOX-bearing chromosomes (chromosomes 2/7/12/17) was evaluated by the phylogenetic analysis of multigene families with triplicated or quadruplicated distribution on these chromosomes. Topology comparison approach categorized the members of 44 families into four distinct co-duplicated groups. Distinct gene families belonging to a particular co-duplicated group, exhibit similar evolutionary history and hence have duplicated simultaneously, whereas genes of two distinct co-duplicated groups do not share their evolutionary history and have not duplicated in concert with each other.ConclusionThe recovery of co-duplicated groups suggests that “ancient segmental duplications and rearrangements” is the most rational model of evolutionary events that have generated the triplicated and quadruplicated paralogy regions seen on the human HOX-bearing chromosomes.     

Multiple Chromosomal Rearrangements Structured the Ancestral Vertebrate Hox-Bearing Protochromosomes

While the proposal that large-scale genome expansions occurred early in vertebrate evolution is widely accepted, the exact mechanisms of the expansion—such as a single or multiple rounds of whole genome duplication, bloc chromosome duplications, large-scale individual gene duplications, or some combination of these—is unclear. Gene families with a single invertebrate member but four vertebrate members, such as the Hox clusters, provided early support for Ohno's hypothesis that two rounds of genome duplication (the 2R-model) occurred in the stem lineage of extant vertebrates. However, despite extensive study, the duplication history of the Hox clusters has remained unclear, calling into question its usefulness in resolving the role of large-scale gene or genome duplications in early vertebrates. Here, we present a phylogenetic analysis of the vertebrate Hox clusters and several linked genes (the Hox “paralogon”) and show that different phylogenies are obtained for Dlx and Col genes than for Hox and ErbB genes. We show that these results are robust to errors in phylogenetic inference and suggest that these competing phylogenies can be resolved if two chromosomal crossover events occurred in the ancestral vertebrate. These results resolve conflicting data on the order of Hox gene duplications and the role of genome duplication in vertebrate evolution and suggest that a period of genome reorganization occurred after genome duplications in early vertebrates.     

Systematic phylogenomic evidence of en bloc duplication of the ancestral 8p11.21-8p21.3-like region

The genomes of many higher organisms, including plants and bony fish, frequently undergo polyploidization, and it has long been hypothesized that these, and other, large-scale genomic duplications have played an important role in the major evolutionary transitions of our past. Here we build upon an early work to show that the human genomic region 8p11.21-8p21.3 has three paralogous regions on chromosomes 4, 5, and 10 that were produced by two rounds of duplications after the protostomian-deuterostomian split and before the actinopterygian-sarcopterygian split. We base our analysis on the phylogenetic reconstruction of the evolutionary history of 38 gene families located in these regions. Using an alignment centered on protein domains, three different phylogenetic methods, and divergence time estimation, this analysis gives more support in favor of two ancient polyploidization events in the vertebrate ancestral genome.     

Exonic remnants of whole-genome duplication reveal cis-regulatory function of coding exons

Using a comparative genomics approach to reconstruct the fate of genomic regulatory blocks (GRBs) and identify exonic remnants that have survived the disappearance of their host genes after whole-genome duplication (WGD) in teleosts, we discover a set of 38 candidate cis-regulatory coding exons (RCEs) with predicted target genes. These elements demonstrate evolutionary separation of overlapping protein-coding and regulatory information after WGD in teleosts. We present evidence that the corresponding mammalian exons are still under both coding and non-coding selection pressure, are more conserved than other protein coding exons in the host gene and several control sets, and share key characteristics with highly conserved non-coding elements in the same regions. Their dual function is corroborated by existing experimental data. Additionally, we show examples of human exon remnants stemming from the vertebrate 2R WGD. Our findings suggest that long-range cis-regulatory inputs for developmental genes are not limited to non-coding regions, but can also overlap the coding sequence of unrelated genes. Thus, exonic regulatory elements in GRBs might be functionally equivalent to those in non-coding regions, calling for a re-evaluation of the sequence space in which to look for long-range regulatory elements and experimentally test their activity.     

The Timing of Timezyme Diversification in Vertebrates

All biological functions in vertebrates are synchronized with daily and seasonal changes in the environment by the time keeping hormone melatonin. Its nocturnal surge is primarily due to the rhythmic activity of the arylalkylamine N-acetyl transferase AANAT, which thus became the focus of many investigations regarding its evolution and function. Various vertebrate isoforms have been reported from cartilaginous fish to mammals but their origin has not been clearly established. Using phylogeny and synteny, we took advantage of the increasing number of available genomes in order to test whether the various rounds of vertebrate whole genome duplications were responsible for the diversification of AANAT. We highlight a gene secondary loss of the AANAT2 in the Sarcopterygii, revealing for the first time that the AAANAT1/2 duplication occurred before the divergence between Actinopterygii (bony fish) and Sarcopterygii (tetrapods, lobe-finned fish, and lungfish). We hypothesize the teleost-specific whole genome duplication (WDG) generated the appearance of the AANAT1a/1b and the AANAT2/2′paralogs, the 2′ isoform being rapidly lost in the teleost common ancestor (ray-finned fish). We also demonstrate the secondary loss of the AANAT1a in a Paracantopterygii (Atlantic cod) and of the 1b in some Ostariophysi (zebrafish and cave fish). Salmonids present an even more diverse set of AANATs that may be due to their specific WGD followed by secondary losses. We propose that vertebrate AANAT diversity resulted from 3 rounds of WGD followed by previously uncharacterized secondary losses. Extant isoforms show subfunctionalized localizations, enzyme activities and affinities that have increased with time since their emergence.     

Gene duplications and the early evolution of neural crest development

Neural crest cells are an important cell type present in all vertebrates, and elaboration of the neural crest is thought to have been a key factor in their evolutionary success. Genomic comparisons suggest there were two major genome duplications in early vertebrate evolution, raising the possibility that evolution of neural crest was facilitated by gene duplications. Here, we review the process of early neural crest formation and its underlying gene regulatory network (GRN) as well as the evolution of important neural crest derivatives. In this context, we assess the likelihood that gene and genome duplications capacitated neural crest evolution, particularly in light of novel data arising from invertebrate chordates.     

Phylogenetic history of paralogous gene quartets on human chromosomes 1, 2, 8 and 20 provides no evidence in favor of the vertebrate octoploidy hypothesis

Fourfold paralogy regions in the human genome have been considered historical remnants of whole-genome duplication events predicted to have occurred early in vertebrate evolution. Taking advantage of the well-annotated and high-quality human genomic sequence map as well as the ever-increasing accessibility of large-scale genomic sequence data from a diverse range of animal species, we investigated the prediction that the ancestral vertebrate genome was shaped by two rapid rounds of whole-genome duplication within a period of 10 million years. Both the map self-comparison approach and a phylogenetic analysis revealed that gene families identified as tetralogous on human chromosomes 1/2/8/20 arose by small-scale duplication events that occurred at widely different time points in animal evolution. Furthermore, the data discount the likelihood that tree topologies of the form ((A,B)(C,D)) are best explained by the octoploidy hypothesis. We instead propose that such symmetrical tree patterns are also consistent with local duplications and rearrangement events.     

ZBED Evolution: Repeated Utilization of DNA Transposons as Regulators of Diverse Host Functions

ZBED genes originate from domesticated hAT DNA transposons and encode regulatory proteins of diverse function in vertebrates. Here we reveal the evolutionary relationship between ZBED genes and demonstrate that they are derived from at least two independent domestication events in jawed vertebrate ancestors. We show that ZBEDs form two monophyletic clades, one of which has expanded through several independent duplications in host lineages. Subsequent diversification of ZBED genes has facilitated regulation of multiple diverse fundamental functions. In contrast to known examples of transposable element exaptation, our results demonstrate a novel unprecedented capacity for the repeated utilization of a family of transposable element-derived protein domains sequestered as regulators during the evolution of diverse host gene functions in vertebrates. Specifically, ZBEDs have contributed to vertebrate regulatory innovation through the donation of modular DNA and protein interacting domains. We identify that C7ORF29, ZBED2, 3, 4, and ZBEDX form a monophyletic group together with ZBED6, that is distinct from ZBED1 genes. Furthermore, we show that ZBED5 is related to Buster DNA transposons and is phylogenetically separate from other ZBEDs. Our results offer new insights into the evolution of regulatory pathways, and suggest that DNA transposons have contributed to regulatory complexity during genome evolution in vertebrates.     

Minor change, major difference: divergent functions of highly conserved cis-regulatory elements subsequent to whole genome duplication events

Within the vertebrate lineage, a high proportion of duplicate genes have been retained after whole genome duplication (WGD) events. It has been proposed that many of these duplicate genes became indispensable because the ancestral gene function was divided between them. In addition, novel functions may have evolved, owing to changes in cis-regulatory elements. Functional analysis of the PAX2/5/8 gene subfamily appears to support at least the first part of this hypothesis. The collective role of these genes has been widely retained, but sub-functions have been differentially partitioned between the genes in different vertebrates. Conserved non-coding elements (CNEs) represent an interesting and readily identifiable class of putative cis-regulatory elements that have been conserved from fish to mammals, an evolutionary distance of 450 million years. Within the PAX2/5/8 gene subfamily, PAX2 is associated with the highest number of CNEs. An additional WGD experienced in the teleost lineage led to two copies of pax2, each of which retained a large proportion of these CNEs. Using a reporter gene assay in zebrafish embryos, we have exploited this rich collection of regulatory elements in order to determine whether duplicate CNEs have evolved different functions. Remarkably, we find that even highly conserved sequences exhibit more functional differences than similarities. We also discover that short flanking sequences can have a profound impact on CNE function. Therefore, if CNEs are to be used as candidate enhancers for transgenic studies or for multi-species comparative analyses, it is paramount that the CNEs are accurately delineated.