Search

Background Studying genomic variation in rapidly evolving pathogens potentially enables identification of genes supporting their “core biology”, being present, functional and expressed by all strains or “flexible biology”, varying between strains. Genes supporting flexible biology may be considered to be “accessory”, whilst the “core” gene set is likely to be important for common features of a pathogen species biology, including virulence on all host genotypes. The wheat-pathogenic fungus Zymoseptoria tritici represents one of the most rapidly evolving threats to global food security and was the focus of this study. Results We constructed a pangenome of 18 European field isolates, with 12 also subjected to RNAseq transcription profiling during infection. Combining this data, we predicted a “core” gene set comprising 9807 sequences which were (1) present in all isolates, (2) lacking inactivating polymorphisms and (3) expressed by all isolates. A large accessory genome, consisting of 45% of the total genes, was also defined. We classified genetic and genomic polymorphism at both chromosomal and individual gene scales. Proteins required for essential functions including virulence had lower-than average sequence variability amongst core genes. Both core and accessory genomes encoded many small, secreted candidate effector proteins that likely interact with plant immunity. Viral vector-mediated transient in planta overexpression of 88 candidates failed to identify any which induced leaf necrosis characteristic of disease. However, functional complementation of a non-pathogenic deletion mutant lacking five core genes demonstrated that full virulence was restored by re-introduction of the single gene exhibiting least sequence polymorphism and highest expression. Conclusions These data support the combined use of pangenomics and transcriptomics for defining genes which represent core, and potentially exploitable, weaknesses in rapidly evolving pathogens.

Zymoseptoria tritici causes Septoria tritici blotch (STB) on wheat. The disease interaction is characterized by clearly defined temporal phases of infection, ultimately resulting in the death of host tissue. Zymoseptoria tritici is a highly polymorphic species with significant intraspecific variation in virulence profiles. We generated a deep transcriptomic sequencing dataset spanning the entire time course of an infection using a previously uncharacterized, highly virulent Z. tritici strain isolated from a Swiss wheat field. We found that seven clusters of gene transcription profiles explained the progression of the infection. The earliest highly up-regulated genes included chloroperoxidases, which may help the fungus cope with plant defences. The onset of necrotrophy was characterized by a concerted up-regulation of proteases, plant cell wall-degrading enzymes and lipases. Functions related to nutrition and growth characterized late necrotrophy and the transition to saprotrophic growth on dead plant tissue. We found that the peak up-regulation of genes essential for mating coincided with the necrotrophic phase. We performed an intraspecies comparative transcriptomics analysis using a comparable time course infection experiment of the genome reference isolate IPO323. Major components of the fungal infection transcriptome were conserved between the two strains. However, individual small, secreted proteins, proteases and cell wall-degrading enzymes showed strongly differentiated transcriptional profiles between isolates. Our analyses illustrate that successful STB infections involve complex transcriptomic remodelling to up-regulate distinct gene functions. Heterogeneity in transcriptomes among isolates may explain some of the considerable variation in virulence and host specialization found within the species.