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

Association mapping with a diverse population of Puccinia graminis f. sp. tritici identified avirulence loci interacting with the barley Rpg1 stem rust resistance gene

BACKGROUND

Wheat stem rust, caused by Puccinia graminis f. sp. tritici (Pgt), is an important disease of barley and wheat. A diverse sexual Pgt population from the Pacific Northwest (PNW) region of the US contains a high proportion of individuals with virulence on the barley stem rust resistance (R) gene, Rpg1. However, the evolutionary mechanisms of this virulence on Rpg1 are mysterious considering that Rpg1 had not been deployed in the region and the gene had remained remarkably durable in the Midwestern US and prairie provinces of Canada.

METHODS AND RESULTS

To identify AvrRpg1 effectors, genome wide association studies (GWAS) were performed using 113 Pgt isolates collected from the PNW (n = 89 isolates) and Midwest (n = 24 isolates) regions of the US. Disease phenotype data were generated on two barley lines Morex and the Golden Promise transgenic (H228.2c) that carry the Rpg1 gene. Genotype data was generated by whole genome sequencing (WGS) of 96 isolates (PNW = 89 isolates and Midwest = 7 isolates) and RNA sequencing (RNAseq) data from 17 Midwestern isolates. Utilizing ~1.2 million SNPs generated from WGS and phenotype data (n = 96 isolates) on the transgenic line H228.2c, 53 marker trait associations (MTAs) were identified. Utilizing ~140 K common SNPs generated from combined analysis of WGS and RNAseq data, two significant MTAs were identified using the cv Morex phenotyping data. The 55 MTAs defined two distinct avirulence loci, on supercontig 2.30 and supercontig 2.11 of the Pgt reference genome of Pgt isolate CRL 75-36-700-3. The major avirulence locus designated AvrRpg1A was identified with the GWAS using both barley lines and was delimited to a 35 kb interval on supercontig 2.30 containing four candidate genes (PGTG_10878, PGTG_10884, PGTG_10885, and PGTG_10886). The minor avirulence locus designated AvrRpg1B identified with cv Morex contained a single candidate gene (PGTG_05433). AvrRpg1A haplotype analysis provided strong evidence that a dominant avirulence gene underlies the locus.

CONCLUSIONS

The association analysis identified strong candidate AvrRpg1 genes. Further analysis to validate the AvrRpg1 genes will fill knowledge gaps in our understanding of rust effector biology and the evolution and mechanism/s of Pgt virulence on Rpg1.

Identification of Candidate Avirulence and Virulence Genes Corresponding to Stem Rust (Puccinia graminis f. sp. tritici) Resistance Genes in Wheat

Stem rust, caused by the biotrophic fungal pathogen Puccinia graminis f. sp. tritici (Pgt), is an important disease of wheat. However, the majority of Pgt virulence/avirulence loci and underlying genes remain uncharacterized due to the constraints of developing bi-parental populations with this obligate biotroph. Genome-wide association studies (GWAS) using a sexual Pgt population mainly collected from the Pacific Northwestern United States were used to identify candidate virulence/avirulence effector genes corresponding to the six wheat Sr genes: Sr5, Sr21, Sr8a, Sr17, Sr9a, and Sr9d. The Pgt isolates were genotyped using whole-genome shotgun sequencing that identified approximately 1.2 million single nucleotide polymorphisms (SNPs) and were phenotyped at the seedling stage on six Sr gene differential lines. Association mapping analyses identified 17 Pgt loci associated with virulence or avirulence phenotypes on six Pgt resistance genes. Among these loci, 16 interacted with a specific Sr gene, indicating Sr-gene specific interactions. However, one avirulence locus interacted with two separate Sr genes (Sr9a and Sr17), suggesting two distinct Sr genes identifying a single avirulence effector. A total of 24 unique effector gene candidates were identified, and haplotype analysis suggests that within this population, AvrSr5, AvrSr21, AvrSr8a, AvrSr17, and AvrSr9a are dominant avirulence genes, while avrSr9d is a dominant virulence gene. The putative effector genes will be fundamental for future effector gene cloning efforts, allowing for further understanding of rust effector biology and the mechanisms underlying virulence evolution in Pgt with respect to race-specific R-genes. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genome-wide atlas of rust resistance loci in wheat

Rust diseases, including leaf rust, stripe/yellow rust, and stem rust, significantly impact wheat (Triticum aestivum L.) yields, causing substantial economic losses every year. Breeding and deployment of cultivars with genetic resistance is the most effective and sustainable approach to control these diseases. The genetic toolkit for wheat breeders to select for rust resistance has rapidly expanded with a multitude of genetic loci identified using the latest advances in genomics, mapping and cloning strategies. The goal of this review was to establish a wheat genome atlas that provides a comprehensive summary of reported loci associated with rust resistance. Our atlas provides a summary of mapped quantitative trait loci (QTL) and characterised genes for the three rusts from 170 publications over the past two decades. A total of 920 QTL or resistance genes were positioned across the 21 chromosomes of wheat based on the latest wheat reference genome (IWGSC RefSeq v2.1). Interestingly, 26 genomic regions contained multiple rust loci suggesting they could have pleiotropic effects on two or more rust diseases. We discuss a range of strategies to exploit this wealth of genetic information to efficiently utilise sources of resistance, including genomic information to stack desirable and multiple QTL to develop wheat cultivars with enhanced resistance to rust disease.

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

AvrSr27 is a zinc-bound effector with a modular structure important for immune recognition

Effector proteins are central to the success of plant pathogens, while immunity in host plants is driven by receptor-mediated recognition of these effectors. Understanding the molecular details of effector-receptor interactions is key for the engineering of novel immune receptors. Here, we experimentally determined the crystal structure of the Puccinia graminis f. sp. tritici (Pgt) effector AvrSr27, which was not accurately predicted using AlphaFold2. We characterised the role of the conserved cysteine residues in AvrSr27 using in vitro biochemical assays and examined Sr27-mediated recognition using transient expression in Nicotiana spp. and wheat protoplasts. The AvrSr27 structure contains a novel β-strand rich modular fold consisting of two structurally similar domains that bind to Zn2+ ions. The N-terminal domain of AvrSr27 is sufficient for interaction with Sr27 and triggering cell death. We identified two Pgt proteins structurally related to AvrSr27 but with low sequence identity that can also associate with Sr27, albeit more weakly. Though only the full-length proteins, trigger Sr27-dependent cell death in transient expression systems. Collectively, our findings have important implications for utilising protein prediction platforms for effector proteins, and those embarking on bespoke engineering of immunity receptors as solutions to plant disease.

A fully haplotype-resolved and nearly gap-free genome assembly of wheat stripe rust fungus

Stripe rust fungus Puccinia striiformis f. sp. tritici (Pst) is a destructive pathogen of wheat worldwide. Pst has a macrocyclic-heteroecious lifecycle, in which one-celled urediniospores are dikaryotic, each nucleus containing one haploid genome. We successfully generated the first fully haplotype-resolved and nearly gap-free chromosome-scale genome assembly of Pst by combining PacBio HiFi sequencing and trio-binning strategy. The genome size of the two haploid assemblies was 75.59 Mb and 75.91 Mb with contig N50 of 4.17 Mb and 4.60 Mb, and both had 18 pseudochromosomes. The high consensus quality values of 55.57 and 59.02 for both haplotypes confirmed the correctness of the assembly. Of the total 18 chromosomes, 15 and 16 were gapless while there were only five and two gaps for the remaining chromosomes of the two haplotypes, respectively. In total, 15,046 and 15,050 protein-coding genes were predicted for the two haplotypes, and the complete BUSCO scores achieved 97.7% and 97.9%, respectively. The genome will lay the foundation for further research on genetic variations and the evolution of rust fungi.

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.

Recent advances in the biosynthesis of fungal glucan structural diversity

Glucans are the most abundant class of macromolecule polymers in fungi, which are commonly found in Ascomycota and Basidiomycota. Fungal glucans are not only essential for cell integrity and function but also crucial for the immense industrial interest in high value applications. They present a variety of structural characteristics at the nanoscale due to the high regulation of genes and the involvement of stochastic processes in synthesis. However, although recent findings have demonstrated the genes of glucans synthesis are relatively conserved across diverse fungi, the formation and organization of diverse glucan structures is still unclear in fungi. Here, we summarize the structural features of fungal glucans and the recent developments in the mechanisms of glucans biosynthesis. Furthermore, we propose the engineering strategies of targeted glucan synthesis and point out the remaining challenges in the synthetic process. Understanding the synthesis process of diverse glucans is necessary for tailoring high value glucan towards specific applications. This engineering strategy contributes to enable the sustainable and efficient production of glucan diversity.

Comparative Analysis of the Avirulence Effectors Produced by the Fungal Stem Rust Pathogen of Wheat

Crops are constantly exposed to pathogenic microbes. Rust fungi are examples of these harmful microorganisms, which have a major economic impact on wheat production. To protect themselves from pathogens like rust fungi, plants employ a multilayered immune system that includes immunoreceptors encoded by resistance genes. Significant efforts have led to the isolation of numerous resistance genes against rust fungi in cereals, especially in wheat. However, the evolution of virulence of rust fungi hinders the durability of resistance genes as a strategy for crop protection. Rust fungi, like other biotrophic pathogens, secrete an arsenal of effectors to facilitate infection, and these are the molecules that plant immunoreceptors target for pathogen recognition and mounting defense responses. When recognized, these effector proteins are referred to as avirulence (Avr) effectors. Despite the many predicted effectors in wheat rust fungi, only five Avr genes have been identified, all from wheat stem rust. Knowledge of the Avr genes and their variation in the fungal population will inform deployment of the most appropriate wheat disease-resistance genes for breeding and farming. The review provides an overview of methodologies as well as the validation techniques that have been used to characterize Avr effectors from wheat stem rust. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Genome-Enabled Analysis of Population Dynamics and Virulence-Associated Loci in the Oat Crown Rust Fungus Puccinia coronata f. sp. avenae

Puccinia coronata f. sp. avenae (Pca) is an important fungal pathogen causing crown rust that impacts oat production worldwide. Genetic resistance for crop protection against Pca is often overcome by the rapid virulence evolution of the pathogen. This study investigated the factors shaping adaptive evolution of Pca using pathogen populations from distinct geographic regions within the United States and South Africa. Phenotypic and genome-wide sequencing data of these diverse Pca collections, including 217 isolates, uncovered phylogenetic relationships and established distinct genetic composition between populations from northern and southern regions from the United States and South Africa. The population dynamics of Pca involve a bidirectional movement of inoculum between northern and southern regions of the United States and contributions from clonality and sexuality. The population from South Africa is solely clonal. A genome-wide association study (GWAS) employing a haplotype-resolved Pca reference genome was used to define 11 virulence-associated loci corresponding to 25 oat differential lines. These regions were screened to determine candidate Avr effector genes. Overall, the GWAS results allowed us to identify the underlying genetic factors controlling pathogen recognition in an oat differential set used in the United States to assign pathogen races (pathotypes). Key GWAS findings support complex genetic interactions in several oat lines, suggesting allelism among resistance genes or redundancy of genes included in the differential set, multiple resistance genes recognizing genetically linked Avr effector genes, or potentially epistatic relationships. A careful evaluation of the composition of the oat differential set accompanied by the development or implementation of molecular markers is recommended. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Soybean-Phakopsora pachyrhizi interactions: towards the development of next-generation disease-resistant plants

Soybean rust (SBR), caused by the obligate biotrophic fungus Phakopsora pachyrhizi, is a devastating foliar disease threatening soybean production. To date, no commercial cultivars conferring durable resistance to SBR are available. The development of long-lasting SBR resistance has been hindered by the lack of understanding of this complex pathosystem, encompassing challenges posed by intricate genetic structures in both the host and pathogen, leading to a gap in the knowledge of gene-for-gene interactions between soybean and P. pachyrhizi. In this review, we focus on recent advancements and emerging technologies that can be used to improve our understanding of the P. pachyrhizi-soybean molecular interactions. We further explore approaches used to combat SBR, including conventional breeding, transgenic approaches and RNA interference, and how advances in our understanding of plant immune networks, the availability of new molecular tools, and the recent sequencing of the P. pachyrhizi genome could be used to aid in the development of better genetic resistance against SBR. Lastly, we discuss the research gaps of this pathosystem and how new technologies can be used to shed light on these questions and to develop durable next-generation SBR-resistant soybean plants.

Virulence Dynamics of the Barley Leaf Rust Pathogen (Puccinia hordei) in the United States from 1989 to 2020

Barley leaf rust, caused by Puccinia hordei, is an important disease of barley worldwide. The pathogen can develop new races that overcome resistance genes, emphasizing the need for monitoring its virulence. This study characterized 519 P. hordei isolates collected in the United States from the 1989 to 2000 and 2010 to 2020 survey periods on 15 Rph (Reaction to Puccinia hordei) genes. We analyzed linearized infection type data to detect virulence patterns across the United States and in five geographical regions: Pacific/West (PW), Southwest (SW), Midwest (MW), Northeast (NE), and Southeast (SE). Over 32 years, we observed high mean infection scores for Rph1.a, Rph4.d, and Rph8.h; intermediate scores for Rph2.b, Rph9.i, Rph10.o, Rph11.p, and Rph13.x; and low scores for Rph3.c, Rph5.e, Rph5.f, Rph7.g, Rph9.z, Rph14.ab, and Rph15.ad. Virulence for Rph2.b, Rph3.c, Rph5.e, Rph9.z, Rph10.o, Rph11.p, and Rph13.x significantly differed between the two survey periods. From 1989 to 2020, regional patterns of virulence were found for Rph5.e, Rph5.f, Rph7.g, and Rph14.ab, while regionalities of virulence for Rph3.c, Rph9.i, Rph9.z were only observed in the 2010 to 2020 survey period. Virulence associations were also detected in the P. hordei population. Notably, isolates that were virulent to Rph5.e and Rph6.f were more likely to be avirulent to Rph7.g and Rph13.x, and vice versa. In decreasing order of effectiveness, Rph15.ad, Rph5.e, Rph3.c, Rph9.z, Rph7.g, Rph5.f, and Rph14.ab were the most effective Rph genes in the United States from 1989 to 2020. Pyramiding Rph15.ad with other widely effective Rph and adult plant resistance genes may provide long-lasting resistance against P. hordei.

Recent co-evolution of two pandemic plant diseases in a multi-hybrid swarm

Most plant pathogens exhibit host specificity but when former barriers to infection break down, new diseases can rapidly emerge. For a number of fungal diseases, there is increasing evidence that hybridization plays a major role in driving host jumps. However, the relative contributions of existing variation versus new mutations in adapting to new host(s) is unclear. Here we reconstruct the evolutionary history of two recently emerged populations of the fungus Pyricularia oryzae that are responsible for two new plant diseases: wheat blast and grey leaf spot of ryegrasses. We provide evidence that wheat blast/grey leaf spot evolved through two distinct mating episodes: the first occurred ~60 years ago, when a fungal individual adapted to Eleusine mated with another individual from Urochloa. Then, about 10 years later, a single progeny from this cross underwent a series of matings with a small number of individuals from three additional host-specialized populations. These matings introduced non-functional alleles of two key host-specificity factors, whose recombination in a multi-hybrid swarm probably facilitated the host jump. We show that very few mutations have arisen since the founding event and a majority are private to individual isolates. Thus, adaptation to the wheat or Lolium hosts appears to have been instantaneous, and driven entirely by selection on repartitioned standing variation, with no obvious role for newly formed mutations.

Arbuscular mycorrhizal fungi heterokaryons have two nuclear populations with distinct roles in host-plant interactions

Arbuscular mycorrhizal fungi (AMF) are prominent root symbionts that can carry thousands of nuclei deriving from two parental strains in a large syncytium. These co-existing genomes can also vary in abundance with changing environmental conditions. Here we assemble the nuclear genomes of all four publicly available AMF heterokaryons using PacBio high-fidelity and Hi-C sequencing. We find that the two co-existing genomes of these strains are phylogenetically related but differ in structure, content and epigenetics. We confirm that AMF heterokaryon genomes vary in relative abundance across conditions and show this can lead to nucleus-specific differences in expression during interactions with plants. Population analyses also reveal signatures of genetic exchange indicative of past events of sexual reproduction in these strains. This work uncovers the origin and contribution of two nuclear genomes in AMF heterokaryons and opens avenues for the improvement and environmental application of these strains.

Nuclear exchange generates population diversity in the wheat leaf rust pathogen Puccinia triticina

In clonally reproducing dikaryotic rust fungi, non-sexual processes such as somatic nuclear exchange are postulated to play a role in diversity but have been difficult to detect due to the lack of genome resolution between the two haploid nuclei. We examined three nuclear-phased genome assemblies of Puccinia triticina, which causes wheat leaf rust disease. We found that the most recently emerged Australian lineage was derived by nuclear exchange between two pre-existing lineages, which originated in Europe and North America. Haplotype-specific phylogenetic analysis reveals that repeated somatic exchange events have shuffled haploid nuclei between long-term clonal lineages, leading to a global P. triticina population representing different combinations of a limited number of haploid genomes. Thus, nuclear exchange seems to be the predominant mechanism generating diversity and the emergence of new strains in this otherwise clonal pathogen. Such genomics-accelerated surveillance of pathogen evolution paves the way for more accurate global disease monitoring.

Uncovering the history of recombination and population structure in western Canadian stripe rust populations through mating type alleles

BACKGROUND

The population structure of crop pathogens such as Puccinia striiformis f. sp. tritici (Pst), the cause of wheat stripe rust, is of interest to researchers looking to understand these pathogens on a molecular level as well as those with an applied focus such as disease epidemiology. Cereal rusts can reproduce sexually or asexually, and the emergence of novel lineages has the potential to cause serious epidemics such as the one caused by the 'Warrior' lineage in Europe. In a global context, Pst lineages in Canada were not well-characterized and the origin of foreign incursions was not known. Additionally, while some Pst mating type genes have been identified in published genomes, there has been no rigorous assessment of mating type diversity and distribution across the species.

RESULTS

We used a whole-genome/transcriptome sequencing approach for the Canadian Pst population to identify lineages in their global context and evidence tracing foreign incursions. More importantly: for the first time ever, we identified nine alleles of the homeodomain mating type locus in the worldwide Pst population and show that previously identified lineages exhibit a single pair of these alleles. Consistently with the literature, we find only two pheromone receptor mating type alleles. We show that the recent population shift from the 'PstS1' lineage to the 'PstS1-related' lineage is also associated with the introduction of a novel mating type allele (Pst-b3-HD) to the Canadian population. We also show evidence for high levels of mating type diversity in samples associated with the Himalayan center of diversity for Pst, including a single Canadian race previously identified as 'PstPr' (probable recombinant) which we identify as a foreign incursion, most closely related to isolates sampled from China circa 2015.

CONCLUSIONS

These data describe a recent shift in the population of Canadian Pst field isolates and characterize homeodomain-locus mating type alleles in the global Pst population which can now be utilized in testing several research questions and hypotheses around sexuality and hybridization in rust fungi.

Resistance that stacks up: engineering rust and mildew disease control in the cereal crops wheat and barley

Staying ahead of the arms race against rust and mildew diseases in cereal crops is essential to maintain and preserve food security. The methodological challenges associated with conventional resistance breeding are major bottlenecks for deploying resistance (R) genes in high-yielding crop varieties. Advancements in our knowledge of plant genomes, structural mechanisms, innovations in bioinformatics, and improved plant transformation techniques have alleviated this bottleneck by permitting rapid gene isolation, functional studies, directed engineering of synthetic resistance and precise genome manipulation in elite crop cultivars. Most cloned cereal R genes encode canonical immune receptors which, on their own, are prone to being overcome through selection for resistance-evading pathogenic strains. However, the increasingly large repertoire of cloned R genes permits multi-gene stacking that, in principle, should provide longer-lasting resistance. This review discusses how these genomics-enabled developments are leading to new breeding and biotechnological opportunities to achieve durable rust and powdery mildew control in cereals.

Plant and pathogen genomics: essential approaches for stem rust resistance gene stacks in wheat

The deployment of disease resistance genes is currently the most economical and environmentally sustainable method of crop protection. However, disease resistance genes can rapidly break down because of constant pathogen evolution, particularly when they are deployed singularly. Polygenic resistance is, therefore, considered the most durable, but combining and maintaining these genes by breeding is a laborious process as effective genes are usually unlinked. The deployment of polygenic resistance with single-locus inheritance is a promising innovation that overcomes these difficulties while enhancing resistance durability. Because of major advances in genomic technologies, increasing numbers of plant resistance genes have been cloned, enabling the development of resistance transgene stacks (RTGSs) that encode multiple genes all located at a single genetic locus. Gene stacks encoding five stem rust resistance genes have now been developed in transgenic wheat and offer both breeding simplicity and potential resistance durability. The development of similar genomic resources in phytopathogens has advanced effector gene isolation and, in some instances, enabled functional validation of individual resistance genes in RTGS. Here, the wheat stem rust pathosystem is used as an illustrative example of how host and pathogen genomic advances have been instrumental in the development of RTGS, which is a strategy applicable to many other agricultural crop species.

From Gene-for-Gene to Resistosomes: Flor's Enduring Legacy

The gene-for-gene model proposed by H. H. Flor has been one of the fundamental precepts of plant-pathogen interactions that has underpinned decades of research towards our current concepts of plant immunity. The broad validity of this model as an elegant and accurate genetic description of specific recognition events between the products of plant resistance (R) and pathogen avirulence (Avr) genes has been demonstrated many times over in a wide variety of plant disease systems. In recent years detailed molecular and structural analyses have provided a deep understanding of the principles by which plant immune receptors recognize pathogen effectors, including providing molecular descriptions of many of the genetic loci in flax and flax rust characterized by Flor. Recent advances in molecular and structural understanding of immune receptor recognition and activation mechanisms have brought the field to a new level, where rational design of novel receptors through engineering approaches is becoming a realizable goal. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Insights into Diversity, Distribution, and Systematics of Rust Genus Puccinia

Puccinia, which comprises 4000 species, is the largest genus of rust fungi and one of the destructive plant pathogenic rust genera that are reported to infect both agricultural and nonagricultural plants with severe illnesses. The presence of bi-celled teliospores is one of the major features of these rust fungi that differentiated them from Uromyces, which is another largest genus of rust fungi. In the present study, an overview of the current knowledge on the general taxonomy and ecology of the rust genus Puccinia is presented. The status of the molecular identification of this genus along with updated species numbers and their current statuses in the 21st century are also presented, in addition to their threats to both agricultural and nonagricultural plants. Furthermore, a phylogenetic analysis based on ITS and LSU DNA sequence data available in GenBank and the published literature was performed to examine the intergeneric relationships of Puccinia. The obtained results revealed the worldwide distribution of Puccinia. Compared with other nations, a reasonable increase in research publications over the current century was demonstrated in Asian countries. The plant families Asteraceae and Poaceae were observed as the most infected in the 21st century. The phylogenetic studies of the LSU and ITS sequence data revealed the polyphyletic nature of Puccinia. In addition, the presences of too short, too lengthy, and incomplete sequences in the NCBI database demonstrate the need for extensive DNA-based analyses for a better understanding of the taxonomic placement of Puccinia.

Major proliferation of transposable elements shaped the genome of the soybean rust pathogen Phakopsora pachyrhizi

With >7000 species the order of rust fungi has a disproportionately large impact on agriculture, horticulture, forestry and foreign ecosystems. The infectious spores are typically dikaryotic, a feature unique to fungi in which two haploid nuclei reside in the same cell. A key example is Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease, one of the world's most economically damaging agricultural diseases. Despite P. pachyrhizi's impact, the exceptional size and complexity of its genome prevented generation of an accurate genome assembly. Here, we sequence three independent P. pachyrhizi genomes and uncover a genome up to 1.25 Gb comprising two haplotypes with a transposable element (TE) content of ~93%. We study the incursion and dominant impact of these TEs on the genome and show how they have a key impact on various processes such as host range adaptation, stress responses and genetic plasticity.

Gapless Genome Assembly of Puccinia triticina Provides Insights into Chromosome Evolution in Pucciniales

Chromosome evolution drives species evolution, speciation, and adaptive radiation. Accurate genome assembly is crucial to understanding chromosome evolution of species, such as dikaryotic fungi. Rust fungi (Pucciniales) in dikaryons represent the largest group of plant pathogens, but the evolutionary process of adaptive radiation in Pucciniales remains poorly understood. Here, we report a gapless genome for the wheat leaf rust fungus Puccinia triticina determined using PacBio high-fidelity (HiFi) sequencing. This gapless assembly contains two sets of chromosomes, showing that one contig represents one chromosome. Comparisons of homologous chromosomes between the phased haplotypes revealed that highly frequent small-scale sequence divergence shapes haplotypic variation. Genome analyses of Puccinia triticina along with other rusts revealed that recent transposable element bursts and extensive segmental gene duplications synergistically highlight the evolution of chromosome structures. Comparative analysis of chromosomes indicated that frequent chromosomal rearrangements may act as a major contributor to rapid radiation of Pucciniales. This study presents the first gapless, phased assembly for a dikaryotic rust fungus and provides insights into adaptive evolution and species radiation in Pucciniales. IMPORTANCE Rust fungi (Pucciniales) are the largest group of plant pathogens. Adaptive radiation is a predominant feature in Pucciniales evolution. Chromosome evolution plays an important role in adaptive evolution. Accurate chromosome-scale assembly is required to understand the role of chromosome evolution in Pucciniales. We took advantage of HiFi sequencing to construct a gapless, phased genome for Puccinia triticina. Further analyses revealed that the evolution of chromosome structures in rust lineage is shaped by the combination of transposable element bursts and segmental gene duplications. Chromosome comparisons of Puccinia triticina and other rusts suggested that frequent chromosomal arrangements may make remarkable contributions to high species diversity of rust fungi. Our results present the first gapless genome for Pucciniales and shed light on the feature of chromosome evolution in Pucciniales.

The broad use of the Pm8 resistance gene in wheat resulted in hypermutation of the AvrPm8 gene in the powdery mildew pathogen

BACKGROUND

Worldwide wheat production is under constant threat by fast-evolving fungal pathogens. In the last decades, wheat breeding for disease resistance heavily relied on the introgression of chromosomal segments from related species as genetic sources of new resistance. The Pm8 resistance gene against the powdery mildew disease has been introgressed from rye into wheat as part of a large 1BL.1RS chromosomal translocation encompassing multiple disease resistance genes and yield components. Due to its high agronomic value, this translocation has seen continuous global use since the 1960s on large growth areas, even after Pm8 resistance was overcome by the powdery mildew pathogen. The long-term use of Pm8 at a global scale provided the unique opportunity to study the consequences of such extensive resistance gene application on pathogen evolution.

RESULTS

Using genome-wide association studies in a population of wheat mildew isolates, we identified the avirulence effector AvrPm8 specifically recognized by Pm8. Haplovariant mining in a global mildew population covering all major wheat growing areas of the world revealed 17 virulent haplotypes of the AvrPm8 gene that grouped into two functional categories. The first one comprised amino acid polymorphisms at a single position along the AvrPm8 protein, which we confirmed to be crucial for the recognition by Pm8. The second category consisted of numerous destructive mutations to the AvrPm8 open reading frame such as disruptions of the start codon, gene truncations, gene deletions, and interference with mRNA splicing. With the exception of a single, likely ancient, gain-of-virulence mutation found in mildew isolates around the world, all AvrPm8 virulence haplotypes were found in geographically restricted regions, indicating that they occurred recently as a consequence of the frequent Pm8 use.

CONCLUSIONS

In this study, we show that the broad and prolonged use of the Pm8 gene in wheat production worldwide resulted in a multitude of gain-of-virulence mechanisms affecting the AvrPm8 gene in the wheat powdery mildew pathogen. Based on our findings, we conclude that both standing genetic variation as well as locally occurring new mutations contributed to the global breakdown of the Pm8 resistance gene introgression.

Global dispersal and diversity of rust fungi in the context of plant health

Long-distance dispersal of plant pathogens at the continental scale may have strong implications on plant health, in particular when incursions result in spread of disease to new territories where the disease was previously absent or insignificant. These dispersions may be caused by airborne transmission of spores or accidental spread via human travel and trade. Recent surveillance efforts of cereal rust fungi have demonstrated that incursion of new strains with superior fitness into areas where the disease is already established may have similar implications on plant health. Since dispersal events are highly stochastic, irrespective of transmission mechanism, critical mitigation efforts include preparedness by coordinated pathogen surveillance activities, host crop diversification, and breeding for disease resistance with low vulnerability to sudden changes in the pathogen population.

Prediction of effector protein structures from fungal phytopathogens enables evolutionary analyses

Elucidating the similarity and diversity of pathogen effectors is critical to understand their evolution across fungal phytopathogens. However, rapid divergence that diminishes sequence similarities between putatively homologous effectors has largely concealed the roots of effector evolution. Here we modelled the structures of 26,653 secreted proteins from 14 agriculturally important fungal phytopathogens, six non-pathogenic fungi and one oomycete with AlphaFold 2. With 18,000 successfully predicted folds, we performed structure-guided comparative analyses on two aspects of effector evolution: uniquely expanded sequence-unrelated structurally similar (SUSS) effector families and common folds present across the fungal species. Extreme expansion of lineage-specific SUSS effector families was found only in several obligate biotrophs, Blumeria graminis and Puccinia graminis. The highly expanded effector families were the source of conserved sequence motifs, such as the Y/F/WxC motif. We identified new classes of SUSS effector families that include known virulence factors, such as AvrSr35, AvrSr50 and Tin2. Structural comparisons revealed that the expanded structural folds further diversify through domain duplications and fusion with disordered stretches. Putatively sub- and neo-functionalized SUSS effectors could reconverge on regulation, expanding the functional pools of effectors in the pathogen infection cycle. We also found evidence that many effector families could have originated from ancestral folds conserved across fungi. Collectively, our study highlights diverse effector evolution mechanisms and supports divergent evolution as a major force in driving SUSS effector evolution from ancestral proteins.

Folding Features and Dynamics of 3D Genome Architecture in Plant Fungal Pathogens

The folding and dynamics of three-dimensional (3D) genome organization are fundamental for eukaryotes executing genome functions but have been largely unexplored in nonmodel fungi. Using high-throughput sequencing coupled with chromosome conformation capture (Hi-C) data, we generated two chromosome-level assemblies for Puccinia striiformis f. sp. tritici, a fungus causing stripe rust disease on wheat, for studying 3D genome architectures of plant pathogenic fungi. The chromatin organization of the fungus followed a combination of the fractal globule model and the equilibrium globule model. Surprisingly, chromosome compartmentalization was not detected. Dynamics of 3D genome organization during two developmental stages of P. striiformis f. sp. tritici indicated that regulation of gene activities might be independent of the changes of genome organization. In addition, chromatin conformation conservation was found to be independent of genome sequence synteny conservation among different fungi. These results highlighted the distinct folding principles of fungal 3D genomes. Our findings should be an important step toward a holistic understanding of the principles and functions of genome architecture across different eukaryotic kingdoms. IMPORTANCE Previously, our understanding of 3D genome architecture has mainly come from model mammals, insects, and plants. However, the organization and regulatory functions of 3D genomes in fungi are largely unknown. In this study, we comprehensively investigated P. striiformis f. sp. tritici, a plant fungal pathogen, and revealed distinct features of the 3D genome, comparing it with the universal folding feature of 3D genomes in higher eukaryotic organisms. We further suggested that there might be distinct regulatory mechanisms of gene expression that are independent of chromatin organization changes during the developmental stages of this rust fungus. Moreover, we showed that the evolutionary pattern of 3D genomes in this fungus is also different from the cases in mammalian genomes. In addition, the genome assembly pipeline and the generated two chromosome-level genomes will be valuable resources. These results highlighted the unexplored distinct features of 3D genome organization in fungi. Therefore, our study provided complementary knowledge to holistically understand the organization and functions of 3D genomes across different eukaryotes.

TRustDB: A comprehensive bioinformatics resource for understanding the complete Wheat-Stem rust host-pathogen interactome

The increasing infectious diseases in wheat immensely reduce crop yield and quality, thus affecting global wheat production. The evolution in phytopathogens hinders the understanding of the disease infection mechanisms. TRustDB is an open-access, comprehensive database that is specifically focused on the disease stem rust (also known as black rust) in Triticum aestivum, which is caused by the fungal pathogen Puccinia graminis (Pgt), strains 'Ug99' and '21-0'. The database aims at a broader focus of providing the researchers with comprehensive tools to predict the protein-protein interactions and avail the functional annotations of the proteins involved in the interactions that cause the disease. The network of the predicted interactome can also be visualized on the browser. Various modules for the functional annotations of the host and pathogen proteins such as subcellular localization, functional domains, gene ontology annotations, pathogen orthologs and effector proteins have been implemented. The host proteins that serve as transcription factors, along with the respective Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways are also available, which further enhance the understanding of the disease infection mechanisms and the defense responses of the host. The database is also linked with several other databases such as InterPro, KEGG pathways, Ensembl and National Center for Biotechnology Information (NCBI). TRustDB has a user-friendly web interface, which can be accessed through . Database URL http://bioinfo.usu.edu/trustdb/.

Changes of Barley Stripe Rust Populations in the United States from 1993 to 2017

Barley stripe rust is a relatively new disease in the United States. The pathogen, Puccinia striiformis f. sp. hordei (Psh), was first observed in Texas in 1991 and has spread north and westwards and mainly caused epidemics in the western United States. A total of 447 isolates collected from 1993 to 2017 were identified as 382 multilocus genotypes (MLGs) using 14 simple sequence repeat markers. The MLGs were clustered into six molecular groups (MGs) using the discriminant analysis of principal components and the hierarchical cluster analysis, and the MGs had significant differences in frequency in different years. MG1 was present in the population prior to the year 2000. MG2, MG3, and MG4 became predominate after 2000. MG5 was detected in all 24 years but more frequent from 2010 to 2017. MG6 was the most recent group detected mainly from 2011 to 2017 and had the highest correlation coefficient with the virulence phenotypes among the MGs. The heterozygosity and genotypic diversity of the Psh populations increased from 2000 to 2017, even more from 2010 to 2017. The results indicate rapid genetic changes from year to year, with major molecular group changes around 2000 and 2010. The possible mechanisms underlying the population changes are discussed.

Population genomics of Puccinia graminis f.sp. tritici highlights the role of admixture in the origin of virulent wheat rust races

Puccinia graminis f.sp. tritici (Pgt) causes stem rust disease in wheat that can result in severe yield losses. The factors driving the evolution of its virulence and adaptation remain poorly characterized. We utilize long-read sequencing to develop a haplotype-resolved genome assembly of a U.S. isolate of Pgt. Using Pgt haplotypes as a reference, we characterize the structural variants (SVs) and single nucleotide polymorphisms in a diverse panel of isolates. SVs impact the repertoire of predicted effectors, secreted proteins involved in host-pathogen interaction, and show evidence of purifying selection. By analyzing global and local genomic ancestry we demonstrate that the origin of 8 out of 12 Pgt clades is linked with either somatic hybridization or sexual recombination between the diverged donor populations. Our study shows that SVs and admixture events appear to play an important role in broadening Pgt virulence and the origin of highly virulent races, creating a resource for studying the evolution of Pgt virulence and preventing future epidemic outbreaks.

Concerted expansion and contraction of immune receptor gene repertoires in plant genomes

Recent reports suggest that cell-surface and intracellular immune receptors function synergistically to activate robust defence against pathogens, but whether they co-evolve is unclear. Here we determined the numbers of cell-surface and intracellular immune receptors in 350 species. Surprisingly, the number of receptor genes that are predicted to encode cell-surface and intracellular immune receptors is strongly correlated. We suggest this is consistent with mutual potentiation of immunity initiated by cell-surface and intracellular receptors being reflected in the concerted co-evolution of the size of their repertoires across plant species.

Applications of In Vitro Tissue Culture Technologies in Breeding and Genetic Improvement of Wheat

Sources of new genetic variability have been limited to existing germplasm in the past. Wheat has been studied extensively for various agronomic traits located throughout the genome. The large size of the chromosomes and the ability of its polyploid genome to tolerate the addition or loss of chromosomes facilitated rapid progress in the early study of wheat genetics using cytogenetic techniques. At the same time, its large genome size has limited the progress in genetic characterization studies focused on diploid species, with a small genome and genetic engineering procedures already developed. Today, the genetic transformation and gene editing procedures offer attractive alternatives to conventional techniques for breeding wheat because they allow one or more of the genes to be introduced or altered into an elite cultivar without affecting its genetic background. Recently, significant advances have been made in regenerating various plant tissues, providing the essential basis for regenerating transgenic plants. In addition, Agrobacterium-mediated, biolistic, and in planta particle bombardment (iPB) gene delivery procedures have been developed for wheat transformation and advanced transgenic wheat development. As a result, several useful genes are now available that have been transferred or would be helpful to be transferred to wheat in addition to the current traditional effort to improve trait values, such as resistance to abiotic and biotic factors, grain quality, and plant architecture. Furthermore, the in planta genome editing method will significantly contribute to the social implementation of genome-edited crops to innovate the breeding pipeline and leverage unique climate adaptations.

Current Status and Future Perspectives of Genomics Research in the Rust Fungi

Rust fungi in Pucciniales have caused destructive plant epidemics, have become more aggressive with new virulence, rapidly adapt to new environments, and continually threaten global agriculture. With the rapid advancement of genome sequencing technologies and data analysis tools, genomics research on many of the devastating rust fungi has generated unprecedented insights into various aspects of rust biology. In this review, we first present a summary of the main findings in the genomics of rust fungi related to variations in genome size and gene composition between and within species. Then we show how the genomics of rust fungi has promoted our understanding of the pathogen virulence and population dynamics. Even with great progress, many questions still need to be answered. Therefore, we introduce important perspectives with emphasis on the genome evolution and host adaptation of rust fungi. We believe that the comparative genomics and population genomics of rust fungi will provide a further understanding of the rapid evolution of virulence and will contribute to monitoring the population dynamics for disease management.

A chromosome-level, fully phased genome assembly of the oat crown rust fungus Puccinia coronata f. sp. avenae: a resource to enable comparative genomics in the cereal rusts

Advances in sequencing technologies as well as development of algorithms and workflows have made it possible to generate fully phased genome references for organisms with nonhaploid genomes such as dikaryotic rust fungi. To enable discovery of pathogen effectors and further our understanding of virulence evolution, we generated a chromosome-scale assembly for each of the 2 nuclear genomes of the oat crown rust pathogen, Puccinia coronata f. sp. avenae (Pca). This resource complements 2 previously released partially phased genome references of Pca, which display virulence traits absent in the isolate of historic race 203 (isolate Pca203) which was selected for this genome project. A fully phased, chromosome-level reference for Pca203 was generated using PacBio reads and Hi-C data and a recently developed pipeline named NuclearPhaser for phase assignment of contigs and phase switch correction. With 18 chromosomes in each haplotype and a total size of 208.10 Mbp, Pca203 has the same number of chromosomes as other cereal rust fungi such as Puccinia graminis f. sp. tritici and Puccinia triticina, the causal agents of wheat stem rust and wheat leaf rust, respectively. The Pca203 reference marks the third fully phased chromosome-level assembly of a cereal rust to date. Here, we demonstrate that the chromosomes of these 3 Puccinia species are syntenous and that chromosomal size variations are primarily due to differences in repeat element content.

Comparative Genome Analyses of Plant Rust Pathogen Genomes Reveal a Confluence of Pathogenicity Factors to Quell Host Plant Defense Responses

Switchgrass rust caused by Puccinia novopanici (P. novopanici) has the ability to significantly affect the biomass yield of switchgrass, an important biofuel crop in the United States. A comparative genome analysis of P. novopanici with rust pathogen genomes infecting monocot cereal crops wheat, barley, oats, maize and sorghum revealed the presence of larger structural variations contributing to their genome sizes. A comparative alignment of the rust pathogen genomes resulted in the identification of collinear and syntenic relationships between P. novopanici and P. sorghi; P. graminis tritici 21–0 (Pgt 21) and P. graminis tritici Ug99 (Pgt Ug99) and between Pgt 21 and P. triticina (Pt). Repeat element analysis indicated a strong presence of retro elements among different Puccinia genomes, contributing to the genome size variation between ~1 and 3%. A comparative look at the enriched protein families of Puccinia spp. revealed a predominant role of restriction of telomere capping proteins (RTC), disulfide isomerases, polysaccharide deacetylases, glycoside hydrolases, superoxide dismutases and multi-copper oxidases (MCOs). All the proteomes of Puccinia spp. share in common a repertoire of 75 secretory and 24 effector proteins, including glycoside hydrolases cellobiohydrolases, peptidyl-propyl isomerases, polysaccharide deacetylases and protein disulfide-isomerases, that remain central to their pathogenicity. Comparison of the predicted effector proteins from Puccinia spp. genomes to the validated proteins from the Pathogen–Host Interactions database (PHI-base) resulted in the identification of validated effector proteins PgtSR1 (PGTG_09586) from P. graminis and Mlp124478 from Melampsora laricis across all the rust pathogen genomes.

Maintenance of divergent lineages of the Rice Blast Fungus Pyricularia oryzae through niche separation, loss of sex and post-mating genetic incompatibilities

Many species of fungal plant pathogens coexist as multiple lineages on the same host, but the factors underlying the origin and maintenance of population structure remain largely unknown. The rice blast fungus Pyricularia oryzae is a widespread model plant pathogen displaying population subdivision. However, most studies of natural variation in P. oryzae have been limited in genomic or geographic resolution, and host adaptation is the only factor that has been investigated extensively as a contributor to population subdivision. In an effort to complement previous studies, we analyzed genetic and phenotypic diversity in isolates of the rice blast fungus covering a broad geographical range. Using single-nucleotide polymorphism genotyping data for 886 isolates sampled from 152 sites in 51 countries, we showed that population subdivision of P. oryzae in one recombining and three clonal lineages with broad distributions persisted with deeper sampling. We also extended previous findings by showing further population subdivision of the recombining lineage into one international and three Asian clusters, and by providing evidence that the three clonal lineages of P. oryzae were found in areas with different prevailing environmental conditions, indicating niche separation. Pathogenicity tests and bioinformatic analyses using an extended set of isolates and rice varieties indicated that partial specialization to rice subgroups contributed to niche separation between lineages, and differences in repertoires of putative virulence effectors were consistent with differences in host range. Experimental crosses revealed that female sterility and early post-mating genetic incompatibilities acted as strong additional barriers to gene flow between clonal lineages. Our results demonstrate that the spread of a fungal pathogen across heterogeneous habitats and divergent populations of a crop species can lead to niche separation and reproductive isolation between distinct, widely distributed, lineages.

Reciprocal recombination genomic signatures in the symbiotic arbuscular mycorrhizal fungi Rhizophagus irregularis

Arbuscular mycorrhizal fungi (AMF) are part of the most widespread fungal-plant symbiosis. They colonize at least 80% of plant species, promote plant growth and plant diversity. These fungi are multinucleated and contain either one or two haploid nuclear genotypes (monokaryon and dikaryon) identified by the alleles at a putative mating-type locus. This taxon has been considered as an ancient asexual scandal because of the lack of observable sexual structures. Despite identification of a putative mating-type locus and functional activation of genes related to mating when two isolates co-exist, it remains unknown if the AMF life cycle involves a sexual or parasexual stage. We used publicly available genome sequences to test if Rhizophagus irregularis dikaryon genomes display signatures of sexual reproduction in the form of reciprocal recombination patterns, or if they display exclusively signatures of parasexual reproduction involving gene conversion. We used short-read and long-read sequence data to identify nucleus-specific alleles within dikaryons and then compared them to orthologous gene sequences from related monokaryon isolates displaying the same putative MAT-types as the dikaryon. We observed that the two nucleus-specific alleles of the dikaryon A5 are more related to the homolog sequences of monokaryon isolates displaying the same putative MAT-type than between each other. We also observed that these nucleus-specific alleles displayed reciprocal recombination signatures. These results confirm that dikaryon and monokaryon isolates displaying the same putative MAT-type are related in their life-cycle. These results suggest that a genetic exchange mechanism, involving reciprocal recombination in dikaryon genomes, allows AMF to generate genetic diversity.

Wheat stem rust recorded for the first time in decades in Ireland

Wheat stem rust, caused by the fungus Puccinia graminis f. sp. tritici (Pgt), occurs in most wheat-growing areas worldwide, and, in western Europe since 2013, has started to re-emerge after many decades of absence. Following this trend across western Europe, in 2020, we also detected and recorded wheat stem rust for the first time in five decades in experimental plots across five locations in Ireland. To examine the potential origin of the Irish Pgt infection in 2020, we carried out transcriptome sequencing on 12 Pgt-infected wheat samples collected across Ireland and compared these to 76 global P. graminis isolates. This analysis identified a close genetic relationship between the Irish Pgt isolates and those from Ethiopia collected in 2015 after a severe stem rust epidemic caused by the TKTTF Pgt race, and with the UK-01 Pgt isolate that was previously assigned to the TKTTF race. Subsequent pathology-based race profiling designated two Irish isolates and recent UK and French Pgt isolates to the TKTTF Pgt race group. This suggests that the Irish Pgt occurrence most probably originated from recent long-distance windborne dispersal of Pgt urediniospores from neighbouring countries in Europe where we confirmed the Pgt TKTTF race continues to be prevalent. The identification of wheat stem rust in Ireland at multiple locations in 2020 illustrates that the disease can occur in Ireland and emphasizes the need to re-initiate local monitoring for this re-emergent threat to wheat production across western Europe.

Sexual reproduction is the null hypothesis for life cycles of rust fungi

Sexual reproduction, mutation, and reassortment of nuclei increase genotypic diversity in rust fungi. Sexual reproduction is inherent to rust fungi, coupled with their coevolved plant hosts in native pathosystems. Rust fungi are hypothesised to exchange nuclei by somatic hybridisation with an outcome of increased genotypic diversity, independent of sexual reproduction. We provide criteria to demonstrate whether somatic exchange has occurred, including knowledge of parental haplotypes and rejection of fertilisation in normal rust life cycles.

Role of Sexual Reproduction in the Evolution of the Wheat Stripe Rust Fungus Races in China

Experimental and population genetic approaches have reshaped our view of how fungal pathogens reproduce, with consequences for our understanding of fungal invasions. Puccinia striiformis f. sp. tritici, the causal agent of stripe rust, poses a severe threat to wheat production worldwide. The sexual stage of P. striiformis f. sp. tritici was discovered >10 years ago, but how it affects the evolution of the pathogen, especially the emergence of the new virulent races, remains largely unknown. Here, using population genetic analyses, we demonstrate that sexual reproduction plays an important role in the evolution of P. striiformis f. sp. tritici races in China, specifically the newly emerged and devastating race virulent to resistance gene Yr26, which is widely used in China and exerts strong selective pressure on the pathogen population. Association analysis identified six genes encoding secreted proteins as candidates for virulence on wheat cultivars carrying the Yr26 resistance gene. Our results highlight the important role of sexual reproduction and selection exerted by hosts in the emergence of new virulent races in China.

Modelling the consequences of the dikaryotic life cycle of mushroom-forming fungi on genomic conflict

Generally, sexual organisms contain two haploid genomes, one from each parent, united in a single diploid nucleus of the zygote which links their fate during growth. A fascinating exception to this is Basidiomycete fungi, where the two haploid genomes remain separate in a dikaryon, retaining the option to fertilize subsequent monokaryons encountered. How the ensuing nuclear competition influences the balance of selection within and between individuals is largely unexplored. We test the consequences of the dikaryotic life cycle for mating success and mycelium-level fitness components. We assume a trade-off between mating fitness at the level of the haploid nucleus and fitness of the fungal mycelium. We show that the maintenance of fertilization potential by dikaryons leads to a higher proportion of fertilized monokaryons, but that the ensuing intradikaryon selection for increased nuclear mating fitness leads to reduced mycelium fitness relative to a diploid life cycle. However, this fitness reduction is lower compared to a hypothetical life cycle where dikaryons can also exchange nuclei. Prohibition of fusion between dikaryons therefore reduces the level of nuclear parasitism. The number of loci influencing fitness is an important determinant of the degree to which average mycelium-level fitness is reduced. The results of this study crucially hinge upon a trade-off between nucleus and mycelium-level fitness. We discuss the evidence for this assumption and the implications of an alternative that there is a positive relationship between nucleus and mycelium-level fitness.

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.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02658-2.

Deciphering the Crosstalk Mechanisms of Wheat-Stem Rust Pathosystem: Genome-Scale Prediction Unravels Novel Host Targets

Triticum aestivum (wheat), a major staple food grain, is affected by various biotic stresses. Among these, fungal diseases cause about 15-20% of yield loss, worldwide. In this study, we performed a comparative analysis of protein-protein interactions between two Puccinia graminis races (Pgt 21-0 and Pgt Ug99) that cause stem (black) rust in wheat. The available molecular techniques to study the host-pathogen interaction mechanisms are expensive and labor-intensive. We implemented two computational approaches (interolog and domain-based) for the prediction of PPIs and performed various functional analysis to determine the significant differences between the two pathogen races. The analysis revealed that T. aestivum-Pgt 21-0 and T. aestivum-Pgt Ug99 interactomes consisted of ∼90M and ∼56M putative PPIs, respectively. In the predicted PPIs, we identified 115 Pgt 21-0 and 34 Pgt Ug99 potential effectors that were highly involved in pathogen virulence and development. Functional enrichment analysis of the host proteins revealed significant GO terms and KEGG pathways such as O-methyltransferase activity (GO:0008171), regulation of signal transduction (GO:0009966), lignin metabolic process (GO:0009808), plastid envelope (GO:0009526), plant-pathogen interaction pathway (ko04626), and MAPK pathway (ko04016) that are actively involved in plant defense and immune signaling against the biotic stresses. Subcellular localization analysis anticipated the host plastid as a primary target for pathogen attack. The highly connected host hubs in the protein interaction network belonged to protein kinase domain including Ser/Thr protein kinase, MAPK, and cyclin-dependent kinase. We also identified 5,577 transcription factors in the interactions, associated with plant defense during biotic stress conditions. Additionally, novel host targets that are resistant to stem rust disease were also identified. The present study elucidates the functional differences between Pgt 21-0 and Pgt Ug99, thus providing the researchers with strain-specific information for further experimental validation of the interactions, and the development of durable, disease-resistant crop lines.

Discovery of a Novel Leaf Rust (Puccinia recondita) Resistance Gene in Rye (Secale cereale L.) Using Association Genomics

Leaf rust constitutes one of the most important foliar diseases in rye (Secale cereale L.). To discover new sources of resistance, we phenotyped 180 lines belonging to a less well-characterized Gülzow germplasm at three field trial locations in Denmark and Northern Germany in 2018 and 2019. We observed lines with high leaf rust resistance efficacy at all locations in both years. A genome-wide association study using 261,406 informative single-nucleotide polymorphisms revealed two genomic regions associated with resistance on chromosome arms 1RS and 7RS, respectively. The most resistance-associated marker on chromosome arm 1RS physically co-localized with molecular markers delimiting Pr3. In the reference genomes Lo7 and Weining, the genomic region associated with resistance on chromosome arm 7RS contained a large number of nucleotide-binding leucine-rich repeat (NLR) genes. Residing in close proximity to the most resistance-associated marker, we identified a cluster of NLRs exhibiting close protein sequence similarity with the wheat leaf rust Lr1 gene situated on chromosome arm 5DL in wheat, which is syntenic to chromosome arm 7RS in rye. Due to the close proximity to the most resistance-associated marker, our findings suggest that the considered leaf rust R gene, provisionally denoted Pr6, could be a Lr1 ortholog in rye.

An automated and combinative method for the predictive ranking of candidate effector proteins of fungal plant pathogens

Fungal plant-pathogens promote infection of their hosts through the release of 'effectors'-a broad class of cytotoxic or virulence-promoting molecules. Effectors may be recognised by resistance or sensitivity receptors in the host, which can determine disease outcomes. Accurate prediction of effectors remains a major challenge in plant pathology, but if achieved will facilitate rapid improvements to host disease resistance. This study presents a novel tool and pipeline for the ranking of predicted effector candidates-Predector-which interfaces with multiple software tools and methods, aggregates disparate features that are relevant to fungal effector proteins, and applies a pairwise learning to rank approach. Predector outperformed a typical combination of secretion and effector prediction methods in terms of ranking performance when applied to a curated set of confirmed effectors derived from multiple species. We present Predector ( https://github.com/ccdmb/predector ) as a useful tool for the ranking of predicted effector candidates, which also aggregates and reports additional supporting information relevant to effector and secretome prediction in a simple, efficient, and reproducible manner.

Genotyping Puccinia striiformis f. sp. tritici Isolates with SSR and SP-SNP Markers Reveals Dynamics of the Wheat Stripe Rust Pathogen in the United States from 1968 to 2009 and Identifies Avirulence-Associated Markers

Stripe rust, caused by Puccinia striiformis f. sp. tritici, is a devastating disease of wheat (Triticum aestivum) in the United States. The fungal pathogen can rapidly evolve, producing new virulent races infecting previously resistant cultivars and genotypes adapting to different environments. The objective of this study was to investigate the long-term population dynamics of P. striiformis f. sp. tritici in the United States. Through genotyping 1,083 isolates taken from 1968 to 2009, using 14 simple sequence repeat (SSR) markers and 92 secreted protein single nucleotide polymorphism (SP-SNP) markers, 614 and 945 genotypes were detected, respectively. In general, the two types of markers produced consistent genetic relationships among the P. striiformis f. sp. tritici populations over the 40-year period. The prior-to-2000 and the 2000-to-2009 populations were significantly different, with the latter showing higher genotypic diversity and higher heterozygosity than the earlier populations. Clustering analyses using genotypes of either SSR or SP-SNP markers revealed three molecular groups (MGs), MG1, MG2, and MG3. The prior-to-2000 and the 2000-to-2009 groups both had evidence of MG1 and MG2; however, MG3 was only found in the 2000-to-2009 population. Some of the isolates in the period of 2000 to 2009 formed individual clusters, suggesting exotic incursions. Other isolates of the same period were clustered with prior-to-2000 isolates, indicating that they were developed from the previously established populations. The data suggest the coexistence of newly introduced populations alongside established populations in the United States. Twenty SP-SNP markers were significantly associated to individual avirulence genes. These results are useful for developing more accurate monitoring systems and provide guidance for disease management.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01123-z.

Whole-genome sequencing of multiple isolates of Puccinia triticina reveals asexual lineages evolving by recurrent mutations

The wheat leaf rust fungus, Puccinia triticina Erikss., is a worldwide pathogen of tetraploid durum and hexaploid wheat. Many races of P. triticina differ for virulence to specific leaf rust resistance genes and are found in most wheat-growing regions of the world. Wheat cultivars with effective leaf rust resistance exert selection pressure on P. triticina populations for virulent race types. The objectives of this study were to examine whole-genome sequence data of 121 P. triticina isolates and to gain insight into race evolution. The collection included isolates comprising of many different race phenotypes collected worldwide from common and durum wheat. One isolate from wild wheat relative Aegilops speltoides and two from Ae. cylindrica were also included for comparison. Based on 121,907 informative variants identified relative to the reference Race 1-1 genome, isolates were clustered into 11 major lineages with 100% bootstrap support. The isolates were also grouped based on variation in 1311 predicted secreted protein genes. In gene-coding regions, all groups had high ratios of nonsynonymous to synonymous mutations and nonsense to readthrough mutations. Grouping of isolates based on two main variation principle components for either genome-wide variation or variation just within the secreted protein genes, indicated similar groupings. Variants were distributed across the entire genome, not just within the secreted protein genes. Our results suggest that recurrent mutation and selection play a major role in differentiation within the clonal lineages.

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 evolved^1,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 diseases³. A number of rust resistance genes have been characterized in wheat and most encode immune receptors of the nucleotide-binding leucine-rich repeat (NLR) class⁴, which recognize pathogen effector proteins known as avirulence (Avr) proteins⁵. However, only two Avr genes have been identified in Pgt so far, AvrSr35 and AvrSr50 (refs. ^6,7), and none in other cereal rusts^8,9. The Sr27 resistance gene was first identified in a wheat line carrying an introgression of the 3R chromosome from Imperial rye¹⁰. 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 Pgt^11,12. Sr27 is effective against Ug99 (ref. ¹³) and other recent Pgt strains^14,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.

Historical Development of the Puccinia triticina Population in South Africa

In contrast to many other countries, the virulence and genetic diversity of the South African Puccinia triticina population before 1980 is unknown, because of the absence of regular and systematic race analysis data and viable rust cultures. Herbarium specimens housed at the National Collection of Fungi, Biosystematics Unit, Plant Health and Protection, Agricultural Research Council, Pretoria, South Africa (SA), provided the opportunity to investigate the genetic development of the population using isolates collected between 1906 and 2010. Five subpopulations that survived between 21 and 82 years in the field were found. While three of these could represent the original races that entered SA during European settlement, two appear to be recent exotic introductions into SA, most probably from other African countries. The demise of the three oldest subpopulations might be from the release of resistant wheat cultivars. The population is clonal, where new virulence develops through single step mutations and selection for virulence. Although a possible case of somatic hybridization was found, sexual reproduction appears to be absent in SA. This study confirmed the importance of annual surveys in SA and its neighboring countries for the timely detection of new virulent races that could threaten wheat production in SA.