Plant Pathogenomics

Plants, like animals, use several lines of defense against pathogen attack. Prominent among genes that confer disease resistance are those encoding nucleotide binding site-leucine rich repeat (NB-LRR) proteins. Likely due to selection pressures caused by pathogens, NB-LRR genes are the most variable gene family in plants, but there appear to be species-specific limits to the number of NB-LRR genes in a genome. Allelic diversity within an individual is also increased by obligatory outcrossing, which leads to genome-wide heterozygosity. In the present study, we therefore compared the NB-LRR gene complement of the selfer Arabidopsis thaliana with its outcrossing close relative A. lyrata. We then complemented and contrasted the interspecific patterns with studies of NB-LRR diversity within A. thaliana. Three important insights are (i) that both species have similar numbers of NB-LRR genes; (ii) that loci with single NB-LRR genes are less variable than tandem arrays; and (iii) that presence-absence polymorphisms within A. thaliana are not strongly correlated with the presence or absence of orthologs in A. lyrata. Although A. thaliana individuals are mostly homozygous and thus potentially less likely to suffer from aberrant interaction of NB-LRR proteins with newly introduced alleles, the number of NB-LRR genes is similar to that in A. lyrata. In intra-and interspecific comparisons, NB-LRR genes are also more variable than Receptor-Like Protein (RLPs) genes. Finally, in contrast to Drosophila, there is a clearly positive relationship between interspecific divergence and intraspecific polymorphisms.

Multilocus sequence analysis (MLSA) and type III effector (T3E) repertoire mining were performed to gain new insights into the genetic relatedness of Xanthomonas oryzae pv. oryzae (Xoo) and Xanthomonas oryzae pv. oryzicola (Xoc), two major bacterial pathogens of rice. Based on a collection of 45 African and Asian strains, we first sequenced and analysed three housekeeping genes by MLSA, Bayesian clustering and a median-joining network approach. Second, we investigated the distribution of 32 T3E genes, which are known to be major virulence factors of plant pathogenic bacteria, in all selected strains, by polymerase chain reaction and dot-blot hybridization methods. The diversity observed within housekeeping genes, as well as within T3E repertoires, clearly showed that both pathogens belong to closely related, but distinct, phylogenetic groups. Interestingly, these evolutionary groups are differentiated according to the geographical origin of the strains, suggesting that populations of Xoo and Xoc might be endemic in Africa and Asia, and thus have evolved separately. We further revealed that T3E gene repertoires of both pathogens comprise core and variable gene suites that probably have distinct roles in pathogenicity and different evolutionary histories. In this study, we carried out a functional analysis of xopO, a differential T3E gene between Xoo and Xoc, to determine the involvement of this gene in tissue specificity. Altogether, our data contribute to a better understanding of the evolutionary history of Xoo and Xoc in Africa and Asia, and provide clues for functional studies aiming to understand the virulence, host and tissue specificity of both rice pathogens.

Sclerotinia sclerotiorum and Botrytis cinerea are closely related necrotrophic plant pathogenic fungi notable for their wide host ranges and environmental persistence. These attributes have made these species models for understanding the complexity of necrotrophic, broad host-range pathogenicity. Despite their similarities, the two species differ in mating behaviour and the ability to produce asexual spores. We have sequenced the genomes of one strain of S. sclerotiorum and two strains of B. cinerea. The comparative analysis of these genomes relative to one another and to other sequenced fungal genomes is provided here. Their 38–39 Mb genomes include 11,860–14,270 predicted genes, which share 83% amino acid identity on average between the two species. We have mapped the S. sclerotiorum assembly to 16 chromosomes and found large-scale co-linearity with the B. cinerea genomes. Seven percent of the S. sclerotiorum genome comprises transposable elements compared to <1% of B. cinerea. The arsenal of genes associated with necrotrophic processes is similar between the species, including genes involved in plant cell wall degradation and oxalic acid production. Analysis of secondary metabolism gene clusters revealed an expansion in number and diversity of B. cinerea–specific secondary metabolites relative to S. sclerotiorum. The potential diversity in secondary metabolism might be involved in adaptation to specific ecological niches. Comparative genome analysis revealed the basis of differing sexual mating compatibility systems between S. sclerotiorum and B. cinerea. The organization of the mating-type loci differs, and their structures provide evidence for the evolution of heterothallism from homothallism. These data shed light on the evolutionary and mechanistic bases of the genetically complex traits of necrotrophic pathogenicity and sexual mating. This resource should facilitate the functional studies designed to better understand what makes these fungi such successful and persistent pathogens of agronomic crops.

Horizontal gene transfer (HGT) can radically alter the genomes of microorganisms, providing the capacity to adapt to new lifestyles, environments, and hosts. However, the extent of HGT between eukaryotes is unclear. Using whole-genome, gene-by-gene phylogenetic analysis we demonstrate an extensive pattern of cross-kingdom HGT between fungi and oomycetes. Comparative genomics, including the de novo genome sequence of Hyphochytrium catenoides, a free-living sister of the oomycetes, shows that these transfers largely converge within the radiation of oomycetes that colonize plant tissues. The repertoire of HGTs includes a large number of putatively secreted proteins; for example, 7.6% of the secreted proteome of the sudden oak death parasite Phytophthora ramorum has been acquired from fungi by HGT. Transfers include gene products with the capacity to break down plant cell walls and acquire sugars, nucleic acids, nitrogen, and phosphate sources from the environment. Predicted HGTs also include proteins implicated in resisting plant defense mechanisms and effector proteins for attacking plant cells. These data are consistent with the hypothesis that some oomycetes became successful plant parasites by multiple acquisitions of genes from fungi.

Emerging plant pathogens have largely been a consequence of the movement of pathogens to new geographic regions. Another documented mechanism for the emergence of plant pathogens is hybridization between individuals of different species or subspecies, which may allow rapid evolution and adaptation to new hosts or environments. Hybrid plant pathogens have traditionally been difficult to detect or confirm, but the increasing ease of cloning and sequencing PCR products now makes the identification of species that consistently have genes or alleles with phylogenetically divergent origins relatively straightforward. We investigated the genetic origin of Phytophthora andina, an increasingly common pathogen of Andean crops Solanum betaceum, S. muricatum, S. quitoense, and several wild Solanum spp. It has been hypothesized that P. andina is a hybrid between the potato late blight pathogen P. infestans and another Phytophthora species. We tested this hypothesis by cloning four nuclear loci to obtain haplotypes and using these loci to infer the phylogenetic relationships of P. andina to P. infestans and other related species. Sequencing of cloned PCR products in every case revealed two distinct haplotypes for each locus in P. andina, such that each isolate had one allele derived from a P. infestans parent and a second divergent allele derived from an unknown species that is closely related but distinct from P. infestans, P. mirabilis, and P. ipomoeae. To the best of our knowledge, the unknown parent has not yet been collected. We also observed sequence polymorphism among P. andina isolates at three of the four loci, many of which segregate between previously described P. andina clonal lineages. These results provide strong support that P. andina emerged via hybridization between P. infestans and another unknown Phytophthora species also belonging to Phytophthora clade 1c.

Recently, genome sequencing of many isolates of genetically monomorphic bacterial human pathogens has given new insights into pathogen microevolution and phylogeography. Here, we report a genome-based micro-evolutionary study of a bacterial plant pathogen, Pseudomonas syringae pv. tomato. Only 267 mutations were identified between five sequenced isolates in 3,543,009 nt of analyzed genome sequence, which suggests a recent evolutionary origin of this pathogen. Further analysis with genome-derived markers of 89 world-wide isolates showed that several genotypes exist in North America and in Europe indicating frequent pathogen movement between these world regions. Genome-derived markers and molecular analyses of key pathogen loci important for virulence and motility both suggest ongoing adaptation to the tomato host. A mutational hotspot was found in the type III-secreted effector gene hopM1. These mutations abolish the cell death triggering activity of the full-length protein indicating strong selection for loss of function of this effector, which was previously considered a virulence factor. Two non-synonymous mutations in the flagellin-encoding gene fliC allowed identifying a new microbe associated molecular pattern (MAMP) in a region distinct from the known MAMP flg22. Interestingly, the ancestral allele of this MAMP induces a stronger tomato immune response than the derived alleles. The ancestral allele has largely disappeared from today's Pto populations suggesting that flagellin-triggered immunity limits pathogen fitness even in highly virulent pathogens. An additional non-synonymous mutation was identified in flg22 in South American isolates. Therefore, MAMPs are more variable than expected differing even between otherwise almost identical isolates of the same pathogen strain.

Magnaporthe oryzae is the causal agent of rice blast disease, a devastating problem worldwide. This fungus has caused breakdown of resistance conferred by newly developed commercial cultivars. To address how the rice blast fungus adapts itself to new resistance genes so quickly, we examined chromosomal locations of AVR-Pita, a subtelomeric gene family corresponding to the Pita resistance gene, in various isolates of M. oryzae (including wheat and millet pathogens) and its related species. We found that AVR-Pita (AVR-Pita1 and AVR-Pita2) is highly variable in its genome location, occurring in chromosomes 1, 3, 4, 5, 6, 7, and supernumerary chromosomes, particularly in rice-infecting isolates. When expressed in M. oryzae, most of the AVR-Pita homologs could elicit Pita-mediated resistance, even those from non-rice isolates. AVR-Pita was flanked by a retrotransposon, which presumably contributed to its multiple translocation across the genome. On the other hand, family member AVR-Pita3, which lacks avirulence activity, was stably located on chromosome 7 in a vast majority of isolates. These results suggest that the diversification in genome location of AVR-Pita in the rice isolates is a consequence of recognition by Pita in rice. We propose a model that the multiple translocation of AVR-Pita may be associated with its frequent loss and recovery mediated by its transfer among individuals in asexual populations. This model implies that the high mobility of AVR-Pita is a key mechanism accounting for the rapid adaptation toward Pita. Dynamic adaptation of some fungal plant pathogens may be achieved by deletion and recovery of avirulence genes using a population as a unit of adaptation.

The pace of technical advancement in microbial genomics has been breathtaking. Since 1995, when the first complete genome sequence of a free-living organism, Haemophilus influenzae, was published, 1554 complete bacterial genome sequences (the majority of which are from pathogens) and 112 complete archaeal genome sequences have been determined, and more than 4800 and 90, respectively, are in progress. A total of 41 complete eukaryotic genome sequences have been determined (19 from fungi), and more than 1100 are in progress. Complete reference genome sequences are available for 2675 viral species, and for some of these species, a large number of strains have been completely sequenced. Nearly 40,000 strains of influenza virus and more than 300,000 strains of human immunodeficiency virus (HIV) type 1 have been partially sequenced. However, the selection of microbes and viruses for genome sequencing is heavily biased toward the tiny minority that are amenable to cultivation in the laboratory, numerically dominant in particular habitats of interest (e.g., the human body), and associated with disease.