Megan C. McDonald

Comparative Pathogenomics Reveals Horizontally Acquired Novel Virulence Genes in Fungi Infecting Cereal Hosts

Comparative analyses of pathogen genomes provide new insights into how pathogens have evolved common and divergent virulence strategies to invade related plant species. Fusarium crown and root rots are important diseases of wheat and barley world-wide. In Australia, these diseases are primarily caused by the fungal pathogen Fusarium pseudograminearum. Comparative genomic analyses showed that the F. pseudograminearum genome encodes proteins that are present in other fungal pathogens of cereals but absent in non-cereal pathogens. In some cases, these cereal pathogen specific genes were also found in bacteria associated with plants. Phylogenetic analysis of selected F. pseudograminearum genes supported the hypothesis of horizontal gene transfer into diverse cereal pathogens. Two horizontally acquired genes with no previously known role in fungal pathogenesis were studied functionally via gene knockout methods and shown to significantly affect virulence of F. pseudograminearum on the cereal hosts wheat and barley. Our results indicate using comparative genomics to identify genes specific to pathogens of related hosts reveals novel virulence genes and illustrates the importance of horizontal gene transfer in the evolution of plant infecting fungal pathogens.

SnTox3 Acts in Effector Triggered Susceptibility to Induce Disease on Wheat Carrying the Snn3 Gene

The necrotrophic fungus Stagonospora nodorum produces multiple proteinaceous host-selective toxins (HSTs) which act in effector triggered susceptibility. Here, we report the molecular cloning and functional characterization of the SnTox3-encoding gene, designated SnTox3, as well as the initial characterization of the SnTox3 protein. SnTox3 is a 693 bp intron-free gene with little obvious homology to other known genes. The predicted immature SnTox3 protein is 25.8 kDa in size. A 20 amino acid signal sequence as well as a possible pro sequence are predicted. Six cysteine residues are predicted to form disulfide bonds and are shown to be important for SnTox3 activity. Using heterologous expression in Pichia pastoris and transformation into an avirulent S. nodorum isolate, we show that SnTox3 encodes the SnTox3 protein and that SnTox3 interacts with the wheat susceptibility gene Snn3. In addition, the avirulent S. nodorum isolate transformed with SnTox3 was virulent on host lines expressing the Snn3 gene. SnTox3-disrupted mutants were deficient in the production of SnTox3 and avirulent on the Snn3 differential wheat line BG220. An analysis of genetic diversity revealed that SnTox3 is present in 60.1% of a worldwide collection of 923 isolates and occurs as eleven nucleotide haplotypes resulting in four amino acid haplotypes. The cloning of SnTox3 provides a fundamental tool for the investigation of the S. nodorum–wheat interaction, as well as vital information for the general characterization of necrotroph–plant interactions.

The cysteine rich necrotrophic effector SnTox1 produced by Stagonospora nodorum triggers susceptibility of wheat lines harboring Snn1.

The wheat pathogen Stagonospora nodorum produces multiple necrotrophic effectors (also called host-selective toxins) that promote disease by interacting with corresponding host sensitivity gene products. SnTox1 was the first necrotrophic effector identified in S. nodorum, and was shown to induce necrosis on wheat lines carrying Snn1. Here, we report the molecular cloning and validation of SnTox1 as well as the preliminary characterization of the mechanism underlying the SnTox1-Snn1 interaction which leads to susceptibility. SnTox1 was identified using bioinformatics tools and verified by heterologous expression in Pichia pastoris. SnTox1 encodes a 117 amino acid protein with the first 17 amino acids predicted as a signal peptide, and strikingly, the mature protein contains 16 cysteine residues, a common feature for some avirulence effectors. The transformation of SnTox1 into an avirulent S. nodorum isolate was sufficient to make the strain pathogenic. Additionally, the deletion of SnTox1 in virulent isolates rendered the SnTox1 mutated strains avirulent on the Snn1 differential wheat line. SnTox1 was present in 85% of a global collection of S. nodorum isolates. We identified a total of 11 protein isoforms and found evidence for strong diversifying selection operating on SnTox1. The SnTox1-Snn1 interaction results in an oxidative burst, DNA laddering, and pathogenesis related (PR) gene expression, all hallmarks of a defense response. In the absence of light, the development of SnTox1-induced necrosis and disease symptoms were completely blocked. By comparing the infection processes of a GFP-tagged avirulent isolate and the same isolate transformed with SnTox1, we conclude that SnTox1 may play a critical role during fungal penetration. This research further demonstrates that necrotrophic fungal pathogens utilize small effector proteins to exploit plant resistance pathways for their colonization, which provides important insights into the molecular basis of the wheat-S. nodorum interaction, an emerging model for necrotrophic pathosystems.

Global diversity and distribution of three necrotrophic effectors in Phaeosphaeria nodorum and related species

Population genetic and phylogenetic studies have shown that Phaeosphaeria nodorum is a member of a species complex that probably shares its center of origin with wheat (Triticum aestivum and Triticum durum). We examined the evolutionary histories of three known necrotrophic effectors (NEs) produced by P. nodorum and compared them with neutral loci. We screened over 1000 individuals for the presence/absence of each effector and assigned each individual to a multi-effector genotype. Diversity at each NE locus was assessed by sequencing c. 200 individuals for each locus. We found significant differences in effector frequency among populations. We propose that these differences reflect the presence/absence of the corresponding susceptibility gene in wheat cultivars. The population harboring the highest sequence diversity was different for each effector locus and never coincided with populations harboring the highest diversity at neutral loci. Coalescent and phylogenetic analyses showed a discontinuous presence of all three NEs among nine closely related Phaeosphaeria species. Only two of the nine species were found to harbor NEs. We present evidence that the three described NEs of P. nodorum were transmitted to its sister species, Phaeosphaeria avenaria tritici 1, via interspecific hybridization.

Is Zymoseptoria tritici a hemibiotroph

The growth of microorganisms in planta is often categorized based on their methods of nutrient acquisition and the physical appearance of symptoms on the host. For example, biotrophs thrive on living tissue while necrotrophic pathogens often quickly lyse cells to access nutrients. Hemibiotrophs are pathogens that initially feed on living host tissue without causing visible symptoms prior to switching to necrotrophy. During infection of wheat, the pathogen Zymoseptoria tritici undergoes a prolonged and asymptomatic phase during which it grows slowly and protects itself from host defenses prior to eliciting a strong necrotic response. However careful analyses of the asymptomatic phase indicate that the pathogen does not alter host growth, casting doubt on the biotrophic nature of this asymptomatic period. Consequently, we question whether Z. tritici is correctly defined as a hemibiotroph.

Fungal phytopathogens encode functional homologues of plant rapid alkalinization factor (RALF) peptides.

In this article, we describe the presence of genes encoding close homologues of an endogenous plant peptide, rapid alkalinization factor (RALF), within the genomes of 26 species of phytopathogenic fungi. Members of the RALF family are key growth factors in plants, and the sequence of the RALF active region is well conserved between plant and fungal proteins. RALF1-like sequences were observed in most cases; however, RALF27-like sequences were present in the Sphaerulina musiva and Septoria populicola genomes. These two species are pathogens of poplar and, interestingly, the closest relative to their respective RALF genes is a poplar RALF27-like sequence. RALF peptides control cellular expansion during plant development, but were originally defined on the basis of their ability to induce rapid alkalinization in tobacco cell cultures. To test whether the fungal RALF peptides were biologically active in plants, we synthesized RALF peptides corresponding to those encoded by two sequenced genomes of the tomato pathogen Fusarium oxysporum f. sp. lycopersici. One of these peptides inhibited the growth of tomato seedlings and elicited responses in tomato and Nicotiana benthamiana typical of endogenous plant RALF peptides (reactive oxygen species burst, induced alkalinization and mitogen-activated protein kinase activation). Gene expression analysis confirmed that a RALF-encoding gene in F. oxysporum f. sp. lycopersici was expressed during infection on tomato. However, a subsequent reverse genetics approach revealed that the RALF peptide was not required by F. oxysporum f. sp. lycopersici for infection on tomato roots. This study has demonstrated the presence of functionally active RALF peptides encoded within phytopathogens that harbour an as yet undetermined role in plant-pathogen interactions.

Phytopathogen emergence in the genomics era.

Phytopathogens are a global threat to plant agriculture and biodiversity. The genomics era has lead to an exponential rise in comparative gene and genome studies of both economically significant and insignificant microorganisms. In this review we highlight some recent comparisons and discuss how they identify shared genes or genomic regions associated with host virulence. The two major mechanisms of rapid genome adaptation - horizontal gene transfer and hybridisation - are reviewed and we consider how intra-specific pan-genome sequences encode alternative host specificity. We also discuss the power that access to expansive gene databases provides in aiding the study of phytopathogen emergence. These databases can rapidly enable the identification of an unknown pathogen and its origin, as well as genomic adaptations required for emergence.

Phylogenetic and population genetic analyses of Phaeosphaeria nodorum and its close relatives indicate cryptic species and an origin in the Fertile Crescent

The origin of the fungal wheat pathogen Phaeosphaeria nodorum remains unclear despite earlier intensive global population genetic and phylogeographical studies. We sequenced 1683 bp distributed across three loci in 355 globally distributed Phaeosphaeria isolates, including 74 collected in Iran near the center of origin of wheat. We identified nine phylogenetically distinct clades, including two previously unknown species tentatively named P1 and P2 collected in Iran. Coalescent analysis indicates that P1 and P2 are sister species of P. nodorum and the other Phaeosphaeria species identified in our analysis. Two species, P. nodorum and P. avenaria f. sp. tritici 1 (Pat1), comprised ~85% of the sampled isolates, making them the dominant wheat-infecting pathogens within the species complex. We designed a PCR-RFLP assay to distinguish P. nodorum from Pat1. Approximately 4% of P. nodorum and Pat1 isolates showed evidence of hybridization. Measures of private allelic richness at SSR and sequence loci suggest that the center of origin of P. nodorum coincides with its host in the Fertile Crescent. We hypothesize that the origin of this species complex is also in the Fertile Crescent, with four species out of nine found exclusively in the Iranian collections.

The discovery of the virulence gene ToxA in the wheat and barley pathogen Bipolaris sorokiniana

Bipolaris sorokiniana is the causal agent of multiple diseases on wheat and barley and is the primary constraint to cereal production throughout South Asia. Despite its significance, the molecular basis of disease is poorly understood. To address this, the genomes of three Australian isolates of B. sorokiniana were sequenced and screened for known pathogenicity genes. Sequence analysis revealed that the isolate BRIP10943 harboured the ToxA gene, which has been associated previously with disease in the wheat pathogens Parastagonospora nodorum and Pyrenophora tritici-repentis. Analysis of the regions flanking ToxA within B. sorokiniana revealed that it was embedded within a 12-kb genomic element nearly identical to the corresponding regions in P. nodorum and P. tritici-repentis. A screen of 35 Australian B. sorokiniana isolates confirmed that ToxA was present in 12 isolates. Sequencing of the ToxA genes within these isolates revealed two haplotypes, which differed by a single non-synonymous nucleotide substitution. Pathogenicity assays showed that a B. sorokiniana isolate harbouring ToxA was more virulent on wheat lines that contained the sensitivity gene when compared with a non-ToxA isolate. This work demonstrates that proteins that confer host-specific virulence can be horizontally acquired across multiple species. This acquisition can dramatically increase the virulence of pathogenic strains on susceptible cultivars, which, in an agricultural setting, can have devastating economic and social impacts.

Utilizing Gene Tree Variation to Identify Candidate Effector Genes in Zymoseptoria tritici

Zymoseptoria tritici is a host-specific, necrotrophic pathogen of wheat. Infection by Z. tritici is characterized by its extended latent period, which typically lasts 2 wks, and is followed by extensive host cell death, and rapid proliferation of fungal biomass. This work characterizes the level of genomic variation in 13 isolates, for which we have measured virulence on 11 wheat cultivars with differential resistance genes. Between the reference isolate, IPO323, and the 13 Australian isolates we identified over 800,000 single nucleotide polymorphisms, of which ∼10% had an effect on the coding regions of the genome. Furthermore, we identified over 1700 probable presence/absence polymorphisms in genes across the Australian isolates using de novo assembly. Finally, we developed a gene tree sorting method that quickly identifies groups of isolates within a single gene alignment whose sequence haplotypes correspond with virulence scores on a single wheat cultivar. Using this method, we have identified < 100 candidate effector genes whose gene sequence correlates with virulence toward a wheat cultivar carrying a major resistance gene.