Molecular Evolution of Cytochrome c Oxidase in High-Performance Fish (Teleostei: Scombroidei)

The 13 peptides encoded by vertebrate mitochondrial DNA (mtDNA) are essential subunits of oxidative phosphorylation (OXPHOS) enzymes. These genes normally experience purifying selection and also coevolve with nuclear-encoded subunits of OXPHOS complexes. However, the role of positive selection on mtDNA evolution is still unclear, as most examples of intergenomic coevolution appear to be the result of compensation by nuclear-encoded genes for mildly deleterious mtDNA mutations, and not simultaneous positive selection in both genomes. Organisms that have experienced strong selective pressures to increase aerobic capacity or adapt to changes in thermal environment may be better candidates in which to examine the impact of positively selected changes on mtDNA evolution. The tuna (suborder Scombroidei, family Scombridae) and billfish (suborder Scombroidei, families Xiphiidae and Istiophoridae) are highly aerobic fish with multiple specializations in muscle energetics, including a high mitochondrial content and regional endothermy. We examined the role of positively selected mtDNA substitutions in the production of these unique phenotypes. Focusing on a catalytic subunit of cytochrome c oxidase (COX II), we found that the rate ratio of nonsynonymous (d(N); amino acid changing)-to-synonymous (d(S); silent) substitutions was not increased in lineages leading to the tuna but was significantly increased in the lineage preceding the billfish. Furthermore, there are a number of individual positively selected sites that, when mapped onto the COX crystal structure, appear to interact with other COX subunits and may affect OXPHOS function and regulation in billfish.     

Positive Darwinian selection operating on the immunoglobulin heavy chain of Antarctic fishes

The cooling of the Southern Ocean to the freezing point of seawater (-1.9 degrees C) over the past 25 million years played a dominant selective role in the evolution of the Antarctic fish fauna. During this period, the perciform suborder Notothenioidei, which is largely endemic to the Antarctic, diversified and developed numerous cold-adapted characters. In this report, we provide compelling evidence that the immunoglobulin heavy chain (IgH) of the notothenioid fishes has undergone adaptive selection. Two and four IgH clones were isolated, respectively, from spleen cDNA libraries prepared from the Antarctic icefish Chaenocephalus aceratus and the yellowbelly rockcod Notothenia coriiceps. The transmembrane region of the membrane form of the rockcod IgM heavy chain was located at the end of the second constant (C(H)) domain, in contrast to other teleost IgMs in which the transmembrane region is located at the end of the third constant domain. Phylogenetic analyses of C(H) regions revealed that rates of nonsynonymous nucleotide substitution were higher than rates of synonymous nucleotide substitution. Many of the nonsynonymous substitutions introduced charge changes, consistent with positive Darwinian selection acting to adapt the structure of the notothenioid immunoglobulins. The rates of nonsynonymous nucleotide substitutions were higher than the rates of synonymous nucleotide substitutions in complementarity determining regions of variable regions, suggesting that diversity at antigen binding sites is enhanced by genomic and/or somatic selection. Results of Southern blot hybridization experiments were consistent with a translocon type of IgH gene organization reminiscent of bony fishes and tetrapods.     

Evolution of glucagon genes

Statistical analyses of DNA sequences of the preproglucagon genes from bovine, human, hamster, and anglerfish suggest that a gene duplication creating two anglerfish genes (AF I and II) occurred about 160 Myr ago, long after the separation of fish and mammals. The analyses further suggest that the internal duplication producing the glucagon and glucagon-like peptide II (GLP-II) regions occurred about 1.2 billion years ago, which would indicate that the GLP-II region was present in the ancestral anglerfish sequence but was silenced or deleted before the gene duplication separating AF I and II. The glucagon-like peptide I (GLP-I) was derived from a duplication of the ancestral glucagon region about 800 Myr ago. The rate of synonymous substitution in these genes is approximately 4.3 x 10(-9) substitutions per year per synonymous site. The rate of nonsynonymous substitution in the signal peptide region is about 1.1 x 10(-9) substitutions per year per nonsynonymous site, a high rate comparable to that in the C-peptide region of preproinsulin. The rate of nonsynonymous substitution in the glicentin-related pancreatic polypeptide (GRPP) region is 0.63 x 10(-9) for the comparisons between mammalian species and 1.8 x 10(-9) for the comparisons between fish and mammals; the moderate rate in mammals suggests a physiological role for GRPP. The glucagon region is extremely conservative; no nonsynonymous substitution is observed in the mammalian genes, and a nonsynonymous rate of 0.18 x 10(-9) was obtained from the comparisons between fish and mammals. In the GLP-I region, the rate of nonsynonymous substitution was estimated to be 0.08 x 10(-9) for the comparisons between mammalian species and 0.30 x 10(- 9) for the comparisons between fish and mammals. In the GLP-II region, the rate was estimated to be 0.25 x 10(-9) for the comparisons between mammalian species. Thus, GLP-I and II are also very conservative, which suggests an important physiological role for these peptides.      

The "fish-specific" Hox cluster duplication is coincident with the origin of teleosts

The Hox gene complement of zebrafish, medaka, and fugu differs from that of other gnathostome vertebrates. These fishes have seven to eight Hox clusters compared to the four Hox clusters described in sarcopterygians and shark. The clusters in different teleost lineages are orthologous, implying that a "fish-specific" Hox cluster duplication has occurred in the stem lineage leading to the most recent common ancestor of zebrafish and fugu. The timing of this event, however, is unknown. To address this question, we sequenced four Hox genes from taxa representing basal actinopterygian and teleost lineages and compared them to known sequences from shark, coelacanth, zebrafish, and other teleosts. The resulting gene genealogies suggest that the fish-specific Hox cluster duplication occurred coincident with the origin of crown group teleosts. In addition, we obtained evidence for an independent Hox cluster duplication in the sturgeon lineage (Acipenseriformes). Finally, results from HoxA11 suggest that duplicated Hox genes have experienced diversifying selection immediately after the duplication event. Taken together, these results support the notion that the duplicated Hox genes of teleosts were causally relevant to adaptive evolution during the initial teleost radiation.     

Neuroglobins from the zebrafish Danio rerio and the pufferfish Tetraodon nigroviridis

Neuroglobin is a recently discovered respiratory, porphyrin-containing protein that is expressed in the brain of mouse and man. Here we show that neuroglobin is also present in the teleost fish. Complete cDNA sequences are reported from the pufferfish Tetraodon nigroviridis and the zebrafish Danio rerio. In addition, the neuroglobin gene of T. nigroviridis was sequenced, demonstrating the conservation of the B12.2, E11.0 and G7.0 introns plus the presence of an additional intron in the 5' noncoding region. The fish neuroglobins each comprise 159 amino acids and are 84.3% identical. Phylogenetic analyses show a basal position of the neuroglobins within the metazoan globin tree. An enhanced amino acid substitution rate was estimated for the fish neuroglobins ( approximately 0.93 x 10(-9) amino acid substitutions per site and year) compared with their mammalian proteins ( approximately 0.39 x 10(-9) replacements per site and year).     

Evidence for a period of directional selection following gene duplication in a neurally expressed locus of triosephosphate isomerase

A striking correlation between neural expression and high net negative charge in some teleost isozymes led to the interesting, yet untested, suggestion that negative charge represents an adaptation (via natural selection) to the neural environment. We examine the evolution of the triosephosphate isomerase (TPI) gene family in fishes for periods of positive selection. Teleost fish express two TPI proteins, including a generally expressed, neutrally charged isozyme and a neurally expressed, negatively charged isozyme; more primitive fish express only a single, generally expressed TPI isozyme. The TPI gene phylogeny constructed from sequences isolated from two teleosts, a single acipenseriform, and other TPI sequences from the databases, supports a single gene duplication event early in the evolution of bony fishes. Comparisons between inferred ancestral TPI sequences indicate that the neural TPI isozyme evolved through a period of positive selection resulting in the biased accumulation of negatively charged amino acids. Further, the number of nucleotide changes required for the observed amino acid substitutions suggests that selection acted on the overall charge of the protein and not on specific key amino acids.     

Atlantic Cod Piscidin and Its Diversification through Positive Selection

Piscidins constitute a family of cationic antimicrobial peptides that are thought to play an important role in the innate immune response of teleosts. On the one hand they show a remarkable diversity, which indicates that they are shaped by positive selection, but on the other hand they are ancient and have specific targets, suggesting that they are constrained by purifying selection. Until now piscidins had only been found in fish species from the superorder Acanthopterygii but we have recently identified a piscidin gene in Atlantic cod (Gadus morhua), thus showing that these antimicrobial peptides are not restricted to evolutionarily modern teleosts. Nucleotide diversity was much higher in the regions of the piscidin gene that code for the mature peptide and its pro domain than in the signal peptide. Maximum likelihood analyses with different evolution models revealed that the piscidin gene is under positive selection. Charge or hydrophobicity-changing amino acid substitutions observed in positively selected sites within the mature peptide influence its amphipathic structure and can have a marked effect on its function. This diversification might be associated with adaptation to new habitats or rapidly evolving pathogens.     

Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene.

Several codon-based models for the evolution of protein-coding DNA sequences are developed that account for varying selection intensity among amino acid sites. The "neutral model" assumes two categories of sites at which amino acid replacements are either neutral or deleterious. The "positive-selection model" assumes an additional category of positively selected sites at which nonsynonymous substitutions occur at a higher rate than synonymous ones. This model is also used to identify target sites for positive selection. The models are applied to a data set of the V3 region of the HIV-1 envelope gene, sequenced at different years after the infection of one patient. The results provide strong support for variable selection intensity among amino acid sites The neutral model is rejected in favor of the positive-selection model, indicating the operation of positive selection in the region. Positively selected sites are found in both the V3 region and the flanking regions.     

The DNA sequence of medaka chromosome LG22

We report the genomic DNA sequence of a single chromosome (linkage group 22; LG22) of the small teleost fish medaka (Oryzias latipes) as a first whole chromosome sequence from a non-mammalian vertebrate. The order and orientation of 633 protein-coding genes were deduced from 18,803,338 bp of DNA sequence, providing the opportunity to analyze chromosome evolution of vertebrate genomes by direct comparison with the human genome. The average number of genes in the "conserved gene cluster" (CGC), a strict definition of "synteny" at the sequence basis, between medaka and human was 1.6. These and other data suggest that approximately 38.8% of pair-wise gene relationships would have been broken from their common ancestor in the human and medaka lineages and further imply that approx 20,000 (15,520-23,280) breaks would have occurred from the entire genome of the common ancestor. These breaks were generated mainly by intra-chromosomal shufflings at a specific era in the vertebrate lineage. These precise comparative genomics allowed us to identify the pieces of ancient chromosomes of the common vertebrate ancestor and estimate chromosomal evolution in the vertebrate lineage.     

Duplication of phospholipase C-delta gene family in fish genomes

Fishes possess more genes than other vertebrates, possibly because of a genome duplication event during the evolution of the teleost (ray-finned) fish lineage. To further explore this idea, we cloned five genes encoding phosphoinositide-specific phospholipase C-delta (PLC-delta), designated respectively PoPLC-deltas, from olive flounder (Paralichthys olivaceus), and we performed phylogenetic analysis and sequence comparison to compare our putative gene products (PoPLC-deltas) with the sequences of known human PLC isoforms. The deduced amino acid sequences shared high sequence identity with human PLC-delta1, -delta3, and -delta4 isozymes and exhibited similar primary structures. In phylogenetic analysis of PoPLC-deltas with PLC-deltas of five teleost fishes (zebrafish, stickleback, medaka, Tetraodon, and Takifugu), three tetrapods (human, chicken, and frog), and two tunicates (sea squirt and pacific sea squirt), whose putative sequences of PLC-delta are available in Ensembl genome browser, the result also indicated that the two paralogous genes corresponding to each PLC-delta isoform originated from fish-specific genome duplication prior to the divergence of teleost fish. Our analyses suggest that an ancestral PLC-delta gene underwent three rounds of genome duplication during the evolution of vertebrates, leading to the six genes of three PLC-delta isoforms in teleost fish.     

Evidence for a divergence in function between two glucocorticoid receptors from a basal teleost

ABSTRACT: BACKGROUND: Duplicated glucocorticoid receptors (GR) are present in most teleost fish. The evolutionary advantage of retaining two GRs is unclear, as no subtype specific functional traits or physiological roles have been defined. To identify factors driving the retention of duplicate GRs in teleosts, the current study examined GRs in representatives of two basal ray-finned fish taxa that emerged either side of the teleost lineage whole genome duplication event (WGD) event, the acipenseriform, Acipenser ruthenus, (pre-WGD) and the osteoglossimorph, Pantodon buchholzi, (post-WGD). RESULTS: The study identified a single GR in A. ruthenus (ArGR) and two GRs in P. buchholzi (PbGR1 and PbGR2). Phylogenetic analyses showed that ArGR formed a distinct branch separate from the teleosts GRs. The teleost GR lineage was subdivded into two sublineages, each of which contained one of the two P. buchholzi GRs. ArGR, PbGR1 and PbGR2 all possess the unique 9 amino acid insert between the zinc-fingers of the DNA-binding domain that is present in one of the teleost GR lineages (GR1), but not the other (GR2). A splice variant of PbGR2 produces an isoform that lacked these 9 amino acids (PbGR2b). Cortisol stimulated transactivation activity of ArGR, PbGR2b and PbGR1 in vitro; with PbGR2b and PbGR1, the glucocorticoid 11-deoxycortisol was a more potent agonist than cortisol. The hormone sensitivity of PbGR2b and PbGR1 differed in the transactivation assay, with PbGR2b having lower EC50 values and greater fold induction. CONCLUSIONS: The difference in transactivation activity sensitivity between duplicated GRs of P. buchholzi suggests potential functional differences between the paralogs emerged early in the teleost lineage. Given the pleiotropic nature of GR function in vertebrates, this finding is in accordance with the hypothesis that duplicated GRs were potentially retained through subfunctionalisation followed by gene sharing. A 9 amino acid insert in the DNA-binding domain emerged in basal ray-finned fish GRs. However, the presence of a PbGR2 splice variant that lacks this insert, as well as the loss of the exon encoding these amino acids in the genes encoding for other teleost GR2 suggests the selection of two receptors with different DNA-binding domain structures in teleosts.     

Genome evolution and biodiversity in teleost fish

Teleost fish, which roughly make up half of the extant vertebrate species, exhibit an amazing level of biodiversity affecting their morphology, ecology and behaviour as well as many other aspects of their biology. This huge variability makes fish extremely attractive for the study of many biological questions, particularly of those related to evolution. New insights gained from different teleost species and sequencing projects have recently revealed several peculiar features of fish genomes that might have played a role in fish evolution and speciation. There is now substantial evidence that a round of tetraploidization/rediploidization has taken place during the early evolution of the ray-finned fish lineage, and that hundreds of duplicate pairs generated by this event have been maintained over hundreds of millions of years of evolution. Differential loss or subfunction partitioning of such gene duplicates might have been involved in the generation of fish variability. In contrast to mammalian genomes, teleost genomes also contain multiple families of active transposable elements, which might have played a role in speciation by affecting hybrid sterility and viability. Finally, the amazing diversity of sex determination systems and the plasticity of sex chromosomes observed in teleost might have been involved in both pre- and postmating reproductive isolation. Comparison of data generated by current and future genome projects as well as complementary studies in other species will allow one to approach the molecular and evolutionary mechanisms underlying genome diversity in fish, and will certainly significantly contribute to our understanding of gene evolution and function in humans and other vertebrates.     

Degree of Selective Constraint as an Explanation of the Different Rates of Evolution of Gender-Specific Mitochondrial DNA Lineages in the Mussel Mytilus

Mussels of the genus Mytilus segregate for a maternally transmitted F lineage and a paternally transmitted M lineage of mitochondrial DNA. Previous studies demonstrated that these lineages are older than the species of the M. edulis complex and that the M lineage evolves faster than the F lineage. Here we show that the latter observation also applies to a region of the molecule with no assigned function. Sequence data for the mitochondrial COIII gene and the ``unassigned' region of the F and M lineages of M. edulis and M. trossulus are used to evaluate various hypotheses that may account for the faster rate of evolution of the M lineage. Tests based on the proportion of synonymous and nonsynonymous substitutions suggest that the M lineage experiences relatively relaxed selection. Further support for this hypothesis comes from an examination of COIII amino acid substitutions at sites defined as either conserved or variable based on the pattern of variation in other mollusks and Drosophila. Most substitutions in the M lineage occur in regions that are also variable among non-Mytilus taxa. We suggest that these differences in selection pressure are a consequence of doubly uniparental mitochondrial DNA transmission in Mytilus.     

Variable evolutionary rates in the molecular evolution of mammalian growth hormones

In mammals pituitary growth hormone (GH) shows a slow basal rate of evolution (0.22 ± 0.03 × 10^–9 substitutions/amino acid site/year) which appears to have increased by at least 25–50-fold on two occasions, during the evolution of primates (to at least 10.8 ± 1.3 X 10^–9 substitutions/amino acid site/year) and artiodactyl ruminants (to at least 5.6 ± 1.3 X 10^–9 substitutions/amino acid site/year). That these rate increases are real, and not due to inadvertent comparison of nonorthologous genes, was established by showing that features of the GH gene sequences that are not expressed as mature hormone do not show corresponding changes in evolutionary rate. Thus, analysis of nonsynonymous substitutions in the coding sequence for the mature protein confirmed the rate increases seen in the primate and ruminant GHs, but analysis of nonsynonymous substitutions in the signal peptide sequence, synonymous substitutions in the coding sequence for signal peptide or mature protein, and 5 and 3 untranslated sequences showed no statistically significant changes in evolutionary rate. Evidence that the increases in evolutionary rate are probably due to positive selection is provided by the observation that in the cases of both ruminant and primate GHs the periods of rapid evolution were followed by a return to a slow rate similar to the basal rate seen in other mammalian GHs. Differences between the biological properties of GHs have been identified which may relate to these periods of rapid adaptive molecular evolution. On the basis of sequence data currently available (but excluding rodent GHs which show an intermediate rate, the basis of which is not clear) for most (90%) of evolutionary time mammalian GHs have been in the slow phase of evolution, with possibly most of the few amino acid substitutions that have occurred being neutral in nature. But most (80%) of the amino acid substitutions that have been introduced into GH during the course of mammalian evolution have been accepted during the rapid phases and were adaptive in nature.     

Effects of gender on exercise-induced growth hormone release.

We examined gender differences in growth hormone (GH) secretion during rest and exercise. Eighteen subjects (9 women and 9 men) were tested on two occasions each [resting condition (R) and exercise condition (Ex)]. Blood was sampled at 10-min intervals from 0600 to 1200 and was assayed for GH by chemiluminescence. At R, women had a 3.69-fold greater mean calculated mass of GH secreted per burst compared with men (5.4 +/- 1.0 vs. 1.7 +/- 0.4 microg/l, respectively) and higher basal (interpulse) GH secretion rates, which resulted in greater GH production rates and serum GH area under the curve (AUC; 1,107 +/- 194 vs. 595 +/- 146 microg x l(-1) x min, women vs. men; P = 0.04). Compared with R, Ex resulted in greater mean mass of GH secreted per burst, greater mean GH secretory burst amplitude, and greater GH AUC (1,196 +/- 211 vs. 506 +/- 90 microg x l(-1) x min, Ex vs. R, respectively; P < 0.001). During Ex, women attained maximal serum GH concentrations significantly earlier than men (24 vs. 32 min after initiation of Ex, respectively; P = 0.004). Despite this temporal disparity, both genders had similar maximal serum GH concentrations. The change in AUC (adjusted for unequal baselines) was similar for men and women (593 +/- 201 vs. 811 +/- 268 microg x l(-1) x min), but there were significant gender-by-condition interactive effects on GH secretory burst mass, pulsatile GH production rate, and maximal serum GH concentration. We conclude that, although women exhibit greater absolute GH secretion rates than men both at rest and during exercise, exercise evokes a similar incremental GH response in men and women. Thus the magnitude of the incremental secretory GH response is not gender dependent.     

Role of selection in fixation of gene duplications

New genes commonly appear through complete or partial duplications of pre-existing genes. Duplications of long DNA segments are constantly produced by rare mutations, may become fixed in a population by selection or random drift, and are subject to divergent evolution of the paralogous sequences after fixation, although gene conversion can impede this process. New data shed some light on each of these processes. Mutations which involve duplications can occur through at least two different mechanisms, backward strand slippage during DNA replication and unequal crossing-over. The background rate of duplication of a complete gene in humans is 10(-9)-10(-10) per generation, although many genes located within hot-spots of large-scale mutation are duplicated much more often. Many gene duplications affect fitness strongly, and are responsible, through gene dosage effects, for a number of genetic diseases. However, high levels of intrapopulation polymorphism caused by presence or absence of long, gene-containing DNA segments imply that some duplications are not under strong selection. The polymorphism to fixation ratios appear to be approximately the same for gene duplications and for presumably selectively neutral nucleotide substitutions, which, according to the McDonald-Kreitman test, is consistent with selective neutrality of duplications. However, this pattern can also be due to negative selection against most of segregating duplications and positive selection for at least some duplications which become fixed. Patterns in post-fixation evolution of duplicated genes do not easily reveal the causes of fixations. Many gene duplications which became fixed recently in a variety of organisms were positively selected because the increased expression of the corresponding genes was beneficial. The effects of gene dosage provide a unified framework for studying all phases of the life history of a gene duplication. Application of well-known methods of evolutionary genetics to accumulating data on new, polymorphic, and fixed duplication will enhance our understanding of the role of natural selection in the evolution by gene duplication.     

Phylogeny and substitution rates of angiosperm actin genes

Forty-four actin genes from five angiosperm species were PCR-amplified, cloned, and sequenced. Phylogenetic analysis of 34 of these actins, along with those previously published, indicates that angiosperm actin genes are monophyletic and underwent several duplications during evolution. Orthologues have been identified between Solanaceae species, as well as between Solanaceae species and soybean. These sequences were used to calculate nucleotide substitution rates. The synonymous rate (6.96 x 10(-9) substitutions/site/year) is similar to that of other nuclear protein-coding genes, but the nonsynonymous rate (0.19 x 10(-9) substitutions/site/year) is 6-19 times higher than that of mammalian actin genes. Relative rate tests indicate that actin genes are evolving at similar rates in monocots and in dicots. Evidence is also presented that some members of the maize actin multigene family have been involved in gene conversion events, that the potato genome contains 24 +/- 12 actin genes, and that potato and tomato diverged 11.6 +/- 3.6 MYA.