Evolution of Linked Avirulence Effectors in Leptosphaeria maculans Is Affected by Genomic Environment and Exposure to Resistance Genes in Host Plants

Brassica napus (canola) cultivars and isolates of the blackleg fungus, Leptosphaeria maculans interact in a ‘gene for gene’ manner whereby plant resistance (R) genes are complementary to pathogen avirulence (Avr) genes. Avirulence genes encode proteins that belong to a class of pathogen molecules known as effectors, which includes small secreted proteins that play a role in disease. In Australia in 2003 canola cultivars with the Rlm1 resistance gene suffered a breakdown of disease resistance, resulting in severe yield losses. This was associated with a large increase in the frequency of virulence alleles of the complementary avirulence gene, AvrLm1, in fungal populations. Surprisingly, the frequency of virulence alleles of AvrLm6 (complementary to Rlm6) also increased dramatically, even though the cultivars did not contain Rlm6. In the L. maculans genome, AvrLm1 and AvrLm6 are linked along with five other genes in a region interspersed with transposable elements that have been degenerated by Repeat-Induced Point (RIP) mutations. Analyses of 295 Australian isolates showed deletions, RIP mutations and/or non-RIP derived amino acid substitutions in the predicted proteins encoded by these seven genes. The degree of RIP mutations within single copy sequences in this region was proportional to their proximity to the degenerated transposable elements. The RIP alleles were monophyletic and were present only in isolates collected after resistance conferred by Rlm1 broke down, whereas deletion alleles belonged to several polyphyletic lineages and were present before and after the resistance breakdown. Thus, genomic environment and exposure to resistance genes in B. napus has affected the evolution of these linked avirulence genes in L. maculans

Phoma stem canker disease on oilseed rape (Brassica napus) in China is caused by Leptosphaeria biglobosa ‘brassicae’

This document is the Accepted Manuscript version of the following article: Ze Liu, Akinwunmi O. Latunde-Dada, Avice M. Hall, Bruce D. L. Fitt, ‘Phoma stem canker disease on oilseed rape (Brassica napus) in China is caused by Leptosphaeria biglobosa ‘brassicae’’, European Journal of Plant Pathology, Vol. 140(4): 841-857, December 2014. The final publication is available at Springer via: http://dx.doi.org/10.1007/s10658-014-0513-7 © Koninklijke Nederlandse Planteziektenkundige Vereniging 2014Phoma stem canker of oilseed rape (Brassica napus) is a globally important disease that is caused by the sibling ascomycete species Leptosphaeria maculans and L. biglobosa. Sixty fungal isolates obtained from oilseed rape stems with phoma stem canker disease symptoms collected from four provinces in China in 1999, 2005 and 2006 were all identified as Leptosphaeria biglobosa, not L. maculans, by PCR diagnostics based on species-specific primers. There were no differences in cultural characteristics (e.g. pigmentation and in vitro growth) between these L. biglobosa isolates from China and those of 37 proven L. biglobosa isolates from Europe or Canada. In studies using amplified fragment length polymorphism (AFLP) markers, Chinese L. biglobosa populations were genetically more similar to European L. biglobosa populations than to the more diverse Canadian L. biglobosa populations. Sequencing of gene fragments of β-tubulin, actin and the internal transcribed spacer (ITS) region of rDNA from L. biglobosa isolates from China, Europe, Australia and Canada showed a closer taxonomic similarity of Chinese L. biglobosa to the European L. biglobosa ‘brassicae’ than to Canadian L. biglobosa ‘canadensis’ or to the Australian L. biglobosa ‘occiaustralensis’ or ‘australensis’ subclades. These results suggest that the Chinese L. biglobosa population in this study is in the same subclade as European L. biglobosa ‘brassicae’ populationsPeer reviewe

Regulation of proteinaceous effector expression in phytopathogenic fungi

Effectors are molecules used by microbial pathogens to facilitate infection via effector-triggered susceptibility or tissue necrosis in their host. Much research has been focussed on the identification and elucidating the function of fungal effectors during plant pathogenesis. By comparison, knowledge of how phytopathogenic fungi regulate the expression of effector genes has been lagging. Several recent studies have illustrated the role of various transcription factors, chromosome-based control, effector epistasis, and mobilisation of endosomes within the fungal hyphae in regulating effector expression and virulence on the host plant. Improved knowledge of effector regulation is likely to assist in improving novel crop protection strategies

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

Multiple Translocation of the AVR-Pita Effector Gene among Chromosomes of the Rice Blast Fungus Magnaporthe oryzae and Related Species

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

Genome sequencing and comparative genomics of the broad host-range pathogen Rhizoctonia solani AG8

Rhizoctonia solani is a soil-borne basidiomycete fungus with a necrotrophic lifestyle which is classified into fourteen reproductively incompatible anastomosis groups (AGs). One of these, AG8, is a devastating pathogen causing bare patch of cereals, brassicas and legumes. R. solani is a multinucleate heterokaryon containing significant heterozygosity within a single cell. This complexity posed significant challenges for the assembly of its genome. We present a high quality genome assembly of R. solani AG8 and a manually curated set of 13,964 genes supported by RNA-seq. The AG8 genome assembly used novel methods to produce a haploid representation of its heterokaryotic state. The whole-genomes of AG8, the rice pathogen AG1-IA and the potato pathogen AG3 were observed to be syntenic and co-linear. Genes and functions putatively relevant to pathogenicity were highlighted by comparing AG8 to known pathogenicity genes, orthology databases spanning 197 phytopathogenic taxa and AG1-IA.We also observed SNP-level “hypermutation” of CpG dinucleotides to TpG between AG8 nuclei, with similarities to repeat-induced point mutation (RIP). Interestingly, gene-coding regions were widely affected along with repetitive DNA, which has not been previously observed for RIP in mononuclear fungi of the Pezizomycotina. The rate of heterozygous SNP mutations within this single isolate of AG8 was observed to be higher than SNP mutation rates observed across populations of most fungal species compared. Comparative analyses were combined to predict biological processes relevant to AG8 and 308 proteins with effector-like characteristics, forming a valuable resource for further study of this pathosystem. Predicted effector-like proteins had elevated levels of non-synonymous point mutations relative to synonymous mutations (dN/dS), suggesting that they may be under diversifying selection pressures. In addition, the distant relationship to sequenced necrotrophs of the Ascomycota suggests the R. solani genome sequence may prove to be a useful resource in future comparative analysis of plant pathogens