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.

Stochastic processes strongly influence HIV-1 evolution during suboptimal protease-inhibitor therapy

It has long been assumed that HIV-1 evolution is best described by deterministic evolutionary models because of the large population size. Recently, however, it was suggested that the effective population size (Ne) may be rather small, thereby allowing chance to influence evolution, a situation best described by a stochastic evolutionary model. To gain experimental evidence supporting one of the evolutionary models, we investigated whether the development of resistance to the protease inhibitor ritonavir affected the evolution of the env gene. Sequential serum samples from five patients treated with ritonavir were used for analysis of the protease gene and the V3 domain of the env gene. Multiple reverse transcription–PCR products were cloned, sequenced, and used to construct phylogenetic trees and to calculate the genetic variation and Ne. Genotypic resistance to ritonavir developed in all five patients, but each patient displayed a unique combination of mutations, indicating a stochastic element in the development of ritonavir resistance. Furthermore, development of resistance induced clear bottleneck effects in the env gene. The mean intrasample genetic variation, which ranged from 1.2% to 5.7% before treatment, decreased significantly (P < 0.025) during treatment. In agreement with these findings, Ne was estimated to be very small (500–15,000) compared with the total HIV-1 RNA copy number. This study combines three independent observations, strong population bottlenecking, small Ne, and selection of different combinations of protease-resistance mutations, all of which indicate that HIV-1 evolution is best described by a stochastic evolutionary model.

Photoassembly of the manganese cluster and oxygen evolution from monomeric and dimeric CP47 reaction center photosystem II complexes

Isolated subcomplexes of photosystem II from spinach (CP47RC), composed of D1, D2, cytochrome b559, CP47, and a number of hydrophobic small subunits but devoid of CP43 and the extrinsic proteins of the oxygen-evolving complex, were shown to reconstitute the Mn4Ca1Clx cluster of the water-splitting system and to evolve oxygen. The photoactivation process in CP47RC dimers proceeds by the same two-step mechanism as observed in PSII membranes and exhibits the same stoichiometry for Mn2+, but with a 10-fold lower affinity for Ca2+ and an increased susceptibility to photodamage. After the lower Ca2+ affinity and the 10-fold smaller absorption cross-section for photons in CP47 dimers is taken into account, the intrinsic rate constant for the rate-limiting calcium-dependent dark step is indistinguishable for the two systems. The monomeric form of CP47RC also showed capacity to photoactivate and catalyze water oxidation, but with lower activity than the dimeric form and increased susceptibility to photodamage. After optimization of the various parameters affecting the photoactivation process in dimeric CP47RC subcores, 18% of the complexes were functionally reconstituted and the quantum efficiency for oxygen production by reactivated centers approached 96% of that observed for reconstituted photosystem II-enriched membranes.

Chloroplast small heat shock proteins: Evidence for atypical evolution of an organelle-localized protein

Knowledge of the origin and evolution of gene families is critical to our understanding of the evolution of protein function. To gain a detailed understanding of the evolution of the small heat shock proteins (sHSPs) in plants, we have examined the evolutionary history of the chloroplast (CP)-localized sHSPs. Previously, these nuclear-encoded CP proteins had been identified only from angiosperms. This study reveals the presence of the CP sHSPs in a moss, Funaria hygrometrica. Two clones for CP sHSPs were isolated from a F. hygrometrica heat shock cDNA library that represent two distinct CP sHSP genes. Our analysis of the CP sHSPs reveals unexpected evolutionary relationships and patterns of sequence conservation. Phylogenetic analysis of the CP sHSPs with other plant CP sHSPs and eukaryotic, archaeal, and bacterial sHSPs shows that the CP sHSPs are not closely related to the cyanobacterial sHSPs. Thus, they most likely evolved via gene duplication from a nuclear-encoded cytosolic sHSP and not via gene transfer from the CP endosymbiont. Previous sequence analysis had shown that all angiosperm CP sHSPs possess a methionine-rich region in the N-terminal domain. The primary sequence of this region is not highly conserved in the F. hygrometrica CP sHSPs. This lack of sequence conservation indicates that sometime in land plant evolution, after the divergence of mosses from the common ancestor of angiosperms but before the monocot–dicot divergence, there was a change in the selective constraints acting on the CP sHSPs.

Tracking molecular evolution of photosynthesis by characterization of a major photosynthesis gene cluster from Heliobacillus mobilis

A DNA sequence has been obtained for a 35.6-kb genomic segment from Heliobacillus mobilis that contains a major cluster of photosynthesis genes. A total of 30 ORFs were identified, 20 of which encode enzymes for bacteriochlorophyll and carotenoid biosynthesis, reaction-center (RC) apoprotein, and cytochromes for cyclic electron transport. Donor side electron-transfer components to the RC include a putative RC-associated cytochrome c553 and a unique four-large-subunit cytochrome bc complex consisting of Rieske Fe-S protein (encoded by petC), cytochrome b6 (petB), subunit IV (petD), and a diheme cytochrome c (petX). Phylogenetic analysis of various photosynthesis gene products indicates a consistent grouping of oxygenic lineages that are distinct and descendent from anoxygenic lineages. In addition, H. mobilis was placed as the closest relative to cyanobacteria, which form a monophyletic origin to chloroplast-based photosynthetic lineages. The consensus of the photosynthesis gene trees also indicates that purple bacteria are the earliest emerging photosynthetic lineage. Our analysis also indicates that an ancient gene-duplication event giving rise to the paralogous bchI and bchD genes predates the divergence of all photosynthetic groups. In addition, our analysis of gene duplication of the photosystem I and photosystem II core polypeptides supports a “heterologous fusion model” for the origin and evolution of oxygenic photosynthesis.

In vitro selection and evolution of functional proteins by using ribosome display

We report here a system with which a correctly folded complete protein and its encoding mRNA both remain attached to the ribosome and can be enriched for the ligand-binding properties of the native protein. We have selected a single-chain fragment (scFv) of an antibody 108-fold by five cycles of transcription, translation, antigen-affinity selection, and PCR. The selected scFv fragments all mutated in vitro by acquiring up to four unrelated amino acid exchanges over the five generations, but they remained fully compatible with antigen binding. Libraries of native folded proteins can now be screened and made to evolve in a cell-free system without any transformation or constraints imposed by the host cell.

Evolution of paired domains: Isolation and sequencing of jellyfish and hydra Pax genes related to Pax-5 and Pax-6

Pax proteins are a family of transcription factors with a highly conserved paired domain; many members also contain a paired-type homeodomain and/or an octapeptide. Nine mammalian Pax genes are known and classified into four subgroups: Pax-1/9, Pax-2/5/8, Pax-3/7, and Pax-4/6. Most of these genes are involved in nervous system development. In particular, Pax-6 is a key regulator that controls eye development in vertebrates and Drosophila. Although the Pax-4/6 subgroup seems to be more closely related to Pax-2/5/8 than to Pax-3/7 or Pax-1/9, its evolutionary origin is unknown. We therefore searched for a Pax-6 homolog and related genes in Cnidaria, which is the lowest phylum of animals that possess a nervous system and eyes. A sea nettle (a jellyfish) genomic library was constructed and two pax genes (Pax-A and -B) were isolated and partially sequenced. Surprisingly, unlike most known Pax genes, the paired box in these two genes contains no intron. In addition, the complete cDNA sequences of hydra Pax-A and -B were obtained. Hydra Pax-B contains both the homeodomain and the octapeptide, whereas hydra Pax-A contains neither. DNA binding assays showed that sea nettle Pax-A and -B and hydra Pax-A paired domains bound to a Pax-5/6 site and a Pax-5 site, although hydra Pax-B paired domain bound neither. An alignment of all available paired domain sequences revealed two highly conserved regions, which cover the DNA binding contact positions. Phylogenetic analysis showed that Pax-A and especially Pax-B were more closely related to Pax-2/5/8 and Pax-4/6 than to Pax-1/9 or Pax-3/7 and that the Pax genes can be classified into two supergroups: Pax-A/Pax-B/Pax-2/5/8/4/6 and Pax-1/9/3/7. From this analysis and the gene structure, we propose that modern Pax-4/6 and Pax-2/5/8 genes evolved from an ancestral gene similar to cnidarian Pax-B, having both the homeodomain and the octapeptide.

The origin and evolution of animal appendages

Animals have evolved diverse appendages adapted for locomotion, feeding and other functions. The genetics underlying appendage formation are best understood in insects and vertebrates. The expression of the Distal-less (Dll) homeoprotein during arthropod limb outgrowth and of Dll orthologs (Dlx) in fish fin and tetrapod limb buds led us to examine whether expression of this regulatory gene may be a general feature of appendage formation in protostomes and deuterostomes. We find that Dll is expressed along the proximodistal axis of developing polychaete annelid parapodia, onychophoran lobopodia, ascidian ampullae, and even echinoderm tube feet. Dll/Dlx expression in such diverse appendages in these six coelomate phyla could be convergent, but this would have required the independent co-option of Dll/Dlx several times in evolution. It appears more likely that ectodermal Dll/Dlx expression along proximodistal axes originated once in a common ancestor and has been used subsequently to pattern body wall outgrowths in a variety of organisms. We suggest that this pre-Cambrian ancestor of most protostomes and the deuterostomes possessed elements of the genetic machinery for and may have even borne appendages.

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.

Empirical laws of survival and evolution: Their universality and implications

Presented analysis of human and fly life tables proves that with the specified accuracy their entire survival and mortality curves are uniquely determined by a single point (e.g., by the birth mortality q0), according to the law, which is universal for species as remote as humans and flies. Mortality at any age decreases with the birth mortality q0. According to life tables, in the narrow vicinity of a certain q0 value (which is the same for all animals of a given species, independent of their living conditions), the curves change very rapidly and nearly simultaneously for an entire population of different ages. The change is the largest in old age. Because probability to survive to the mean reproductive age quantifies biological fitness and evolution, its universal rapid change with q0 (which changes with living conditions) manifests a new kind of an evolutionary spurt of an entire population. Agreement between theoretical and life table data is explicitly seen in the figures. Analysis of the data on basic metabolism reduces it to the maximal mean lifespan (for animals from invertebrates to mammals), or to the maximal mean fission time (for bacteria), and universally scales them with the total number of body atoms only. Phenomenological origin of this unification and universality of metabolism, survival, and evolution is suggested. Their implications and challenges are discussed.

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.

Attenuated T2 relaxation by mutual cancellation of dipole–dipole coupling and chemical shift anisotropy indicates an avenue to NMR structures of very large biological macromolecules in solution

Fast transverse relaxation of 1H, 15N, and 13C by dipole-dipole coupling (DD) and chemical shift anisotropy (CSA) modulated by rotational molecular motions has a dominant impact on the size limit for biomacromolecular structures that can be studied by NMR spectroscopy in solution. Transverse relaxation-optimized spectroscopy (TROSY) is an approach for suppression of transverse relaxation in multidimensional NMR experiments, which is based on constructive use of interference between DD coupling and CSA. For example, a TROSY-type two-dimensional 1H,15N-correlation experiment with a uniformly 15N-labeled protein in a DNA complex of molecular mass 17 kDa at a 1H frequency of 750 MHz showed that 15N relaxation during 15N chemical shift evolution and 1HN relaxation during signal acquisition both are significantly reduced by mutual compensation of the DD and CSA interactions. The reduction of the linewidths when compared with a conventional two-dimensional 1H,15N-correlation experiment was 60% and 40%, respectively, and the residual linewidths were 5 Hz for 15N and 15 Hz for 1HN at 4°C. Because the ratio of the DD and CSA relaxation rates is nearly independent of the molecular size, a similar percentagewise reduction of the overall transverse relaxation rates is expected for larger proteins. For a 15N-labeled protein of 150 kDa at 750 MHz and 20°C one predicts residual linewidths of 10 Hz for 15N and 45 Hz for 1HN, and for the corresponding uniformly 15N,2H-labeled protein the residual linewidths are predicted to be smaller than 5 Hz and 15 Hz, respectively. The TROSY principle should benefit a variety of multidimensional solution NMR experiments, especially with future use of yet somewhat higher polarizing magnetic fields than are presently available, and thus largely eliminate one of the key factors that limit work with larger molecules.

Asian genotypes of JC virus in Native Americans and in a Pacific Island population: Markers of viral evolution and human migration

The human polyomavirus JC (JCV) causes the central nervous system demyelinating disease progressive multifocal leukoencephalopathy. Previously, we showed that 40% of Caucasians in the United States excrete JCV in the urine as detected by PCR. We have now studied 68 Navaho from New Mexico, 25 Flathead from Montana, and 29 Chamorro from Guam. By using PCR amplification of a fragment of the VP1 gene, JCV DNA was detected in the urine of 45 (66%) Navaho, 14 (56%) Flathead, and 20 (69%) Chamorro. Genotyping of viral DNAs in these cohorts by cycle sequencing showed predominantly type 2 (Asian), rather than type 1 (European). Type 1 is the major type in the United States and Hungary. Type 2 can be further subdivided into 2A, 2B, and 2C. Type 2A is found in China and Japan. Type 2B is a subtype related to the East Asian type, and is now found in Europe and the United States. The large majority (56–89%) of strains excreted by Native Americans and Pacific Islanders were the type 2A subtype, consistent with the origin of these strains in Asia. These findings indicate that JCV infection of Native Americans predates contact with Europeans, and likely predates migration of Amerind ancestors across the Bering land bridge around 12,000–30,000 years ago. If JCV had already differentiated into stable modern genotypes and subtypes prior to first settlement, the origin of JCV in humans may date from 50,000 to 100,000 years ago or more. We conclude that JCV may have coevolved with the human species, and that it provides a convenient marker for human migrations in both prehistoric and modern times.

Somatic polyploidization and cellular proliferation drive body size evolution in nematodes

Most of the hypodermis of a rhabditid nematode such as Caenorhabditis elegans is a single syncytium. The size of this syncytium (as measured by body size) has evolved repeatedly in the rhabditid nematodes. Two cellular mechanisms are important in the evolution of body size: changes in the numbers of cells that fuse with the syncytium, and the extent of its acellular growth. Thus nematodes differ from mammals and other invertebrates in which body size evolution is caused by changes in cell number alone. The evolution of acellular syncytial growth in nematodes is also associated with changes in the ploidy of hypodermal nuclei. These nuclei are polyploid as a consequence of iterative rounds of endoreduplication, and this endocycle has evolved repeatedly. The association between acellular growth and endoreduplication is also seen in C. elegans mutations that interrupt transforming growth factor-β signaling and that result in dwarfism and deficiencies in hypodermal ploidy. The transforming growth factor-β pathway is a candidate for being involved in nematode body size evolution.