Distinct Life Histories Impact Dikaryotic Genome Evolution in the Rust Fungus Puccinia striiformis Causing Stripe Rust in Wheat

Use of Field pathogenomics approach for Puccinia striiformis f. sp. tritici race identification and phylogenomic delineation in North India

Stripe rust of wheat caused by Puccinia striiformis f. sp. tritici (Pst) is among the top 10 most important fungal phytopathogens in the world, threatening the global wheat production. Continuous pathogen evolution and air borne nature of the pathogen, increased the importance of pathogen population structure studies in case of wheat-Puccinia pathosystem. The current study aimed to characterize the Pst pathotype distribution in North India. A total of 61 stripe rust infected samples were collected from sub-mountainous areas of Punjab (50), Himachal Pradesh (HP) (8) and Jammu & Kashmir (J&K) (3) during 2021-24. Virulence profiling results of 58 samples using Indian stripe rust differentials revealed the prevalence of five pathotypes (238S119, 110S119, 46S119, 47S103 and 6S0) in North India. Among which 238S119 was the most predominant (51.72%) one followed by 110S119 (27.59%) and 46S119 (17.24%). No virulence was found against stripe rust resistant Yr5, Yr10, Yr15, Yr24 and YrSP genes indicating these genes are the most effective resistant genes in North India providing complete resistance till date. A total of 23 field infected samples were used for transcriptome analysis (RNA-seq) generating an average of 31.77 million raw reads. The clustered assembly generated mean assembly size of ∼ 45.10 Mb containing 7,587-2,49,571 contigs with 47.99% GC content, N50 value of 749.87 with 7.48-94.90% BUSCO score. Phylogenomic analysis of farmer's field samples revealed the existence of three different Pst lineages in North India i.e., Punjab lineage, Himachal lineage and Kashmir lineage, where Himachal lineage exhibited highest genetic diversity. A total of three (3) SSR and 14 KASP markers developed in previous studies were tested on farmer's field isolates for pathotype identification. Eight KASP markers showed reproducible polymorphism for pathotype specific detection. Correlating the KASP assay with results of virulence profiling revealed that these markers have the ability to detect the virulence change from five established pathotypes under field conditions. The findings of present study provided better understanding in Pst pathotype distribution in North India and pathotype detection using consistent polymorphic markers along with deciphering the behaviour of Pst under field conditions and deducing their evolutionary relationship in North India.

Developing a crop- wild-reservoir pathogen system to understand pathogen evolution and emergence

Crop pathogens reduce yield and contribute to global malnourishment. Surveillance not only detects presence/absence but also reveals genetic diversity, which can inform our understanding of rapid adaptation and control measures. An often neglected aspect is that pathogens may also use crop wild relatives as alternative hosts. This study develops the beet (Beta vulgaris) rust (Uromyces beticola) system to explore how crop pathogens evolve to evade resistance using a wild reservoir. We test predictions that crop selection will drive virulence gene differentiation and affect rates of sex between crop- and wild-host rust populations. We sequenced, assembled, and annotated the 588 Mb beet rust genome, developed a novel leaf peel pathogen DNA extraction protocol, and analysed genetic diversity in 42 wild and crop isolates. We found evidence for two populations: one containing exclusively wild-host isolates; the other containing all crop-host isolates, plus five wild isolates. Effectors showed greater diversity in the exclusively wild population and greater differentiation between populations. Preliminary evidence suggests the rates of sexual reproduction may differ between populations. This study highlights how differences in pathogen populations might be used to identify genes important for survival on crops and how reproduction might impact adaptation. These findings are relevant to all crop-reservoir systems and will remain unnoticed without comparison to wild reservoirs.

Genomics Research on the Road of Studying Biology and Virulence of Cereal Rust Fungi

Rust fungi are highly destructive pathogens that pose a significant threat to crop production worldwide, especially cereals. Obligate biotrophy and, in many cases, complex life cycles make rust fungi particularly challenging to study. However, recent rapid advances in sequencing technologies and genomic analysis tools have revolutionised rust fungal research. It is anticipated that the increasing availability and ongoing substantial improvements in genome assemblies will propel the field of rust biology into the post-genomic era, instigating a cascade of research endeavours encompassing multi-omics and gene discoveries. This is especially the case for many cereal rust pathogens, for which continental-scale studies of virulence have been conducted over many years and historical collections of viable isolates have been sequenced and assembled. Genomic analysis plays a crucial role in uncovering the underlying causes of the high variability of virulence and the complexity of population dynamics in rust fungi. Here, we provide an overview of progress in rust genomics, discuss the strategies employed in genomic analysis, and elucidate the strides that will drive cereal rust biology into the post-genomic era.

Migration of wheat stripe rust from the primary oversummering region to neighboring regions in China

Changing climate and changes in cropping systems have greatly affected outbreaks of plant diseases. Wheat stripe rust is a disease posing a threat to global wheat production, caused by Puccinia striiformis f. sp. tritici (Pst). Pst oversummering regions play a crucial role in the emergence of new races in China. To unveil the migration pattern of oversummering to adjacent regions, we develop a set of KASP-SNP marker from 28 Pst whole-genome sequences to investigate the population structure in the oversummering and its adjacent regions. A set of 19 Chinese wheat differentials is used to characterize the virulence patterns of 308 sampled Pst isolates. By integrating virulence characterization, population genetic analysis, air trajectory simulation and field disease monitoring, two main Pst dispersal routes are identified. Inocula from Eastern Qinghai are dispersed to Western and Eastern Liupan Mountain, and reach Guanzhong Plain. The second route originates from Middle Gansu, then through Longnan, and reaches the Guanzhong Plain via Eastern Liupan Mountain. Both dispersal routes result in Pst inoculum spreading to the Huang-Huai-Hai region, the main wheat-growing region in China. The proposed migration routes can be used to develop disease management strategies at a regional and national scale.

From Lesions to Lessons: Two Decades of Filamentous Plant Pathogen Genomics

Many filamentous microorganisms, such as fungi and oomycetes, have evolved the ability to colonize plants and cause devastating crop diseases. Coevolutionary conflicts with their hosts have shaped the genomes of these plant pathogens. Over the past 20 years, genomics and genomics-enabled technologies have revealed remarkable diversity in genome size, architecture, and gene regulatory mechanisms. Technical and conceptual advances continue to provide novel insights into evolutionary dynamics, diversification of distinct genomic compartments, and facilitated molecular disease diagnostics. In this review, we discuss how genomics has advanced our understanding of genome organization and plant-pathogen coevolution and provide a perspective on future developments in the field. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Beyond Asexual: Genomics-Driven Progress in Unveiling Sexual Reproduction in Cereal Rust Fungi

Recent advances in genomics technologies have revolutionized our understanding of cereal rust fungi, providing unprecedented insights into the complexities of their sexual life cycle. Genomic approaches, including long-read sequencing, genome assembly, and haplotype phasing technologies, have revealed critical insights into mating systems, genetic diversity, virulence evolution, and host adaptation. Population genomics studies have uncovered diverse reproductive strategies across different cereal rust species and geographic regions, highlighting the interplay between sexual recombination and asexual reproduction. Transcriptomics have begun to unravel the gene expression networks driving sexual reproduction, and complementary omics approaches such as proteomics and metabolomics offer potential insights into the underlying molecular processes. Despite this progress, many aspects of cereal rust sexual reproduction remain elusive. Integrating multiple omics approaches with advanced cell biology techniques can help address these knowledge gaps, particularly in understanding sexual reproduction and its role in pathogen evolution. This comprehensive approach will be crucial for developing more targeted and resilient crop protection strategies, ultimately contributing to global food security. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

From Natural Hosts to Agricultural Threats: The Evolutionary Journey of Phytopathogenic Fungi

Since the domestication of plants, pathogenic fungi have consistently threatened crop production, evolving genetically to develop increased virulence under various selection pressures. Understanding their evolutionary trends is crucial for predicting and designing control measures against future disease outbreaks. This paper reviews the evolution of fungal pathogens from natural habitats to agricultural settings, focusing on eight significant phytopathogens: Pyricularia oryzae, Botrytis cinerea, Puccinia spp., Fusarium graminearum, F. oxysporum, Blumeria graminis, Zymoseptoria tritici, and Colletotrichum spp. Also, we explore the mechanism used to understand evolutionary trends in these fungi. The studied pathogens have evolved in agroecosystems through either (1) introduction from elsewhere; or (2) local origins involving co-evolution with host plants, host shifts, or genetic variations within existing strains. Genetic variation, generated via sexual recombination and various asexual mechanisms, often drives pathogen evolution. While sexual recombination is rare and mainly occurs at the center of origin of the pathogen, asexual mechanisms such as mutations, parasexual recombination, horizontal gene or chromosome transfer, and chromosomal structural variations are predominant. Farming practices like mono-cropping resistant cultivars and prolonged use of fungicides with the same mode of action can drive the emergence of new pathotypes. Furthermore, host range does not necessarily impact pathogen adaptation and evolution. Although halting pathogen evolution is impractical, its pace can be slowed by managing selective pressures, optimizing farming practices, and enforcing quarantine regulations. The study of pathogen evolution has been transformed by advancements in molecular biology, genomics, and bioinformatics, utilizing methods like next-generation sequencing, comparative genomics, transcriptomics and population genomics. However, continuous research remains essential to monitor how pathogens evolve over time and to develop proactive strategies that mitigate their impact on agriculture.

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.

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.

Genome biology and evolution of mating-type loci in four cereal rust fungi

Permanent heterozygous loci, such as sex- or mating-compatibility regions, often display suppression of recombination and signals of genomic degeneration. In Basidiomycota, two distinct loci confer mating compatibility. These loci encode homeodomain (HD) transcription factors and pheromone receptor (Pra)-ligand allele pairs. To date, an analysis of genome level mating-type (MAT) loci is lacking for obligate biotrophic basidiomycetes in the Pucciniales, an order containing serious agricultural plant pathogens. Here, we focus on four species of Puccinia that infect oat and wheat, including P. coronata f. sp. avenae, P. graminis f. sp. tritici, P. triticina and P. striiformis f. sp. tritici. MAT loci are located on two separate chromosomes supporting previous hypotheses of a tetrapolar mating compatibility system in the Pucciniales. The HD genes are multiallelic in all four species while the PR locus appears biallelic, except for P. graminis f. sp. tritici, which potentially has multiple alleles. HD loci are largely conserved in their macrosynteny, both within and between species, without strong signals of recombination suppression. Regions proximal to the PR locus, however, displayed signs of recombination suppression and genomic degeneration in the three species with a biallelic PR locus. Our observations support a link between recombination suppression, genomic degeneration, and allele diversity of MAT loci that is consistent with recent mathematical modelling and simulations. Finally, we confirm that MAT genes are expressed during the asexual infection cycle, and we propose that this may support regulating nuclear maintenance and pairing during infection and spore formation. Our study provides insights into the evolution of MAT loci of key pathogenic Puccinia species. Understanding mating compatibility can help predict possible combinations of nuclear pairs, generated by sexual reproduction or somatic recombination, and the potential evolution of new virulent isolates of these important plant pathogens.

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.

In Planta Study Localizes an Effector Candidate from Austropuccinia psidii Strain MF-1 to the Nucleus and Demonstrates In Vitro Cuticular Wax-Dependent Differential Expression

Austropuccinia psidii is a biotrophic fungus that causes myrtle rust. First described in Brazil, it has since spread to become a globally important pathogen that infects more than 480 myrtaceous species. One of the most important commercial crops affected by A. psidii is eucalypt, a widely grown forestry tree. The A. psidii-Eucalyptus spp. interaction is poorly understood, but pathogenesis is likely driven by pathogen-secreted effector molecules. Here, we identified and characterized a total of 255 virulence effector candidates using a genome assembly of A. psidii strain MF-1, which was recovered from Eucalyptus grandis in Brazil. We show that the expression of seven effector candidate genes is modulated by cell wax from leaves sourced from resistant and susceptible hosts. Two effector candidates with different subcellular localization predictions, and with specific gene expression profiles, were transiently expressed with GFP-fusions in Nicotiana benthamiana leaves. Interestingly, we observed the accumulation of an effector candidate, Ap28303, which was upregulated under cell wax from rust susceptible E. grandis and described as a peptidase inhibitor I9 domain-containing protein in the nucleus. This was in accordance with in silico analyses. Few studies have characterized nuclear effectors. Our findings open new perspectives on the study of A. psidii-Eucalyptus interactions by providing a potential entry point to understand how the pathogen manipulates its hosts in modulating physiology, structure, or function with effector proteins.

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.

Genomic analysis, trajectory tracking, and field surveys reveal sources and long-distance dispersal routes of wheat stripe rust pathogen in China

Identifying sources of phytopathogen inoculum and determining their contributions to disease outbreaks are essential for predicting disease development and establishing control strategies. Puccinia striiformis f. sp. tritici (Pst), the causal agent of wheat stripe rust, is an airborne fungal pathogen with rapid virulence variation that threatens wheat production through its long-distance migration. Because of wide variation in geographic features, climatic conditions, and wheat production systems, Pst sources and related dispersal routes in China are largely unclear. In the present study, we performed genomic analyses of 154 Pst isolates from all major wheat-growing regions in China to determine Pst population structure and diversity. Through trajectory tracking, historical migration studies, genetic introgression analyses, and field surveys, we investigated Pst sources and their contributions to wheat stripe rust epidemics. We identified Longnan, the Himalayan region, and the Guizhou Plateau, which contain the highest population genetic diversities, as the Pst sources in China. Pst from Longnan disseminates mainly to eastern Liupan Mountain, the Sichuan Basin, and eastern Qinghai; that from the Himalayan region spreads mainly to the Sichuan Basin and eastern Qinghai; and that from the Guizhou Plateau migrates mainly to the Sichuan Basin and the Central Plain. These findings improve our current understanding of wheat stripe rust epidemics in China and emphasize the need for managing stripe rust on a national scale.

Puccinia striiformis f. sp. tritici effectors in wheat immune responses

The obligate biotrophic fungus Puccinia striiformis f. sp. tritici, which causes yellow (stripe) rust disease, is among the leading biological agents resulting in tremendous yield losses on global wheat productions per annum. The combatting strategies include, but are not limited to, fungicide applications and the development of resistant cultivars. However, evolutionary pressure drives rapid changes, especially in its "effectorome" repertoire, thus allowing pathogens to evade and breach resistance. The extracellular and intracellular effectors, predominantly secreted proteins, are tactical arsenals aiming for many defense processes of plants. Hence, the identity of the effectors and the molecular mechanisms of the interactions between the effectors and the plant immune system have long been targeted in research. The obligate biotrophic nature of P. striiformis f. sp. tritici and the challenging nature of its host, the wheat, impede research on this topic. Next-generation sequencing and novel prediction algorithms in bioinformatics, which are accompanied by in vitro and in vivo validation approaches, offer a speedy pace for the discovery of new effectors and investigations of their biological functions. Here, we briefly review recent findings exploring the roles of P. striiformis f. sp. tritici effectors together with their cellular/subcellular localizations, host responses, and interactors. The current status and the challenges will be discussed. We hope that the overall work will provide a broader view of where we stand and a reference point to compare and evaluate new findings.

Whole genome resequencing and comparative genome analysis of three Puccinia striiformis f. sp. tritici pathotypes prevalent in India

Stripe rust disease of wheat, caused by Puccinia striiformis f. sp. tritici, (Pst) is one of the most serious diseases of wheat worldwide. In India, virulent stripe rust races have been constantly evolving in the North-Western Plains Zone leading to the failure of some of the most widely grown resistant varieties in the region. With the goal of studying the recent evolution of virulent races in this region, we conducted whole-genome re-sequencing of three prevalent Indian Pst pathotypes Pst46S119, Pst78S84 and Pst110S119. We assembled 58.62, 58.33 and 55.78 Mb of Pst110S119, Pst46S119 and Pst78S84 genome, respectively and found that pathotypes were highly heterozygous. Comparative phylogenetic analysis indicated the recent evolution of pathotypes Pst110S119 and Pst78S84 from Pst46S119. Pathogenicity-related genes classes (CAZyme, proteases, effectors, and secretome proteins) were identified and found to be under positive selection. Higher rate of gene families expansion were also observed in the three pathotypes. A strong association between the effector genes and transposable elements may be the source of the rapid evolution of these strains. Phylogenetic analysis differentiated the Indian races in this study from other known United States, European, African, and Asian races. Diagnostic markers developed for the identification of three Pst pathotypes will help tracking of yellow rust at farmers field and strategizing resistance gene deployment.

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.

The evolving battle between yellow rust and wheat: implications for global food security

Wheat (Triticum aestivum L.) is a global commodity, and its production is a key component underpinning worldwide food security. Yellow rust, also known as stripe rust, is a wheat disease caused by the fungus Puccinia striiformis Westend f. sp. tritici (Pst), and results in yield losses in most wheat growing areas. Recently, the rapid global spread of genetically diverse sexually derived Pst races, which have now largely replaced the previous clonally propagated slowly evolving endemic populations, has resulted in further challenges for the protection of global wheat yields. However, advances in the application of genomics approaches, in both the host and pathogen, combined with classical genetic approaches, pathogen and disease monitoring, provide resources to help increase the rate of genetic gain for yellow rust resistance via wheat breeding while reducing the carbon footprint of the crop. Here we review key elements in the evolving battle between the pathogen and host, with a focus on solutions to help protect future wheat production from this globally important disease.

Incursions of divergent genotypes, evolution of virulence and host jumps shape a continental clonal population of the stripe rust pathogen Puccinia striiformis

Long-distance migration and host adaptation by transboundary plant pathogens often brings detrimental effects to important agroecosystems. Efficient surveillance as a basis for responding to the dynamics of such pathogens is often hampered by a lack of information on incursion origin, evolutionary pathways and the genetic basis of rapidly evolving virulence across larger timescales. Here, we studied these genetic features by using historical isolates of the obligate biotrophic pathogen Puccinia striiformis f. sp. tritici (Pst), which causes one of the most widespread and devastating diseases, stripe (yellow) rust, of wheat. Through a combination of genotypic, phenotypic and genomic analyses, we assigned eight Pst isolates representing putative exotic Pst incursions into Australia to four previously defined genetic groups, PstS0, PstS1, PstS10 and PstS13. We showed that isolates of an additional incursion of P. striiformis, known locally as P. striiformis f. sp. pseudo-hordei, had a new and unique multilocus SSR genotype (MLG). We provide results of overall genomic variation of representative Pst isolates from each genetic group by comparative genomic analyses. We showed that isolates within the PstS1 and PstS13 genetic groups are most distinct at the whole-genome variant level from isolates belonging to genetic group PstS0, whereas the isolate from the PstS10 genetic group is intermediate. We further explored variable gene content, including putative effectors, representing both shared but also unique genetic changes that have occurred following introduction, some of which may additionally account for local adaptation of these isolates to triticale. Our genotypic and genomic data revealed new genetic insights into the evolution of diverse phenotypes of rust pathogens following incursion into a geographically isolated continental region.

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.

A Chromosome-Scale Assembly of the Wheat Leaf Rust Pathogen Puccinia triticina Provides Insights Into Structural Variations and Genetic Relationships With Haplotype Resolution

Despite the global economic importance of the wheat leaf rust pathogen Puccinia triticina (Pt), genomic resources for Pt are limited and chromosome-level assemblies of Pt are lacking. Here, we present a complete haplotype-resolved genome assembly at a chromosome-scale for Pt using the Australian pathotype 64-(6),(7),(10),11 (Pt64; North American race LBBQB) built upon the newly developed technologies of PacBio and Hi-C sequencing. PacBio reads with ∼200-fold coverage (29.8 Gb data) were assembled by Falcon and Falcon-unzip and subsequently scaffolded with Hi-C data using Falcon-phase and Proximo. This approach allowed us to construct 18 chromosome pseudomolecules ranging from 3.5 to 12.3 Mb in size for each haplotype of the dikaryotic genome of Pt64. Each haplotype had a total length of ∼147 Mb, scaffold N 50 of ∼9.4 Mb, and was ∼93% complete for BUSCOs. Each haplotype had ∼29,800 predicted genes, of which ∼2,000 were predicted as secreted proteins (SPs). The investigation of structural variants (SVs) between haplotypes A and B revealed that 10% of the total genome was spanned by SVs, highlighting variations previously undetected by short-read based assemblies. For the first time, the mating type (MAT) genes on each haplotype of Pt64 were identified, which showed that MAT loci a and b are located on two chromosomes (chromosomes 7 and 14), representing a tetrapolar type. Furthermore, the Pt64 assembly enabled haplotype-based evolutionary analyses for 21 Australian Pt isolates, which highlighted the importance of a haplotype resolved reference when inferring genetic relationships using whole genome SNPs. This Pt64 assembly at chromosome-scale with full phase information provides an invaluable resource for genomic and evolutionary research, which will accelerate the understanding of molecular mechanisms underlying Pt-wheat interactions and facilitate the development of durable resistance to leaf rust in wheat and sustainable control of rust disease.

Austropuccinia psidii, causing myrtle rust, has a gigabase-sized genome shaped by transposable elements

Austropuccinia psidii, originating in South America, is a globally invasive fungal plant pathogen that causes rust disease on Myrtaceae. Several biotypes are recognized, with the most widely distributed pandemic biotype spreading throughout the Asia-Pacific and Oceania regions over the last decade. Austropuccinia psidii has a broad host range with more than 480 myrtaceous species. Since first detected in Australia in 2010, the pathogen has caused the near extinction of at least three species and negatively affected commercial production of several Myrtaceae. To enable molecular and evolutionary studies into A. psidii pathogenicity, we assembled a highly contiguous genome for the pandemic biotype. With an estimated haploid genome size of just over 1 Gb (gigabases), it is the largest assembled fungal genome to date. The genome has undergone massive expansion via distinct transposable element (TE) bursts. Over 90% of the genome is covered by TEs predominantly belonging to the Gypsy superfamily. These TE bursts have likely been followed by deamination events of methylated cytosines to silence the repetitive elements. This in turn led to the depletion of CpG sites in TEs and a very low overall GC content of 33.8%. Compared to other Pucciniales, the intergenic distances are increased by an order of magnitude indicating a general insertion of TEs between genes. Overall, we show how TEs shaped the genome evolution of A. psidii and provide a greatly needed resource for strategic approaches to combat disease spread.

Evolution and Adaptation of Forest and Crop Pathogens in the Anthropocene

Anthropocene marks the era when human activity is making a significant impact on earth, its ecological and biogeographical systems. The domestication and intensification of agricultural and forest production systems have had a large impact on plant and tree health. Some pathogens benefitted from these human activities and have evolved and adapted in response to the expansion of crop and forest systems, resulting in global outbreaks. Global pathogen genomics data including population genomics and high-quality reference assemblies are crucial for understanding the evolution and adaptation of pathogens. Crops and forest trees have remarkably different characteristics, such as reproductive time and the level of domestication. They also have different production systems for disease management with more intensive management in crops than forest trees. By comparing and contrasting results from pathogen population genomic studies done on widely different agricultural and forest production systems, we can improve our understanding of pathogen evolution and adaptation to different selection pressures. We find that in spite of these differences, similar processes such as hybridization, host jumps, selection, specialization, and clonal expansion are shaping the pathogen populations in both crops and forest trees. We propose some solutions to reduce these impacts and lower the probability of global pathogen outbreaks so that we can envision better management strategies to sustain global food production as well as ecosystem services.

Increased virulence of Puccinia coronata f. sp.avenae populations through allele frequency changes at multiple putative Avr loci

Pathogen populations are expected to evolve virulence traits in response to resistance deployed in agricultural settings. However, few temporal datasets have been available to characterize this process at the population level. Here, we examined two temporally separated populations of Puccinia coronata f. sp. avenae (Pca), which causes crown rust disease in oat (Avena sativa) sampled from 1990 to 2015. We show that a substantial increase in virulence occurred from 1990 to 2015 and this was associated with a genetic differentiation between populations detected by genome-wide sequencing. We found strong evidence for genetic recombination in these populations, showing the importance of the alternate host in generating genotypic variation through sexual reproduction. However, asexual expansion of some clonal lineages was also observed within years. Genome-wide association analysis identified seven Avr loci associated with virulence towards fifteen Pc resistance genes in oat and suggests that some groups of Pc genes recognize the same pathogen effectors. The temporal shift in virulence patterns in the Pca populations between 1990 and 2015 is associated with changes in allele frequency in these genomic regions. Nucleotide diversity patterns at a single Avr locus corresponding to Pc38, Pc39, Pc55, Pc63, Pc70, and Pc71 showed evidence of a selective sweep associated with the shift to virulence towards these resistance genes in all 2015 collected isolates.

A haplotype-phased genome of wheat stripe rust pathogen Puccinia striiformis f. sp. tritici, race PST-130 from the Western USA

More virulent and aggressive races of Puccinia striiformis f. sp. tritici (Pst), the pathogen causing wheat stripe rust, have been spreading around the world since 2000 causing large grain yield losses. A better understanding of the genome and genetic diversity of these new Pst races will be useful to develop new strategies to ameliorate these losses. In this study, we generated an improved genome assembly of a post-2000 virulent race from the Western USA designated as PST-130. We implemented a haplotype phasing strategy using the diploid-aware assembler, Falcon-Unzip and new long-read technology from PacBio to phase the two genomes of this dikaryotic organism. The combination of these new technologies resulted in an improved PST-130 assembly with only 151 contigs (85.4 Mb, N50 of 1.44 Mb), and a complementary assembly (haplotigs) with 458 contigs (65.9 Mb, N50 of 0.235 Mb, PRJNA650506). This new assembly improved gene predictions resulting in 228 more predicted complete genes than in the initial Illumina assembly (29,178 contigs, N50 of 5 kb). The alignment of the non-repetitive primary and haplotig contigs revealed and average of 5.22 SNP/kb, with 39.1% showing <2 SNP/kb and 15.9% >10 SNP/kb. This large divergent regions may represent introgressions of chromosome segments from more divergent Pst races in regions where a complete sexual cycle and recombination are possible. We hypothesize that some of the divergent regions in PST-130 may be related to the European "Warrior" race PST-DK0911 because this genome is more similar to PST-130 (3.18 SNP/kb) than to the older European race PST-104E (3.75 SNP/kb). Complete phasing of additional Pst genomes or sequencing individual nuclei will facilitate the tracing of the haploid genomes introduced by the new Pst races into local populations.