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

Bioinformatics and evolution of vertebrate nociceptin and opioid receptors

G protein-coupled receptors (GPCRs) are ancestrally related membrane proteins on cells that mediate the pharmacological effect of most drugs and neurotransmitters. GPCRs are the largest group of membrane receptor proteins encoded in the human genome. One of the most famous types of GPCRs is the opioid receptors. Opioid family receptors consist of four closely related proteins expressed in all vertebrate brains and spinal cords examined to date. The three classical types of opioid receptors shown unequivocally to mediate analgesia in animal models and in humans are the mu- (MOR), delta- (DOR), and kappa-(KOR) opioid receptor proteins. The fourth and most recent member of the opioid receptor family discovered is the nociceptin or orphanin FQ receptor (ORL). The role of ORL and its ligands in producing analgesia is not as clear, with both analgesic and hyperalgesic effects reported. All four opioid family receptor genes were cloned from expressed mRNA in a number of vertebrate species, and there are enough sequences presently available to carry out bioinformatic analysis. This chapter presents the results of a comparative analysis of vertebrate opioid receptors using pharmacological studies, bioinformatics, and the latest data from human whole-genome studies. Results confirm our initial hypotheses that the four opioid receptor genes most likely arose by whole-genome duplication, that there is an evolutionary vector of opioid receptor type divergence in sequence and function, and that the hMOR gene shows evidence of positive selection or adaptive evolution in Homo sapiens.

Directional Darwinian Selection in proteins

Background

Molecular evolution is a very active field of research, with several complementary approaches, including dN/dS, HON90, MM01, and others. Each has documented strengths and weaknesses, and no one approach provides a clear picture of how natural selection works at the molecular level. The purpose of this work is to present a simple new method that uses quantitative amino acid properties to identify and characterize directional selection in proteins.

Methods

Inferred amino acid replacements are viewed through the prism of a single physicochemical property to determine the amount and direction of change caused by each replacement. This allows the calculation of the probability that the mean change in the single property associated with the amino acid replacements is equal to zero (H₀: μ = 0; i.e., no net change) using a simple two-tailed t-test.

Results

Example data from calanoid and cyclopoid copepod cytochrome oxidase subunit I sequence pairs are presented to demonstrate how directional selection may be linked to major shifts in adaptive zones, and that convergent evolution at the whole organism level may be the result of convergent protein adaptations.

Conclusions

Rather than replace previous methods, this new method further complements existing methods to provide a holistic glimpse of how natural selection shapes protein structure and function over evolutionary time.

Positive selection in glycolysis among Australasian stick insects

BACKGROUND

The glycolytic pathway is central to cellular energy production. Selection on individual enzymes within glycolysis, particularly phosphoglucose isomerase (Pgi), has been associated with metabolic performance in numerous organisms. Nonetheless, how whole energy-producing pathways evolve to allow organisms to thrive in different environments and adopt new lifestyles remains little explored. The Lanceocercata radiation of Australasian stick insects includes transitions from tropical to temperate climates, lowland to alpine habitats, and winged to wingless forms. This permits a broad investigation to determine which steps within glycolysis and what sites within enzymes are the targets of positive selection. To address these questions we obtained transcript sequences from seven core glycolysis enzymes, including two Pgi paralogues, from 29 Lanceocercata species.

RESULTS

Using maximum likelihood methods a signature of positive selection was inferred in two core glycolysis enzymes. Pgi and Glyceraldehyde 3-phosphate dehydrogenase (Gaphd) genes both encode enzymes linking glycolysis to the pentose phosphate pathway. Positive selection among Pgi paralogues and orthologues predominately targets amino acids with residues exposed to the protein's surface, where changes in physical properties may alter enzyme performance.

CONCLUSION

Our results suggest that, for Lancerocercata stick insects, adaptation to new stressful lifestyles requires a balance between maintaining cellular energy production, efficiently exploiting different energy storage pools and compensating for stress-induced oxidative damage.

Cellular and biochemical characterization of two closely related triosephosphate isomerases from Trichomonas vaginalis

The glycolytic enzyme triosephosphate isomerase catalyses the isomerization between glyceraldehyde 3-phosphate and dihydroxyacetone phosphate. Here we report that Trichomonas vaginalis contains 2 fully functional tpi genes. Both genes are located in separated chromosomal context with different promoter regulatory elements and encode ORFs of 254 amino acids; the only differences between them are the character of 4 amino acids located in α-helices 1, 2 and 8. Semi-quantitative RT-PCR assays showed that tpi2 transcript is approximately 3·3-fold more abundant than tpi1. Using an anti-TvTIM2 polyclonal antibody it was demonstrated that TIM proteins have a cytoplasmic localization and both enzymes are able to complement an Escherichia coli strain carrying a deletion of its endogenous tpi gene. Both TIM proteins assemble as dimers and their secondary structure assessment is essentially identical to TIM from Saccharomyces cerevisiae. The kinetic catalytic constants of the recombinant enzymes using glyceraldehyde-3-phosphate as substrate are similar to the catalytic constants of TIMs from other organisms including parasitic protozoa. As T. vaginalis depends on glycolysis for ATP production, we speculate 2 possible reasons to maintain a duplicated tpi copy on its genome: an increase in gene dosage or an early event of neofunctionalization of TIM as a moonlighting protein.

Dynamic programming procedure for searching optimal models to estimate substitution rates based on the maximum-likelihood method

The substitution rate in a gene can provide valuable information for understanding its functionality and evolution. A widely used method to estimate substitution rates is the maximum-likelihood method implemented in the CODEML program in the PAML package. A limited number of branch models, chosen based on a priori information or an interest in a particular lineage(s), are tested, whereas a large number of potential models are neglected. A complementary approach is also needed to test all or a large number of possible models to search for the globally optional model(s) of maximum likelihood. However, the computational time for this search even in a small number of sequences becomes impractically long. Thus, it is desirable to explore the most probable spaces to search for the optimal models. Using dynamic programming techniques, we developed a simple computational method for searching the most probable optimal branch-specific models in a practically feasible computational time. We propose three search methods to find the optimal models, which explored O(n) (method 1) to O(n(2)) (method 2 and method 3) models when the given phylogeny has n branches. In addition, we derived a formula to calculate the number of all possible models, revealing the complexity of finding the optimal branch-specific model. We show that in a reanalysis of over 50 previously published studies, the vast majority obtained better models with significantly higher likelihoods than the conventional hypothesis model methods.

A coleopteran triosephosphate isomerase: X-ray structure and phylogenetic impact of insect sequences

A coleopteran triosephosphate isomerase (TIM) from Tenebrio molitor (yellow mealworm beetle) was recombinantly expressed in Escherichia coli and characterized with respect to thermal stability, kinetic parameters and oligomeric state. The enzyme was successfully crystallized and the structure determined by X-ray analysis to 2.0 A resolution. This is the first example of an invertebrate TIM. We compare structural features with known structures of TIMs from microorganisms, plants and vertebrates, and discuss the utility of the Tenebrio TIM sequence, together with several newly sequenced insect TIMs, for molecular phylogenetic analysis.

BD SIMULATIONS OF THE IONIC STRENGTH DEPENDENCE OF THE INTERACTIONS BETWEEN TRIOSE PHOSPHATE ISOMERASE AND F-ACTIN

Functional protein-protein interactions are essential for many physiological processes. For example, the association of glycolytic enzymes to F-actin is proposed to be one mechanism through which glycolytic enzymes are compartmentalized, and as a result, play essential roles such as regulation of the glycolytic pathway and increasing glycolytic flux. Many glycolytic enzymes including fructose-1,6-bisphophate aldolase, glyceraldedhye-3-phosphate dehydrogenase, and lactate dehydrogenase bind F-actin strongly. Other glycolytic enzymes including triose phosphate isomerase (TPI) do not interact with F-actin significantly. Herein, Brownian dynamics (BD) simulations determine the energetics of the association of F-actin with the glycolytic enzyme triose phosphate isomerase as a function of ionic strength. This is the first thorough control study examining how well BD reproduces the experimental observations that the binding of TPI to F-actin is very weak and falls off rapidly as ionic strength increases. The BD results confirm experimental observations that the degree of association diminishes as ionic strength increases and that the interaction of TPI with F-actin is weakly nonspecific to nonexistent.

Adaptive evolution of young gene duplicates in mammals

Duplicate genes act as a source of genetic material from which new functions arise. They exist in large numbers in every sequenced eukaryotic genome and may be responsible for many differences in phenotypes between species. However, recent work searching for the targets of positive selection in humans has largely ignored duplicated genes due to complications in orthology assignment. Here we find that a high proportion of young gene duplicates in the human, macaque, mouse, and rat genomes have experienced adaptive natural selection. Approximately 10% of all lineage-specific duplicates show evidence for positive selection on their protein sequences, larger than any reported amount of selection among single-copy genes in these lineages using similar methods. We also find that newly duplicated genes that have been transposed to new chromosomal locations are significantly more likely to have undergone positive selection than the ancestral copy. Human-specific duplicates evolving under adaptive natural selection include a surprising number of genes involved in neuronal and cognitive functions. Our results imply that genome scans for selection that ignore duplicated loci are missing a large fraction of all adaptive substitutions. The results are also in agreement with the classical model of evolution by gene duplication, supporting a common role for neofunctionalization in the long-term maintenance of gene duplicates.

Cross-species analysis of the glycolytic pathway by comparison of molecular interaction fields

The electrostatic potential of an enzyme is a key determinant of its substrate interactions and catalytic turnover. Here we invoke comparative analysis of protein electrostatic potentials, along with sequence and structural analysis, to classify and characterize all the enzymes in an entire pathway across a set of different organisms. The electrostatic potentials of the enzymes from the glycolytic pathway of 11 eukaryotes were analyzed by qPIPSA (quantitative protein interaction property similarity analysis). The comparison allows the functional assignment of neuron-specific isoforms of triosephosphate isomerase from zebrafish, the identification of unusual protein surface interaction properties of the mosquito glucose-6-phosphate isomerase and the functional annotation of ATP-dependent phosphofructokinases and cofactor-dependent phosphoglycerate mutases from plants. We here show that plants possess two parallel pathways to convert glucose. One is similar to glycolysis in humans, the other is specialized to let plants adapt to their environmental conditions. We use differences in electrostatic potentials to estimate kinetic parameters for the triosephosphate isomerases from nine species for which published parameters are not available. Along the core glycolytic pathway, phosphoglycerate mutase displays the most conserved electrostatic potential. The largest cross-species variations are found for glucose-6-phosphate isomerase, enolase and fructose-1,6-bisphosphate aldolase. The extent of conservation of electrostatic potentials along the pathway is consistent with the absence of a single rate-limiting step in glycolysis.

Molecular evolution of the betagamma lens crystallin superfamily: evidence for a retained ancestral function in gamma N crystallins?

Within the vertebrate eye, betagamma crystallins are extremely stable lens proteins that are uniquely adapted to increase refractory power while maintaining transparency. Unlike alpha crystallins, which are well-characterized, multifunctional proteins that have important functions both in and out of the lens, betagamma lens crystallins are a diverse group of proteins with no clear ancestral or contemporary nonlens role. We carried out phylogenetic and molecular evolutionary analyses of the betagamma-crystallin superfamily in order to study the evolutionary history of the gamma N crystallins, a recently discovered, biochemically atypical family suggested to possess a divergent or ancestral function. By including nonlens, betagamma-motif-containing sequences in our analysis as outgroups, we confirmed the phylogenetic position of the gamma N family as sister to other gamma crystallins. Using maximum likelihood codon models to estimate lineage-specific nonsynonymous-to-synonymous rate ratios revealed strong positive selection in all of the early lineages within the betagamma family, with the striking exception of the lineage leading to the gamma N crystallins which was characterized by strong purifying selection. Branch-site analysis, used to identify candidate sites involved in functional divergence between gamma N crystallins and its sister clade containing all other gamma crystallins, identified several positively selected changes at sites of known functional importance in the betagamma crystallin protein structure. Further analyses of a fish-specific gamma N crystallin gene duplication revealed a more recent episode of positive selection in only one of the two descendant lineages (gamma N2). Finally, from the guppy, Poecilia reticulata, we isolated complete gamma N1 and gamma N2 coding sequence data from cDNA and partial coding sequence data from genomic DNA in order to confirm the presence of a novel gamma N2 intron, discovered through data mining of two pufferfish genomes. We conclude that the function of the gamma N family likely resembles the ancestral vertebrate betagamma crystallin more than other betagamma families. Furthermore, owing to the presence of an additional intron in some fish gamma N2 crystallins, and the inferred action of positive selection following the fish-specific gamma N duplication, we suggest that further study of fish gamma N crystallins will be critical in further elucidating possible ancestral functions of gamma N crystallins and any nonstructural role they may have.

The evolution of vertebrate opioid receptors

The proteins that mediate the analgesic and other effects of opioid drugs and endogenous opioid peptides are known as opioid receptors. Opioid receptors consist of a family of four closely-related proteins belonging to the large superfamily of G-protein coupled receptors. The three types of opioid receptors shown unequivocally to mediate analgesia in animal models are the mu (MOR), delta (DOR), and kappa (KOR) opioid receptor proteins. The role of the fourth member of the opioid receptor family, the nociceptin or orphanin FQ receptor (ORL), is not as clear as hyperalgesia, analgesia, and no effect was reported after administration of ORL agonists. There are now cDNA sequences for all four types of opioid receptors that are expressed in the brain of six species from three different classes of vertebrates. This review presents a comparative analysis of vertebrate opioid receptors using bioinformatics and data from recent human genome studies. Results indicate that opioid receptors arose by gene duplication, that there is a vector of opioid receptor divergence, and that MOR shows evidence of rapid evolution.

Electric charge divergence in proteins: insights into the evolution of their three-dimensional properties

Previous studies of protein evolution have identified important mutations in various proteins that affect a small number of residues, but dramatically alter protein function. However, the evolutionary process underlying the three-dimensional protein properties, which are determined by a much larger number of residues, remains unclear. Based on a comparative evolutionary analysis of teleost phosphoglucose isomerases (PGIs; EC 5.3.1.9), we previously demonstrated that the relatively weak selection on many amino acid sites has played an important role in the evolution of protein electric charge as a model of three-dimensional protein properties. To ascertain the generality of this finding, we sought further evidence of this type of protein evolution. For this purpose, we analyzed the vertebrate isoforms of fructose-1,6-bisphosphate aldolase (ALD; EC 4.1.2.13), for which electric charges are known to have diverged after gene duplication. The results showed that the divergence in electric charge between the ALD isoforms was also driven by weak selection on many amino acid sites, as in PGI, confirming the generality of earlier findings. To obtain further insights, ALD and PGI were compared to the proteins pancreatic ribonuclease (EC 3.1.27.5) and triose-phosphate isomerase (EC 5.3.1.1), for which electric charges likely evolved through a well-defined mode of molecular evolution; namely, strong selection on specific amino acid sites. Comparison of the number and composition of amino acids on the protein surface suggested that the absolute number of evolutionarily changeable amino acids in a protein affects the strength of selection pressure acting on individual amino acid sites.

Highly divergent duplicate mannose-6-phosphate isomerase (Mpi) genes in the blue mussel, Mytilus edulis

Gene duplication has been hypothesized to play a major role in the evolution of genes and genomes and in generating phenotypic diversity among the proteins those genes encode. We have identified duplicate genes for the glycolytic enzyme mannose-6-phosphate isomerase (MPI: EC 5.3.1.8) in marine mussels in the genus Mytilus. Overall, there was only 52% sequence identity (72% similarity) between the proteins encoded by these two genes, designated as Mpi-A and Mpi-B. Based on a comparison of the rate of non-synonymous substitution between orthologous and paralogous Mpi coding sequences obtained from Mytilus edulis and the congener M. trossulus we estimate that the duplication of Mpi in mussels occurred ~170 MYA. We detected paralog-specific differences in the ratio of non-synonymous to synonymous substitutions (k(a)/k(s)) and in the predicted net charge for MPI-A and MPI-B. Using a real-time quantitative RT-PCR assay we observed substantial changes in Mpi-A and Mpi-B transcript levels between tissue types; the strongest expression of Mpi-A was observed in mantle tissue while Mpi-B expression exceeded that of Mpi-A in gill and hepatopancreas tissues. Taken together, these observations suggest that different functional roles have evolved for these two Mytilus Mpi genes subsequent to gene duplication.

Post-duplication charge evolution of phosphoglucose isomerases in teleost fishes through weak selection on many amino acid sites

Background

The partitioning of ancestral functions among duplicated genes by neutral evolution, or subfunctionalization, has been considered the primary process for the evolution of novel proteins (neofunctionalization). Nonetheless, how a subfunctionalized protein can evolve into a more adaptive protein is poorly understood, mainly due to the limitations of current analytical methods, which can detect only strong selection for amino acid substitutions involved in adaptive molecular evolution. In this study, we employed a comparative evolutionary approach to this question, focusing on differences in the structural properties of a protein, specifically the electric charge, encoded by fish-specific duplicated phosphoglucose isomerase (Pgi) genes.

Results

Full-length cDNA cloning, RT-PCR based gene expression analyses, and comparative sequence analyses showed that after subfunctionalization with respect to the expression organ of duplicate Pgi genes, the net electric charge of the PGI-1 protein expressed mainly in internal tissues became more negative, and that of PGI-2 expressed mainly in muscular tissues became more positive. The difference in net protein charge was attributable not to specific amino acid sites but to the sum of various amino acid sites located on the surface of the PGI molecule.

Conclusion

This finding suggests that the surface charge evolution of PGI proteins was not driven by strong selection on individual amino acid sites leading to permanent fixation of a particular residue, but rather was driven by weak selection on a large number of amino acid sites and consequently by steady directional and/or purifying selection on the overall structural properties of the protein, which is derived from many modifiable sites. The mode of molecular evolution presented here may be relevant to various cases of adaptive modification in proteins, such as hydrophobic properties, molecular size, and electric charge.

Cloning and bioinformatics of amphibian mu, delta, kappa, and nociceptin opioid receptors expressed in brain tissue: evidence for opioid receptor divergence in mammals

Opioid agonists produce analgesia in humans and other mammals by binding to three distinct types of G protein-coupled receptors; mu (MOR), delta (DOR), and kappa (KOR) opioid receptors. A fourth member of the opioid receptor family is the nociceptin or orphanin FQ receptor (ORL), however the role of the ORL receptor in analgesia is less clear. In the Northern grass frog, Rana pipiens, systemic and central administration of morphine and selective MOR, DOR, and KOR agonists produced dose-dependent antinociceptive effects blocked by the general opioid antagonist, naltrexone. The present study reports on the sequence, expression, and bioinformatics of four opioid receptor cDNAs cloned from Rana pipiens; rpMOR, rpDOR, rpKOR, and rpORL. These were the first opioid receptors cloned from a species of Class Amphibia, are selectively expressed in brain tissue, and show 70-84% identity to their homologous mammalian opioid receptors. Comparisons within species showed that MOR, DOR, and KOR proteins are significantly less divergent in earlier-evolved vertebrates compared to humans and other mammals. Among the four types of opioid receptors, MOR proteins show the least sequence variation among the six vertebrate species. Additionally, phylogenetic analysis supports the hypothesis that the family of opioid receptor proteins are coded by four genes that arose from two gene duplications of a single ancestral opioid receptor gene.

Three rounds (1R/2R/3R) of genome duplications and the evolution of the glycolytic pathway in vertebrates

Background

Evolution of the deuterostome lineage was accompanied by an increase in systematic complexity especially with regard to highly specialized tissues and organs. Based on the observation of an increased number of paralogous genes in vertebrates compared with invertebrates, two entire genome duplications (2R) were proposed during the early evolution of vertebrates. Most glycolytic enzymes occur as several copies in vertebrate genomes, which are specifically expressed in certain tissues. Therefore, the glycolytic pathway is particularly suitable for testing theories of the involvement of gene/genome duplications in enzyme evolution.

Results

We assembled datasets from genomic databases of at least nine vertebrate species and at least three outgroups (one deuterostome and two protostomes), and used maximum likelihood and Bayesian methods to construct phylogenies of the 10 enzymes of the glycolytic pathway. Through this approach, we intended to gain insights into the vertebrate specific evolution of enzymes of the glycolytic pathway. Many of the obtained gene trees generally reflect the history of two rounds of duplication during vertebrate evolution, and were in agreement with the hypothesis of an additional duplication event within the lineage of teleost fish. The retention of paralogs differed greatly between genes, and no direct link to the multimeric structure of the active enzyme was found.

Conclusion

The glycolytic pathway has subsequently evolved by gene duplication and divergence of each constituent enzyme with taxon-specific individual gene losses or lineage-specific duplications. The tissue-specific expression might have led to an increased retention for some genes since paralogs can subdivide the ancestral expression domain or find new functions, which are not necessarily related to the original function.

Varying signals of the effects of natural selection during teleost growth hormone gene evolution

The growth hormone (GH) gene of teleost fish exhibits a higher degree of variability compared with other vertebrate groups. However, the different selective constraints at the sequence level are not well understood. In this study, maximum-likelihood (ML) models of codon substitutions were used to investigate Darwinian adaptive evolution of the GH gene in teleost fishes. Complete GH gene sequences of 54 fish species were classified into 4 orders, and the variable nature of GH was examined by determining the dNand dSrate variation and the rates of molecular evolution for each teleost order. The results indicate that although the overall evolution rate for teleost GH is high ((1.15 ± 0.01) × 10–9substitutions/(aa site·y)) compared with the "slow phases" in mammals ((0.21 to 0.28 ± 0.05) × 10–9), the vital structure of this gene has been retained. While the majority of the amino acid changes appear to be due to relaxation of purifying selection, some positively selected sites were detected in regions with no specifically identified role in protein function. The positively selected regions observed in salmoniformes lineage suggests a possible role for positive selection driving functional divergence in paralogous forms of the GH gene after whole-genome duplication in this lineage.Key words: teleost fish, growth hormone, positive selection, synonymous substitution, non-synonymous substitution, molecular evolution.

Influence of hydrogeographic history and hybridization on the distribution of genetic variation in the pupfishes Cyprinodon atrorus and C. bifasciatus

The evolutionary importance of hybridization in animals has been subject of much debate. In this study, we examined the influence of hydrogeographic history and hybridization on the present distribution of nuclear and mitochondrial DNA variation in two pupfish species, Cyprinodon atrorus and Cyprinodon bifasciatus. Results presented here indicate that there has been limited introgression of nuclear genes; however, mtDNA introgression has been substantial, with complete replacement of the C. bifasciatus mitochondrial genome by that of C. atrorus. Subsequent to this replacement, there has been diversification of mitochondrial haplotypes along major geographic regions in the basin. Evidence was also found that mitochondrial replacement follows a predictable, cyclical pattern in this system, with isolation and diversification followed by re-contact and replacement of C. bifasciatus mitochondrial haplotypes by those of C. atrorus. This pattern is best explained by a combination of a numeric bias towards C. atrorus and mating site selection rather than selection for C. atrorus mitochondrial genome. These results demonstrate the important role hybridization can play in evolution.

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.

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.

A microsatellite genetic linkage map for Xiphophorus

Interspecies hybrids between distinct species of the genus Xiphophorus are often used in varied research investigations to identify genomic regions associated with the inheritance of complex traits. There are 24 described Xiphophorus species and a greater number of pedigreed strains; thus, the number of potential interspecies hybrid cross combinations is quite large. Previously, select Xiphophorus experimental crosses have been shown to exhibit differing characteristics between parental species and among the hybrid fishes derived from crossing them, such as widely differing susceptibilities to chemical or physical agents. For instance, genomic regions harboring tumor suppressor and oncogenes have been identified via linkage association of these loci with a small set of established genetic markers. The power of this experimental strategy is related to the number of genetic markers available in the Xiphophorus interspecies cross of interest. Thus, we have undertaken the task of expanding the suite of easily scored markers by characterization of Xiphophorus microsatellite sequences. Using a cross between Xiphophorus maculatus and X. andersi, we report a linkage map predominantly composed of microsatellite markers. All 24 acrocentric chromosome sets of Xiphophorus are represented in the assembled linkage map with an average intergenomic distance of 7.5 cM. Since both male and female F1 hybrids were used to produce backcross progeny, these recombination rates were compared between "male" and "female" maps. Although several genomic regions exhibit differences in map length, male- and female-derived maps are similar. Thus Xiphophorus, in contrast to zebrafish, Danio rerio, and several other vertebrate species, does not show sex-specific differences in recombination. The microsatellite markers we report can be easily adapted to any Xiphophorus interspecies and some intraspecies crosses, and thus provide a means to directly compare results derived from independent experiments.

Analytical Methods for Studying the Evolution of Paralogs Using Duplicate Gene Datasets

Gene duplication is widely viewed as an important source of raw material for functional innovation in proteins because at least some duplicate copies will evolve new or slightly modified functions. The study of the molecular processes by which functional innovation occurs interests both evolutionary biologists and protein chemists, and the development of methods to investigate these processes has led to a productive meeting of disciplines and an availability of complementary approaches for exploring datasets. This has resulted in insights into past events, prediction of current function, and prediction of future change. The methods fall broadly into two categories: those that rely on detection of shifts in selective constraints and those that rely on detection of correlations between molecular changes and functional shifts. Strengths and limitations of the methods are evaluated here in the context of the question being addressed, the input required, and the specific metric that is evaluated in each test.

Using models of nucleotide evolution to build phylogenetic trees

Molecular phylogenetics and its applications are popular and useful tools for making comparative investigations in genetics; however, estimating phylogenetic trees is not always straightforward. Some phylogenetic estimators use an explicit model of nucleotide evolution to estimate evolutionary parameters such as branch lengths and tree topology. There are many models to choose from, and use of the optimal model for a particular data set is important to avoid a loss of power and accuracy in phylogenetic estimations. Here, we review some molecular evolutionary forces and the parameters included in some common models of evolution used to interpret resulting patterns of molecular variation. We present some statistical methods of selecting a particular model of nucleotide evolution, and provide an empirical example of model selection. Statistical model selection strikes a balance between the bias introduced by some models and the increased variance of parameter estimates that results from using other models.

Neglected Taxonomy of Rare Desert Fishes: Congruent Evidence for Two Species of Leatherside Chub

Conservation biologists rely heavily on taxonomy to set the scope for biological monitoring and recovery planning of rare or threatened species. Yet, taxonomic boundaries are seldom evaluated as falsifiable hypotheses that can be statistically tested. Here, we examine species boundaries in leatherside chub (Teleostei, Cyprinidae), an imperiled desert fish native to the Bonneville Basin and upper Snake River drainages of western North America. Recent molecular data hint that this fish could be composed of two distinct taxa that are geographically separated into northern and southern species. To formally test this hypothesis, we evaluated leatherside chub using several different categories of species concepts, including criteria dependent on phylogenetic, morphological, and ecological data. We found that leatherside chub is composed of two reciprocally monophyletic clades (candidate species) characterized by numerous fixed genetic differences for both mitochondrial and nuclear DNA markers; mtDNA sequence divergence between the two clades approached 8%. The candidate species also showed significant differences in cranial shape, revealed by morphometric analysis. Finally, controlled growth and foraging experiments using representative populations from each clade show that candidate species appear to be locally adapted to the thermal environments where they now occur. Combined, these three lines of evidence support the hypothesis that leatherside chub is composed of two species. Moreover, all lines of evidence place these two species within the genusLepidomeda, a group consisting of three additional species of endangered spinedace fishes, and one extinct species, all native to the Colorado River system. Hence, we elevate the two clades of leatherside chub to distinct species status (Lepidomeda copeiin the north andL. aliciaein the south), and argue that each warrants independent conservation and recovery action.

Reconstructing patterns of reticulate evolution in plants

Until recently, rigorously reconstructing the many hybrid speciation events in plants has not been practical because of the limited number of molecular markers available for plant phylogenetic reconstruction and the lack of good, biologically based methods for inferring reticulation (network) events. This situation should change rapidly with the development of multiple nuclear markers for phylogenetic reconstruction and new methods for reconstructing reticulate evolution. These developments will necessitate a much greater incorporation of population genetics into phylogenetic reconstruction than has been common. Population genetic events such as gene duplication coupled with lineage sorting and meiotic and sexual recombination have always had the potential to affect phylogenetic inference. For tree reconstruction, these problems are usually minimized by using uniparental markers and nuclear markers that undergo rapid concerted evolution. Because reconstruction of reticulate speciation events will require nuclear markers that lack these characteristics, effects of population genetics on phylogenetic inference will need to be addressed directly. Current models and methods that allow hybrid speciation to be detected and reconstructed are discussed, with a focus on how lineage sorting and meiotic and sexual recombination affect network reconstruction. Approaches that would allow inference of phylogenetic networks in their presence are suggested.

Reconstructing patterns of reticulate evolution in plants

Until recently, rigorously reconstructing the many hybrid speciation events in plants has not been practical because of the limited number of molecular markers available for plant phylogenetic reconstruction and the lack of good, biologically based methods for inferring reticulation (network) events. This situation should change rapidly with the development of multiple nuclear markers for phylogenetic reconstruction and new methods for reconstructing reticulate evolution. These developments will necessitate a much greater incorporation of population genetics into phylogenetic reconstruction than has been common. Population genetic events such as gene duplication coupled with lineage sorting and meiotic and sexual recombination have always had the potential to affect phylogenetic inference. For tree reconstruction, these problems are usually minimized by using uniparental markers and nuclear markers that undergo rapid concerted evolution. Because reconstruction of reticulate speciation events will require nuclear markers that lack these characteristics, effects of population genetics on phylogenetic inference will need to be addressed directly. Current models and methods that allow hybrid speciation to be detected and reconstructed are discussed, with a focus on how lineage sorting and meiotic and sexual recombination affect network reconstruction. Approaches that would allow inference of phylogenetic networks in their presence are suggested.

Expression, Purification and the 1.8Å Resolution Crystal Structure of Human Neuron Specific Enolase

Human neuron-specific enolase (NSE) or isozyme gamma has been expressed with a C-terminal His-tag in Escherichia coli. The enzyme has been purified, crystallized and its crystal structure determined. In the crystals the enzyme forms the asymmetric complex NSE x Mg2 x SO4/NSE x Mg x Cl, where "/" separates the dimer subunits. The subunit that contains the sulfate (or phosphate) ion and two magnesium ions is in the closed conformation observed in enolase complexes with the substrate or its analogues; the other subunit is in the open conformation observed in enolase subunits without bound substrate or analogues. This indicates negative cooperativity for ligand binding between subunits. Electrostatic charge differences between isozymes alpha and gamma, -19 at physiological pH, are concentrated in the regions of the molecular surface that are negatively charged in alpha, i.e. surface areas negatively charged in alpha are more negatively charged in gamma, while areas that are neutral or positively charged tend to be charge-conserved.

Basal jawed vertebrate phylogeny inferred from multiple nuclear DNA-coded genes

Background

Phylogenetic analyses of jawed vertebrates based on mitochondrial sequences often result in confusing inferences which are obviously inconsistent with generally accepted trees. In particular, in a hypothesis by Rasmussen and Arnason based on mitochondrial trees, cartilaginous fishes have a terminal position in a paraphyletic cluster of bony fishes. No previous analysis based on nuclear DNA-coded genes could significantly reject the mitochondrial trees of jawed vertebrates.

Results

We have cloned and sequenced seven nuclear DNA-coded genes from 13 vertebrate species. These sequences, together with sequences available from databases including 13 jawed vertebrates from eight major groups (cartilaginous fishes, bichir, chondrosteans, gar, bowfin, teleost fishes, lungfishes and tetrapods) and an outgroup (a cyclostome and a lancelet), have been subjected to phylogenetic analyses based on the maximum likelihood method.

Conclusion

Cartilaginous fishes have been inferred to be basal to other jawed vertebrates, which is consistent with the generally accepted view. The minimum log-likelihood difference between the maximum likelihood tree and trees not supporting the basal position of cartilaginous fishes is 18.3 ± 13.1. The hypothesis by Rasmussen and Arnason has been significantly rejected with the minimum log-likelihood difference of 123 ± 23.3. Our tree has also shown that living holosteans, comprising bowfin and gar, form a monophyletic group which is the sister group to teleost fishes. This is consistent with a formerly prevalent view of vertebrate classification, although inconsistent with both of the current morphology-based and mitochondrial sequence-based trees. Furthermore, the bichir has been shown to be the basal ray-finned fish. Tetrapods and lungfish have formed a monophyletic cluster in the tree inferred from the concatenated alignment, being consistent with the currently prevalent view. It also remains possible that tetrapods are more closely related to ray-finned fishes than to lungfishes.

Selection and gene duplication: a view from the genome

Immediately after a gene duplication event, the duplicate genes have redundant functions. Is natural selection therefore completely relaxed after duplication? Does one gene evolve more rapidly than the other? Several recent genome-wide studies have suggested that duplicate genes are always under purifying selection and do not always evolve at the same rate.