The wheat Pm3 resistance gene against the powdery mildew pathogen occurs as an allelic series encoding functionally different immune receptors which induce resistance upon recognition of isolate-specific avirulence (AVR) effectors from the pathogen. Here, we describe the identification of five effector proteins from the mildew pathogens of wheat, rye, and the wild grass Dactylis glomerata , specifically recognized by the PM3B, PM3C and PM3D receptors. Together with the earlier identified AVRPM3 A2/F2 , the recognized AVRs of PM3B/C, (AVRPM3 B2/C2 ), and PM3D (AVRPM3 D3 ) belong to a large group of proteins with low sequence homology but predicted structural similarities. AvrPm3 b2/c2 and AvrPm3 d3 are conserved in all tested isolates of wheat and rye mildew, and non-host infection assays demonstrate that Pm3b , Pm3c , and Pm3d are also restricting the growth of rye mildew on wheat. Furthermore, divergent AVR homologues from non-adapted rye and Dactylis mildews are recognized by PM3B, PM3C, or PM3D, demonstrating their involvement in host specificity.

Summary Blumeria graminis f. sp. tritici ( B.g. tritici ) is the causal agent of the wheat powdery mildew disease. The highly fragmented B.g. tritici genome available so far has prevented a systematic analysis of effector genes that are known to be involved in host adaptation. To study the diversity and evolution of effector genes we produced a chromosome‐scale assembly of the B.g. tritici genome. The genome assembly and annotation was achieved by combining long‐read sequencing with high‐density genetic mapping, bacterial artificial chromosome fingerprinting and transcriptomics. We found that the 166.6 Mb B.g. tritici genome encodes 844 candidate effector genes, over 40% more than previously reported. Candidate effector genes have characteristic local genomic organization such as gene clustering and enrichment for recombination‐active regions and certain transposable element families. A large group of 412 candidate effector genes shows high plasticity in terms of copy number variation in a global set of 36 isolates and of transcription levels. Our data suggest that copy number variation and transcriptional flexibility are the main drivers for adaptation in B.g. tritici . The high repeat content may play a role in providing a genomic environment that allows rapid evolution of effector genes with selection as the driving force.

 Pm1a, the first powdery mildew resistance gene described in wheat, is part of a complex resistance (R) gene cluster located in a distal region of chromosome 7AL that has suppressed genetic recombination.  A nucleotide-binding, leucine-rich repeat (NLR) immune receptor gene was isolated using mutagenesis and R gene enrichment sequencing (MutRenSeq). Stable transformation confirmed Pm1a identity which induced a strong resistance phenotype in transgenic plants upon challenge with avirulent Blumeria graminis (wheat powdery mildew) pathogens.  A high-density genetic map of a B. graminis family segregating for Pm1a avirulence combined with pathogen genome resequencing and RNAseq identified AvrPm1a effector gene candidates. In planta expression identified an effector, with an N terminal Y/FxC motif, that induced a strong hypersensitive response when co-expressed with Pm1a in Nicotiana benthamiana.  Single chromosome enrichment sequencing (ChromSeq) and assembly of chromosome 7A suggested that suppressed recombination around the Pm1a region was due to a rearrangement involving chromosomes 7A, 7B and 7D. The cloning of Pm1a and its identification in a highly rearranged region of chromosome 7A provides insight into the role of chromosomal rearrangements in the evolution of this complex resistance cluster.

Introgressions of chromosomal segments from related species into wheat are important sources of resistance against fungal diseases. The durability and effectiveness of introgressed resistance genes upon agricultural deployment is highly variable—a phenomenon that remains poorly understood, as the corresponding fungal avirulence genes are largely unknown. Until its breakdown, the Pm17 resistance gene introgressed from rye to wheat provided broad resistance against powdery mildew ( Blumeria graminis ). Here, we used quantitative trait locus (QTL) mapping to identify the corresponding wheat mildew avirulence effector AvrPm17 . It is encoded by two paralogous genes that exhibit signatures of reoccurring gene conversion events and are members of a mildew sublineage specific effector cluster. Extensive haplovariant mining in wheat mildew and related sublineages identified several ancient virulent AvrPm17 variants that were present as standing genetic variation in wheat powdery mildew prior to the Pm17 introgression, thereby paving the way for the rapid breakdown of the Pm17 resistance. QTL mapping in mildew identified a second genetic component likely corresponding to an additional resistance gene present on the 1AL.1RS translocation carrying Pm17. This gene remained previously undetected due to suppressed recombination within the introgressed rye chromosomal segment. We conclude that the initial effectiveness of 1AL.1RS was based on simultaneous introgression of two genetically linked resistance genes. Our results demonstrate the relevance of pathogen-based genetic approaches to disentangling complex resistance loci in wheat. We propose that identification and monitoring of avirulence gene diversity in pathogen populations become an integral part of introgression breeding to ensure effective and durable resistance in wheat.

Summary Recognition of the AVRPM3A2/F2 avirulence protein from powdery mildew by the wheat PM3A/F immune receptor induces a hypersensitive response after co‐expression in Nicotiana benthamiana. The molecular determinants of this interaction and how they shape natural AvrPm3 a2/f2 allelic diversity are unknown.We sequenced the AvrPm3 a2/f2 gene in a worldwide collection of 272 mildew isolates. Using the natural polymorphisms of AvrPm3 a2/f2 as well as sequence information from related gene family members, we tested 85 single‐residue‐altered AVRPM3A2/F2 variants with PM3A, PM3F and PM3FL 456P/Y458H (modified for improved signaling) in Nicotiana benthamiana for effects on recognition.An intact AvrPm3 a2/f2 gene was found in all analyzed isolates and the protein variant recognized by PM3A/F occurred globally at high frequencies. Single‐residue alterations in AVRPM3A2/F2 mostly disrupted, but occasionally enhanced, the recognition response by PM3A, PM3F and PM3FL 456P/Y458H. Residues enhancing hypersensitive responses constituted a protein domain separate from both naturally occurring polymorphisms and positively selected residues of the gene family.These results demonstrate the utility of using gene family sequence diversity to screen residues for their role in recognition. This approach identified a putative interaction surface in AVRPM3A2/F2 not polymorphic in natural alleles. We conclude that molecular mechanisms besides recognition drive AvrPm3 a2/f2 diversification.