Isoprenoid biosynthesis: The evolution of two ancient and distinct pathways across genomes

Isopentenyl diphosphate (IPP) is the central intermediate in the biosynthesis of isoprenoids, the most ancient and diverse class of natural products. Two distinct routes of IPP biosynthesis occur in nature: the mevalonate pathway and the recently discovered deoxyxylulose 5-phosphate (DXP) pathway. The evolutionary history of the enzymes involved in both routes and the phylogenetic distribution of their genes across genomes suggest that the mevalonate pathway is germane to archaebacteria, that the DXP pathway is germane to eubacteria, and that eukaryotes have inherited their genes for IPP biosynthesis from prokaryotes. The occurrence of genes specific to the DXP pathway is restricted to plastid-bearing eukaryotes, indicating that these genes were acquired from the cyanobacterial ancestor of plastids. However, the individual phylogenies of these genes, with only one exception, do not provide evidence for a specific affinity between the plant genes and their cyanobacterial homologues. The results suggest that lateral gene transfer between eubacteria subsequent to the origin of plastids has played a major role in the evolution of this pathway.

On roads not taken in the evolution of protein catalysts: Antibody steroid isomerases that use an enamine mechanism

Reactive immunization has emerged as a new tool for the study of biological catalysis. A powerful application resulted in catalytic antibodies that use an enamine mechanism akin to that used by the class I aldolases. With regard to the evolution of enzyme mechanisms, we investigated the utility of an enamine pathway for the allylic rearrangement exemplified by Δ5-3-ketosteroid isomerase (KSI; EC 5.3.3.1). Our aldolase antibodies were found to catalyze the isomerization of both steroid model compounds and steroids. The kinetic and chemical studies showed that the antibodies afforded rate accelerations up to a factor of 104 by means of an enamine mechanism in which imine formation was the rate-determining step. In light of our observations and the enzyme studies by other workers, we suggest that an enamine pathway could have been an early, viable KSI mechanism. Although this pathway is amenable to optimization for increased catalytic power, it appears that certain factors precluded its evolution in known KSI enzymes.

Ecological factors rather than temporal factors dominate the evolution of vesicular stomatitis virus

Vesicular stomatitis New Jersey virus (VSV-NJ) is a rhabdovirus that causes economically important disease in cattle and other domestic animals in endemic areas from southeastern United States to northern South America. Its negatively stranded RNA genome is capable of undergoing rapid evolution, which allows phylogenetic analysis and molecular epidemiology studies to be performed. Previous epidemiological studies in Costa Rica showed the existence of at least two distinct ecological zones of high VSV-NJ activity, one located in the highlands (premontane tropical moist forest) and the other in the lowlands (tropical dry forest). We wanted to test the hypothesis that the viruses circulating in these ecological zones were genetically distinct. For this purpose, we sequenced the hypervariable region of the phosphoprotein gene for 50 VSV-NJ isolates from these areas. Phylogenetic analysis showed that viruses from each ecological zone had distinct genotypes. These genotypes were maintained in each area for periods of up to 8 years. This evolutionary pattern of VSV-NJ suggests an adaptation to ecological factors that could exert selective pressure on the virus. As previous data indicated an absence of virus adaptation to factors related to the bovine host (including immunological pressure), it appears that VSV genetic divergence represents positive selection to adapt to specific vectors and/or reservoirs at each ecological zone.

The carbamate kinase-like carbamoyl phosphate synthetase of the hyperthermophilic archaeon Pyrococcus furiosus, a missing link in the evolution of carbamoyl phosphate biosynthesis

Microbial carbamoyl phosphate synthetases (CPS) use glutamine as nitrogen donor and are composed of two subunits (or domains), one exhibiting glutaminase activity, the other able to synthesize carbamoyl phosphate (CP) from bicarbonate, ATP, and ammonia. The pseudodimeric organization of this synthetase suggested that it has evolved by duplication of a smaller kinase, possibly a carbamate kinase (CK). In contrast to other prokaryotes the hyperthermophilic archaeon Pyrococcus furiosus was found to synthesize CP by using ammonia and not glutamine. We have purified the cognate enzyme and found it to be a dimer of two identical subunits of Mr 32,000. Its thermostability is considerable, 50% activity being retained after 1 h at 100°C or 3 h at 95°C. The corresponding gene was cloned by PCR and found to present about 50% amino acid identity with known CKs. The stoichiometry of the reaction (two ATP consumed per CP synthesized) and the ability of the enzyme to catalyze at high rate a bicarbonate-dependent ATPase reaction however clearly distinguish P. furiosus CPS from ordinary CKs. Thus the CPS of P. furiosus could represent a primeval step in the evolution of CPS from CK. Our results suggest that the first event in this evolution was the emergence of a primeval synthetase composed of subunits able to synthesize both carboxyphosphate and CP; this step would have preceded the duplication assumed to have generated the two subdomains of modern CPSs. The gene coding for this CK-like CPS was called cpkA.

Population dynamics, demographic stochasticity, and the evolution of cooperation

A basic evolutionary problem posed by the Iterated Prisoner’s Dilemma game is to understand when the paradigmatic cooperative strategy Tit-for-Tat can invade a population of pure defectors. Deterministically, this is impossible. We consider the role of demographic stochasticity by embedding the Iterated Prisoner’s Dilemma into a population dynamic framework. Tit-for-Tat can invade a population of defectors when their dynamics exhibit short episodes of high population densities with subsequent crashes and long low density periods with strong genetic drift. Such dynamics tend to have reddened power spectra and temporal distributions of population size that are asymmetric and skewed toward low densities. The results indicate that ecological dynamics are important for evolutionary shifts between adaptive peaks.

Aminoacylation of tRNA in the evolution of an aminoacyl-tRNA synthetase

Aminoacyl-tRNA synthetases catalyze aminoacylation of tRNAs by joining an amino acid to its cognate tRNA. The selection of the cognate tRNA is jointly determined by separate structural domains that examine different regions of the tRNA. The cysteine-tRNA synthetase of Escherichia coli has domains that select for tRNAs containing U73, the GCA anticodon, and a specific tertiary structure at the corner of the tRNA L shape. The E. coli enzyme does not efficiently recognize the yeast or human tRNACys, indicating the evolution of determinants for tRNA aminoacylation from E. coli to yeast to human and the coevolution of synthetase domains that interact with these determinants. By successively modifying the yeast and human tRNACys to ones that are efficiently aminoacylated by the E. coli enzyme, we have identified determinants of the tRNA that are important for aminoacylation but that have diverged in the course of evolution. These determinants provide clues to the divergence of synthetase domains. We propose that the domain for selecting U73 is conserved in evolution. In contrast, we propose that the domain for selecting the corner of the tRNA L shape diverged early, after the separation between E. coli and yeast, while that for selecting the GCA-containing anticodon loop diverged late, after the separation between yeast and human.

The evolution of begging: Signaling and sibling competition

In many species, young solicit food from their parents, which respond by feeding them. Because of the difference in genetic make-up between parents and their offspring and the consequent conflict, this interaction is often studied as a paradigm for the evolution of communication. Existent theoretical models demonstrate that chick signaling and parent responding can be stable if solicitation is a costly signal. The marginal cost of producing stronger signals allows the system to converge to an equilibrium: young beg with intensity that reflects their need, and parents use this information to maximize their own inclusive fitness. However, we show that there is another equilibrium where chicks do not beg and parents’ provisioning effort is optimal with respect to the statistically probable distribution of chicks’ states. Expected fitness for parents and offspring at the nonsignaling equilibrium is higher than at the signaling equilibrium. Because nonsignaling is stable and it is likely to be the ancestral condition, we would like to know how natural systems evolved from nonsignaling to signaling. We suggest that begging may have evolved through direct sibling fighting before the establishment of a parental response, that is, that nonsignaling squabbling leads to signaling. In multiple-offspring broods, young following a condition-dependent strategy in the contest for resources provide information about their condition. Parents can use this information even though it is not an adaptation for communication, and evolution will lead the system to the signaling equilibrium. This interpretation implies that signaling evolved in multiple-offspring broods, but given that signaling is evolutionarily stable, it would also be favored in species which secondarily evolved single-chick broods.

The human sex-reversing ATRX gene has a homologue on the marsupial Y chromosome, ATRY: Implications for the evolution of mammalian sex determination

Mutations in the ATRX gene on the human X chromosome cause X-linked α-thalassemia and mental retardation. XY patients with deletions or mutations in this gene display varying degrees of sex reversal, implicating ATRX in the development of the human testis. To explore further the role of ATRX in mammalian sex differentiation, the homologous gene was cloned and characterized in a marsupial. Surprisingly, active homologues of ATRX were detected on the marsupial Y as well as the X chromosome. The Y-borne copy (ATRY) displays testis-specific expression. This, as well as the sex reversal of ATRX patients, suggests that ATRY is involved in testis development in marsupials and may represent an ancestral testis-determining mechanism that predated the evolution of SRY as the primary mammalian male sex-determining gene. There is no evidence for a Y-borne ATRX homologue in mouse or human, implying that this gene has been lost in eutherians and its role supplanted by the evolution of SRY from SOX3 as the dominant determiner of male differentiation.

Positive Darwinian selection drives the evolution of several female reproductive proteins in mammals

Rapid evolution driven by positive Darwinian selection is a recurrent theme in male reproductive protein evolution. In contrast, positive selection has never been demonstrated for female reproductive proteins. Here, we perform phylogeny-based tests on three female mammalian fertilization proteins and demonstrate positive selection promoting their divergence. Two of these female fertilization proteins, the zona pellucida glycoproteins ZP2 and ZP3, are part of the mammalian egg coat. Several sites identified in ZP3 as likely to be under positive selection are located in a region previously demonstrated to be involved in species-specific sperm-egg interaction, suggesting the selective pressure is related to male-female interaction. The results provide long-sought evidence for two evolutionary hypotheses: sperm competition and sexual conflict.

The evolution of coloration and toxicity in the poison frog family (Dendrobatidae)

The poison frogs (family Dendrobatidae) are terrestrial anuran amphibians displaying a wide range of coloration and toxicity. These frogs generally have been considered to be aposematic, but relatively little research has been carried out to test the predictions of this hypothesis. Here we use a comparative approach to test one prediction of the hypothesis of aposematism: that coloration will evolve in tandem with toxicity. Recently, we developed a phylogenetic hypothesis of the evolutionary relationships among representative species of poison frogs, using sequences from three regions of mitochondrial DNA. In our analysis, we use that DNA-based phylogeny and comparative analysis of independent contrasts to investigate the correlation between coloration and toxicity in the poison frog family (Dendrobatidae). Information on the toxicity of different species was obtained from the literature. Two different measures of the brightness and extent of coloration were used. (i) Twenty-four human observers were asked to rank different photos of each different species in the analysis in terms of contrast to a leaf-littered background. (ii) Color photos of each species were scanned into a computer and a computer program was used to obtain a measure of the contrast of the colors of each species relative to a leaf-littered background. Comparative analyses of the results were carried out with two different models of character evolution: gradual change, with branch lengths proportional to the amount of genetic change, and punctuational change, with all change being associated with speciation events. Comparative analysis using either method or model indicated a significant correlation between the evolution of toxicity and coloration across this family. These results are consistent with the hypothesis that coloration in this group is aposematic.

Long term trends in the evolution of H(3) HA1 human influenza type A

We have studied the HA1 domain of 254 human influenza A(H3N2) virus genes for clues that might help identify characteristics of hemagglutinins (HAs) of circulating strains that are predictive of that strain’s epidemic potential. Our preliminary findings include the following. (i) The most parsimonious tree found requires 1,260 substitutions of which 712 are silent and 548 are replacement substitutions. (ii) The HA1 portion of the HA gene is evolving at a rate of 5.7 nucleotide substitutions/year or 5.7 × 10−3 substitutions/site per year. (iii) The replacement substitutions are distributed randomly across the three positions of the codon when allowance is made for the number of ways each codon can change the encoded amino acid. (iv) The replacement substitutions are not distributed randomly over the branches of the tree, there being 2.2 times more changes per tip branch than for non-tip branches. This result is independent of how the virus was amplified (egg grown or kidney cell grown) prior to sequencing or if sequencing was carried out directly on the original clinical specimen by PCR. (v) These excess changes on the tip branches are probably the result of a bias in the choice of strains to sequence and the detection of deleterious mutations that had not yet been removed by negative selection. (vi) There are six hypervariable codons accumulating replacement substitutions at an average rate that is 7.2 times that of the other varied codons. (vii) The number of variable codons in the trunk branches (the winners of the competitive race against the immune system) is 47 ± 5, significantly fewer than in the twigs (90 ± 7), which in turn is significantly fewer variable codons than in tip branches (175 ± 8). (viii) A minimum of one of every 12 branches has nodes at opposite ends representing viruses that reside on different continents. This is, however, no more than would be expected if one were to randomly reassign the continent of origin of the isolates. (ix) Of 99 codons with at least four mutations, 31 have ratios of non-silent to silent changes with probabilities less than 0.05 of occurring by chance, and 14 of those have probabilities <0.005. These observations strongly support positive Darwinian selection. We suggest that the small number of variable positions along the successful trunk lineage, together with knowledge of the codons that have shown positive selection, may provide clues that permit an improved prediction of which strains will cause epidemics and therefore should be used for vaccine production.

The evolution of plant nuclear genes

We analyze the evolutionary dynamics of three of the best-studied plant nuclear multigene families. The data analyzed derive from the genes that encode the small subunit of ribulose-1,5-bisphosphate carboxylase (rbcS), the gene family that encodes the enzyme chalcone synthase (Chs), and the gene family that encodes alcohol dehydrogenases (Adh). In addition, we consider the limited evolutionary data available on plant transposable elements. New Chs and rbcS genes appear to be recruited at about 10 times the rate estimated for Adh genes, and this is correlated with a much smaller average gene family size for Adh genes. In addition, duplication and divergence in function appears to be relatively common for Chs genes in flowering plant evolution. Analyses of synonymous nucleotide substitution rates for Adh genes in monocots reject a linear relationship with clock time. Replacement substitution rates vary with time in a complex fashion, which suggests that adaptive evolution has played an important role in driving divergence following gene duplication events. Molecular population genetic studies of Adh and Chs genes reveal high levels of molecular diversity within species. These studies also reveal that inter- and intralocus recombination are important forces in the generation allelic novelties. Moreover, illegitimate recombination events appear to be an important factor in transposable element loss in plants. When we consider the recruitment and loss of new gene copies, the generation of allelic diversity within plant species, and ectopic exchange among transposable elements, we conclude that recombination is a pervasive force at all levels of plant evolution.