Cereal rust fungi genomics and the pursuit of virulence and avirulence factors

Countrywide inter-epidemic region migration pattern suggests the role of southwestern population in wheat stripe rust epidemics in China

Understanding countrywide pathogen population structure and inter-epidemic region spread is crucial for deciphering crop potential losses. Wheat stripe rust caused by Puccinia striiformis f. sp. tritici is a destructive disease that affects worldwide wheat production, widespread in China, representing largest epidemic region globally. This study aimed to understand the population structure and migration route of P. striiformis f. sp. tritici across China based on sampling from 15 provinces representing six epidemic zones, viz., over-summering, over-wintering, eastern, Yun-Gui, Xinjiang and Tibet epidemic regions. High genotypic diversity was recorded in over-summering, Tibet and over-wintering epidemic regions. Epidemic regions partly explain population subdivision with variable divergence (F_ST  = 0.005-0.344). Xinjiang and Tibet epidemic regions were independent epidemic zones with least sharing of genotypes. Among other epidemic zones, i.e. over-summering, over-wintering, eastern and Yun-Gui epidemic zones, re-sampling MLGs, clustering-based structure, DAPC analyses, relative migration and low divergence (F_ST from 0.006 to 0.073) revealed frequent geneflow. Yun-Gui epidemic regions, with a potential for both over-summering and over-wintering, could play an important role in causing epidemics in main wheat-cultivating areas of China. High diversity, recombination signatures and inter-epidemic region migration patterns need to be considered in host-resistant cultivar development in China and neighbouring countries, considering risk of long-distance migration capacity of pathogen.

An Avirulence Gene Cluster in the Wheat Stripe Rust Pathogen (Puccinia striiformis f. sp. tritici) Identified through Genetic Mapping and Whole-Genome Sequencing of a Sexual Population

Puccinia striiformis f. sp. tritici, the causal agent of wheat stripe (yellow) rust, is an obligate, biotrophic fungus. It was difficult to study the genetics of the pathogen due to the lack of sexual reproduction. The recent discovery of alternate hosts for P. striiformis f. sp. tritici makes it possible to study inheritance and map genes involved in its interaction with plant hosts. To identify avirulence (Avr) genes in P. striiformis f. sp. tritici, we developed a segregating population by selfing isolate 12-368 on barberry (Berberis vulgaris) plants under controlled conditions. The dikaryotic sexual population segregated for avirulent/virulent phenotypes on nine Yr single-gene lines. The parental and progeny isolates were whole-genome sequenced at >30× coverage using Illumina HiSeq PE150 technology. A total of 2,637 high-quality markers were discovered by mapping the whole-genome sequencing (WGS) reads to the reference genome of strain 93-210 and used to construct a genetic map, consisting of 41 linkage groups, spanning 7,715.0 centimorgans (cM) and covering 68 Mb of the reference genome. The recombination rate was estimated to be 1.81 ± 2.32 cM/10 kb. Quantitative trait locus analysis mapped six Avr gene loci to the genetic map, including an Avr cluster harboring four Avr genes, AvYr7, AvYr43, AvYr44, and AvYrExp2 Aligning the genetic map to the reference genome identified Avr candidates and narrowed them to a small genomic region (<200 kb). The discovery of the Avr gene cluster is useful for understanding pathogen evolution, and the identification of candidate genes is an important step toward cloning Avr genes for studying molecular mechanisms of pathogen-host interactions.IMPORTANCE Stripe rust is a destructive disease of wheat worldwide. Growing resistant cultivars is the most effective, easy-to-use, economical, and environmentally friendly strategy for the control of the disease. However, P. striiformis f. sp. tritici can produce new virulent races that may circumvent race-specific resistance. Therefore, understanding the genetic basis of the interactions between wheat genes for resistance and P. striiformis f. sp. tritici genes for avirulence is useful for improving cultivar resistance for more effective control of the disease. This study developed a high-quality map that facilitates genomic and genetic studies of important traits related to pathogen pathogenicity and adaptation to different environments and crop cultivars carrying different resistance genes. The information on avirulence/virulence genes identified in this study can be used for guiding breeding programs to select combinations of genes for developing new cultivars with effective resistance to mitigate this devastating disease.

The intrinsic and regulated proteomes of barley seeds in response to fungal infection

Barley is an important cereal grain used for beer brewing, animal feed, and human food consumption. Fungal disease can impact barley production, as it causes substantial yield loss and lowers seed quality. We used sequential window acquisition of all theoretical ions mass spectrometry (SWATH-MS) to measure and quantify the relative abundance of proteins within seeds of different barley varieties under various fungal pathogen burdens (ProteomeXchange Datasets PXD011303 and PXD014093). Fungal burden in the leaves and stems of barley resulted in changes to the seed proteome. However, these changes were minimal and showed substantial variation among barley samples infected with different pathogens. The limited effect of intrinsic disease resistance on the seed proteome is consistent with the main mediators of disease resistance being present in the leaves and stems of the plant. The seeds of barley varieties accredited for use as malt had higher levels of proteins associated with starch synthesis and beer quality. The proteomic workflows developed and implemented here have potential application in quality control, breeding and processing of barley, and other agricultural products.

Host-induced silencing of essential genes in Puccinia triticina through transgenic expression of RNAi sequences reduces severity of leaf rust infection in wheat

Leaf rust, caused by the pathogenic fungus Puccinia triticina (Pt), is one of the most serious biotic threats to sustainable wheat production worldwide. This obligate biotrophic pathogen is prevalent worldwide and is known for rapid adaptive evolution to overcome resistant wheat varieties. Novel disease control approaches are therefore required to minimize the yield losses caused by Pt. Having shown previously the potential of host-delivered RNA interference (HD-RNAi) in functional screening of Pt genes involved in pathogenesis, we here evaluated the use of this technology in transgenic wheat plants as a method to achieve protection against wheat leaf rust (WLR) infection. Stable expression of hairpin RNAi constructs with sequence homology to Pt MAP-kinase (PtMAPK1) or a cyclophilin (PtCYC1) encoding gene in susceptible wheat plants showed efficient silencing of the corresponding genes in the interacting fungus resulting in disease resistance throughout the T2 generation. Inhibition of Pt proliferation in transgenic lines by in planta-induced RNAi was associated with significant reduction in target fungal transcript abundance and reduced fungal biomass accumulation in highly resistant plants. Disease protection was correlated with the presence of siRNA molecules specific to targeted fungal genes in the transgenic lines harbouring the complementary HD-RNAi construct. This work demonstrates that generating transgenic wheat plants expressing RNAi-inducing transgenes to silence essential genes in rust fungi can provide effective disease resistance, thus opening an alternative way for developing rust-resistant crops.

Divergent and convergent modes of interaction between wheat and Puccinia graminis f. sp. tritici isolates revealed by the comparative gene co-expression network and genome analyses

BACKGROUND

Two opposing evolutionary constraints exert pressure on plant pathogens: one to diversify virulence factors in order to evade plant defenses, and the other to retain virulence factors critical for maintaining a compatible interaction with the plant host. To better understand how the diversified arsenals of fungal genes promote interaction with the same compatible wheat line, we performed a comparative genomic analysis of two North American isolates of Puccinia graminis f. sp. tritici (Pgt).

RESULTS

The patterns of inter-isolate divergence in the secreted candidate effector genes were compared with the levels of conservation and divergence of plant-pathogen gene co-expression networks (GCN) developed for each isolate. Comprative genomic analyses revealed substantial level of interisolate divergence in effector gene complement and sequence divergence. Gene Ontology (GO) analyses of the conserved and unique parts of the isolate-specific GCNs identified a number of conserved host pathways targeted by both isolates. Interestingly, the degree of inter-isolate sub-network conservation varied widely for the different host pathways and was positively associated with the proportion of conserved effector candidates associated with each sub-network. While different Pgt isolates tended to exploit similar wheat pathways for infection, the mode of plant-pathogen interaction varied for different pathways with some pathways being associated with the conserved set of effectors and others being linked with the diverged or isolate-specific effectors.

CONCLUSIONS

Our data suggest that at the intra-species level pathogen populations likely maintain divergent sets of effectors capable of targeting the same plant host pathways. This functional redundancy may play an important role in the dynamic of the "arms-race" between host and pathogen serving as the basis for diverse virulence strategies and creating conditions where mutations in certain effector groups will not have a major effect on the pathogen's ability to infect the host.

Comparative Genomics Integrated with Association Analysis Identifies Candidate Effector Genes Corresponding to Lr20 in Phenotype-Paired Puccinia triticina Isolates from Australia

Leaf rust is one of the most common and damaging diseases of wheat, and is caused by an obligate biotrophic basidiomycete, Puccinia triticina (Pt). In the present study, 20 Pt isolates from Australia, comprising 10 phenotype-matched pairs with contrasting pathogenicity for Lr20, were analyzed using whole genome sequencing. Compared to the reference genome of the American Pt isolate 1-1 BBBD Race 1, an average of 404,690 single nucleotide polymorphisms (SNPs) per isolate was found and the proportion of heterozygous SNPs was above 87% in the majority of the isolates, demonstrating a high level of polymorphism and a high rate of heterozygosity. From the genome-wide SNPs, a phylogenetic tree was inferred, which consisted of a large clade of 15 isolates representing diverse presumed clonal lineages including 14 closely related isolates and the more diverged isolate 670028, and a small clade of five isolates characterized by lower heterozygosity level. Principle component analysis detected three distinct clusters, corresponding exactly to the two major subsets of the small clade and the large clade comprising all 15 isolates without further separation of isolate 670028. While genome-wide association analysis identified 302 genes harboring at least one SNP associated with Lr20 virulence (p < 0.05), a Wilcoxon rank sum test revealed that 36 and 68 genes had significant (p < 0.05) and marginally significant (p < 0.1) differences in the counts of non-synonymous mutations between Lr20 avirulent and virulent groups, respectively. Twenty of these genes were predicted to have a signal peptide without a transmembrane segment, and hence identified as candidate effector genes corresponding to Lr20. SNP analysis also implicated the potential involvement of epigenetics and small RNA in Pt pathogenicity. Future studies are thus warranted to investigate the biological functions of the candidate effectors as well as the gene regulation mechanisms at epigenetic and post-transcription levels. Our study is the first to integrate phenotype-genotype association with effector prediction in Pt genomes, an approach that may circumvent some of the technical difficulties in working with obligate rust fungi and accelerate avirulence gene identification.

Antifungal peptides: To be or not to be membrane active

Most antifungal peptides (AFPs), if not all, have membrane activity, while some also have alternative targets. Fungal membranes share many characteristics with mammalian membranes with only a few differences, such as differences in sphingolipids, phosphatidylinositol (PI) content and the main sterol is ergosterol. Fungal membranes are also more negative and a better target for cationic AFPs. Targeting just the fungal membrane lipids such as phosphatidylinositol and/or ergosterol by AFPs often translates into mammalian cell toxicity. Conversely, a specific AFP target in the fungal pathogen, such as glucosylceramide, mannosyldiinositol phosphorylceramide or a fungal protein target translates into high pathogen selectivity. However, a lower target concentration, absence or change in the specific fungal target can naturally lead to resistance, although such resistance in turn could result in reduced pathogen virulence. The question is then to be or not to be membrane active - what is the best choice for a successful AFP? In this review we deliberate on this question by focusing on the recent advances in our knowledge on how natural AFPs target fungi.

High genome heterozygosity and endemic genetic recombination in the wheat stripe rust fungus

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most destructive diseases of wheat. Here we report a 110-Mb draft sequence of Pst isolate CY32, obtained using a 'fosmid-to-fosmid' strategy, to better understand its race evolution and pathogenesis. The Pst genome is highly heterozygous and contains 25,288 protein-coding genes. Compared with non-obligate fungal pathogens, Pst has a more diverse gene composition and more genes encoding secreted proteins. Re-sequencing analysis indicates significant genetic variation among six isolates collected from different continents. Approximately 35% of SNPs are in the coding sequence regions, and half of them are non-synonymous. High genetic diversity in Pst suggests that sexual reproduction has an important role in the origin of different regional races. Our results show the effectiveness of the 'fosmid-to-fosmid' strategy for sequencing dikaryotic genomes and the feasibility of genome analysis to understand race evolution in Pst and other obligate pathogens.

Proteomic characterization of the Rph15 barley resistance gene-mediated defence responses to leaf rust

Background

Leaf rust, caused by the biotrophic fungal pathogen Puccinia hordei, is one of the most important foliar disease of barley (Hordeum vulgare) and represents a serious threat in many production regions of the world. The leaf rust resistance gene Rph15 is of outstanding interest for resistance breeding because it confers resistance to over 350 Puccinia hordei isolates collected from around the world. Molecular and biochemical mechanisms responsible for the Rph15 effectiveness are currently not investigated. The aim of the present work was to study the Rph15-based defence responses using a proteomic approach.

Results

Protein pattern changes in response to the leaf rust pathogen infection were investigated in two barley near isogenic lines (NILs), Bowman (leaf rust susceptible) and Bowman-Rph15 (leaf rust resistant), differing for the introgression of the leaf rust resistance gene Rph15. Two infection time points, 24 hours and four days post inoculation (dpi), were analysed. No statistically significant differences were identified at the early time point, while at 4 dpi eighteen protein spots were significantly up or down regulated with a fold-change equal or higher than two in response to pathogen infection. Almost all the pathogen-responsive proteins were identified in the Bowman-Rph15 resistant NIL. Protein spots were characterized by LC-MS/MS analysis and found to be involved in photosynthesis and energy metabolism, carbohydrate metabolism, protein degradation and defence. Proteomic data were complemented by transcriptional analysis of the respective genes. The identified proteins can be related to modulation of the photosynthetic apparatus components, re-direction of the metabolism to sustain defence responses and deployment of defence proteins.

Conclusions

The identification of leaf rust infection-modulated defence responses restricted to the resistant NIL support the hypothesis that basal defence responses of Bowman, but not the Rph15 resistance gene-based ones, are suppressed or delayed by pathogen effectors to levels below the detection power of the adopted proteomic approach. Additionally, Rph15-mediated resistance processes identified mainly resides on a modulation of primary metabolism, affecting photosyntesis and carbohydrate pool.

Lr34-mediated leaf rust resistance in wheat: transcript profiling reveals a high energetic demand supported by transient recruitment of multiple metabolic pathways

The wheat gene Lr34 confers partial resistance to all races of Puccinia triticina, the causal agent of wheat leaf rust. However, the biological basis for the exceptional durability of Lr34 is unclear. We used the Affymetrix GeneChip Wheat Genome Array to compare transcriptional changes of near-isogenic lines of Thatcher wheat in a compatible interaction, an incompatible interaction conferred by the resistance gene Lr1, and the race-nonspecific response conditioned by Lr34 3 and 7 days postinoculation (dpi) with P. triticina. No differentially expressed genes were detected in Lr1 plants at either timepoint whereas, in the compatible Thatcher interaction, differentially expressed genes were detected only at 7 dpi. In contrast, differentially expressed genes were identified at both timepoints in P. triticina-inoculated Lr34 plants. At 3 dpi, upregulated genes associated with Lr34-mediated resistance encoded various defense and stress-related proteins, secondary metabolism enzymes, and transcriptional regulation and cellular-signaling proteins. Further, coordinated upregulation of key genes in several metabolic pathways that can contribute to increased carbon flux through the tricarboxylic cycle was detected. This indicates that Lr34-mediated resistance imposes a high energetic demand that leads to the induction of multiple metabolic responses to support cellular energy requirements. These metabolic responses were not sustained through 7 dpi, and may explain why Lr34 fails to inhibit the pathogen fully but does increase the latent period.

Plant-pathogen interactions: what is proteomics telling us?

Over the years, several studies have been performed to analyse plant-pathogen interactions. Recently, functional genomic strategies, including proteomics and transcriptomics, have contributed to the effort of defining gene and protein function and expression profiles. Using these 'omic' approaches, pathogenicity- and defence-related genes and proteins expressed during phytopathogen infections have been identified and enormous datasets have been accumulated. However, the understanding of molecular plant-pathogen interactions is still an intriguing area of investigation. Proteomics has dramatically evolved in the pursuit of large-scale functional assignment of candidate proteins and, by using this approach, several proteins expressed during phytopathogenic interactions have been identified. In this review, we highlight the proteins expressed during plant-virus, plant-bacterium, plant-fungus and plant-nematode interactions reported in proteomic studies, and discuss these findings considering the advantages and limitations of current proteomic tools.

Leaf rust resistance gene Lr1, isolated from bread wheat (Triticum aestivum L.) is a member of the large psr567 gene family

In hexaploid wheat, leaf rust resistance gene Lr1 is located at the distal end of the long arm of chromosome 5D. To clone this gene, an F(1)-derived doubled haploid population and a recombinant inbred line population from a cross between the susceptible cultivar AC Karma and the resistant line 87E03-S2B1 were phenotyped for resistance to Puccinia triticina race 1-1 BBB that carries the avirulence gene Avr1. A high-resolution genetic map of the Lr1 locus was constructed using microsatellite, resistance gene analog (RGA), BAC end (BE), and low pass (LP) markers. A physical map of the locus was constructed by screening a hexaploid wheat BAC library from cultivar Glenlea that is known to have Lr1. The locus comprised three RGAs from a gene family related to RFLP marker Xpsr567. Markers specific to each paralog were developed. Lr1 segregated with RGA567-5 while recombinants were observed for the other two RGAs. Transformation of the susceptible cultivar Fielder with RGA567-5 demonstrated that it corresponds to the Lr1 resistance gene. In addition, the candidate gene was also confirmed by virus-induced gene silencing. Twenty T (1) lines from resistant transgenic line T (0)-938 segregated for resistance, partial resistance and susceptibility to Avr1 corresponding to a 1:2:1 ratio for a single hemizygous insertion. Transgene presence and expression correlated with the phenotype. The resistance phenotype expressed by Lr1 seemed therefore to be dependant on the zygosity status. T (3)-938 sister lines with and without the transgene were further tested with 16 virulent and avirulent rust isolates. Rust reactions were all as expected for Lr1 thereby providing additional evidence toward the Lr1 identity of RGA567-5. Sequence analysis of Lr1 indicated that it is not related to the previously isolated Lr10 and Lr21 genes and unlike these genes, it is part of a large gene family.