Comparative genomics of duplicate γ-glutamyl transferase genes in teleosts: medaka (Oryzias latipes), stickleback (Gasterosteus aculeatus), green spotted pufferfish (Tetraodon nigroviridis), fugu (Takifugu rubripes), and zebrafish (Danio rerio)

The availability of multiple teleost (bony fish) genomes is providing unprecedented opportunities to understand the diversity and function of gene duplication events using comparative genomics. Here we examine multiple paralogous genes of γ-glutamyl transferase (GGT) in several distantly related teleost species including medaka, stickleback, green spotted pufferfish, fugu, and zebrafish. Through mining genome databases, we have identified multiple GGT orthologs. Duplicate (paralogous) GGT sequences for GGT1 (GGT1 a and b), GGTL1 (GGTL1 a and b), and GGTL3 (GGTL3 a and b) were identified for each species. Phylogenetic analysis suggests that GGTs are ancient proteins conserved across most metazoan phyla and those paralogous GGTs in teleosts likely arose from the serial 3R genome duplication events. A third GGTL1 gene (GGTL1c) was found in green spotted pufferfish; however, this gene is not present in medaka, stickleback, or fugu. Similarly, one or both paralogs of GGTL3 appear to have been lost in green spotted pufferfish, fugu, and zebrafish. Syntenic relationships were highly maintained between duplicated teleost chromosomes, among teleosts and across ray-finned (Actinopterygii) and lobe-finned (Sarcopterygii) species. To assess subfunction partitioning, six medaka GGT genes were cloned and assessed for developmental and tissue-specific expression. On the basis of these data, we propose a modification of the "duplication-degeneration-complementation" model of subfunction partitioning where quantitative differences rather than absolute differences in gene expression are observed between gene paralogs. Our results demonstrate that multiple GGT genes have been retained within teleost genomes. Questions remain, however, regarding the functional roles of multiple GGTs in these species.     

Phylogenetic Timing of the Fish-Specific Genome Duplication Correlates with the Diversification of Teleost Fish

For many genes, ray-finned fish (Actinopterygii) have two paralogous copies, where only one ortholog is present in tetrapods. The discovery of an additional, almost-complete set of Hox clusters in teleosts (zebrafish, pufferfish, medaka, and cichlid) but not in basal actinopterygian lineages (Polypterus) led to the formulation of the fish-specific genome duplication hypothesis. The phylogenetic timing of this genome duplication during the evolution of ray-finned fish is unknown, since only a few species of basal fish lineages have been investigated so far. In this study, three nuclear genes (fzd8, sox11, tyrosinase) were sequenced from sturgeons (Acipenseriformes), gars (Semionotiformes), bony tongues (Osteoglossomorpha), and a tenpounder (Elopomorpha). For these three genes, two copies have been described previously teleosts (e.g., zebrafish, pufferfish), but only one orthologous copy is found in tetrapods. Individual gene trees for these three genes and a concatenated dataset support the hypothesis that the fish-specific genome duplication event took place after the split of the Acipenseriformes and the Semionotiformes from the lineage leading to teleost fish but before the divergence of Osteoglossiformes. If these three genes were duplicated during the proposed fish-specific genome duplication event, then this event separates the species-poor early-branching lineages from the species-rich teleost lineage. The additional number of genes resulting from this event might have facilitated the evolutionary radiation and the phenotypic diversification of the teleost fish.[Reviewing Editor: Martin Kreitman]     

9种硬骨鱼类转座组注释

通过生物信息学手段对9种硬骨鱼转座组进行注释。结果表明9种硬骨鱼类转座组大小和构成差异显著,其转座组含量从高到低分别为斑马鱼、矛尾鱼、青鳉鱼、罗非鱼、花斑剑尾鱼、大西洋鳕鱼、三刺鱼、金娃娃和红鳍东方鲀,转座子含量和基因组大小呈正相关。DNA转座子在硬骨鱼类中具有多样性高和含量差异大的特点(0.50%–38.37%),是硬骨鱼类转座组差异的主要决定因素,其中h AT和Tc/Mariner超家族是硬骨鱼类主要的DNA转座子。RNA转座子在硬骨鱼类中也具有多样性高的特点,其中LINE转座子占硬骨鱼类基因组的0.53%–5.75%,共检测到14个超家族分布,其中L1、L2、RTE和Rex转座子扩增较为明显,LTR转座子除了在斑马鱼和三刺鱼中含量达到5.58%和2.51%,在大多硬骨鱼类基因组中的含量低于2%,在硬骨鱼类中共检测到6个LTR转座子(Copia、DIRS、ERV、Gypsy、Ngaro和Pao)超家族分布,其中扩增最为明显的是Gypsy。而SINE转座子在硬骨鱼类中扩增最弱,仅在斑马鱼和矛尾鱼中分别达到3.28%和5.64%,在其他7个物种中低于1%。SINE中t RNA、5S和MIR三个超家族在部分硬骨鱼类中有一定程度扩增。本研究表明硬骨鱼类转座组具有多样性丰富、差异大的特点,转座组差异与硬骨鱼基因组大小有很强的相关性,转座组是决定硬骨鱼基因组大小的重要因素。     

Ancient vertebrate conserved noncoding elements have been evolving rapidly in teleost fishes

Vertebrate genomes contain thousands of conserved noncoding elements (CNEs) that often function as tissue-specific enhancers. In this study, we have identified CNEs in human, dog, chicken, Xenopus, and four teleost fishes (zebrafish, stickleback, medaka, and fugu) using elephant shark, a cartilaginous vertebrate, as the base genome and investigated the evolution of these ancient vertebrate CNEs (aCNEs) in bony vertebrate lineages. Our analysis shows that aCNEs have been evolving at different rates in different bony vertebrate lineages. Although 78-83% of CNEs have diverged beyond recognition ("lost") in different teleost fishes, only 24% and 40% have been lost in the chicken and mammalian lineages, respectively. Relative rate tests of substitution rates in CNEs revealed that the teleost fish CNEs have been evolving at a significantly higher rate than those in other bony vertebrates. In the ray-finned fish lineage, 68% of aCNEs were lost before the divergence of the four teleosts. This implicates the "fish-specific" whole-genome duplication in the accelerated evolution and the loss of a large number of both copies of duplicated CNEs in teleost fishes. The aCNEs are rich in tissue-specific enhancers and thus many of them are likely to be evolutionarily constrained cis-regulatory elements. The rapid evolution of aCNEs might have affected the expression patterns driven by them. Transgenic zebrafish assay of some human CNE enhancers that have been lost in teleosts has indicated instances of conservation or changes in trans-acting factors between mammals and fishes.     

Temporal pattern of loss/persistence of duplicate genes involved in signal transduction and metabolic pathways after teleost-specific genome duplication

Background  

Recent genomic studies have revealed a teleost-specific third-round whole genome duplication (3R-WGD) event occurred in a common ancestor of teleost fishes. However, it is unclear how the genes duplicated in this event were lost or persisted during the diversification of teleosts, and therefore, how many of the duplicated genes contribute to the genetic differences among teleosts. This subject is also important for understanding the process of vertebrate evolution through WGD events. We applied a comparative evolutionary approach to this question by focusing on the genes involved in long-term potentiation, taste and olfactory transduction, and the tricarboxylic acid cycle, based on the whole genome sequences of four teleosts; zebrafish, medaka, stickleback, and green spotted puffer fish.     

Length and GC Content Variability of Introns among Teleostean Genomes in the Light of the Metabolic Rate Hypothesis

A comparative analysis of five teleostean genomes, namely zebrafish, medaka, three-spine stickleback, fugu and pufferfish was performed with the aim to highlight the nature of the forces driving both length and base composition of introns (i.e., bpi and GCi). An inter-genome approach using orthologous intronic sequences was carried out, analyzing independently both variables in pairwise comparisons. An average length shortening of introns was observed at increasing average GCi values. The result was not affected by masking transposable and repetitive elements harbored in the intronic sequences. The routine metabolic rate (mass specific temperature-corrected using the Boltzmann''s factor) was measured for each species. A significant correlation held between average differences of metabolic rate, length and GC content, while environmental temperature of fish habitat was not correlated with bpi and GCi. Analyzing the concomitant effect of both variables, i.e., bpi and GCi, at increasing genomic GC content, a decrease of bpi and an increase of GCi was observed for the significant majority of the intronic sequences (from ∼40% to ∼90%, in each pairwise comparison). The opposite event, concomitant increase of bpi and decrease of GCi, was counter selected (from <1% to ∼10%, in each pairwise comparison). The results further support the hypothesis that the metabolic rate plays a key role in shaping genome architecture and evolution of vertebrate genomes.     

Genomic organization and characterization of two vomeronasal 1 receptor-like genes (ora1 and ora2) in Atlantic salmon Salmo salar

Olfactory receptors are encoded by three large multigene superfamilies (OR, V1R and V2R) in mammals. Fish do not possess a vomeronasal system; therefore, it has been proposed that their V1R-like genes be classified as olfactory receptors related to class A G protein-coupled receptors (ora). Unlike mammalian genomes, which contain more than a hundred V1R genes, the five species of teleost fish that have been investigated to date appear to have six ora genes (ora1-6) except for pufferfish that have lost ora1. The common ancestor of salmonid fishes is purported to have undergone a whole genome duplication. As salmonids have a life history that requires the use of olfactory cues to navigate back to their natal habitats to spawn, we set out to determine if ora1 or ora2 is duplicated in a representative species, Atlantic salmon (Salmo salar). We used an oligonucleotide probe designed from a conserved sequence of several teleost ora2 genes to screen an Atlantic salmon BAC library (CHORI-214). Hybridization-positive BACs belonged to a single fingerprint contig of the Atlantic salmon physical map. All were also positive for ora2 by PCR. One of these BACs was chosen for further study, and shotgun sequencing of this BAC identified two V1R-like genes, ora1 and ora2, that are in a head-to-head conformation as is seen in some other teleosts. The gene products, ora1 and ora2, are highly conserved among teleosts. We only found evidence for a single ora1-2 locus in the Atlantic salmon genome, which was mapped to linkage group 6. Fluorescent in situ hybridization (FISH) analysis placed ora1-2 on chromosome 12. Conserved synteny was found surrounding the ora1 and ora2 genes in Atlantic salmon, medaka and three-spined stickleback, but not zebrafish.     

Linked by Ancestral Bonds: Multiple Whole-Genome Duplications and Reticulate Evolution in a Brassicaceae Tribe

Pervasive hybridization and whole-genome duplications (WGDs) influenced genome evolution in several eukaryotic lineages. Although frequent and recurrent hybridizations may result in reticulate phylogenies, the evolutionary events underlying these reticulations, including detailed structure of the ancestral diploid and polyploid genomes, were only rarely reconstructed. Here, we elucidate the complex genomic history of a monophyletic clade from the mustard family (Brassicaceae), showing contentious relationships to the early-diverging clades of this model plant family. Genome evolution in the crucifer tribe Biscutelleae (∼60 species, 5 genera) was dominated by pervasive hybridizations and subsequent genome duplications. Diversification of an ancestral diploid genome into several divergent but crossable genomes was followed by hybridizations between these genomes. Whereas a single genus (Megadenia) remained diploid, the four remaining genera originated by allopolyploidy (Biscutella, Lunaria, Ricotia) or autopolyploidy (Heldreichia). The contentious relationships among the Biscutelleae genera, and between the tribe and other early diverged crucifer lineages, are best explained by close genomic relatedness among the recurrently hybridizing ancestral genomes. By using complementary cytogenomics and phylogenomics approaches, we demonstrate that the origin of a monophyletic plant clade can be more complex than a parsimonious assumption of a single WGD spurring postpolyploid cladogenesis. Instead, recurrent hybridization among the same and/or closely related parental genomes may phylogenetically interlink diploid and polyploid genomes despite the incidence of multiple independent WGDs. Our results provide new insights into evolution of early-diverging Brassicaceae lineages and elucidate challenges in resolving the contentious relationships within and between land plant lineages with pervasive hybridization and WGDs.     

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.     

Comprehensive analysis of medaka major histocompatibility complex (MHC) class II genes: Implications for evolution in teleosts

The major histocompatibility complex (MHC) class II molecules play central roles in adaptive immunity by regulating immune response via the activation of CD4 T cells. The full complement of the MHC class II genes has been elucidated only in mammalian species to date. To understand the evolution of these genes, we performed their first comprehensive analysis in nonmammalian species using a teleost, medaka (Oryzias latipes). Based on a database search, cDNA cloning, and genomic PCR, medaka was shown to possess five pairs of expressed class II genes, comprising one IIA and one IIB gene. Each pair was located on a different chromosome and was not linked to the class I genes. Only one pair showed a high degree of polymorphism and was considered to be classical class II genes, whereas the other four pairs were nonclassical. Phylogenetic analysis of all medaka class II genes and most reported teleost class II genes revealed that the IIA and IIB genes formed separate clades, each containing three well-corresponding lineages. One lineage contained three medaka genes and all known classical class II genes of Ostariophysi and Euteleostei and was presumed to be an original lineage of the teleost MHC class II genes. The other two lineages contained one nonclassical medaka gene each and some Euteleostei genes. These results indicate that multiple lineages of the teleost MHC class II genes have been conserved for hundreds of millions of years and that the tightly linked IIA and IIB genes have undergone concerted evolution.     

Search for enhancers: teleost models in comparative genomic and transgenic analysis of cis regulatory elements

Homology searches between DNA sequences of evolutionary distant species (phylogenetic footprinting) offer a fast detection method for regulatory sequences. Because of the small size of their genomes, tetraodontid species such as the Japanese pufferfish and green spotted pufferfish have become attractive models for comparative genomics. A disadvantage of the tetraodontid species is, however, that they cannot be bred and manipulated routinely under laboratory conditions, so these species are less attractive for developmental and genetic analysis. In contrast, an increasing arsenal of transgene techniques with the developmental model species zebrafish and medaka are being used for functional analysis of cis regulatory sequences. The main disadvantage is the much larger genome. While comparison between many loci proved the suitability of phylogenetic footprinting using fish and mammalian sequences, fast rate of change in enhancer structure and gene duplication within teleosts may obscure detection of homologies. Here we discuss the contribution and potentials provided by different teleost models for the detection and functional analysis of conserved cis-regulatory elements.     

Genomic organization of repetitive DNAs in the cichlid fish <Emphasis Type="Italic">Astronotus ocellatus</Emphasis>

To contribute to the knowledge of fish genomes, we identified and characterized by means of nucleotide sequencing and physical chromosome mapping, three classes of repetitive DNAs in the genome of the South American cichlid fish Astronotus ocellatus. The first class corresponds to a satellite DNA family (AoSat) that shares similarity with a centromeric satellite DNA of the pufferfish Tetraodon nigroviridis. The second repetitive DNA class (AoRex3) is related to the retrotransposon Rex3, which is widely distributed among teleost fishes. The last repetitive element (AoLINE) shows a high similarity to the CR1-like LINE element of other teleosts. The three isolated repetitive elements are clustered in the centromeric heterochromatin of all chromosomes of the complement. The repetitive sequences are not randomly distributed in the genome, suggesting a pattern of compartmentalization on chromosomes.     

Comparative analysis of teleost fish genomes reveals preservation of different ancient clock duplicates in different fishes

Clock (Circadian locomotor output cycle kaput) was the first vertebrate circadian clock gene identified in a mouse forward genetics mutagenesis screen. It encodes a bHLH-PAS protein that is highly conserved throughout evolution. Tetrapods also have the second Clock gene, Clock2 or Npas2 (Neuronal PAS domain protein 2). Conversely, the fruit fly, an invertebrate, has only one clock gene. Interrogation of the five teleost fish genome databases revealed that the zebrafish and the Japanese pufferfish (fugu) each have three clock genes, whereas the green spotted pufferfish (tetraodon), the Japanese medaka fish and the three-spine stickleback each have two clock genes. Phylogenetic and splice site analyses indicated that zebrafish and fugu each have two clock1 genes, clock1a and clock1b and one clock2; tetraodon also have clock1a and clock1b but do not have clock2; and medaka and stickleback each have clock1b and one clock2. Genome neighborhood analysis further showed that clock1a/clock1b in zebrafish, fugu and tetraodon is an ancient duplicate. While the dN/dS ratios of these three fish clock duplicates are all <1, indicating that purifying selection has acted upon them; the Tajima relative rate test showed that all three fish clock duplicates have asymmetric evolutionary rates, implicating that one of these duplicates have been under positive selection or relaxed functional constraint. These results support the view that teleost fish clock genes were generated from an ancient genome-wide duplication, and differential gene loss after the duplication resulted in retention of different ancient duplicates in different teleost fishes, which could have contributed to the evolution of the distinct fish circadian clock mechanisms.     

Comparative genomic analysis of teleost fish <Emphasis Type="Italic">bmal</Emphasis> genes

Bmal1 (Brain and muscle ARNT like 1) gene is a key circadian clock gene. Tetrapods also have the second Bmal gene, Bmal2. Fruit fly has only one bmal1/cycle gene. Interrogation of the five teleost fish genome sequences coupled with phylogenetic and splice site analyses found that zebrafish have two bmal1 genes, bmal1a and bmal1b, and bmal2a; Japanese pufferfish (fugu), green spotted pufferfish (tetraodon) and Japanese medaka fish each have two bmal2 genes, bmal2a and bmal2b, and bmal1a; and three-spine stickleback have bmal1a and bmal2b. Syntenic analysis further indicated that zebrafish bmal1a/bmal1b, and fugu, tetraodon and medaka bmal2a/bmal2b are ancient duplicates. Although the dN/dS ratios of these four fish bmal duplicates are all <1, implicating they have been under purifying selection, the Tajima relative rate test showed that fugu, tetraodon and medaka bmal2a/bmal2b have asymmetric evolutionary rates, suggesting that one of these duplicates have been subject to positive selection or relaxed functional constraint. These results support the notion that teleost fish bmal genes were derived from the fish-specific genome duplication (FSGD), divergent resolution following the duplication led to retaining different ancient bmal duplicates in different fishes, which could have shaped the evolution of the complex teleost fish timekeeping mechanisms. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Survey sequencing and comparative analysis of the elephant shark (Callorhinchus milii) genome

Owing to their phylogenetic position, cartilaginous fishes (sharks, rays, skates, and chimaeras) provide a critical reference for our understanding of vertebrate genome evolution. The relatively small genome of the elephant shark, Callorhinchus milii, a chimaera, makes it an attractive model cartilaginous fish genome for whole-genome sequencing and comparative analysis. Here, the authors describe survey sequencing (1.4× coverage) and comparative analysis of the elephant shark genome, one of the first cartilaginous fish genomes to be sequenced to this depth. Repetitive sequences, represented mainly by a novel family of short interspersed element–like and long interspersed element–like sequences, account for about 28% of the elephant shark genome. Fragments of approximately 15,000 elephant shark genes reveal specific examples of genes that have been lost differentially during the evolution of tetrapod and teleost fish lineages. Interestingly, the degree of conserved synteny and conserved sequences between the human and elephant shark genomes are higher than that between human and teleost fish genomes. Elephant shark contains putative four Hox clusters indicating that, unlike teleost fish genomes, the elephant shark genome has not experienced an additional whole-genome duplication. These findings underscore the importance of the elephant shark as a critical reference vertebrate genome for comparative analysis of the human and other vertebrate genomes. This study also demonstrates that a survey-sequencing approach can be applied productively for comparative analysis of distantly related vertebrate genomes.     

Evolution of genome size and chromosome number in the carnivorous plant genus Genlisea (Lentibulariaceae), with a new estimate of the minimum genome size in angiosperms

MethodsNuclear genome sizes were measured from cultivated plant material for a comprehensive sampling of taxa, including nearly half of all species of Genlisea and representing all major lineages. Flow cytometric measurements were conducted in parallel in two laboratories in order to compare the consistency of different methods and controls. Chromosome counts were performed for the majority of taxa, comparing different staining techniques for the ultrasmall chromosomes.ConclusionsGenlisea is an ideal candidate model organism for the understanding of genome reduction as the genus includes species with both relatively large (∼1700 Mbp) and ultrasmall (∼61 Mbp) genomes. This comparative, phylogeny-based analysis of genome sizes and karyotypes in Genlisea provides essential data for selection of suitable species for comparative whole-genome analyses, as well as for further studies on both the molecular and cytogenetic basis of genome reduction in plants.     

Examining Phylogenetic Relationships Among Gibbon Genera Using Whole Genome Sequence Data Using an Approximate Bayesian Computation Approach

Gibbons are believed to have diverged from the larger great apes ∼16.8 MYA and today reside in the rainforests of Southeast Asia. Based on their diploid chromosome number, the family Hylobatidae is divided into four genera, Nomascus, Symphalangus, Hoolock, and Hylobates. Genetic studies attempting to elucidate the phylogenetic relationships among gibbons using karyotypes, mitochondrial DNA (mtDNA), the Y chromosome, and short autosomal sequences have been inconclusive . To examine the relationships among gibbon genera in more depth, we performed second-generation whole genome sequencing (WGS) to a mean of ∼15× coverage in two individuals from each genus. We developed a coalescent-based approximate Bayesian computation (ABC) method incorporating a model of sequencing error generated by high coverage exome validation to infer the branching order, divergence times, and effective population sizes of gibbon taxa. Although Hoolock and Symphalangus are likely sister taxa, we could not confidently resolve a single bifurcating tree despite the large amount of data analyzed. Instead, our results support the hypothesis that all four gibbon genera diverged at approximately the same time. Assuming an autosomal mutation rate of 1 × 10⁻⁹/site/year this speciation process occurred ∼5 MYA during a period in the Early Pliocene characterized by climatic shifts and fragmentation of the Sunda shelf forests. Whole genome sequencing of additional individuals will be vital for inferring the extent of gene flow among species after the separation of the gibbon genera.     

Evolutionary Patterns and Selective Pressures of Odorant/Pheromone Receptor Gene Families in Teleost Fishes

Background

Teleost fishes do not have a vomeronasal organ (VNO), and their vomeronasal receptors (V1Rs, V2Rs) are expressed in the main olfactory epithelium (MOE), as are odorant receptors (ORs) and trace amine-associated receptors (TAARs). In this study, to obtain insights into the functional distinction among the four chemosensory receptor families in teleost fishes, their evolutionary patterns were examined in zebrafish, medaka, stickleback, fugu, and spotted green pufferfish.

Methodology/Principal Findings

Phylogenetic analysis revealed that many lineage-specific gene gains and losses occurred in the teleost fish TAARs, whereas only a few gene gains and losses have taken place in the teleost fish vomeronasal receptors. In addition, synonymous and nonsynonymous nucleotide substitution rate ratios (K_A/K_S) in TAARs tended to be higher than those in ORs and V2Rs.

Conclusions/Significance

Frequent gene gains/losses and high K_A/K_S in teleost TAARs suggest that receptors in this family are used for detecting some species-specific chemicals such as pheromones. Conversely, conserved repertoires of V1R and V2R families in teleost fishes may imply that receptors in these families perceive common odorants for teleosts, such as amino acids. Teleost ORs showed intermediate evolutionary pattern between TAARs and vomeronasal receptors. Many teleost ORs seem to be used for common odorants, but some ORs may have evolved to recognize lineage-specific odors.