Modular Evolution of PGC-1α in Vertebrates

In mammals, the peroxisome proliferator activated receptor (PPAR)γ coactivator-1α (PGC-1α) is a central regulator of mitochondrial gene expression, acting in concert with nuclear respiratory factor-1 (NRF-1) and the PPARs. Its role as a “master regulator” of oxidative capacity is clear in mammals, but its role in other vertebrates is ambiguous. In lower vertebrates, although PGC-1α seems to play a role in coordinating the PPARα axis as in mammals, it does not appear to be involved in NRF-1 regulation of mitochondrial content. To evaluate the evolutionary patterns of this coactivator in fish and mammals, we investigated the evolutionary trajectories of PGC-1α homologs in representative vertebrate lineages. A phylogeny of the PGC-1 paralogs suggested that the family diversified through repeated genome duplication events early in vertebrate evolution. Bayesian and maximum likelihood phylogenetic reconstructions of PGC-1α in representative vertebrate species revealed divergent evolutionary dynamics across the different functional domains of the protein. Specifically, PGC-1α exhibited strong conservation of the activation/PPAR interaction domain across vertebrates, whereas the NRF-1 and MEF2c interaction domains experienced accelerated rates of evolution in actinopterygian (fish lineages) compared to sarcopterygians (tetrapod lineages). Furthermore, analysis of the amino acid sequence of these variable domains revealed successive serine- and glutamine-rich insertions within the teleost lineages, with important ramifications for PGC-1α function in these lineages. Collectively, these results suggest modular evolution of the PGC-1α protein in vertebrates that could allow for lineage-specific divergences in the coactivating capabilities of this regulator.     

Evolutionary crossroads in developmental biology: cyclostomes (lamprey and hagfish)

Lampreys and hagfish, which together are known as the cyclostomes or 'agnathans', are the only surviving lineages of jawless fish. They diverged early in vertebrate evolution, before the origin of the hinged jaws that are characteristic of gnathostome (jawed) vertebrates and before the evolution of paired appendages. However, they do share numerous characteristics with jawed vertebrates. Studies of cyclostome development can thus help us to understand when, and how, key aspects of the vertebrate body evolved. Here, we summarise the development of cyclostomes, highlighting the key species studied and experimental methods available. We then discuss how studies of cyclostomes have provided important insight into the evolution of fins, jaws, skeleton and neural crest.     

The complete mitochondrial genome of the nudibranch Roboastra europaea (Mollusca: Gastropoda) supports the monophyly of opisthobranchs

The complete nucleotide sequence (14,472 bp) of the mitochondrial genome of the nudibranch Roboastra europaea (Gastropoda: Opisthobranchia) was determined. This highly compact mitochondrial genome is nearly identical in gene organization to that found in opisthobranchs and pulmonates (Euthyneura) but not to that in prosobranchs (a paraphyletic group including the most basal lineages of gastropods). The newly determined mitochondrial genome differs only in the relative position of the trnC gene when compared with the mitochondrial genome of Pupa strigosa, the only opisthobranch mitochondrial genome sequenced so far. Pupa and Roboastra represent the most basal and derived lineages of opisthobranchs, respectively, and their mitochondrial genomes are more similar in sequence when compared with those of pulmonates. All phylogenetic analyses (maximum parsimony, minimum evolution, maximum likelihood, and Bayesian) based on the deduced amino acid sequences of all mitochondrial protein-coding genes supported the monophyly of opisthobranchs. These results are in agreement with the classical view that recognizes Opisthobranchia as a natural group and contradict recent phylogenetic studies of the group based on shorter sequence data sets. The monophyly of opisthobranchs was further confirmed when a fragment of 2,500 nucleotides including the mitochondrial cox1, rrnL, nad6, and nad5 genes was analyzed in several species representing five different orders of opisthobranchs with all common methods of phylogenetic inference. Within opisthobranchs, the polyphyly of cephalaspideans and the monophyly of nudibranchs were recovered. The evolution of mitochondrial tRNA rearrangements was analyzed using the cox1+rrnL+nad6+nad5 gene phylogeny. The relative position of the trnP gene between the trnA and nad6 genes was found to be a synapomorphy of opisthobranchs that supports their monophyly.     

Replication-Associated Mutational Pressure (RMP) Governs Strand-Biased Compositional Asymmetry (SCA) and Gene Organization in Animal Mitochondrial Genomes

The nucleotide composition of the light (L-) and heavy (H-) strands of animal mitochondrial genomes is known to exhibit strand-biased compositional asymmetry (SCA). One of the possibilities is the existence of a replication-associated mutational pressure (RMP) that may introduce characteristic nucleotide changes among mitochondrial genomes of different animal lineages. Here, we discuss the influence of RMP on nucleotide and amino acid compositions as well as gene organization. Among animal mitochondrial genomes, RMP may represent the major force that compels the evolution of mitochondrial protein-coding genes, coupled with other process-based selective pressures, such as on components of translation machinery- tRNAs and their anticodons. Through comparative analyses of sequenced mitochondrial genomes among diverse animal lineages and literature reviews, we suggest a strong RMP effect, observed among invertebrate mitochondrial genes as compared to those of vertebrates, that is either a result of positive selection on the invertebrate or a relaxed selective pressure on the vertebrate mitochondrial genes.     

Characterization of novel GPCR gene coding locus in amphioxus genome: gene structure, expression, and phylogenetic analysis with implications for its involvement in chemoreception

Chemosensation is the primary sensory modality in almost all metazoans. The vertebrate olfactory receptor genes exist as tandem clusters in the genome, so that identifying their evolutionary origin would be useful for understanding the expansion of the sensory world in relation to a large-scale genomic duplication event in a lineage leading to the vertebrates. In this study, I characterized a novel GPCR (G-protein-coupled receptor) gene-coding locus from the amphioxus genome. The genomic DNA contains an intronless ORF whose deduced amino acid sequence encodes a seven-transmembrane protein with some amino acid residues characteristic of vertebrate olfactory receptors (ORs). Surveying counterparts in the Ciona intestinalis (Asidiacea, Urochordata) genome by querying BLAST programs against the Ciona genomic DNA sequence database resulted in the identification of a remotely related gene. In situ hybridization analysis labeled primary sensory neurons in the rostral epithelium of amphioxus adults. Based on these findings, together with comparison of the developmental gene expression between amphioxus and vertebrates, I postulate that chemoreceptive primary sensory neurons in the rostrum are an ancient cell population traceable at least as far back in phylogeny as the common ancestor of amphioxus and vertebrates.     

Evolution of constrained gonadotropin-releasing hormone ligand conformation and receptor selectivity

Gonadotropin-releasing hormone (GnRH) is the central regulator of reproduction in vertebrates. GnRHs have recently been identified in protochordates and retain the conserved N- and C-terminal domains involved in receptor binding and activation. GnRHs of the jawed vertebrates have a central achiral amino acid (glycine) that favors a type II' beta-turn such that the N- and C-terminal domains are closely apposed in binding the GnRH receptor. However, protochordate GnRHs have a chiral amino acid in this position, suggesting that they bind their receptors in a more extended form. We demonstrate here that a protochordate GnRH receptor does not distinguish GnRHs with achiral or chiral amino acids, whereas GnRH receptors of jawed vertebrates are highly selective for GnRHs with the central achiral glycine. The poor activity of the protochordate GnRH was increased >10-fold at vertebrate receptors by replacement of the chiral amino acid with glycine or a d-amino acid, which favor the type II' beta-turn. Structural analysis of the GnRHs using ion mobility-mass spectrometry and molecular modeling showed a greater propensity for a type II' beta-turn in GnRHs with glycine or a d-amino acid, which correlates with binding affinity at vertebrate receptors. These findings indicate that the substitution of glycine for a chiral amino acid in GnRH during evolution allows a more constrained conformation for receptor binding and that this subtle single amino acid substitution in a site remote from the ligand functional domains has marked effects on its structure and activity.     

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.     

Multiple origins and rapid evolution of duplicated mitochondrial genes in parthenogenetic geckos (Heteronotia binoei; Squamata, Gekkonidae)

Accumulating evidence for alternative gene orders demonstrates that vertebrate mitochondrial genomes are more evolutionarily dynamic than previously thought. Several lineages of parthenogenetic lizards contain large, tandem duplications that include rRNA, tRNA, and protein-coding genes, as well as the control region. Such duplications are hypothesized as intermediate stages in gene rearrangement, but the early stages of their evolution have not been previously studied. To better understand the evolutionary dynamics of duplicated segments of mitochondrial DNA, we sequenced 10 mitochondrial genomes from recently formed ( approximately 300,000 years ago) hybrid parthenogenetic geckos of the Heteronotia binoei complex and 1 from a sexual form. These genomes included some with an arrangement typical of vertebrates and others with tandem duplications varying in size from 5.7 to 9.4 kb, each with different gene contents and duplication endpoints. These results, together with phylogenetic analyses, indicate independent and frequent origins of the duplications. Small, direct repeats at the duplication endpoints imply slipped-strand error as a mechanism generating the duplications as opposed to a false initiation/termination of DNA replication mechanism that has been invoked to explain duplications in other lizard mitochondrial systems. Despite their recent origin, there is evidence for nonfunctionalization of genes due primarily to deletions, and the observed pattern of gene disruption supports the duplication-deletion model for rearrangement of mtDNA gene order. Conversely, the accumulation of mutations between these recent duplicates provides no evidence for gene conversion, as has been reported in some other systems. These results demonstrate that, despite their long-term stasis in gene content and arrangement in some lineages, vertebrate mitochondrial genomes can be evolutionary dynamic even at short timescales.     

Evolution of the WANCY region in amniote mitochondrial DNA

In most vertebrate mitochondrial genomes, the site for initiation of light-strand replication, OL, is found within a cluster of five transfer RNA (tRNA) genes (tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), and tRNA(Tyr)). This region and part of the adjacent cytochrome c oxydase subunit I (COI) gene were sequenced for two crocodilian, two turtle, and one snake species and for Sphenodon punctatus; part of the adjacent nicotinamide adenine dinucleotide dehydrogenase subunit 2 (ND2) gene was also sequenced for the crocodilian and turtle species. All had the typical vertebrate gene order. The turtles and the snake have a lengthy noncoding sequence between the tRNA(Asn) and tRNA(Cys) genes that we assumed to be homologous to the mammalian OL. The crocodilians and Sphenodon lack such a sequence, a condition they share with birds. Most proposed phylogenies for the amniotes require that OL at this position was lost at least twice during their diversification or was evolved independently more than once. Within the five tRNA genes, frequencies of substitutions are much higher in loops than in stems. Many loops vary dramatically in size among the species; in the most extreme case, the D-arm of the Sphenodon tRNA(Cys) is a "D-arm replacement" loop of seven nucleotides. Frequency of transitions in stems is relatively uniform across tRNAs, but frequency of transversions varies greatly. Mismatches in stems are infrequent, and their relative frequency in a specific tRNA is unrelated to the frequency of substitution in the corresponding gene. Several features of mammalian mitochondrial tRNAs are conserved in WANCY tRNAs throughout amniotes. The inferred initiation codon for COI is GTG in crocodilians, turtles, and the snake, a condition they share with fishes, certain amphibians, and birds. TTG appears to be the initiation codon for COI in Sphenodon; if correct, this would be a novel initiation codon for vertebrate mitochondrial DNA. Phylogenetic analyses of the inferred amino acid sequences of ND2 and COI support the sister-group relationship of birds and crocodilians and suggest that mammals are an early derived lineage within the amniotes.      

Sequence and gene organization of the chicken mitochondrial genome. A novel gene order in higher vertebrates

The 16,775 base-pair mitochondrial genome of the white Leghorn chicken has been cloned and sequenced. The avian genome encodes the same set of genes (13 proteins, 2 rRNAs and 22 tRNAs) as do other vertebrate mitochondrial DNAs and is organized in a very similar economical fashion. There are very few intergenic nucleotides and several instances of overlaps between protein or tRNA genes. The protein genes are highly similar to their mammalian and amphibian counterparts and are translated according to the same variant genetic code. Despite these highly conserved features, the chicken mitochondrial genome displays two distinctive characteristics. First, it exhibits a novel gene order, the contiguous tRNA(Glu) and ND6 genes are located immediately adjacent to the displacement loop region of the molecule, just ahead of the contiguous tRNA(Pro), tRNA(Thr) and cytochrome b genes, which border the displacement loop region in other vertebrate mitochondrial genomes. This unusual gene order is conserved among the galliform birds. Second, a light-strand replication origin, equivalent to the conserved sequence found between the tRNA(Cys) and tRNA(Asn) genes in all vertebrate mitochondrial genomes sequenced thus far, is absent in the chicken genome. These observations indicate that galliform mitochondrial genomes departed from their mammalian and amphibian counterparts during the course of evolution of vertebrate species. These unexpected characteristics represent useful markers for investigating phylogenetic relationships at a higher taxonomic level.     

A comparative and phylogenetic analysis of the alpha-actinin rod domain

Alpha-actinin is a ubiquitous actin-binding protein, composed of 3 domains; an actin-binding domain and a calcium-binding domain at the termini, connected by a rod domain composed by 1, 2, or 4 spectrin repeats (SRs). To understand how the rod domain has evolved during evolution, we have analyzed and compared the amino acid residue heterogeneity and phylogeny of the SRs of alpha-actinins of vertebrates, invertebrates, fungi, and several protozoa. The repeats of vertebrate alpha-actinins show a high degree of similarity, whereas repeats of invertebrates, fungi, and, in particular, of protozoa are more divergent. In the phylogeny, SR1 of all species were clustered together, independent of the number of repeats in the protein. It was also obvious that the second and last repeat in fungi (SR2) grouped with the fourth and last repeat of vertebrates and invertebrates (SR4). Therefore, the phylogeny implied that the rod domain of the cenancestral alpha-actinin only contained one SR. It was also obvious that SR2 of fungi are related to SR4 of vertebrates and invertebrates, implying that in the second intragenic duplication 2 repeats (i.e., what become SR2 and SR3) were inserted between the initial 2 repeats that become SR1 and SR4.     

Molecular evolution of the vertebrate hexokinase gene family: Identification of a conserved fifth vertebrate hexokinase gene

Hexokinases (HK) phosphorylate sugar immediately upon its entry into cells allowing these sugars to be metabolized. A total of four hexokinases have been characterized in a diversity of vertebrates—HKI, HKII, HKIII, and HKIV. HKIV is often called glucokinase (GCK) and has half the molecular weight of the other hexokinases, as it only has one hexokinase domain, while other vertebrate HKs have two. Differing hypothesis has been proposed to explain the diversification of the hexokinase gene family. We used a genomic approach to characterize hexokinase genes in a diverse array of vertebrate species and close relatives. Surprisingly we identified a fifth hexokinase-like gene, HKDC1 that exists and is expressed in diverse vertebrates. Analysis of the amino acid sequence of HKDC1 suggests that it may function as a hexokinase. To understand the evolution of the vertebrate hexokinase gene family we established a phylogeny of the hexokinase domain in all of the vertebrate hexokinase genes, as well as hexokinase genes from close relatives of the vertebrates. Our phylogeny demonstrates that duplication of the hexokinase domain, yielding a HK with two hexokinase domains, occurred prior to the diversification of the hexokinase gene family. We also establish that GCK evolved from a two hexokinase domain-containing gene, but has lost its N-terminal hexokinase domain. We also show that parallel changes in enzymatic function of HKI and HKIII have occurred.     

Vertebrate serpins: construction of a conflict-free phylogeny by combining exon-intron and diagnostic site analyses

A combination of three independent biological features, genomic organization, diagnostic amino acid sites, and rare indels, was used to elucidate the phylogeny of the vertebrate serpin (serine protease inhibitor) superfamily. A strong correlation between serpin gene families displaying (1) a conserved exon-intron pattern and (2) family-specific combinations of amino acid residues at specific sites suggests that present-day vertebrates encompass six serpin gene families which evolved from primordial genes by massive intron insertion before or during early vertebrate radiation. Introns placed at homologous positions in the gene sequences in combination with diagnostic sequence characters may also constitute a reliable kinship indicator for other protein superfamilies.     

Genome biology of the cyclostomes and insights into the evolutionary biology of vertebrate genomes

The jawless vertebrates (lamprey and hagfish) are the closest extant outgroups to all jawed vertebrates (gnathostomes) and can therefore provide critical insight into the evolution and basic biology of vertebrate genomes. As such, it is notable that the genomes of lamprey and hagfish possess a capacity for rearrangement that is beyond anything known from the gnathostomes. Like the jawed vertebrates, lamprey and hagfish undergo rearrangement of adaptive immune receptors. However, the receptors and the mechanisms for rearrangement that are utilized by jawless vertebrates clearly evolved independently of the gnathostome system. Unlike the jawed vertebrates, lamprey and hagfish also undergo extensive programmed rearrangements of the genome during embryonic development. By considering these fascinating genome biologies in the context of proposed (albeit contentious) phylogenetic relationships among lamprey, hagfish, and gnathostomes, we can begin to understand the evolutionary history of the vertebrate genome. Specifically, the deep shared ancestry and rapid divergence of lampreys, hagfish and gnathostomes is considered evidence that the two versions of programmed rearrangement present in lamprey and hagfish (embryonic and immune receptor) were present in an ancestral lineage that existed more than 400 million years ago and perhaps included the ancestor of the jawed vertebrates. Validating this premise will require better characterization of the genome sequence and mechanisms of rearrangement in lamprey and hagfish.     

The complete mitochondrial genome of the Chinese hook snout carp Opsariichthys bidens (Actinopterygii: Cypriniformes) and an alternative pattern of mitogenomic evolution in vertebrate

The complete mitochondrial genome sequence of the Chinese hook snout carp, Opsariichthys bidens, was newly determined using the long and accurate polymerase chain reaction method. The 16,611-nucleotide mitogenome contains 13 protein-coding genes, two rRNA genes (12S, 16S), 22 tRNA genes, and a noncoding control region. We use these data and homologous sequence data from multiple other ostariophysan fishes in a phylogenetic evaluation to test hypothesis pertaining to codon usage pattern of O. bidens mitochondrial protein genes as well as to re-examine the ostariophysan phylogeny. The mitochondrial genome of O. bidens reveals an alternative pattern of vertebrate mitochondrial evolution. For the mitochondrial protein genes of O. bidens, the most frequently used codon generally ends with either A or C, with C preferred over A for most fourfold degenerate codon families; the relative synonymous codon usage of G-ending codons is greatly elevated in all categories. The codon usage pattern of O. bidens mitochondrial protein genes is remarkably different from the general pattern found previously in the relatively closely related zebrafish and most other vertebrate mitochondria. Nucleotide bias at third codon positions is the main cause of codon bias in the mitochondrial protein genes of O. bidens, as it is biased particularly in favor of C over A. Bayesian analysis of 12 concatenated mitochondrial protein sequences for O. bidens and 46 other teleostean taxa supports the monophyly of Cypriniformes and Otophysi and results in a robust estimate of the otophysan phylogeny.     

Rapid evolutionary emergence of the combinatorial recognition repertoire

Although the capacity of cells to respond to environmental challengessuch as oxidative damage are ancient evolutionary developmentsthat have been carried through to modern higher vertebratesas "innate" immunity, the characteristic immune response ofvertebrates is a relatively recent evolutionary developmentthat is present only in jawed vertebrates. The vertebrate "combinatorial"response is defined by the presence of lymphocytes as specificantigen recognition cells and by the complete panel of antibodies,T cell receptors, and major histocompatibility complex moleculesall of which are members of the immunoglobulin family. Its emergencein evolution was an extremely rapid event (approximately 10million years) that was catalyzed by the horizontal transferof recombinase activator genes (RAG) from microbes to an ancestraljawed vertebrate. RAGs occur in jawed vertebrates, but havenot been found in invertebrates and other intermediate species.We propose that antigen recognition capacity contributed bythis novel combinatorial mechanism gave jawed vertebrates theability to recognize the entire range of potential antigenicmolecular structures, including self components and moleculesof infectious microbes not shared with vertebrates. The contrastwithin the vertebrates is striking because the most ancientextant jawed vertebrates, sharks and their kin, have the completepanoply of T-cell receptors, antibodies, MHC products and RAGgenes, whereas agnathans possess cells resembling lymphocytesbut ostensibly lack all of the molecules definitive of combinatorialimmunity. Another vertebrate innovation may have been the utilizationof nuclear receptor superfamily, in the regulation of lymphocytesand other cells of the immune lineage. Unlike, RAG, however,this superfamily occurs in all metazoans with the exceptionof sponges.     

Molecular evolution of cytochrome c oxidase subunit I in primates: is there coevolution between mitochondrial and nuclear genomes?

Phylogenetic analyses carried out on cytochrome c oxidase (COX) subunit I mitochondrial genes from 14 primates representing the major branches of the order and four outgroup nonprimate eutherians revealed that transversions and amino acid replacements (i.e., the more slowly occurring sequence changes) contained lower levels of homoplasy and thus provided more accurate information on cladistic relationships than transitions (i.e., the more rapidly occurring sequence changes). Several amino acids, each with a high likelihood of functionality involving the binding of cytochrome c or interaction with COX VIII, have changed in Anthropoidea, the primate suborder grouping New World monkey, Old World monkey, ape, and human lineages. They are conserved in other mammalian lineages and in nonanthropoid primates. Maximum-likelihood ancestral COX I nucleotide sequences were determined utilizing a near most parsimonious branching arrangement for the primate sequences that was consistent with previously hypothesized primate cladistic relationships based on larger and more diverse data sets. Relative rate tests of COX I mitochondrial sequences showed an elevated nonsynonymous (N) substitution rate for anthropoid-nonanthropoid comparisons. This finding for the largest mitochondrial (mt) DNA-encoded subunit is consistent with previous observations of elevated nonsynonymous substitution/synonymous substitution (S) rates in primates for mt-encoded COX II and for the nuclear-encoded COX IV and COX VIIa-H. Other COX-related proteins, including cytochrome c and cytochrome b, also show elevated amino acid replacement rates or N/S during similar time frames, suggesting that this group of interacting genes is likely to have coevolved during primate evolution.