Fungi are of primary ecological, biotechnological and economic importance. Many fundamental biological processes that are shared by animals and fungi are studied in fungi due to their experimental tractability. Many fungi are pathogens or mutualists and are model systems to analyse effector genes and their mechanisms of diversification. In this study, we report the genome sequence of the phytopathogenic ascomycete Leptosphaeria maculans and characterize its repertoire of protein effectors. The L. maculans genome has an unusual bipartite structure with alternating distinct guanine and cytosine-equilibrated and adenine and thymine (AT)-rich blocks of homogenous nucleotide composition. The AT-rich blocks comprise one-third of the genome and contain effector genes and families of transposable elements, both of which are affected by repeat-induced point mutation, a fungal-specific genome defence mechanism. This genomic environment for effectors promotes rapid sequence diversification and underpins the evolutionary potential of the fungus to adapt rapidly to novel host-derived constraints.

Throughout the history of agriculture, many new crop species (polyploids or artificial hybrids) have been introduced to diversify products or to increase yield. However, little is known about how these new crops influence the evolution of new pathogens and diseases. Triticale is an artificial hybrid of wheat and rye, and it was resistant to the fungal pathogen powdery mildew (Blumeria graminis) until 2001 (refs. 1,2,3). We sequenced and compared the genomes of 46 powdery mildew isolates covering several formae speciales. We found that B. graminis f. sp. triticale, which grows on triticale and wheat, is a hybrid between wheat powdery mildew (B. graminis f. sp. tritici) and mildew specialized on rye (B. graminis f. sp. secalis). Our data show that the hybrid of the two mildews specialized on two different hosts can infect the hybrid plant species originating from those two hosts. We conclude that hybridization between mildews specialized on different species is a mechanism of adaptation to new crops introduced by agriculture.

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.

In cereals, several mildew resistance genes occur as large allelic series; for example, in wheat (Triticum aestivum and Triticum turgidum), 17 functional Pm3 alleles confer agronomically important race-specific resistance to powdery mildew (Blumeria graminis). The molecular basis of race specificity has been characterized in wheat, but little is known about the corresponding avirulence genes in powdery mildew. Here, we dissected the genetics of avirulence for six Pm3 alleles and found that three major Avr loci affect avirulence, with a common locus_1 involved in all AvrPm3-Pm3 interactions. We cloned the effector gene AvrPm3^a2/f2 from locus_2, which is recognized by the Pm3a and Pm3f alleles. Induction of a Pm3 allele-dependent hypersensitive response in transient assays in Nicotiana benthamiana and in wheat demonstrated specificity. Gene expression analysis of Bcg1 (encoded by locus_1) and AvrPm3 ^a2/f2 revealed significant differences between isolates, indicating that in addition to protein polymorphisms, expression levels play a role in avirulence. We propose a model for race specificity involving three components: an allele-specific avirulence effector, a resistance gene allele, and a pathogen-encoded suppressor of avirulence. Thus, whereas a genetically simple allelic series controls specificity in the plant host, recognition on the pathogen side is more complex, allowing flexible evolutionary responses and adaptation to resistance genes.

Summary: There is a large diversity of genetically defined resistance genes in bread wheat against the powdery mildew pathogen Blumeria graminis (B. g.) f. sp. tritici. Many confer race‐specific resistance to this pathogen, but until now only the mildew avirulence gene AvrPm3^a2/f2 that is recognized by Pm3a/f was known molecularly. We performed map‐based cloning and genome‐wide association studies to isolate a candidate for the mildew avirulence gene AvrPm2. We then used transient expression assays in Nicotiana benthamiana to demonstrate specific and strong recognition of AvrPm2 by Pm2. The virulent AvrPm2 allele arose from a conserved 12 kb deletion, while there is no protein sequence diversity in the gene pool of avirulent B. g. tritici isolates. We found one polymorphic AvrPm2 allele in B. g. triticale and one orthologue in B. g. secalis and both are recognized by Pm2. AvrPm2 belongs to a small gene family encoding structurally conserved RNase‐like effectors, including Avr_a13 from B. g. hordei, the cognate Avr of the barley resistance gene Mla13. These results demonstrate the conservation of functional avirulence genes in two cereal powdery mildews specialized on different hosts, thus providing a possible explanation for successful introgression of resistance genes from rye or other grass relatives to wheat.

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.