Heterogeneous natural selection on oxidative phosphorylation genes among fishes with extreme high and low aerobic performance.

BACKGROUND: Oxidative phosphorylation (OXPHOS) is the primary source of ATP in eukaryotes and serves as a mechanistic link between variation in genotypes and energetic phenotypes. While several physiological and anatomical factors may lead to increased aerobic capacity, variation in OXPHOS proteins may influence OXPHOS efficiency and facilitate adaptation in organisms with varied energy demands. Although there is evidence that natural selection acts on OXPHOS genes, the focus has been on detection of directional (positive) selection on specific phylogenetic branches where traits that increase energetic demands appear to have evolved. We examined patterns of selection in a broader evolutionary context, i.e., on multiple lineages of fishes with extreme high and low aerobic performance. RESULTS: We found that patterns of natural selection on mitochondrial OXPHOS genes are complex among fishes with different swimming performance. Positive selection is not consistently associated with high performance taxa and appears to be strongest on lineages containing low performance taxa. In contrast, within high performance lineages, purifying (negative) selection appears to predominate. CONCLUSIONS: We provide evidence that selection on OXPHOS varies in both form and intensity within and among lineages through evolutionary time. These results provide evidence for fluctuating selection on OXPHOS associated with divergence in aerobic performance. However, in contrast to previous studies, positive selection was strongest on low performance taxa suggesting that adaptation of OXPHOS involves many factors beyond enhancing ATP production in high performance taxa. The broader pattern indicates a complex interplay between organismal adaptations, ATP demand, and OXPHOS function.

Mitochondrial-nuclear interactions: compensatory evolution or variable functional constraint among vertebrate oxidative phosphorylation genes?

Oxidative phosphorylation (OXPHOS), the major energy-producing pathway in aerobic organisms, includes protein subunits encoded by both mitochondrial (mt) and nuclear (nu) genomes. How these independent genomes have coevolved is a long-standing question in evolutionary biology. Although mt genes evolve faster than most nu genes, maintenance of OXPHOS structural stability and functional efficiency may involve correlated evolution of mt and nu OXPHOS genes. The nu OXPHOS genes might be predicted to exhibit accelerated evolutionary rates to accommodate the elevated substitution rates of mt OXPHOS subunits with which they interact. Evolutionary rates of nu OXPHOS genes should, therefore, be higher than that of nu genes that are not involved in OXPHOS (nu non-OXPHOS). We tested the compensatory evolution hypothesis by comparing the evolutionary rates (synonymous substitution rate dS and nonsynonymous substitution rate dN) among 13 mt OXPHOS genes, 60 nu OXPHOS genes, and 77 nu non-OXPHOS genes in vertebrates (7 fish and 40 mammal species). The results from a combined analysis of all OXPHOS subunits fit the predictions of the hypothesis. However, results from two OXPHOS complexes did not fit this pattern when analyzed separately. We found that the d(N) of nu OXPHOS genes for "core" subunits (those involved in the major catalytic activity) was lower than that of "noncore" subunits, whereas there was no significant difference in d(N) between genes for nu non-OXPHOS and core subunits. This latter finding suggests that compensatory changes play a minor role in the evolution of OXPHOS genes and that the observed accelerated nu substitution rates are due largely to reduced functional constraint on noncore subunits.

Multi-locus phylogenetic analysis reveals the pattern and tempo of bony fish evolution.

Over half of all vertebrates are "fishes", which exhibit enormous diversity in morphology, physiology, behavior, reproductive biology, and ecology. Investigation of fundamental areas of vertebrate biology depend critically on a robust phylogeny of fishes, yet evolutionary relationships among the major actinopterygian and sarcopterygian lineages have not been conclusively resolved. Although a consensus phylogeny of teleosts has been emerging recently, it has been based on analyses of various subsets of actinopterygian taxa, but not on a full sample of all bony fishes. Here we conducted a comprehensive phylogenetic study on a broad taxonomic sample of 61 actinopterygian and sarcopterygian lineages (with a chondrichthyan outgroup) using a molecular data set of 21 independent loci. These data yielded a resolved phylogenetic hypothesis for extant Osteichthyes, including 1) reciprocally monophyletic Sarcopterygii and Actinopterygii, as currently understood, with polypteriforms as the first diverging lineage within Actinopterygii; 2) a monophyletic group containing gars and bowfin (= Holostei) as sister group to teleosts; and 3) the earliest diverging lineage among teleosts being Elopomorpha, rather than Osteoglossomorpha. Relaxed-clock dating analysis employing a set of 24 newly applied fossil calibrations reveals divergence times that are more consistent with paleontological estimates than previous studies. Establishing a new phylogenetic pattern with accurate divergence dates for bony fishes illustrates several areas where the fossil record is incomplete and provides critical new insights on diversification of this important vertebrate group.

The tree of life and a new classification of bony fishes.

The tree of life of fishes is in a state of flux because we still lack a comprehensive phylogeny that includes all major groups. The situation is most critical for a large clade of spiny-finned fishes, traditionally referred to as percomorphs, whose uncertain relationships have plagued ichthyologists for over a century. Most of what we know about the higher-level relationships among fish lineages has been based on morphology, but rapid influx of molecular studies is changing many established systematic concepts. We report a comprehensive molecular phylogeny for bony fishes that includes representatives of all major lineages. DNA sequence data for 21 molecular markers (one mitochondrial and 20 nuclear genes) were collected for 1410 bony fish taxa, plus four tetrapod species and two chondrichthyan outgroups (total 1416 terminals). Bony fish diversity is represented by 1093 genera, 369 families, and all traditionally recognized orders. The maximum likelihood tree provides unprecedented resolution and high bootstrap support for most backbone nodes, defining for the first time a global phylogeny of fishes. The general structure of the tree is in agreement with expectations from previous morphological and molecular studies, but significant new clades arise. Most interestingly, the high degree of uncertainty among percomorphs is now resolved into nine well-supported supraordinal groups. The order Perciformes, considered by many a polyphyletic taxonomic waste basket, is defined for the first time as a monophyletic group in the global phylogeny. A new classification that reflects our phylogenetic hypothesis is proposed to facilitate communication about the newly found structure of the tree of life of fishes. Finally, the molecular phylogeny is calibrated using 60 fossil constraints to produce a comprehensive time tree. The new time-calibrated phylogeny will provide the basis for and stimulate new comparative studies to better understand the evolution of the amazing diversity of fishes.

Current and historical hybridization with differential introgression among three species of cyprinid fishes (genus Cyprinella).

Hybridization is common among freshwater fishes, particular among the Cyprinidae. We used two mitochondrial genes and one nuclear gene to characterize hybridization among two species pairs of Cyprinella in southwestern North America. Genalogical patterns revealed that C. lutrensis and C. venusta are currently hybridizing in several localities producing apparent F(1), F(2) and backcross generations, yet there was no evidence for introgression outside of local hybrid zones. Alternatively, mitochondrial haplotypes from C. lutrensis appear to have introgressed into a C. lepida population in the Nueces River completely replacing the native C. lepida haplotype. There was no evidence of introgression of nuclear DNA and there does not appear to be ongoing hybridization. The population of C. lepida from the nearby Frio River exhibits no evidence of hybridization with C. lutrensis. Thus, contact between C. lutrensis and C. venusta results in the formation of localized hybrid swarms, while contact between C. lutrensis and C. lepida has resulted in complete mitochondrial introgression in the Nueces River or no apparent hybridization in the Frio River. The three different outcomes of contact between these species illustrate the variable nature of interspecific reproductive interactions and provide an excellent system in which to better understand the factors influencing hybridization among freshwater fishes.

Diversification on an ecologically constrained adaptive landscape.

We used phylogenetic analysis of body-size ecomorphs in a crustacean species complex to gain insight into how spatial complexity of ecological processes generates and maintains biological diversity. Studies of geographically widespread species of Hyalella amphipods show that phenotypic evolution is tightly constrained in a manner consistent with adaptive responses to alternative predation regimes. A molecular phylogeny indicates that evolution of Hyalella ecomorphs is characterized by parallel evolution and by phenotypic stasis despite substantial levels of underlying molecular change. The phylogeny suggests that species diversification sometimes occurs by niche shifts, and sometimes occurs without a change in niche. Moreover, diversification in the Hyalella ecomorphs has involved the repeated evolution of similar phenotypic forms that exist in similar ecological settings, a hallmark of adaptive evolution. The evolutionary stasis observed in clades separated by substantial genetic divergence, but existing in similar habitats, is also suggestive of stabilizing natural selection acting to constrain phenotypic evolution within narrow bounds. We interpret the observed decoupling of genetic and phenotypic diversification in terms of adaptive radiation on an ecologically constrained adaptive landscape, and suggest that ecological constraints, perhaps acting together with genetic and functional constraints, may explain the parallel evolution and evolutionary stasis inferred by the phylogeny.

Why weight?

Whether phylogenetic data should be differentially or equally weighted is currently debated. Further, if differential weighting is to be explored, there is no consensus among investigators as to which weighting scheme is most appropriate. Mitochondrial genome data offer a powerful tool in assessment of differential weighting schemes because taxa can be selected from which a highly corroborated phylogeny is available (so that accuracy can be assessed), and it can be assumed that different data partitions share the same history (so that gene-sorting issues are not so problematic). Using mitochondrial data from 17 mammalian genomes, we evaluated the most commonly used weighting schemes, such as successive weighting, transversion weighting, codon-based weighting, and amino acid coding, and compared them to more complex weighting schemes including a 6-parameter weighting, pseudoreplicate reweighting, and tri-level weighting. We found that the most commonly used weighting schemes perform the worst with these data. Some of the more complex schemes perform well, however, none of them is consistently superior. These results support ones biases; if one has a predilection to avoid differential weighting, these data support equally weighted parsimony and maximum likelihood. Others might be encouraged by these results to try weighting as a form of data exploration.

Spatial covariation of mutation and nonsynonymous substitution rates in vertebrate mitochondrial genomes.

Mitochondrial genomes encode fundamental subunits of the basic energy producing machinery of eukaryotic cells that are under strong functional constraint. Paradoxically, these genes evolve rapidly in general, and there is substantial variation in evolutionary rates among genes within genomes. In order to investigate spatial variation in selection intensity, we conducted tests of neutrality using ratios of synonymous to nonsynonymous substitutions (dN/dS = omega) on numerous protein gene segments from fishes and mammals. Values of omega were very low for nearly all genomic regions. However, values of both omega and dN varied in a clinal pattern with increasing distance from the light-strand origin of replication. Spatial heterogeneity of nonsynonymous substitution rates exhibits a significantly positive correlation with variation in mutation rates that are related to the mode of mitochondrial DNA replication. The finding that nonsynonymous substitution rates are proportional to mutation rates is expected if a majority of substitutions are selectively neutral or slightly deleterious. Spatial patterns of among-gene variation in nonsynonymous rates were highly similar between fishes and mammals, suggesting that forces governing mitochondrial gene evolution have remained relatively constant over 450 Myr of vertebrate evolution. Conservation of substitution patterns despite major shifts in thermal habit and metabolic demands among taxa implicates a conserved replication mechanism controlling relative mutation rates as a major determinant of mitochondrial protein evolution.

Nuclear gene genealogies reveal historical, demographic and selective factors associated with speciation in field crickets.

Population genetics theory predicts that genetic drift should eliminate shared polymorphism, leading to monophyly or exclusivity of populations, when the elapsed time between lineage-splitting events is large relative to effective population size. We examined patterns of nucleotide variation in introns at four nuclear loci to relate processes affecting the history of genes to patterns of divergence among natural populations and species. Ancestral polymorphisms were shared among three recognized species, Gryllus firmus, G. pennsylvanicus, and G. ovisopis, and genealogical patterns suggest that successive speciation events occurred recently and rapidly relative to effective population size. High levels of shared polymorphism among these morphologically, behaviorally, and ecologically distinct species indicate that only a small fraction of the genome needs to become differentiated for speciation to occur. Among the four nuclear gene loci there was a 10-fold range in nucleotide diversity, and patterns of polymorphism and divergence suggest that natural selection has acted to maintain or eliminate variation at some loci. While nuclear gene genealogies may have limited applications in phylogeography or other approaches dependent on population monophyly, they provide important insights into the historical, demographic, and selective forces that shape speciation.

SIGNIFICANT ROLE FOR HISTORICAL EFFECTS IN THE EVOLUTION OF REPRODUCTIVE ISOLATION: EVIDENCE FROM PATTERNS OF INTROGRESSION BETWEEN THE CYPRINID FISHES, LUXILUS CORNUTUS AND LUXILUS CHRYSOCEPHALUS.

Samples of Luxilus cornutus, Luxilus chrysocephalus, and their hybrids were collected along hypothesized routes of dispersal from Pleistocene refugia to examine the significance of geographic variation in patterns of introgression between these species. Patterns of allozyme and mitochondrial DNA (mtDNA) variation were generally consistent with those from previous studies. Tests of Hardy-Weinberg equilibrium revealed significant deficiencies of heterozygotes in all samples, indicating some form of reproductive isolation. Mitochondrial DNAs of each species were not equally represented in F1 hybrids; however, this bias was eliminated when the two largest samples were excluded from the analysis. Backcross hybrids exhibited biased mtDNA introgression, as samples from Lake Erie (eastern) and Lake Michigan (western) drainages showed significant excesses of mtDNAs from L. chrysocephalus and L. cornutus, respectively, relative to frequencies of diagnostic allozyme markers. The extent and direction of allozyme and mtDNA introgression was quantified by calculating isolation index values from morphologically "pure" individuals of each species from each locality. Analysis of variance of these measures identified limited introgression of allozyme variants with no geographic pattern, but significant differences in direction of mtDNA introgression between drainages (i.e., postglacial dispersal route). Association between patterns of mtDNA introgression and dispersal route across the latitudinal width of the contact zone is best explained by genetic divergence during past isolation of ancestral populations from these drainages. These results identify a significant role for historical effects in the evolution of reproductive isolation and the process of speciation.