Rapid evolution in plant chitinases: Molecular targets of selection in plant-pathogen coevolution

Many pathogen recognition genes, such as plant R-genes, undergo rapid adaptive evolution, providing evidence that these genes play a critical role in plant-pathogen coevolution. Surprisingly, whether rapid adaptive evolution also occurs in genes encoding other kinds of plant defense proteins is unknown. Unlike recognition proteins, plant chitinases attack pathogens directly, conferring disease resistance by degrading chitin, a component of fungal cell walls. Here, we show that nonsynonymous substitution rates in plant class I chitinase often exceed synonymous rates in the plant genus Arabis (Cruciferae) and in other dicots, indicating a succession of adaptively driven amino acid replacements. We identify individual residues that are likely subject to positive selection by using codon substitution models and determine the location of these residues on the three-dimensional structure of class I chitinase. In contrast to primate lysozymes and plant class III chitinases, structural and functional relatives of class I chitinase, the adaptive replacements of class I chitinase occur disproportionately in the active site cleft. This highly unusual pattern of replacements suggests that fungi directly defend against chitinolytic activity through enzymatic inhibition or other forms of chemical resistance and identifies target residues for manipulating chitinolytic activity. These data also provide empirical evidence that plant defense proteins not involved in pathogen recognition also evolve in a manner consistent with rapid coevolutionary interactions.

In vitro evolution of a T cell receptor with high affinity for peptide/MHC

T cell receptors (TCRs) exhibit genetic and structural diversity similar to antibodies, but they have binding affinities that are several orders of magnitude lower. It has been suggested that TCRs undergo selection in vivo to maintain lower affinities. Here, we show that there is not an inherent genetic or structural limitation on higher affinity. Higher-affinity TCR variants were generated in the absence of in vivo selective pressures by using yeast display and selection from a library of Vα CDR3 mutants. Selected mutants had greater than 100-fold higher affinity (KD ≈ 9 nM) for the peptide/MHC ligand while retaining a high degree of peptide specificity. Among the high-affinity TCR mutants, a strong preference was found for CDR3α that contained Pro or Gly residues. Finally, unlike the wild-type TCR, a soluble monomeric form of a high-affinity TCR was capable of directly detecting peptide/MHC complexes on antigen-presenting cells. These findings prove that affinity maturation of TCRs is possible and suggest a strategy for engineering TCRs that can be used in targeting specific peptide/MHC complexes for diagnostic and therapeutic purposes.

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.

Independent evolution of the prochlorophyte and green plant chlorophyll a/b light-harvesting proteins

The prochlorophytes are oxygenic prokaryotes differing from other cyanobacteria by the presence of a light-harvesting system containing both chlorophylls (Chls) a and b and by the absence of phycobilins. We demonstrate here that the Chl a/b binding proteins from all three known prochlorophyte genera are closely related to IsiA, a cyanobacterial Chl a-binding protein induced by iron starvation, and to CP43, a constitutively expressed Chl a antenna protein of photosystem II. The prochlorophyte Chl a/b protein (pcb) genes do not belong to the extended gene family encoding eukaryotic Chl a/b and Chl a/c light-harvesting proteins. Although higher plants and prochlorophytes share common pigment complements, their light-harvesting systems have evolved independently.

Paleoproterozoic snowball Earth: Extreme climatic and geochemical global change and its biological consequences

Geological, geophysical, and geochemical data support a theory that Earth experienced several intervals of intense, global glaciation (“snowball Earth” conditions) during Precambrian time. This snowball model predicts that postglacial, greenhouse-induced warming would lead to the deposition of banded iron formations and cap carbonates. Although global glaciation would have drastically curtailed biological productivity, melting of the oceanic ice would also have induced a cyanobacterial bloom, leading to an oxygen spike in the euphotic zone and to the oxidative precipitation of iron and manganese. A Paleoproterozoic snowball Earth at 2.4 Giga-annum before present (Ga) immediately precedes the Kalahari Manganese Field in southern Africa, suggesting that this rapid and massive change in global climate was responsible for its deposition. As large quantities of O2 are needed to precipitate this Mn, photosystem II and oxygen radical protection mechanisms must have evolved before 2.4 Ga. This geochemical event may have triggered a compensatory evolutionary branching in the Fe/Mn superoxide dismutase enzyme, providing a Paleoproterozoic calibration point for studies of molecular evolution.

Mitochondrial carbonic anhydrase CA VB: Differences in tissue distribution and pattern of evolution from those of CA VA suggest distinct physiological roles

A cDNA for a second mouse mitochondrial carbonic anhydrase (CA) called CA VB was identified by homology to the previously characterized murine CA V, now called CA VA. The full-length cDNA encodes a 317-aa precursor that contains a 33-aa classical mitochondrial leader sequence. Comparison of products expressed from cDNAs for murine CA VB and CA VA in COS cells revealed that both expressed active CAs that localized in mitochondria, and showed comparable activities in crude extracts and in mitochondria isolated from transfected COS cells. Northern blot analyses of total RNAs from mouse tissues and Western blot analyses of mouse tissue homogenates showed differences in tissue-specific expression between CA VB and CA VA. CA VB was readily detected in most tissues, while CA VA expression was limited to liver, skeletal muscle, and kidney. The human orthologue of murine CA VB was recently reported also. Comparison of the CA domain sequence of human CA VB with that reported here shows that the CA domains of CA VB are much more highly conserved between mouse and human (95% identity) than the CA domains of mouse and human CA VAs (78% identity). Analysis of phylogenetic relationships between these and other available human and mouse CA isozyme sequences revealed that mammalian CA VB evolved much more slowly than CA VA, accepting amino acid substitutions at least 4.5 times more slowly since each evolved from its respective human–mouse ancestral gene around 90 million years ago. Both the differences in tissue distribution and the much greater evolutionary constraints on CA VB sequences suggest that CA VB and CA VA have evolved to assume different physiological roles.

Molecular cytogenetic dissection of human chromosomes 3 and 21 evolution

Chromosome painting in placental mammalians illustrates that genome evolution is marked by chromosomal synteny conservation and that the association of chromosomes 3 and 21 may be the largest widely conserved syntenic block known for mammals. We studied intrachromosomal rearrangements of the syntenic block 3/21 by using probes derived from chromosomal subregions with a resolution of up to 10–15 Mbp. We demonstrate that the rearrangements visualized by chromosome painting, mostly translocations, are only a fraction of the actual chromosomal changes that have occurred during evolution. The ancestral segment order for both primates and carnivores is still found in some species in both orders. From the ancestral primate/carnivore condition an inversion is needed to derive the pig homolog, and a fission of chromosome 21 and a pericentric inversion is needed to derive the Bornean orangutan condition. Two overlapping inversions in the chromosome 3 homolog then would lead to the chromosome form found in humans and African apes. This reconstruction of the origin of human chromosome 3 contrasts with the generally accepted scenario derived from chromosome banding in which it was proposed that only one pericentric inversion was needed. From the ancestral form for Old World primates (now found in the Bornean orangutan) a pericentric inversion and centromere shift leads to the chromosome ancestral for all Old World monkeys. Intrachromosomal rearrangements, as shown here, make up a set of potentially plentiful and informative markers that can be used for phylogenetic reconstruction and a more refined comparative mapping of the genome.

Adaptive evolution of color vision of the Comoran coelacanth (Latimeria chalumnae)

The coelacanth, a “living fossil,” lives near the coast of the Comoros archipelago in the Indian Ocean. Living at a depth of about 200 m, the Comoran coelacanth receives only a narrow range of light, at about 480 nm. To detect the entire range of “color” at this depth, the coelacanth appears to use only two closely related paralogous RH1 and RH2 visual pigments with the optimum light sensitivities (λmax) at 478 nm and 485 nm, respectively. The λmax values are shifted about 20 nm toward blue compared with those of the corresponding orthologous pigments. Mutagenesis experiments show that each of these coadapted changes is fully explained by two amino acid replacements.

Convergent evolution of Trichomonas vaginalis lactate dehydrogenase from malate dehydrogenase

Lactate dehydrogenase (LDH) is present in the amitochondriate parasitic protist Trichomonas vaginalis and some but not all other trichomonad species. The derived amino acid sequence of T. vaginalis LDH (TvLDH) was found to be more closely related to the cytosolic malate dehydrogenase (MDH) of the same species than to any other LDH. A key difference between the two T. vaginalis sequences was that Arg91 of MDH, known to be important in coordinating the C-4 carboxyl of oxalacetate/malate, was replaced by Leu91 in LDH. The change Leu91Arg by site-directed mutagenesis converted TvLDH into an MDH. The reverse single amino acid change Arg91Leu in TvMDH, however, gave a product with no measurable LDH activity. Phylogenetic reconstructions indicate that TvLDH arose from an MDH relatively recently.

Multiple independent origins of Shigella clones of Escherichia coli and convergent evolution of many of their characteristics

The evolutionary relationships of 46 Shigella strains representing each of the serotypes belonging to the four traditional Shigella species (subgroups), Dysenteriae, Flexneri, Boydii, and Sonnei, were determined by sequencing of eight housekeeping genes in four regions of the chromosome. Analysis revealed a very similar evolutionary pattern for each region. Three clusters of strains were identified, each including strains from different subgroups. Cluster 1 contains the majority of Boydii and Dysenteriae strains (B1–4, B6, B8, B10, B14, and B18; and D3–7, D9, and D11–13) plus Flexneri 6 and 6A. Cluster 2 contains seven Boydii strains (B5, B7, B9, B11, B15, B16, and B17) and Dysenteriae 2. Cluster 3 contains one Boydii strain (B12) and the Flexneri serotypes 1–5 strains. Sonnei and three Dysenteriae strains (D1, D8, and D10) are outside of the three main clusters but, nonetheless, are clearly within Escherichia coli. Boydii 13 was found to be distantly related to E. coli. Shigella strains, like the other pathogenic forms of E. coli, do not have a single evolutionary origin, indicating convergent evolution of Shigella phenotypic properties. We estimate the three main Shigella clusters to have evolved within the last 35,000 to 270,000 years, suggesting that shigellosis was one of the early infectious diseases of humans.

Directed evolution of antibody fragments with monovalent femtomolar antigen-binding affinity

Single-chain antibody mutants have been evolved in vitro with antigen-binding equilibrium dissociation constant Kd = 48 fM and slower dissociation kinetics (half-time > 5 days) than those for the streptavidin–biotin complex. These mutants possess the highest monovalent ligand-binding affinity yet reported for an engineered protein by over two orders of magnitude. Optimal kinetic screening of randomly mutagenized libraries of 105–107 yeast surface-displayed antibodies enabled a >1,000-fold decrease in the rate of dissociation after four cycles of affinity mutagenesis and screening. The consensus mutations are generally nonconservative by comparison with naturally occurring mouse Fv sequences and with residues that do not contact the fluorescein antigen in the wild-type complex. The existence of these mutants demonstrates that the antibody Fv architecture is not intrinsically responsible for an antigen-binding affinity ceiling during in vivo affinity maturation.

Postsymbiotic plasmid acquisition and evolution of the repA1-replicon in Buchnera aphidicola

Buchnera aphidicola is an obligate, strictly vertically transmitted, bacterial symbiont of aphids. It supplies its host with essential amino acids, nutrients required by aphids but deficient in their diet of plant phloem sap. Several lineages of Buchnera show adaptation to their nutritional role in the form of plasmid-mediated amplification of key-genes involved in the biosynthesis of tryptophan (trpEG) and leucine (leuABCD). Phylogenetic analyses of these plasmid-encoded functions have thus far suggested the absence of horizontal plasmid exchange among lineages of Buchnera. Here, we describe three new Buchnera plasmids, obtained from species of the aphid host families Lachnidae and Pemphigidae. All three plasmids belong to the repA1 family of Buchnera plasmids, which is characterized by the presence of a repA1-replicon responsible for replication initiation. A comprehensive analysis of this family of plasmids unexpectedly revealed significantly incongruent phylogenies for different plasmid and chromosomally encoded loci. We infer from these incongruencies a case of horizontal plasmid transfer in Buchnera. This process may have been mediated by secondary endosymbionts, which occasionally undergo horizontal transmission in aphids.

Adaptive evolution of a candidate gene for aging in Drosophila

Examination of the phenotypic effects of specific mutations has been extensively used to identify candidate genes affecting traits of interest. However, such analyses do not reveal anything about the evolutionary forces acting at these loci, or whether standing allelic variation contributes to phenotypic variance in natural populations. The Drosophila gene methuselah (mth) has been proposed as having major effects on organismal stress response and longevity phenotype. Here, we examine patterns of polymorphism and divergence at mth in population level samples of Drosophila melanogaster, D. simulans, and D. yakuba. Mth has experienced an unusually high level of adaptive amino acid divergence concentrated in the intra- and extracellular loop domains of the receptor protein, suggesting the historical action of positive selection on those regions of the molecule that modulate signal transduction. Further analysis of single nucleotide polymorphisms (SNPs) in D. melanogaster provided evidence for contemporary and spatially variable selection at the mth locus. In ten surveyed populations, the most common mth haplotype exhibited a 40% cline in frequency that coincided with population level differences in multiple life-history traits including lifespan. This clinal pattern was not associated with any particular SNP in the coding region, indicating that selection is operating at a closely linked site that may be involved in gene expression. Together, these consistently nonneutral patterns of inter- and intraspecific variation suggest adaptive evolution of a signal transduction pathway that may modulate lifespan in nature.

Purifying selection and birth-and-death evolution in the ubiquitin gene family

Ubiquitin is a highly conserved protein that is encoded by a multigene family. It is generally believed that this gene family is subject to concerted evolution, which homogenizes the member genes of the family. However, protein homogeneity can be attained also by strong purifying selection. We therefore studied the proportion (pS) of synonymous nucleotide differences between members of the ubiquitin gene family from 28 species of fungi, plants, and animals. The results have shown that pS is generally very high and is often close to the saturation level, although the protein sequence is virtually identical for all ubiquitins from fungi, plants, and animals. A small proportion of species showed a low level of pS values, but these values appeared to be caused by recent gene duplication. It was also found that the number of repeat copies of the gene family varies considerably with species, and some species harbor pseudogenes. These observations suggest that the members of this gene family evolve almost independently by silent nucleotide substitution and are subjected to birth-and-death evolution at the DNA level.

Mutational inactivation of the proapoptotic gene BAX confers selective advantage during tumor clonal evolution

A remarkable instability at simple repeated sequences characterizes gastrointestinal cancer of the microsatellite mutator phenotype (MMP). Mutations in the DNA mismatch repair gene family underlie the MMP, a landmark for hereditary nonpolyposis colorectal cancer. These tumors define a distinctive pathway for carcinogenesis because they display a particular spectrum of mutated cancer genes containing target repeats for mismatch repair deficiency. One such gene is BAX, a proapoptotic member of the Bcl-2 family of proteins, which plays a key role in programmed cell death. More than half of colon and gastric cancers of the MMP contain BAX frameshifts in a (G)8 mononucleotide tract. However, the functional significance of these mutations in tumor progression has not been established. Here we show that inactivation of the wild-type BAX allele by de novo frameshift mutations confers a strong advantage during tumor clonal evolution. Tumor subclones with only mutant alleles frequently appeared after inoculation into nude mice of single-cell clones of colon tumor cell lines with normal alleles. In contrast, no clones of BAX-expressing cells were found after inoculation of homozygous cell clones without wild-type BAX. These results support the interpretation that BAX inactivation contributes to tumor progression by providing a survival advantage. In this context, survival analyses show that BAX mutations are indicators of poor prognosis for both colon and gastric cancer of the MMP.