Haplotype mapping of a diploid non-meiotic organism using existing and induced aneuploidies

Haplotype maps (HapMaps) reveal underlying sequence variation and facilitate the study of recombination and genetic diversity. In general, HapMaps are produced by analysis of Single-Nucleotide Polymorphism (SNP) segregation in large numbers of meiotic progeny. Candida albicans, the most common human fungal pathogen, is an obligate diploid that does not appear to undergo meiosis. Thus, standard methods for haplotype mapping cannot be used. We exploited naturally occurring aneuploid strains to determine the haplotypes of the eight chromosome pairs in the C. albicans laboratory strain SC5314 and in a clinical isolate. Comparison of the maps revealed that the clinical strain had undergone a significant amount of genome rearrangement, consisting primarily of crossover or gene conversion recombination events. SNP map haplotyping revealed that insertion and activation of the UAU1 cassette in essential and non-essential genes can result in whole chromosome aneuploidy. UAU1 is often used to construct homozygous deletions of targeted genes in C. albicans; the exact mechanism (trisomy followed by chromosome loss versus gene conversion) has not been determined. UAU1 insertion into the essential ORC1 gene resulted in a large proportion of trisomic strains, while gene conversion events predominated when UAU1 was inserted into the non-essential LRO1 gene. Therefore, induced aneuploidies can be used to generate HapMaps, which are essential for analyzing genome alterations and mitotic recombination events in this clonal organism.     

The Parasexual Cycle in Candida albicans Provides an Alternative Pathway to Meiosis for the Formation of Recombinant Strains

Candida albicans has an elaborate, yet efficient, mating system that promotes conjugation between diploid a and α strains. The product of mating is a tetraploid a/α cell that must undergo a reductional division to return to the diploid state. Despite the presence of several “meiosis-specific” genes in the C. albicans genome, a meiotic program has not been observed. Instead, tetraploid products of mating can be induced to undergo efficient, random chromosome loss, often producing strains that are diploid, or close to diploid, in ploidy. Using SNP and comparative genome hybridization arrays we have now analyzed the genotypes of products from the C. albicans parasexual cycle. We show that the parasexual cycle generates progeny strains with shuffled combinations of the eight C. albicans chromosomes. In addition, several isolates had undergone extensive genetic recombination between homologous chromosomes, including multiple gene conversion events. Progeny strains exhibited altered colony morphologies on laboratory media, demonstrating that the parasexual cycle generates phenotypic variants of C. albicans. In several fungi, including Saccharomyces cerevisiae and Schizosaccharomyces pombe, the conserved Spo11 protein is integral to meiotic recombination, where it is required for the formation of DNA double-strand breaks. We show that deletion of SPO11 prevented genetic recombination between homologous chromosomes during the C. albicans parasexual cycle. These findings suggest that at least one meiosis-specific gene has been re-programmed to mediate genetic recombination during the alternative parasexual life cycle of C. albicans. We discuss, in light of the long association of C. albicans with warm-blooded animals, the potential advantages of a parasexual cycle over a conventional sexual cycle.     

Als3 is a Candida albicans invasin that binds to cadherins and induces endocytosis by host cells

Candida albicans is the most common cause of hematogenously disseminated and oropharyngeal candidiasis. Both of these diseases are characterized by fungal invasion of host cells. Previously, we have found that C. albicans hyphae invade endothelial cells and oral epithelial cells in vitro by inducing their own endocytosis. Therefore, we set out to identify the fungal surface protein and host cell receptors that mediate this process. We found that the C. albicans Als3 is required for the organism to be endocytosed by human umbilical vein endothelial cells and two different human oral epithelial lines. Affinity purification experiments with wild-type and an als3Δ/als3Δ mutant strain of C. albicans demonstrated that Als3 was required for C. albicans to bind to multiple host cell surface proteins, including N-cadherin on endothelial cells and E-cadherin on oral epithelial cells. Furthermore, latex beads coated with the recombinant N-terminal portion of Als3 were endocytosed by Chinese hamster ovary cells expressing human N-cadherin or E-cadherin, whereas control beads coated with bovine serum albumin were not. Molecular modeling of the interactions of the N-terminal region of Als3 with the ectodomains of N-cadherin and E-cadherin indicated that the binding parameters of Als3 to either cadherin are similar to those of cadherin–cadherin binding. Therefore, Als3 is a fungal invasin that mimics host cell cadherins and induces endocytosis by binding to N-cadherin on endothelial cells and E-cadherin on oral epithelial cells. These results uncover the first known fungal invasin and provide evidence that C. albicans Als3 is a molecular mimic of human cadherins.     

The genome sequence of Caenorhabditis briggsae: a platform for comparative genomics

The soil nematodes Caenorhabditis briggsae and Caenorhabditis elegans diverged from a common ancestor roughly 100 million years ago and yet are almost indistinguishable by eye. They have the same chromosome number and genome sizes, and they occupy the same ecological niche. To explore the basis for this striking conservation of structure and function, we have sequenced the C. briggsae genome to a high-quality draft stage and compared it to the finished C. elegans sequence. We predict approximately 19,500 protein-coding genes in the C. briggsae genome, roughly the same as in C. elegans. Of these, 12,200 have clear C. elegans orthologs, a further 6,500 have one or more clearly detectable C. elegans homologs, and approximately 800 C. briggsae genes have no detectable matches in C. elegans. Almost all of the noncoding RNAs (ncRNAs) known are shared between the two species. The two genomes exhibit extensive colinearity, and the rate of divergence appears to be higher in the chromosomal arms than in the centers. Operons, a distinctive feature of C. elegans, are highly conserved in C. briggsae, with the arrangement of genes being preserved in 96% of cases. The difference in size between the C. briggsae (estimated at approximately 104 Mbp) and C. elegans (100.3 Mbp) genomes is almost entirely due to repetitive sequence, which accounts for 22.4% of the C. briggsae genome in contrast to 16.5% of the C. elegans genome. Few, if any, repeat families are shared, suggesting that most were acquired after the two species diverged or are undergoing rapid evolution. Coclustering the C. elegans and C. briggsae proteins reveals 2,169 protein families of two or more members. Most of these are shared between the two species, but some appear to be expanding or contracting, and there seem to be as many as several hundred novel C. briggsae gene families. The C. briggsae draft sequence will greatly improve the annotation of the C. elegans genome. Based on similarity to C. briggsae, we found strong evidence for 1,300 new C. elegans genes. In addition, comparisons of the two genomes will help to understand the evolutionary forces that mold nematode genomes.     

αADα Hybrids of Cryptococcus neoformans: Evidence of Same-Sex Mating in Nature and Hybrid Fitness

Cryptococcus neoformans is a ubiquitous human fungal pathogen that causes meningoencephalitis in predominantly immunocompromised hosts. The fungus is typically haploid, and sexual reproduction involves two individuals with opposite mating types/sexes, α and a. However, the overwhelming predominance of mating type (MAT) α over a in C. neoformans populations limits α–a mating in nature. Recently it was discovered that C. neoformans can undergo same-sex mating under laboratory conditions, especially between α isolates. Whether same-sex mating occurs in nature and contributes to the current population structure was unknown. In this study, natural αADα hybrids that arose by fusion between two α cells of different serotypes (A and D) were identified and characterized, providing definitive evidence that same-sex mating occurs naturally. A novel truncated allele of the mating-type-specific cell identity determinant SXI1α was also identified as a genetic factor likely involved in this process. In addition, laboratory-constructed αADα strains exhibited hybrid vigor both in vitro and in vivo, providing a plausible explanation for their relative abundance in nature despite the fact that AD hybrids are inefficient in meiosis/sporulation and are trapped in the diploid state. These findings provide insights on the origins, genetic mechanisms, and fitness impact of unisexual hybridization in the Cryptococcus population.     

Genetic Mechanisms of Rare Matings of the Yeast SACCHAROMYCES CEREVISIAE Heterozygous for Mating Type

Yeast heterozygous for mating type lacks the ability to conjugate as judged by the mass-mating technique and accordingly is designated "non-mater". However, the non-mater shows rare mating ability with a frequency of less than 10^-6. In the present study, the RD auxotroph mating method was mainly employed with the intention of examining the rare mating ability of various non-maters, using lactate ethanol minimal medium as a selective medium for hybridization. Crosses of aα x a, aα x a, aaα x a, aαα x a, etc. resulted in the production of respective hybrids with a relatively high frequency of about 10^-6 to 10^-7, whereas crosses of aaα x a, aαα x α, aaαα x a, aaαα x α, etc. resulted in hybrids with an extremely low frequency of about less than 10^-8. Genetic analyses revealed that the rare matings were mostly caused by the presence of cells derived from the non-maters in which mating type had converted to a homozygous genotype. Mitotic recombination was shown to be a likely explanation for most of the conversion, judging from associated exchange of an outside marker, thr4. By successive employment of the RD auxotroph mating method, it was possible to produce a series of polyploid yeasts, triploids to octoploids. The DNA content and the cell volume were observed to increase parallel to the elevated ploidy states.     

Mitotic recombination accelerates adaptation in the fungus Aspergillus nidulans

Understanding the prevalence of sexual reproduction in eukaryotes is a hard problem. At least two aspects still defy a fully satisfactory explanation, the functional significance of genetic recombination and the great variation among taxa in the relative lengths of the haploid and diploid phases in the sexual cycle. We have performed an experimental study to explore the specific advantages of haploidy or diploidy in the fungus Aspergillus nidulans. Comparing the rate of adaptation to a novel environment between haploid and isogenic diploid strains over 3,000 mitotic generations, we demonstrate that diploid strains, which during the experiment have reverted to haploidy following parasexual recombination, reach the highest fitness. This is due to the accumulation of recessive deleterious mutations in diploid nuclei, some of which show their combined beneficial effect in haploid recombinants. Our findings show the adaptive significance of mitotic recombination combined with flexibility in the timing of ploidy level transition if sign epistasis is an important determinant of fitness.     

Asymmetrical reinforcement and Wolbachia infection in Drosophila

Reinforcement refers to the evolution of increased mating discrimination against heterospecific individuals in zones of geographic overlap and can be considered a final stage in the speciation process. One the factors that may affect reinforcement is the degree to which hybrid matings result in the permanent loss of genes from a species' gene pool. Matings between females of Drosophila subquinaria and males of D. recens result in high levels of offspring mortality, due to interspecific cytoplasmic incompatibility caused by Wolbachia infection of D. recens. Such hybrid inviability is not manifested in matings between D. recens females and D. subquinaria males. Here we ask whether the asymmetrical hybrid inviability is associated with a corresponding asymmetry in the level of reinforcement. The geographic ranges of D. recens and D. subquinaria were found to overlap across a broad belt of boreal forest in central Canada. Females of D. subquinaria from the zone of sympatry exhibit much stronger levels of discrimination against males of D. recens than do females from allopatric populations. In contrast, such reproductive character displacement is not evident in D. recens, consistent with the expected effects of unidirectional cytoplasmic incompatibility. Furthermore, there is substantial behavioral isolation within D. subquinaria, because females from populations sympatric with D. recens discriminate against allopatric conspecific males, whereas females from populations allopatric with D. recens show no discrimination against any conspecific males. Patterns of general genetic differentiation among populations are not consistent with patterns of behavioral discrimination, which suggests that the behavioral isolation within D. subquinaria results from selection against mating with Wolbachia-infected D. recens. Interspecific cytoplasmic incompatibility may contribute not only to post-mating isolation, an effect already widely recognized, but also to reinforcement, particularly in the uninfected species. The resulting reproductive character displacement not only increases behavioral isolation from the Wolbachia-infected species, but may also lead to behavioral isolation between populations of the uninfected species. Given the widespread occurrence of Wolbachia among insects, it thus appears that there are multiple ways by which these endosymbionts may directly and indirectly contribute to reproductive isolation and speciation.     

A screen for recessive speciation genes expressed in the gametes of F1 hybrid yeast

Diploid hybrids of Saccharomyces cerevisiae and its closest relative, Saccharomyces paradoxus, are viable, but the sexual gametes they produce are not. One of several possible causes of this gamete inviability is incompatibility between genes from different species—such incompatible genes are usually called “speciation genes.” In diploid F1 hybrids, which contain a complete haploid genome from each species, the presence of compatible alleles can mask the effects of (recessive) incompatible speciation genes. But in the haploid gametes produced by F1 hybrids, recessive speciation genes may be exposed, killing the gametes and thus preventing F1 hybrids from reproducing sexually. Here I present the results of an experiment to detect incompatibilities that kill hybrid gametes. I transferred nine of the 16 S. paradoxus chromosomes individually into S. cerevisiae gametes and tested the ability of each to replace its S. cerevisiae homeolog. All nine chromosomes were compatible, producing nine viable haploid strains, each with 15 S. cerevisiae chromosomes and one S. paradoxus chromosome. Thus, none of these chromosomes contain speciation genes that were capable of killing the hybrid gametes that received them. This is a surprising result that suggests that such speciation genes do not play a major role in yeast speciation.     

Rare-Cell Fusion Events between Diploid and Haploid Strains of the Sexually Agglutinative Yeast HANSENULA WINGEI

Diploids of the yeast Hansenula wingei are nonagglutinative and do not form zygotes in mixed cultures with either sexually agglutinative haploid mating type. However, a low frequency of diploid x haploid cell fusions (about 10^-3) is detectable by prototrophic selection. This frequency of rare diploid x haploid matings is not increased after the diploid culture is induced for sexual agglutination. Therefore, we conclude that genes that repress mating are different from those that repress sexual agglutination.——Six prototrophs isolated from one diploid x haploid cross had an average DNA value (µg DNA per 10⁸ cells) of 6.19, compared to 2.53 and 4.35 for the haploid and diploid strains, respectively. Four prototrophs were clearly cell-fusion products because they contained genes from both the diploid and the haploid partners. However, genetic analysis of the prototrophs yielded results inconsistent with triploid meiosis; all six isolates yielded a 2:2 segregation for the mating-type alleles and linked genes.——Mitotic segregation of monosomic (2n-1) cells lacking one homolog of the chromosome carrying the mating-type locus is proposed to explain the rare production of sexually active cells in the diploid cultures. Fusion between such monosomic cells and normal haploids is thought to have produced 3n-1 cells, disomic for the chromosome carrying the mating-type locus. We conclude that in the diploid strain we studied, the physiological mechanisms repressing sexual agglutination and conjugation function efficiently, but events occuring during mitosis lead to a low frequency of genetically altered cells in the population.     

Widespread discordance of gene trees with species tree in Drosophila: evidence for incomplete lineage sorting

The phylogenetic relationship of the now fully sequenced species Drosophila erecta and D. yakuba with respect to the D. melanogaster species complex has been a subject of controversy. All three possible groupings of the species have been reported in the past, though recent multi-gene studies suggest that D. erecta and D. yakuba are sister species. Using the whole genomes of each of these species as well as the four other fully sequenced species in the subgenus Sophophora, we set out to investigate the placement of D. erecta and D. yakuba in the D. melanogaster species group and to understand the cause of the past incongruence. Though we find that the phylogeny grouping D. erecta and D. yakuba together is the best supported, we also find widespread incongruence in nucleotide and amino acid substitutions, insertions and deletions, and gene trees. The time inferred to span the two key speciation events is short enough that under the coalescent model, the incongruence could be the result of incomplete lineage sorting. Consistent with the lineage-sorting hypothesis, substitutions supporting the same tree were spatially clustered. Support for the different trees was found to be linked to recombination such that adjacent genes support the same tree most often in regions of low recombination and substitutions supporting the same tree are most enriched roughly on the same scale as linkage disequilibrium, also consistent with lineage sorting. The incongruence was found to be statistically significant and robust to model and species choice. No systematic biases were found. We conclude that phylogenetic incongruence in the D. melanogaster species complex is the result, at least in part, of incomplete lineage sorting. Incomplete lineage sorting will likely cause phylogenetic incongruence in many comparative genomics datasets. Methods to infer the correct species tree, the history of every base in the genome, and comparative methods that control for and/or utilize this information will be valuable advancements for the field of comparative genomics.     

A quantitative study of the division cycle of Caulobacter crescentus stalked cells

Progression of a cell through the division cycle is tightly controlled at different steps to ensure the integrity of genome replication and partitioning to daughter cells. From published experimental evidence, we propose a molecular mechanism for control of the cell division cycle in Caulobacter crescentus. The mechanism, which is based on the synthesis and degradation of three “master regulator” proteins (CtrA, GcrA, and DnaA), is converted into a quantitative model, in order to study the temporal dynamics of these and other cell cycle proteins. The model accounts for important details of the physiology, biochemistry, and genetics of cell cycle control in stalked C. crescentus cell. It reproduces protein time courses in wild-type cells, mimics correctly the phenotypes of many mutant strains, and predicts the phenotypes of currently uncharacterized mutants. Since many of the proteins involved in regulating the cell cycle of C. crescentus are conserved among many genera of α-proteobacteria, the proposed mechanism may be applicable to other species of importance in agriculture and medicine.     

Major structural differences and novel potential virulence mechanisms from the genomes of multiple campylobacter species

Sequencing and comparative genome analysis of four strains of Campylobacter including C. lari RM2100, C. upsaliensis RM3195, and C. coli RM2228 has revealed major structural differences that are associated with the insertion of phage- and plasmid-like genomic islands, as well as major variations in the lipooligosaccharide complex. Poly G tracts are longer, are greater in number, and show greater variability in C. upsaliensis than in the other species. Many genes involved in host colonization, including racR/S, cadF, cdt, ciaB, and flagellin genes, are conserved across the species, but variations that appear to be species specific are evident for a lipooligosaccharide locus, a capsular (extracellular) polysaccharide locus, and a novel Campylobacter putative licABCD virulence locus. The strains also vary in their metabolic profiles, as well as their resistance profiles to a range of antibiotics. It is evident that the newly identified hypothetical and conserved hypothetical proteins, as well as uncharacterized two-component regulatory systems and membrane proteins, may hold additional significant information on the major differences in virulence among the species, as well as the specificity of the strains for particular hosts.     

Odorant-binding proteins OBP57d and OBP57e affect taste perception and host-plant preference in Drosophila sechellia

Despite its morphological similarity to the other species in the Drosophila melanogaster species complex, D. sechellia has evolved distinct physiological and behavioral adaptations to its host plant Morinda citrifolia, commonly known as Tahitian Noni. The odor of the ripe fruit of M. citrifolia originates from hexanoic and octanoic acid. D. sechellia is attracted to these two fatty acids, whereas the other species in the complex are repelled. Here, using interspecies hybrids between D. melanogaster deficiency mutants and D. sechellia, we showed that the Odorant-binding protein 57e (Obp57e) gene is involved in the behavioral difference between the species. D. melanogaster knock-out flies for Obp57e and Obp57d showed altered behavioral responses to hexanoic acid and octanoic acid. Furthermore, the introduction of Obp57d and Obp57e from D. simulans and D. sechellia shifted the oviposition site preference of D. melanogaster Obp57d/e^KO flies to that of the original species, confirming the contribution of these genes to D. sechellia's specialization to M. citrifolia. Our finding of the genes involved in host-plant determination may lead to further understanding of mechanisms underlying taste perception, evolution of plant–herbivore interactions, and speciation.     

Adaptive gene expression divergence inferred from population genomics

Detailed studies of individual genes have shown that gene expression divergence often results from adaptive evolution of regulatory sequence. Genome-wide analyses, however, have yet to unite patterns of gene expression with polymorphism and divergence to infer population genetic mechanisms underlying expression evolution. Here, we combined genomic expression data—analyzed in a phylogenetic context—with whole genome light-shotgun sequence data from six Drosophila simulans lines and reference sequences from D. melanogaster and D. yakuba. These data allowed us to use molecular population genetics to test for neutral versus adaptive gene expression divergence on a genomic scale. We identified recent and recurrent adaptive evolution along the D. simulans lineage by contrasting sequence polymorphism within D. simulans to divergence from D. melanogaster and D. yakuba. Genes that evolved higher levels of expression in D. simulans have experienced adaptive evolution of the associated 3′ flanking and amino acid sequence. Concomitantly, these genes are also decelerating in their rates of protein evolution, which is in agreement with the finding that highly expressed genes evolve slowly. Interestingly, adaptive evolution in 5′ cis-regulatory regions did not correspond strongly with expression evolution. Our results provide a genomic view of the intimate link between selection acting on a phenotype and associated genic evolution.     

The genetic system controlling recombination in the silkworm

The nature of recombination modifiers was investigated in Bombyx mori lines selected for high (H) and low (L) recombination rates between the p^S and Y loci in chromosome 2. Since the mean recombination rates for the H x L and L x H F₁ crosses were approximately intermediate between those of high and low lines, the cytoplasmic maternal effect and difference in the activity of recombination modifiers between marked and unmarked second chromosomes were not detected. The H x (L x H), H x (H x L), L x (L x H) and L x (H x L) backcrosses indicated the presence of additive and dominance effects of marked and unmarked second chromosomes and the remaining chromosomes.——Recombination rates between the p^S and Y loci in chromosome 2 and half-nonrecombination rates between the pe and re loci in chromosome 5 of high and low lines indicated that these recombination modifiers caused changes in the recombination frequency between p^S and Y in chromosome 2, but not between pe and re in chromosome 5.——There were no differences in viability between individuals having the second chromosomes of the recombinant types [p^S +, p^Y (H); p^S +, + Y (L)] and those of the nonrecombinant types [p^S Y, p + (H); p^S Y, + + (L)] in both high and low lines. Mean recombination rates measured in cis [p^S Y/p + (H); p^S Y/+ + (L)] and trans [p^S +/p Y (H); p^S +/+ Y (L)] males were the same in the high but not in the low line. No segregation of a single recombination modifier was indicated by the distribution of recombination rates measured in trans males [p^S +/p Y (H); p^S +/+ Y (L)] of high and low lines. Accordingly, the recombination modifiers distributed on chromosome 2 in the heterozygous condition were not gross chromosomal aberrations, but polygenic factors in the low line.     

Distribution and Frequency of Tandem Duplications of the rII Region of Bacteriophage T4d

Genetic analyses of 49 duplications of the rII region of bacteriophage T4D suggests that there is a non-random relationship between the end points of duplicated segments, that relaxed packaging restrictions have little if any effect on the distribution of duplications, that segregation is 3–4 times more frequent than normal recombination for the same interval, and that non-tandem duplications are rare. Extrapolation of the r1231 x rJ101 cross data suggests that the minimum frequency of duplications/genome is 1.7 x 10^-6, but possibly 3.4 x 10^-4.     

Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans

The population genetic perspective is that the processes shaping genomic variation can be revealed only through simultaneous investigation of sequence polymorphism and divergence within and between closely related species. Here we present a population genetic analysis of Drosophila simulans based on whole-genome shotgun sequencing of multiple inbred lines and comparison of the resulting data to genome assemblies of the closely related species, D. melanogaster and D. yakuba. We discovered previously unknown, large-scale fluctuations of polymorphism and divergence along chromosome arms, and significantly less polymorphism and faster divergence on the X chromosome. We generated a comprehensive list of functional elements in the D. simulans genome influenced by adaptive evolution. Finally, we characterized genomic patterns of base composition for coding and noncoding sequence. These results suggest several new hypotheses regarding the genetic and biological mechanisms controlling polymorphism and divergence across the Drosophila genome, and provide a rich resource for the investigation of adaptive evolution and functional variation in D. simulans.     

A Comparison of Mutation Rates for Specific Loci and Chromosome Regions in Dysgenic Hybrid Males of DROSOPHILA MELANOGASTER

The mutation rates of specific loci and chromosome regions were estimated for two types of dysgenic hybrid males. These came from crosses between P or Q males and M females in the P-M system of hybrid dysgenesis. The M x P hybrids were the more mutable for each of the loci and chromosome regions tested. The Beadex locus was highly mutable in these hybrids but did not mutate at all in the sample of gametes from the M x Q hybrids. The singed locus had 75% of the mutability of Beadex in the M x P hybrids; it was also mutable in the M x Q hybrids. The white locus was only slightly mutable in the M x P hybrids and not at all mutable in the M x Q hybrids. The mutations in singed and white probably arose from the insertion of P elements into these loci; the mutations at Beadex probably involved the action of a P element located near this locus on the X chromosome of the P strain that was used in the experiments. Mutations in two chromosome regions, one including the zeste-white loci and the other near the miniature locus, were much more frequent in the M x P hybrids than in the M x Q hybrids. These mutations also probably arose from P element insertions. The implication is that insertion mutations occur infrequently in the M x Q hybrids, possibly because most of the P elements they carry are defective. In M x P hybrids, there is variation among loci with respect to P elements mutagenesis, indicating that P elements possess a degree of insertional specificity.     

Type III effector diversification via both pathoadaptation and horizontal transfer in response to a coevolutionary arms race

The concept of the coevolutionary arms race holds a central position in our understanding of pathogen–host interactions. Here we identify the molecular mechanisms and follow the stepwise progression of an arms race in a natural system. We show how the evolution and function of the HopZ family of type III secreted effector proteins carried by the plant pathogen Pseudomonas syringae are influenced by a coevolutionary arms race between pathogen and host. We surveyed 96 isolates of P. syringae and identified three homologs (HopZ1, HopZ2, and HopZ3) distributed among ~45% of the strains. All alleles were sequenced and their expression was confirmed. Evolutionary analyses determined that the diverse HopZ1 homologs are ancestral to P. syringae, and have diverged via pathoadaptive mutational changes into three functional and two degenerate forms, while HopZ2 and HopZ3 have been brought into P. syringae via horizontal transfer from other ecologically similar bacteria. A PAML selection analysis revealed that the C terminus of HopZ1 is under strong positive selection. Despite the extensive genetic variation observed in this family, all three homologs have cysteine–protease activity, although their substrate specificity may vary. The introduction of the ancestral hopZ1 allele into strains harboring alternate alleles results in a resistance protein-mediated defense response in their respective hosts, which is not observed with the endogenous allele. These data indicate that the P. syringae HopZ family has undergone allelic diversification via both pathoadaptive mutational changes and horizontal transfer in response to selection imposed by the host defense system. This genetic diversity permits the pathogen to avoid host defenses while still maintaining a virulence-associated protease, thereby allowing it to thrive on its current host, while simultaneously impacting its host range.