Pooled effector library screening in protoplasts rapidly identifies novel Avr genes.

Crop breeding for durable disease resistance is challenging due to the rapid evolution of pathogen virulence. While progress in resistance (R) gene cloning and stacking has accelerated in recent years1-3, the identification of corresponding avirulence (Avr) genes in many pathogens is hampered by the lack of high-throughput screening options. To address this technology gap, we developed a platform for pooled library screening in plant protoplasts to allow rapid identification of interacting R-Avr pairs. We validated this platform by isolating known and novel Avr genes from wheat stem rust (Puccinia graminis f. sp. tritici) after screening a designed library of putative effectors against individual R genes. Rapid Avr gene identification provides molecular tools to understand and track pathogen virulence evolution via genotype surveillance, which in turn will lead to optimized R gene stacking and deployment strategies. This platform should be broadly applicable to many crop pathogens and could potentially be adapted for screening genes involved in other protoplast-selectable traits.

Haplotype variants of the stripe rust resistance gene Yr28 in Aegilops tauschii.

KEY MESSAGE: Stripe rust resistance gene YrAet672 from Aegilops tauschii accession CPI110672 encodes a nucleotide-binding and leucine-rich repeat domain containing protein similar to YrAS2388 and both these members were haplotypes of Yr28. New sources of host resistance are required to counter the continued emergence of new pathotypes of the wheat stripe rust pathogen Puccinia striiformis Westend. f. sp. tritici Erikss. (Pst). Here, we show that CPI110672, an Aegilops tauschii accession from Turkmenistan, carries a single Pst resistance gene, YrAet672, that is effective against multiple Pst pathotypes, including the four predominant Pst lineages present in Australia. The YRAet672 locus was fine mapped to the short arm of chromosome 4D, and a nucleotide-binding and leucine-rich repeat gene was identified at the locus. A transgene encoding the YrAet672 genomic sequence, but lacking a copy of a duplicated sequence present in the 3' UTR, was transformed into wheat cultivar Fielder and Avocet S. This transgene conferred a weak resistance response, suggesting that the duplicated 3' UTR region was essential for function. Subsequent analyses demonstrated that YrAet672 is the same as two other Pst resistance genes described in Ae. tauschii, namely YrAS2388 and Yr28. They were identified as haplotypes encoding identical protein sequences but are polymorphic in non-translated regions of the gene. Suppression of resistance conferred by YrAet672 and Yr28 in synthetic hexaploid wheat lines (AABBDD) involving Langdon (AABB) as the tetraploid parent was associated with a reduction in transcript accumulation.

Physical separation of haplotypes in dikaryons allows benchmarking of phasing accuracy in Nanopore and HiFi assemblies with Hi-C data.

BACKGROUND: Most animals and plants have more than one set of chromosomes and package these haplotypes into a single nucleus within each cell. In contrast, many fungal species carry multiple haploid nuclei per cell. Rust fungi are such species with two nuclei (karyons) that contain a full set of haploid chromosomes each. The physical separation of haplotypes in dikaryons means that, unlike in diploids, Hi-C chromatin contacts between haplotypes are false-positive signals. RESULTS: We generate the first chromosome-scale, fully-phased assembly for the dikaryotic leaf rust fungus Puccinia triticina and compare Nanopore MinION and PacBio HiFi sequence-based assemblies. We show that false-positive Hi-C contacts between haplotypes are predominantly caused by phase switches rather than by collapsed regions or Hi-C read mis-mappings. We introduce a method for phasing of dikaryotic genomes into the two haplotypes using Hi-C contact graphs, including a phase switch correction step. In the HiFi assembly, relatively few phase switches occur, and these are predominantly located at haplotig boundaries and can be readily corrected. In contrast, phase switches are widespread throughout the Nanopore assembly. We show that haploid genome read coverage of 30-40 times using HiFi sequencing is required for phasing of the leaf rust genome, with 0.7% heterozygosity, and that HiFi sequencing resolves genomic regions with low heterozygosity that are otherwise collapsed in the Nanopore assembly. CONCLUSIONS: This first Hi-C based phasing pipeline for dikaryons and comparison of long-read sequencing technologies will inform future genome assembly and haplotype phasing projects in other non-haploid organisms.

The stem rust effector protein AvrSr50 escapes Sr50 recognition by a substitution in a single surface-exposed residue.

Pathogen effectors are crucial players during plant colonisation and infection. Plant resistance mostly relies on effector recognition to activate defence responses. Understanding how effector proteins escape from plant surveillance is important for plant breeding and resistance deployment. Here we examined the role of genetic diversity of the stem rust (Puccinia graminis f. sp. tritici (Pgt)) AvrSr50 gene in determining recognition by the corresponding wheat Sr50 resistance gene. We solved the crystal structure of a natural variant of AvrSr50 and used site-directed mutagenesis and transient expression assays to dissect the molecular mechanisms explaining gain of virulence. We report that AvrSr50 can escape recognition by Sr50 through different mechanisms including DNA insertion, stop codon loss or by amino-acid variation involving a single substitution of the AvrSr50 surface-exposed residue Q121. We also report structural homology of AvrSr50 to cupin superfamily members and carbohydrate-binding modules indicating a potential role in binding sugar moieties. This study identifies key polymorphic sites present in AvrSr50 alleles from natural stem rust populations that play important roles to escape from Sr50 recognition. This constitutes an important step to better understand Pgt effector evolution and to monitor AvrSr50 variants in natural rust populations.

The stem rust fungus Puccinia graminis f. sp. tritici induces centromeric small RNAs during late infection that are associated with genome-wide DNA methylation.

BACKGROUND: Silencing of transposable elements (TEs) is essential for maintaining genome stability. Plants use small RNAs (sRNAs) to direct DNA methylation to TEs (RNA-directed DNA methylation; RdDM). Similar mechanisms of epigenetic silencing in the fungal kingdom have remained elusive. RESULTS: We use sRNA sequencing and methylation data to gain insight into epigenetics in the dikaryotic fungus Puccinia graminis f. sp. tritici (Pgt), which causes the devastating stem rust disease on wheat. We use Hi-C data to define the Pgt centromeres and show that they are repeat-rich regions (~250 kb) that are highly diverse in sequence between haplotypes and, like in plants, are enriched for young TEs. DNA cytosine methylation is particularly active at centromeres but also associated with genome-wide control of young TE insertions. Strikingly, over 90% of Pgt sRNAs and several RNAi genes are differentially expressed during infection. Pgt induces waves of functionally diversified sRNAs during infection. The early wave sRNAs are predominantly 21 nts with a 5' uracil derived from genes. In contrast, the late wave sRNAs are mainly 22-nt sRNAs with a 5' adenine and are strongly induced from centromeric regions. TEs that overlap with late wave sRNAs are more likely to be methylated, both inside and outside the centromeres, and methylated TEs exhibit a silencing effect on nearby genes. CONCLUSIONS: We conclude that rust fungi use an epigenetic silencing pathway that might have similarity with RdDM in plants. The Pgt RNAi machinery and sRNAs are under tight temporal control throughout infection and might ensure genome stability during sporulation.

Genomics accelerated isolation of a new stem rust avirulence gene-wheat resistance gene pair.

Stem rust caused by the fungus Puccinia graminis f. sp. tritici (Pgt) is a devastating disease of the global staple crop wheat. Although this disease was largely controlled in the latter half of the twentieth century, new virulent strains of Pgt, such as Ug99, have recently evolved1,2. These strains have caused notable losses worldwide and their continued spread threatens global wheat production. Breeding for disease resistance provides the most cost-effective control of wheat rust diseases3. A number of rust resistance genes have been characterized in wheat and most encode immune receptors of the nucleotide-binding leucine-rich repeat (NLR) class4, which recognize pathogen effector proteins known as avirulence (Avr) proteins5. However, only two Avr genes have been identified in Pgt so far, AvrSr35 and AvrSr50 (refs. 6,7), and none in other cereal rusts8,9. The Sr27 resistance gene was first identified in a wheat line carrying an introgression of the 3R chromosome from Imperial rye10. Although not deployed widely in wheat, Sr27 is widespread in the artificial crop species Triticosecale (triticale), which is a wheat-rye hybrid and is a host for Pgt11,12. Sr27 is effective against Ug99 (ref. 13) and other recent Pgt strains14,15. Here, we identify both the Sr27 gene in wheat and the corresponding AvrSr27 gene in Pgt and show that virulence to Sr27 can arise experimentally and in the field through deletion mutations, copy number variation and expression level polymorphisms at the AvrSr27 locus.

A recombined Sr26 and Sr61 disease resistance gene stack in wheat encodes unrelated NLR genes.

The re-emergence of stem rust on wheat in Europe and Africa is reinforcing the ongoing need for durable resistance gene deployment. Here, we isolate from wheat, Sr26 and Sr61, with both genes independently introduced as alien chromosome introgressions from tall wheat grass (Thinopyrum ponticum). Mutational genomics and targeted exome capture identify Sr26 and Sr61 as separate single genes that encode unrelated (34.8%) nucleotide binding site leucine rich repeat proteins. Sr26 and Sr61 are each validated by transgenic complementation using endogenous and/or heterologous promoter sequences. Sr61 orthologs are absent from current Thinopyrum elongatum and wheat pan genome sequences, contrasting with Sr26 where homologues are present. Using gene-specific markers, we validate the presence of both genes on a single recombinant alien segment developed in wheat. The co-location of these genes on a small non-recombinogenic segment simplifies their deployment as a gene stack and potentially enhances their resistance durability.

A five-transgene cassette confers broad-spectrum resistance to a fungal rust pathogen in wheat.

Breeding wheat with durable resistance to the fungal pathogen Puccinia graminis f. sp. tritici (Pgt), a major threat to cereal production, is challenging due to the rapid evolution of pathogen virulence. Increased durability and broad-spectrum resistance can be achieved by introducing more than one resistance gene, but combining numerous unlinked genes by breeding is laborious. Here we generate polygenic Pgt resistance by introducing a transgene cassette of five resistance genes into bread wheat as a single locus and show that at least four of the five genes are functional. These wheat lines are resistant to aggressive and highly virulent Pgt isolates from around the world and show very high levels of resistance in the field. The simple monogenic inheritance of this multigene locus greatly simplifies its use in breeding. However, a new Pgt isolate with virulence to several genes at this locus suggests gene stacks will need strategic deployment to maintain their effectiveness.

Stem rust resistance in wheat is suppressed by a subunit of the mediator complex.

Stem rust is an important disease of wheat that can be controlled using resistance genes. The gene SuSr-D1 identified in cultivar 'Canthatch' suppresses stem rust resistance. SuSr-D1 mutants are resistant to several races of stem rust that are virulent on wild-type plants. Here we identify SuSr-D1 by sequencing flow-sorted chromosomes, mutagenesis, and map-based cloning. The gene encodes Med15, a subunit of the Mediator Complex, a conserved protein complex in eukaryotes that regulates expression of protein-coding genes. Nonsense mutations in Med15b.D result in expression of stem rust resistance. Time-course RNAseq analysis show a significant reduction or complete loss of differential gene expression at 24 h post inoculation in med15b.D mutants, suggesting that transcriptional reprogramming at this time point is not required for immunity to stem rust. Suppression is a common phenomenon and this study provides novel insight into suppression of rust resistance in wheat.

Emergence of the Ug99 lineage of the wheat stem rust pathogen through somatic hybridisation.

Parasexuality contributes to diversity and adaptive evolution of haploid (monokaryotic) fungi. However, non-sexual genetic exchange mechanisms are not defined in dikaryotic fungi (containing two distinct haploid nuclei). Newly emerged strains of the wheat stem rust pathogen, Puccinia graminis f. sp. tritici (Pgt), such as Ug99, are a major threat to global food security. Here, we provide genomics-based evidence supporting that Ug99 arose by somatic hybridisation and nuclear exchange between dikaryons. Fully haplotype-resolved genome assembly and DNA proximity analysis reveal that Ug99 shares one haploid nucleus genotype with a much older African lineage of Pgt, with no recombination or chromosome reassortment. These findings indicate that nuclear exchange between dikaryotes can generate genetic diversity and facilitate the emergence of new lineages in asexual fungal populations.

De Novo Genome Assembly and Comparative Genomics of the Barley Leaf Rust Pathogen Puccinia hordei Identifies Candidates for Three Avirulence Genes.

Puccinia hordei (Ph) is a damaging pathogen of barley throughout the world. Despite its importance, almost nothing is known about the genomics of this pathogen, and a reference genome is lacking. In this study, the first reference genome was assembled for an Australian isolate of Ph ("Ph560") using long-read SMRT sequencing. A total of 838 contigs were assembled, with a total size of 207 Mbp. This included both haplotype collapsed and separated regions, consistent with an estimated haploid genome size of about 150Mbp. An annotation pipeline that combined RNA-Seq of Ph-infected host tissues and homology to proteins from four other Puccinia species predicted 25,543 gene models of which 1,450 genes were classified as encoding secreted proteins based on the prediction of a signal peptide and no transmembrane domain. Genome resequencing using short-read technology was conducted for four additional Australian strains, Ph612, Ph626, Ph608 and Ph584, which are considered to be simple mutational derivatives of Ph560 with added virulence to one or two of three barley leaf rust resistance genes (viz. Rph3, Rph13 and Rph19). To identify candidate genes for the corresponding avirulence genes AvrRph3, AvrRph13 and AvrRph19, genetic variation in predicted secreted protein genes between the strains was correlated to the virulence profiles of each, identifying 35, 29 and 46 candidates for AvrRph13, AvrRph3 and AvrRph19, respectively. The identification of these candidate genes provides a strong foundation for future efforts to isolate these three avirulence genes, investigate their biological properties, and develop new diagnostic tests for monitoring pathogen virulence.

Components of Brachypodium distachyon resistance to nonadapted wheat stripe rust pathogens are simply inherited.

Phytopathogens have a limited range of host plant species that they can successfully parasitise ie. that they are adapted for. Infection of plants by nonadapted pathogens often results in an active resistance response that is relatively poorly characterised because phenotypic variation in this response often does not exist within a plant species, or is too subtle for genetic dissection. In addition, complex polygenic inheritance often underlies these resistance phenotypes and mutagenesis often does not impact upon this resistance, presumably due to genetic or mechanistic redundancy. Here it is demonstrated that phenotypic differences in the resistance response of Brachypodium distachyon to the nonadapted wheat stripe rust pathogen Puccinia striiformis f. sp. tritici (Pst) are genetically tractable and simply inherited. Two dominant loci were identified on B. distachyon chromosome 4 that each reduce attempted Pst colonisation compared with sib and parent lines without these loci. One locus (Yrr1) is effective against diverse Australian Pst isolates and present in two B. distachyon mapping families as a conserved region that was reduced to 5 candidate genes by fine mapping. A second locus, Yrr2, shows Pst race-specificity and encodes a disease resistance gene family typically associated with host plant resistance. These data indicate that some components of resistance to nonadapted pathogens are genetically tractable in some instances and may mechanistically overlap with host plant resistance to avirulent adapted pathogens.

De Novo Assembly and Phasing of Dikaryotic Genomes from Two Isolates of Puccinia coronata f. sp. avenae, the Causal Agent of Oat Crown Rust.

Oat crown rust, caused by the fungus Pucinnia coronata f. sp. avenae, is a devastating disease that impacts worldwide oat production. For much of its life cycle, P. coronata f. sp. avenae is dikaryotic, with two separate haploid nuclei that may vary in virulence genotype, highlighting the importance of understanding haplotype diversity in this species. We generated highly contiguous de novo genome assemblies of two P. coronata f. sp. avenae isolates, 12SD80 and 12NC29, from long-read sequences. In total, we assembled 603 primary contigs for 12SD80, for a total assembly length of 99.16 Mbp, and 777 primary contigs for 12NC29, for a total length of 105.25 Mbp; approximately 52% of each genome was assembled into alternate haplotypes. This revealed structural variation between haplotypes in each isolate equivalent to more than 2% of the genome size, in addition to about 260,000 and 380,000 heterozygous single-nucleotide polymorphisms in 12SD80 and 12NC29, respectively. Transcript-based annotation identified 26,796 and 28,801 coding sequences for isolates 12SD80 and 12NC29, respectively, including about 7,000 allele pairs in haplotype-phased regions. Furthermore, expression profiling revealed clusters of coexpressed secreted effector candidates, and the majority of orthologous effectors between isolates showed conservation of expression patterns. However, a small subset of orthologs showed divergence in expression, which may contribute to differences in virulence between 12SD80 and 12NC29. This study provides the first haplotype-phased reference genome for a dikaryotic rust fungus as a foundation for future studies into virulence mechanisms in P. coronata f. sp. avenaeIMPORTANCE Disease management strategies for oat crown rust are challenged by the rapid evolution of Puccinia coronata f. sp. avenae, which renders resistance genes in oat varieties ineffective. Despite the economic importance of understanding P. coronata f. sp. avenae, resources to study the molecular mechanisms underpinning pathogenicity and the emergence of new virulence traits are lacking. Such limitations are partly due to the obligate biotrophic lifestyle of P. coronata f. sp. avenae as well as the dikaryotic nature of the genome, features that are also shared with other important rust pathogens. This study reports the first release of a haplotype-phased genome assembly for a dikaryotic fungal species and demonstrates the amenability of using emerging technologies to investigate genetic diversity in populations of P. coronata f. sp. avenae.

Loss of AvrSr50 by somatic exchange in stem rust leads to virulence for Sr50 resistance in wheat.

Race-specific resistance genes protect the global wheat crop from stem rust disease caused by Puccinia graminis f. sp. tritici (Pgt) but are often overcome owing to evolution of new virulent races of the pathogen. To understand virulence evolution in Pgt, we identified the protein ligand (AvrSr50) recognized by the Sr50 resistance protein. A spontaneous mutant of Pgt virulent to Sr50 contained a 2.5 mega-base pair loss-of-heterozygosity event. A haustorial secreted protein from this region triggers Sr50-dependent defense responses in planta and interacts directly with the Sr50 protein. Virulence alleles of AvrSr50 have arisen through DNA insertion and sequence divergence, and our data provide molecular evidence that in addition to sexual recombination, somatic exchange can play a role in the emergence of new virulence traits in Pgt.

Genome analysis and avirulence gene cloning using a high-density RADseq linkage map of the flax rust fungus, Melampsora lini.

BACKGROUND: Rust fungi are an important group of plant pathogens that cause devastating losses in agricultural, silvicultural and natural ecosystems. Plants can be protected from rust disease by resistance genes encoding receptors that trigger a highly effective defence response upon recognition of specific pathogen avirulence proteins. Identifying avirulence genes is crucial for understanding how virulence evolves in the field. RESULTS: To facilitate avirulence gene cloning in the flax rust fungus, Melampsora lini, we constructed a high-density genetic linkage map using single nucleotide polymorphisms detected in restriction site-associated DNA sequencing (RADseq) data. The map comprises 13,412 RADseq markers in 27 linkage groups that together span 5860 cM and contain 2756 recombination bins. The marker sequences were used to anchor 68.9 % of the M. lini genome assembly onto the genetic map. The map and anchored assembly were then used to: 1) show that M. lini has a high overall meiotic recombination rate, but recombination distribution is uneven and large coldspots exist; 2) show that substantial genome rearrangements have occurred in spontaneous loss-of-avirulence mutants; and 3) identify the AvrL2 and AvrM14 avirulence genes by map-based cloning. AvrM14 is a dual-specificity avirulence gene that encodes a predicted nudix hydrolase. AvrL2 is located in the region of the M. lini genome with the lowest recombination rate and encodes a small, highly-charged proline-rich protein. CONCLUSIONS: The M. lini high-density linkage map has greatly advanced our understanding of virulence mechanisms in this pathogen by providing novel insights into genome variability and enabling identification of two new avirulence genes.

Changing the Game: Using Integrative Genomics to Probe Virulence Mechanisms of the Stem Rust Pathogen Puccinia graminis f. sp. tritici.

The recent resurgence of wheat stem rust caused by new virulent races of Puccinia graminis f. sp. tritici (Pgt) poses a threat to food security. These concerns have catalyzed an extensive global effort toward controlling this disease. Substantial research and breeding programs target the identification and introduction of new stem rust resistance (Sr) genes in cultivars for genetic protection against the disease. Such resistance genes typically encode immune receptor proteins that recognize specific components of the pathogen, known as avirulence (Avr) proteins. A significant drawback to deploying cultivars with single Sr genes is that they are often overcome by evolution of the pathogen to escape recognition through alterations in Avr genes. Thus, a key element in achieving durable rust control is the deployment of multiple effective Sr genes in combination, either through conventional breeding or transgenic approaches, to minimize the risk of resistance breakdown. In this situation, evolution of pathogen virulence would require changes in multiple Avr genes in order to bypass recognition. However, choosing the optimal Sr gene combinations to deploy is a challenge that requires detailed knowledge of the pathogen Avr genes with which they interact and the virulence phenotypes of Pgt existing in nature. Identifying specific Avr genes from Pgt will provide screening tools to enhance pathogen virulence monitoring, assess heterozygosity and propensity for mutation in pathogen populations, and confirm individual Sr gene functions in crop varieties carrying multiple effective resistance genes. Toward this goal, much progress has been made in assembling a high quality reference genome sequence for Pgt, as well as a Pan-genome encompassing variation between multiple field isolates with diverse virulence spectra. In turn this has allowed prediction of Pgt effector gene candidates based on known features of Avr genes in other plant pathogens, including the related flax rust fungus. Upregulation of gene expression in haustoria and evidence for diversifying selection are two useful parameters to identify candidate Avr genes. Recently, we have also applied machine learning approaches to agnostically predict candidate effectors. Here, we review progress in stem rust pathogenomics and approaches currently underway to identify Avr genes recognized by wheat Sr genes.

The wheat durable, multipathogen resistance gene Lr34 confers partial blast resistance in rice.

The wheat gene Lr34 confers durable and partial field resistance against the obligate biotrophic, pathogenic rust fungi and powdery mildew in adult wheat plants. The resistant Lr34 allele evolved after wheat domestication through two gain-of-function mutations in an ATP-binding cassette transporter gene. An Lr34-like fungal disease resistance with a similar broad-spectrum specificity and durability has not been described in other cereals. Here, we transformed the resistant Lr34 allele into the japonica rice cultivar Nipponbare. Transgenic rice plants expressing Lr34 showed increased resistance against multiple isolates of the hemibiotrophic pathogen Magnaporthe oryzae, the causal agent of rice blast disease. Host cell invasion during the biotrophic growth phase of rice blast was delayed in Lr34-expressing rice plants, resulting in smaller necrotic lesions on leaves. Lines with Lr34 also developed a typical, senescence-based leaf tip necrosis (LTN) phenotype. Development of LTN during early seedling growth had a negative impact on formation of axillary shoots and spikelets in some transgenic lines. One transgenic line developed LTN only at adult plant stage which was correlated with lower Lr34 expression levels at seedling stage. This line showed normal tiller formation and more importantly, disease resistance in this particular line was not compromised. Interestingly, Lr34 in rice is effective against a hemibiotrophic pathogen with a lifestyle and infection strategy that is different from obligate biotrophic rusts and mildew fungi. Lr34 might therefore be used as a source in rice breeding to improve broad-spectrum disease resistance against the most devastating fungal disease of rice.

SHOEBOX Modulates Root Meristem Size in Rice through Dose-Dependent Effects of Gibberellins on Cell Elongation and Proliferation.

Little is known about how the size of meristem cells is regulated and whether it participates in the control of meristem size in plants. Here, we report our findings on shoebox (shb), a mild gibberellin (GA) deficient rice mutant that has a short root meristem size. Quantitative analysis of cortical cell length and number indicates that shb has shorter, rather than fewer, cells in the root meristem until around the fifth day after sowing, from which the number of cortical cells is also reduced. These defects can be either corrected by exogenous application of bioactive GA or induced in wild-type roots by a dose-dependent inhibitory effect of paclobutrazol on GA biosynthesis, suggesting that GA deficiency is the primary cause of shb mutant phenotypes. SHB encodes an AP2/ERF transcription factor that directly activates transcription of the GA biosynthesis gene KS1. Thus, root meristem size in rice is modulated by SHB-mediated GA biosynthesis that regulates the elongation and proliferation of meristem cells in a developmental stage-specific manner.