Genome expansion and gene loss in powdery mildew fungi reveal tradeoffs in extreme parasitism

Unraveling the pathogenomics of Rhizoctonia solani infecting proso millet (Panicum miliaceum L.): genomic perspective on ruthless virulence and adaptive evolution

Introduction

Banded sheath blight (Bsb), caused by Rhizoctonia solani, is an emerging threat to proso millet cultivation, significantly impacting yield and grain quality. This study on the pathogenomics of R. solani seeks to unravel its genetic mechanisms, identify key virulence factors, decode host-pathogen interactions, and pinpoint molecular targets for effective control strategies.

Methods

R. solani isolates were collected from various regions across India, resulting in six distinct isolates. These isolates were comprehensively characterized through morphological observations, molecular analyses, and virulence assessments to gain comprehensive insights into their diversity and pathogenic potential. The most virulent strain, designated VAP-1, infecting proso millet, was sequenced using the Illumina platform and de novo assembled using the SPAdes assembler, resulting in a highly complete genome. Functional regions of the genome were predicted and annotated using Funannotate. A subsequent comparative genomics study and secretome analysis were conducted to support functional genomic investigations.

Results

The VAP-1 genome assembly resulted in a total size of 47.12 Mb, with approximately 17.62% of the genome consisting of repetitive sequences, predominantly dominated by interspersed elements (around 97.8%). These interspersed elements were primarily classified as retrotransposons (72%), with DNA transposons comprising a smaller proportion (5%), while the remaining interspersed sequences were not fully annotated. Functional analysis of the genome revealed significant enrichment in KEGG pathways, including "Carbohydrate metabolism," "Translation," "Signal transduction," and "Transport and catabolism." In addition, Gene Ontology (GO) terms such as "Proteolysis," "Membrane," and "ATP binding" were notably enriched. The secretory protein profile of the VAP-1 genome from R. solani features key proteins from the major facilitator superfamily (MFS) transporters, (Trans) glycosidases, P-loop containing nucleoside triphosphate hydrolases, and galactose oxidase, all within the central domain superfamily. Glycoside hydrolases represent the largest class of CAZymes in the VAP-1 genome. Comparative genomic analysis of VAP-1 with other R. solani strains infecting Poaceae (e.g., rice) and non-Poaceae (e.g., sugar beet and tobacco) hosts showed that VAP-1 clusters closely with rice-infecting strains at the species level, yet exhibits a greater divergence in genomic similarity from strains infecting sugar beet and tobacco. Notably, variations were observed in important secretory proteins, such as multiple base deletions in MFS proteins across strains infecting proso millet, rice, and sugar beet.

Discussion

Functional analysis of the VAP-1 genome has unveiled a wealth of insights, though we have only begun to scratch the surface. KEGG and GO annotations point to critical proteins that are essential for host infection, providing the pathogen with a potent arsenal for successful penetration, survival, and dissemination within the host. The secretory proteins encoded in the VAP-1 genome play a pivotal role in equipping the pathogen with the necessary tools to degrade plant cell wall polymers, release cell wall-bound saccharides, and break down polysaccharides for energy utilization and host colonization. Notable variations were observed in several secretome superfamily proteins within the VAP-1 strain. These findings underscore the genomic diversity present within R. solani strains and suggest possible adaptations that may contribute to host specificity.

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.

Technical Advances Drive the Molecular Understanding of Effectors from Wheat and Barley Powdery Mildew Fungi

Pathogens manipulate host physiology through the secretion of virulence factors (effectors) to invade and proliferate on the host. The molecular functions of effectors inside plant hosts have been of interest in the field of molecular plant-microbe interactions. Obligate biotrophic pathogens, such as rusts and powdery mildews, cannot proliferate outside of plant hosts. In addition to the inhibition of the plant's immune components, these pathogens are under particular pressure to extract nutrients efficiently from the host. Understanding the molecular basis of infections mediated by obligate biotrophic pathogens is significant because of their impact in modern agriculture. In addition, powdery mildews serve as excellent models for obligate biotrophic cereal pathogens. Here, we summarize the current knowledge on the effectorome of the barley and wheat powdery mildews and putative molecular virulence functions of effectors. We emphasize the availability of comprehensive genomic, transcriptomic, and proteomic resources and discuss the methodological approaches used for identifying candidate effectors, assessing effector virulence traits, and identifying effector targets in the host. We highlight established and more recently employed methodologies, discuss limitations, and suggest additional strategies. We identify open questions and discuss how addressing them with currently available resources will enhance our understanding of Blumeria candidates for secretor effector proteins. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Sugars, Lipids and More: New Insights Into Plant Carbon Sources During Plant-Microbe Interactions

Heterotrophic microbes rely on host-derived carbon sources for their growth and survival. Depriving pathogens of plant carbon is therefore a promising strategy for protecting plants from disease and reducing yield losses. Importantly, this carbon starvation-mediated resistance is expected to be more broad-spectrum and durable than race-specific R-gene-mediated resistance. Although sugars are well characterized as major carbon sources for bacteria, emerging evidence suggests that plant-derived lipids are likely to be an essential carbon source for some fungal microbes, particularly biotrophs. Here, we comprehensively discuss the dual roles of carbon sources (mainly sugars and lipids) and their transport processes in immune signalling and microbial nutrition. We summarize recent findings revealing the crucial roles of lipids as susceptibility factors at all stages of pathogen infection. In particular, we discuss the potential pathways by which lipids and other plant carbon sources are delivered to biotrophs, including protein-mediated transport, vesicle trafficking and autophagy. Finally, we highlight knowledge gaps and offer suggestions for clarifying the mechanisms that underlie nutrient uptake by biotrophs, providing guidance for future research on the application of carbon starvation-mediated resistance.

A barley MLA immune receptor is activated by a fungal nonribosomal peptide effector for disease susceptibility

The barley Mla locus contains functionally diversified genes that encode intracellular nucleotide-binding leucine-rich repeat receptors (NLRs) and confer strain-specific immunity to biotrophic and hemibiotrophic fungal pathogens. In this study, we isolated a barley gene Scs6, which is an allelic variant of Mla genes but confers susceptibility to the isolate ND90Pr (BsND90Pr) of the necrotrophic fungus Bipolaris sorokiniana. We generated Scs6 transgenic barley lines and showed that Scs6 is sufficient to confer susceptibility to BsND90Pr in barley genotypes naturally lacking the receptor. The Scs6-encoded NLR (SCS6) is activated by a nonribosomal peptide (NRP) effector produced by BsND90Pr to induce cell death in barley and Nicotiana benthamiana. Domain swaps between MLAs and SCS6 reveal that the SCS6 leucine-rich repeat domain is a specificity determinant for receptor activation by the NRP effector. Scs6 is maintained in both wild and domesticated barley populations. Our phylogenetic analysis suggests that Scs6 is a Hordeum-specific innovation. We infer that SCS6 is a bona fide immune receptor that is likely directly activated by the nonribosomal peptide effector of BsND90Pr for disease susceptibility in barley. Our study provides a stepping stone for the future development of synthetic NLR receptors in crops that are less vulnerable to modification by necrotrophic pathogens.

Horizontal transfers between fungal Fusarium species contributed to successive outbreaks of coffee wilt disease

Outbreaks of fungal diseases have devastated plants and animals throughout history. Over the past century, the repeated emergence of coffee wilt disease caused by the fungal pathogen Fusarium xylarioides severely impacted coffee production across sub-Saharan Africa. To improve the disease management of such pathogens, it is crucial to understand their genetic structure and evolutionary potential. We compared the genomes of 13 historic strains spanning 6 decades and multiple disease outbreaks to investigate population structure and host specialisation. We found that F. xylarioides comprised at least 4 distinct lineages: 1 host-specific to Coffea arabica, 1 to C. canephora var. robusta, and 2 historic lineages isolated from various Coffea species. The presence/absence of large genomic regions across populations, the higher genetic similarities of these regions between species than expected based on genome-wide divergence and their locations in different loci in genomes across populations showed that horizontal transfers of effector genes from members of the F. oxysporum species complex contributed to host specificity. Multiple transfers into F. xylarioides populations matched different parts of the F. oxysporum mobile pathogenicity chromosome and were enriched in effector genes and transposons. Effector genes in this region and other carbohydrate-active enzymes important in the breakdown of plant cell walls were shown by transcriptomics to be highly expressed during infection of C. arabica by the fungal arabica strains. Widespread sharing of specific transposons between F. xylarioides and F. oxysporum, and the correspondence of a putative horizontally transferred regions to a Starship (large mobile element involved in horizontal gene transfers in fungi), reinforce the inference of horizontal transfers and suggest that mobile elements were involved. Our results support the hypothesis that horizontal gene transfers contributed to the repeated emergence of coffee wilt disease.

To be or not to be a nonhost species: A case study of the Leptosphaeria maculans and Brassica carinata interaction

Leptosphaeria maculans is one of the major fungal pathogens on oilseed rape (Brassica napus), causing stem canker disease. The closely related Brassica species B. nigra, B. juncea, and B. carinata display extreme resistance toward stem canker. In this study, we demonstrate the nonhost status of B. carinata toward L. maculans in France through field experiments and inoculations performed in controlled conditions. A few isolates moderately adapted to B. carinata in controlled conditions were recovered in the field on B. nigra leaves, allowing us to investigate the unusual B. carinata-L. maculans interactions using molecular, macroscopic, and microscopic analyses. A cross between a L. maculans isolate adapted to B. napus and an isolate moderately adapted to B. carinata allowed the generation, in the lab, of recombinant L. maculans strains better adapted to B. carinata than the natural parental isolate obtained from B. nigra, and highlighted the polygenic determinism of the adaptation of L. maculans to B. carinata and B. napus. This biological material will allow further investigation of the molecular determinants of the adaptation of L. maculans to nonhost species and elucidate the genetic resistance basis of B. carinata.

A fungal plant pathogen overcomes mlo-mediated broad-spectrum disease resistance by rapid gene loss

Hosts and pathogens typically engage in a coevolutionary arms race. This also applies to phytopathogenic powdery mildew fungi, which can rapidly overcome plant resistance and perform host jumps. Using experimental evolution, we show that the powdery mildew pathogen Blumeria hordei is capable of breaking the agriculturally important broad-spectrum resistance conditioned by barley loss-of-function mlo mutants. Partial mlo virulence of evolved B. hordei isolates is correlated with a distinctive pattern of adaptive mutations, including small-sized (c. 8-40 kb) deletions, of which one is linked to the de novo insertion of a transposable element. Occurrence of the mutations is associated with a transcriptional induction of effector protein-encoding genes that is absent in mlo-avirulent isolates on mlo mutant plants. The detected mutational spectrum comprises the same loci in at least two independently isolated mlo-virulent isolates, indicating convergent multigenic evolution. The mutational events emerged in part early (within the first five asexual generations) during experimental evolution, likely generating a founder population in which incipient mlo virulence was later stabilized by additional events. This work highlights the rapid dynamic genome evolution of an obligate biotrophic plant pathogen with a transposon-enriched genome.

Diverse epistatic effects in barley-powdery mildew interactions localize to host chromosome hotspots

Barley Mildew locus a (Mla) encodes a multi-allelic series of nucleotide-binding leucine-rich repeat (NLR) receptors that specify recognition to diverse cereal diseases. We exploited time-course transcriptome dynamics of barley and derived immune mutants infected with the powdery mildew fungus, Blumeria hordei (Bh), to infer gene effects governed by Mla6 and two other loci significant to disease development, Blufensin1 (Bln1), and Required for Mla6 resistance3 (rar3 = Sgt1 ΔKL308-309 ). Interactions of Mla6 and Bln1 resulted in diverse epistatic effects on the Bh-induced barley transcriptome, differential immunity to Pseudomonas syringae expressing the effector protease AvrPphB, and reaction to Bh. From a total of 468 barley NLRs, 115 were grouped under different gene effect models; genes classified under these models localized to host chromosome hotspots. The corresponding Bh infection transcriptome was classified into nine co-expressed modules, linking differential expression with pathogen structures, signifying that disease is regulated by an inter-organismal network that diversifies the response.

First genome assembly of the order Strepsiptera using PacBio HiFi reads reveals a miniature genome

Twisted-wing insects (Strepsiptera) are an enigmatic order of parasites with unusual life histories and striking sexual dimorphism. Males emerge from hosts as free-living winged adults, while females from most species remain as endoparasites that retain larval traits. Due to scarce genomic data and phylogenetic controversies, Strepsiptera was only recently placed as the closest living relative to beetles (Coleoptera). Here, we report the first PacBio HiFi genome assembly of the strepsipteran Xenos peckii (Xenidae). This de novo assembly size is 72.1 Mb, with a BUSCO score of 87.4%, N50 of 7.3 Mb, 23.4% GC content, and 38.41% repeat content. We identified 8 contigs that contain >75% of the assembly and reflect the haploid chromosome number reported from karyotypic data, and 3 contigs that exhibit sex chromosome coverage patterns. Additionally, the mitochondrial genome is 16,111 bp long and has 37 genes. This long-read assembly for Strepsiptera reveals a miniature genome and provides a unique tool to understand complex genome evolution associated with a parasitic lifestyle and extreme sexual dimorphism.

Genomic insights into the evolution of secondary metabolism of Escovopsis and its allies, specialized fungal symbionts of fungus-farming ants

The metabolic intimacy of symbiosis often demands the work of specialists. Natural products and defensive secondary metabolites can drive specificity by ensuring infection and propagation across host generations. But in contrast to bacteria, little is known about the diversity and distribution of natural product biosynthetic pathways among fungi and how they evolve to facilitate symbiosis and adaptation to their host environment. In this study, we define the secondary metabolism of Escovopsis and closely related genera, symbionts in the gardens of fungus-farming ants. We ask how the gain and loss of various biosynthetic pathways correspond to divergent lifestyles. Long-read sequencing allowed us to define the chromosomal features of representative Escovopsis strains, revealing highly reduced genomes composed of seven to eight chromosomes. The genomes are highly syntenic with macrosynteny decreasing with increasing phylogenetic distance, while maintaining a high degree of mesosynteny. An ancestral state reconstruction analysis of biosynthetic pathways revealed that, while many secondary metabolites are shared with non-ant-associated Sordariomycetes, 56 pathways are unique to the symbiotic genera. Reflecting adaptation to diverging ant agricultural systems, we observe that the stepwise acquisition of these pathways mirrors the ecological radiations of attine ants and the dynamic recruitment and replacement of their fungal cultivars. As different clades encode characteristic combinations of biosynthetic gene clusters, these delineating profiles provide important insights into the possible mechanisms underlying specificity between these symbionts and their fungal hosts. Collectively, our findings shed light on the evolutionary dynamic nature of secondary metabolism in Escovopsis and its allies, reflecting adaptation of the symbionts to an ancient agricultural system.IMPORTANCEMicrobial symbionts interact with their hosts and competitors through a remarkable array of secondary metabolites and natural products. Here, we highlight the highly streamlined genomic features of attine-associated fungal symbionts. The genomes of Escovopsis species, as well as species from other symbiont genera, many of which are common with the gardens of fungus-growing ants, are defined by seven chromosomes. Despite a high degree of metabolic conservation, we observe some variation in the symbionts' potential to produce secondary metabolites. As the phylogenetic distribution of the encoding biosynthetic gene clusters coincides with attine transitions in agricultural systems, we highlight the likely role of these metabolites in mediating adaptation by a group of highly specialized symbionts.

Retracing the evolution of Pneumocystis species, with a focus on the human pathogen Pneumocystis jirovecii

SUMMARYEvery human being is presumed to be infected by the fungus Pneumocystis jirovecii at least once in his or her lifetime. This fungus belongs to a large group of species that appear to exclusively infect mammals, with P. jirovecii being the only one known to cause disease in humans. The mystery of P. jirovecii origin and speciation is just beginning to unravel. Here, we provide a review of the major steps of P. jirovecii evolution. The Pneumocystis genus likely originated from soil or plant-associated organisms during the period of Cretaceous ~165 million years ago and successfully shifted to mammals. The transition coincided with a substantial loss of genes, many of which are related to the synthesis of nutrients that can be scavenged from hosts or cell wall components that could be targeted by the mammalian immune system. Following the transition, the Pneumocystis genus cospeciated with mammals. Each species specialized at infecting its own host. Host specialization is presumably built at least partially upon surface glycoproteins, whose protogene was acquired prior to the genus formation. P. jirovecii appeared at ~65 million years ago, overlapping with the emergence of the first primates. P. jirovecii and its sister species P. macacae, which infects macaques nowadays, may have had overlapping host ranges in the distant past. Clues from molecular clocks suggest that P. jirovecii did not cospeciate with humans. Molecular evidence suggests that Pneumocystis speciation involved chromosomal rearrangements and the mounting of genetic barriers that inhibit gene flow among species.

Whole-genome sequencing of marine water-derived Curvularia verruculosa KHW-7: a pioneering study

Marine microorganisms are renowned for being a rich source of new secondary metabolites that are significant to humans. The fungi strain KHW-7 was isolated from the seawater collected from the Gulf of Khambhat, India, and identified as Curvularia verruculosa KHW-7. On a next-generation sequencing platform, C. verruculosa KHW-7's whole-genome sequencing (WGS) and gene annotation were carried out using several bioinformatic methods. The 31.59 MB genome size, 52.3% GC, and 158 bp mean read length were discovered using WGS. This genome also contained 9,745 protein-coding genes, including 852 secreted proteins and 2048 transmembrane proteins. The antiSMASH algorithm used to analyze genomes found 25 secondary metabolite biosynthetic gene clusters (BGCs) that are abundant in terpene, non-ribosomal peptide synthetase (NRPS), and polyketides type 1 (T1PKS). To our knowledge, this is the first whole-genome sequence report of C. verruculosa. The WGS analysis of C. verruculosa KHW-7 indicated that this marine-derived fungus could be an efficient generator of bioactive secondary metabolites and an important industrial enzyme, both of which demand further investigation and development.

Genome sequencing and comparative genomics reveal insights into pathogenicity and evolution of Fusarium zanthoxyli, the causal agent of stem canker in prickly ash

BACKGROUND

Fusarium zanthoxyli is a destructive pathogen causing stem canker in prickly ash, an ecologically and economically important forest tree. However, the genome lack of F. zanthoxyli has hindered research on its interaction with prickly ash and the development of precise control strategies for stem canker.

RESULTS

In this study, we sequenced and annotated a relatively high-quality genome of F. zanthoxyli with a size of 43.39 Mb, encoding 11,316 putative genes. Pathogenicity-related factors are predicted, comprising 495 CAZymes, 217 effectors, 156 CYP450s, and 202 enzymes associated with secondary metabolism. Besides, a comparative genomics analysis revealed Fusarium and Colletotrichum diverged from a shared ancestor approximately 141.1 ~ 88.4 million years ago (MYA). Additionally, a phylogenomic investigation of 12 different phytopathogens within Fusarium indicated that F. zanthoxyli originated approximately 34.6 ~ 26.9 MYA, and events of gene expansion and contraction within them were also unveiled. Finally, utilizing conserved domain prediction, the results revealed that among the 59 unique genes, the most enriched domains were PnbA and ULP1. Among the 783 expanded genes, the most enriched domains were PKc_like kinases and those belonging to the APH_ChoK_Like family.

CONCLUSION

This study sheds light on the genetic basis of F. zanthoxyli's pathogenicity and evolution which provides valuable information for future research on its molecular interactions with prickly ash and the development of effective strategies to combat stem canker.

Plant-Entomopathogenic Fungi Interaction: Recent Progress and Future Prospects on Endophytism-Mediated Growth Promotion and Biocontrol

Entomopathogenic fungi, often acknowledged primarily for their insecticidal properties, fulfill diverse roles within ecosystems. These roles encompass endophytism, antagonism against plant diseases, promotion of the growth of plants, and inhabitation of the rhizosphere, occurring both naturally and upon artificial inoculation, as substantiated by a growing body of contemporary research. Numerous studies have highlighted the beneficial aspects of endophytic colonization. This review aims to systematically organize information concerning the direct (nutrient acquisition and production of phytohormones) and indirect (resistance induction, antibiotic and secondary metabolite production, siderophore production, and mitigation of abiotic and biotic stresses) implications of endophytic colonization. Furthermore, a thorough discussion of these mechanisms is provided. Several challenges, including isolation complexities, classification of novel strains, and the impact of terrestrial location, vegetation type, and anthropogenic reluctance to use fungal entomopathogens, have been recognized as hurdles. However, recent advancements in biotechnology within microbial research hold promising solutions to many of these challenges. Ultimately, the current constraints delineate potential future avenues for leveraging endophytic fungal entomopathogens as dual microbial control agents.

Long-term and rapid evolution in powdery mildew fungi

The powdery mildew fungi (Erysiphaceae) are globally distributed plant pathogens with a range of more than 10,000 plant hosts. In this review, we discuss the long- and short-term evolution of these obligate biotrophic fungi and outline their diversity with respect to morphology, lifestyle, and host range. We highlight their remarkable ability to rapidly overcome plant immunity, evolve fungicide resistance, and broaden their host range, for example, through adaptation and hybridization. Recent advances in genomics and proteomics, particularly in cereal powdery mildews (genus Blumeria), provided first insights into mechanisms of genomic adaptation in these fungi. Transposable elements play key roles in shaping their genomes, where even close relatives exhibit diversified patterns of recent and ongoing transposon activity. These transposons are ubiquitously distributed in the powdery mildew genomes, resulting in a highly adaptive genome architecture lacking obvious regions of conserved gene space. Transposons can also be neofunctionalized to encode novel virulence factors, particularly candidate secreted effector proteins, which may undermine the plant immune system. In cereals like barley and wheat, some of these effectors are recognized by plant immune receptors encoded by resistance genes with numerous allelic variants. These effectors determine incompatibility ("avirulence") and evolve rapidly through sequence diversification and copy number variation. Altogether, powdery mildew fungi possess plastic genomes that enable their fast evolutionary adaptation towards overcoming plant immunity, host barriers, and chemical stress such as fungicides, foreshadowing future outbreaks, host range shifts and expansions as well as potential pandemics by these pathogens.

Powdery mildew effectors AVRA1 and BEC1016 target the ER J-domain protein HvERdj3B required for immunity in barley

The barley powdery mildew fungus, Blumeria hordei (Bh), secretes hundreds of candidate secreted effector proteins (CSEPs) to facilitate pathogen infection and colonization. One of these, CSEP0008, is directly recognized by the barley nucleotide-binding leucine-rich-repeat (NLR) receptor MLA1 and therefore is designated AVRA1. Here, we show that AVRA1 and the sequence-unrelated Bh effector BEC1016 (CSEP0491) suppress immunity in barley. We used yeast two-hybrid next-generation interaction screens (Y2H-NGIS), followed by binary Y2H and in planta protein-protein interactions studies, and identified a common barley target of AVRA1 and BEC1016, the endoplasmic reticulum (ER)-localized J-domain protein HvERdj3B. Silencing of this ER quality control (ERQC) protein increased Bh penetration. HvERdj3B is ER luminal, and we showed using split GFP that AVRA1 and BEC1016 translocate into the ER signal peptide-independently. Overexpression of the two effectors impeded trafficking of a vacuolar marker through the ER; silencing of HvERdj3B also exhibited this same cellular phenotype, coinciding with the effectors targeting this ERQC component. Together, these results suggest that the barley innate immunity, preventing Bh entry into epidermal cells, requires ERQC. Here, the J-domain protein HvERdj3B appears to be essential and can be regulated by AVRA1 and BEC1016. Plant disease resistance often occurs upon direct or indirect recognition of pathogen effectors by host NLR receptors. Previous work has shown that AVRA1 is directly recognized in the cytosol by the immune receptor MLA1. We speculate that the AVRA1 J-domain target being inside the ER, where it is inapproachable by NLRs, has forced the plant to evolve this challenging direct recognition.

The massive 340 megabase genome of Anisogramma anomala, a biotrophic ascomycete that causes eastern filbert blight of hazelnut

BACKGROUND

The ascomycete fungus Anisogramma anomala causes Eastern Filbert Blight (EFB) on hazelnut (Corylus spp.) trees. It is a minor disease on its native host, the American hazelnut (C. americana), but is highly destructive on the commercially important European hazelnut (C. avellana). In North America, EFB has historically limited commercial production of hazelnut to west of the Rocky Mountains. A. anomala is an obligately biotrophic fungus that has not been grown in continuous culture, rendering its study challenging. There is a 15-month latency before symptoms appear on infected hazelnut trees, and only a sexual reproductive stage has been observed. Here we report the sequencing, annotation, and characterization of its genome.

RESULTS

The genome of A. anomala was assembled into 108 scaffolds totaling 342,498,352 nt with a GC content of 34.46%. Scaffold N50 was 33.3 Mb and L50 was 5. Nineteen scaffolds with lengths over 1 Mb constituted 99% of the assembly. Telomere sequences were identified on both ends of two scaffolds and on one end of another 10 scaffolds. Flow cytometry estimated the genome size of A. anomala at 370 Mb. The genome exhibits two-speed evolution, with 93% of the assembly as AT-rich regions (32.9% GC) and the other 7% as GC-rich (57.1% GC). The AT-rich regions consist predominantly of repeats with low gene content, while 90% of predicted protein coding genes were identified in GC-rich regions. Copia-like retrotransposons accounted for more than half of the genome. Evidence of repeat-induced point mutation (RIP) was identified throughout the AT-rich regions, and two copies of the rid gene and one of dim-2, the key genes in the RIP mutation pathway, were identified in the genome. Consistent with its homothallic sexual reproduction cycle, both MAT1-1 and MAT1-2 idiomorphs were found. We identified a large suite of genes likely involved in pathogenicity, including 614 carbohydrate active enzymes, 762 secreted proteins and 165 effectors.

CONCLUSIONS

This study reveals the genomic structure, composition, and putative gene function of the important pathogen A. anomala. It provides insight into the molecular basis of the pathogen's life cycle and a solid foundation for studying EFB.

Deciphering the Genomic Landscape and Virulence Mechanisms of the Wheat Powdery Mildew Pathogen Blumeria graminis f. sp. tritici Wtn1: Insights from Integrated Genome Assembly and Conidial Transcriptomics

A high-quality genome sequence from an Indian isolate of Blumeria graminis f. sp. tritici Wtn1, a persistent threat in wheat farming, was obtained using a hybrid method. The assembly of over 9.24 million DNA-sequence reads resulted in 93 contigs, totaling a 140.61 Mb genome size, potentially encoding 8480 genes. Notably, more than 73.80% of the genome, spanning approximately 102.14 Mb, comprises retro-elements, LTR elements, and P elements, influencing evolution and adaptation significantly. The phylogenomic analysis placed B. graminis f. sp. tritici Wtn1 in a distinct monocot-infecting clade. A total of 583 tRNA anticodon sequences were identified from the whole genome of the native virulent strain B. graminis f. sp. tritici, which comprises distinct genome features with high counts of tRNA anticodons for leucine (70), cysteine (61), alanine (58), and arginine (45), with only two stop codons (Opal and Ochre) present and the absence of the Amber stop codon. Comparative InterProScan analysis unveiled "shared and unique" proteins in B. graminis f. sp. tritici Wtn1. Identified were 7707 protein-encoding genes, annotated to different categories such as 805 effectors, 156 CAZymes, 6102 orthologous proteins, and 3180 distinct protein families (PFAMs). Among the effectors, genes like Avra10, Avrk1, Bcg-7, BEC1005, CSEP0105, CSEP0162, BEC1016, BEC1040, and HopI1 closely linked to pathogenesis and virulence were recognized. Transcriptome analysis highlighted abundant proteins associated with RNA processing and modification, post-translational modification, protein turnover, chaperones, and signal transduction. Examining the Environmental Information Processing Pathways in B. graminis f. sp. tritici Wtn1 revealed 393 genes across 33 signal transduction pathways. The key pathways included yeast MAPK signaling (53 genes), mTOR signaling (38 genes), PI3K-Akt signaling (23 genes), and AMPK signaling (21 genes). Additionally, pathways like FoxO, Phosphatidylinositol, the two-component system, and Ras signaling showed significant gene representation, each with 15-16 genes, key SNPs, and Indels in specific chromosomes highlighting their relevance to environmental responses and pathotype evolution. The SNP and InDel analysis resulted in about 3.56 million variants, including 3.45 million SNPs, 5050 insertions, and 5651 deletions within the whole genome of B. graminis f. sp. tritici Wtn1. These comprehensive genome and transcriptome datasets serve as crucial resources for understanding the pathogenicity, virulence effectors, retro-elements, and evolutionary origins of B. graminis f. sp. tritici Wtn1, aiding in developing robust strategies for the effective management of wheat powdery mildew.

Guanosine-specific single-stranded ribonuclease effectors of a phytopathogenic fungus potentiate host immune responses

Plants activate immunity upon recognition of pathogen-associated molecular patterns. Although phytopathogens have evolved a set of effector proteins to counteract plant immunity, some effectors are perceived by hosts and induce immune responses. Here, we show that two secreted ribonuclease effectors, SRN1 and SRN2, encoded in a phytopathogenic fungus, Colletotrichum orbiculare, induce cell death in a signal peptide- and catalytic residue-dependent manner, when transiently expressed in Nicotiana benthamiana. The pervasive presence of SRN genes across Colletotrichum species suggested the conserved roles. Using a transient gene expression system in cucumber (Cucumis sativus), an original host of C. orbiculare, we show that SRN1 and SRN2 potentiate host pattern-triggered immunity responses. Consistent with this, C. orbiculare SRN1 and SRN2 deletion mutants exhibited increased virulence on the host. In vitro analysis revealed that SRN1 specifically cleaves single-stranded RNAs at guanosine, leaving a 3'-end phosphate. Importantly, the potentiation of C. sativus responses by SRN1 and SRN2, present in the apoplast, depends on ribonuclease catalytic residues. We propose that the pathogen-derived apoplastic guanosine-specific single-stranded endoribonucleases lead to immunity potentiation in plants.

Long noncoding RNAs emerge from transposon-derived antisense sequences and may contribute to infection stage-specific transposon regulation in a fungal phytopathogen

BACKGROUND

The genome of the obligate biotrophic phytopathogenic barley powdery mildew fungus Blumeria hordei is inflated due to highly abundant and possibly active transposable elements (TEs). In the absence of the otherwise common repeat-induced point mutation transposon defense mechanism, noncoding RNAs could be key for regulating the activity of TEs and coding genes during the pathogenic life cycle.

RESULTS

We performed time-course whole-transcriptome shotgun sequencing (RNA-seq) of total RNA derived from infected barley leaf epidermis at various stages of fungal pathogenesis and observed significant transcript accumulation and time point-dependent regulation of TEs in B. hordei. Using a manually curated consensus database of 344 TEs, we discovered phased small RNAs mapping to 104 consensus transposons, suggesting that RNA interference contributes significantly to their regulation. Further, we identified 5,127 long noncoding RNAs (lncRNAs) genome-wide in B. hordei, of which 823 originated from the antisense strand of a TE. Co-expression network analysis of lncRNAs, TEs, and coding genes throughout the asexual life cycle of B. hordei points at extensive positive and negative co-regulation of lncRNAs, subsets of TEs and coding genes.

CONCLUSIONS

Our work suggests that similar to mammals and plants, fungal lncRNAs support the dynamic modulation of transcript levels, including TEs, during pivotal stages of host infection. The lncRNAs may support transcriptional diversity and plasticity amid loss of coding genes in powdery mildew fungi and may give rise to novel regulatory elements and virulence peptides, thus representing key drivers of rapid evolutionary adaptation to promote pathogenicity and overcome host defense.

Near-complete de novo assembly of Tricholoma bakamatsutake chromosomes revealed the structural divergence and differentiation of Tricholoma genomes

Tricholoma bakamatsutake, which is an edible ectomycorrhizal fungus associated with Fagaceae trees, may have diverged before the other species in Tricholoma section Caligata. We generated a highly contiguous whole-genome sequence for T. bakamatsutake SF-Tf05 isolated in an Oak (Quercus salicina) forest in Japan. The assembly of high-fidelity long reads, with a median read length of 12.3 kb, resulted in 13 chromosome-sized contigs comprising 142,068,211 bases with an average guanine and cytosine (GC) content of 43.94%. The 13 chromosomes were predicted to encode 11,060 genes. A contig (122,566 bases) presumably containing the whole circular mitochondrial genome was also recovered. The chromosome-wide comparison of T. bakamatsutake and Tricholoma matsutake (TMA_r1.0) indicated that the basic number of chromosomes (13) was conserved, but the structures of the corresponding chromosomes diverged, with multiple inversions and translocations. Gene conservation and cluster analyses revealed at least 3 phylogenetic clades in Tricholoma section Caligata. Specifically, all T. bakamatsutake strains belonged to the "bakamatsutake" clade, which is most proximal to the "caligatum" clade consisting of Tricholoma caligatum and Tricholoma fulvocastaneum. The constructed highly contiguous nearly telomere-to-telomere genome sequence of a T. bakamatsutake isolate will serve as a fundamental resource for future research on the evolution and differentiation of Tricholoma species.

Cell-type-specific responses to fungal infection in plants revealed by single-cell transcriptomics

Pathogen infection is a dynamic process. Here, we employ single-cell transcriptomics to investigate plant response heterogeneity. By generating an Arabidopsis thaliana leaf atlas encompassing 95,040 cells during infection by a fungal pathogen, Colletotrichum higginsianum, we unveil cell-type-specific gene expression, notably an enrichment of intracellular immune receptors in vasculature cells. Trajectory inference identifies cells that had different interactions with the invading fungus. This analysis divulges transcriptional reprogramming of abscisic acid signaling specifically occurring in guard cells, which is consistent with a stomatal closure dependent on direct contact with the fungus. Furthermore, we investigate the transcriptional plasticity of genes involved in glucosinolate biosynthesis in cells at the fungal infection sites, emphasizing the contribution of the epidermis-expressed MYB122 to disease resistance. This work underscores spatially dynamic, cell-type-specific plant responses to a fungal pathogen and provides a valuable resource that supports in-depth investigations of plant-pathogen interactions.

An effector of Erysiphe necator translocates to chloroplasts and plasma membrane to suppress host immunity in grapevine

The powdery mildew (Erysiphe necator) is a prevalent pathogen hampering grapevine growth in the vineyard. An arsenal of candidate secreted effector proteins (CSEPs) was encoded in the E. necator genome, but it is largely unclear what role CSEPs plays during the E. necator infection. In the present study, we identified a secreted effector CSEP080 of E. necator, which was located in plant chloroplasts and plasma membrane. Transient expressing CSEP080 promotes plant photosynthesis and inhibits INF1-induced cell death in tobacco leaves. We found that CSEP080 was a necessary effector for the E. necator pathogenicity, which interacted with grapevine chloroplast protein VviB6f (cytochrome b6-f complex iron-sulfur subunit), affecting plant photosynthesis. Transient silencing VviB6f increased the plant hydrogen peroxide production, and the plant resistance to powdery mildew. In addition, CSEP080 manipulated the VviPE (pectinesterase) to promote pectin degradation. Our results demonstrated the molecular mechanisms that an effector of E. necator translocates to host chloroplasts and plasma membrane, which suppresses with the grapevine immunity system by targeting the chloroplast protein VviB6f to suppress hydrogen peroxide accumulation and manipulating VviPE to promote pectin degradation.

A chromosome-scale genome assembly of the grape powdery mildew pathogen Erysiphe necator reveals its genomic architecture and previously unknown features of its biology

Erysiphe necator is an obligate fungal pathogen that causes grape powdery mildew, globally the most important disease on grapevines. Previous attempts to obtain a quality genome assembly for this pathogen were hindered by its high repetitive DNA content. Here, chromatin conformation capture (Hi-C) with long-read PacBio sequencing was combined to obtain a chromosome-scale assembly and a high-quality annotation for E. necator isolate EnFRAME01. The resulting 81.1 Mb genome assembly is 98% complete and consists of 34 scaffolds, 11 of which represent complete chromosomes. All chromosomes contain large centromeric-like regions and lack synteny to the 11 chromosomes of the cereal PM pathogen Blumeria graminis. Further analysis of their composition showed that repeats and transposable elements (TEs) occupy 62.7% of their content. TEs were almost evenly interspersed outside centromeric and telomeric regions and massively overlapped with regions of annotated genes, suggesting that they could have a significant functional impact. Abundant gene duplicates were observed as well, particularly in genes encoding candidate secreted effector proteins. Moreover, younger in age gene duplicates exhibited more relaxed selection pressure and were more likely to be located physically close in the genome than older duplicates. A total of 122 genes with copy number variations among six isolates of E. necator were also identified and were enriched in genes that were duplicated in EnFRAME01, indicating they may reflect an adaptive variation. Taken together, our study illuminates higher-order genomic architectural features of E. necator and provides a valuable resource for studying genomic structural variations in this pathogen. IMPORTANCE Grape powdery mildew caused by the ascomycete fungus Erysiphe necator is economically the most important and recurrent disease in vineyards across the world. The obligate biotrophic nature of E. necator hinders the use of typical genetic methods to elucidate its pathogenicity and adaptation to adverse conditions, and thus comparative genomics has been a major method to study its genome biology. However, the current reference genome of E. necator isolate C-strain is highly fragmented with many non-coding regions left unassembled. This incompleteness prohibits in-depth comparative genomic analyses and the study of genomic structural variations (SVs) that are known to affect several aspects of microbial life, including fitness, virulence, and host adaptation. By obtaining a chromosome-scale genome assembly and a high-quality gene annotation for E. necator, we reveal the organization of its chromosomal content, unearth previously unknown features of its biology, and provide a reference for studying genomic SVs in this pathogen.

Suppression of Chitin-Triggered Immunity by Plant Fungal Pathogens: A Case Study of the Cucurbit Powdery Mildew Fungus Podosphaera xanthii

Fungal pathogens are significant plant-destroying microorganisms that present an increasing threat to the world's crop production. Chitin is a crucial component of fungal cell walls and a conserved MAMP (microbe-associated molecular pattern) that can be recognized by specific plant receptors, activating chitin-triggered immunity. The molecular mechanisms underlying the perception of chitin by specific receptors are well known in plants such as rice and Arabidopsis thaliana and are believed to function similarly in many other plants. To become a plant pathogen, fungi have to suppress the activation of chitin-triggered immunity. Therefore, fungal pathogens have evolved various strategies, such as prevention of chitin digestion or interference with plant chitin receptors or chitin signaling, which involve the secretion of fungal proteins in most cases. Since chitin immunity is a very effective defensive response, these fungal mechanisms are believed to work in close coordination. In this review, we first provide an overview of the current understanding of chitin-triggered immune signaling and the fungal proteins developed for its suppression. Second, as an example, we discuss the mechanisms operating in fungal biotrophs such as powdery mildew fungi, particularly in the model species Podosphaera xanthii, the main causal agent of powdery mildew in cucurbits. The key role of fungal effector proteins involved in the modification, degradation, or sequestration of immunogenic chitin oligomers is discussed in the context of fungal pathogenesis and the promotion of powdery mildew disease. Finally, the use of this fundamental knowledge for the development of intervention strategies against powdery mildew fungi is also discussed.

Genome-wide characterization of L-aspartate oxidase genes in wheat and their potential roles in the responses to wheat disease and abiotic stresses

L-aspartate oxidase (AO) is the first enzyme in NAD⁺ biosynthesis and is widely distributed in plants, animals, and microorganisms. Recently, AO family members have been reported in several plants, including Arabidopsis thaliana and Zea mays. Research on AO in these plants has revealed that AO plays important roles in plant growth, development, and biotic stresses; however, the nature and functions of AO proteins in wheat are still unclear. In this study, nine AO genes were identified in the wheat genome via sequence alignment and conserved protein domain analysis. These nine wheat AO genes (TaAOs) were distributed on chromosomes 2, 5, and 6 of sub-genomes A, B, and D. Analysis of the phylogenetic relationships, conserved motifs, and gene structure showed that the nine TaAOs were clustered into three groups, and the TaAOs in each group had similar conserved motifs and gene structure. Meanwhile, the subcellular localization analysis of transient expression mediated by Agrobacterium tumetioniens indicated that TaAO3-6D was localized to chloroplasts. Prediction of cis-elements indicated that a large number of cis-elements involved in responses to ABA, SA, and antioxidants/electrophiles, as well as photoregulatory responses, were found in TaAO promoters, which suggests that the expression of TaAOs may be regulated by these factors. Finally, transcriptome and real-time PCR analysis showed that the expression of TaAOs belonging to Group III was strongly induced in wheat infected by F. graminearum during anthesis, while the expression of TaAOs belonging to Group I was heavily suppressed. Additionally, the inducible expression of TaAOs belonging to Group III during anthesis in wheat spikelets infected by F. graminearum was repressed by ABA. Finally, expression of almost all TaAOs was induced by exposure to cold treatment. These results indicate that TaAOs may participate in the response of wheat to F. graminearum infection and cold stress, and ABA may play a negative role in this process. This study lays a foundation for further investigation of TaAO genes and provides novel insights into their biological functions.

Structural and functional analysis of the cerato-platanin-like protein Cpl1 suggests diverging functions in smut fungi

Plant-pathogenic fungi are causative agents of the majority of plant diseases and can lead to severe crop loss in infected populations. Fungal colonization is achieved by combining different strategies, such as avoiding and counteracting the plant immune system and manipulating the host metabolome. Of major importance are virulence factors secreted by fungi, which fulfil diverse functions to support the infection process. Most of these proteins are highly specialized, with structural and biochemical information often absent. Here, we present the atomic structures of the cerato-platanin-like protein Cpl1 from Ustilago maydis and its homologue Uvi2 from Ustilago hordei. Both proteins adopt a double-Ψβ-barrel architecture reminiscent of cerato-platanin proteins, a class so far not described in smut fungi. Our structure-function analysis shows that Cpl1 binds to soluble chitin fragments via two extended grooves at the dimer interface of the two monomer molecules. This carbohydrate-binding mode has not been observed previously and expands the repertoire of chitin-binding proteins. Cpl1 localizes to the cell wall of U. maydis and might synergize with cell wall-degrading and decorating proteins during maize infection. The architecture of Cpl1 harbouring four surface-exposed loop regions supports the idea that it might play a role in the spatial coordination of these proteins. While deletion of cpl1 has only mild effects on the virulence of U. maydis, a recent study showed that deletion of uvi2 strongly impairs U. hordei virulence. Our structural comparison between Cpl1 and Uvi2 reveals sequence variations in the loop regions that might explain a diverging function.

Site-specific analysis reveals candidate cross-kingdom small RNAs, tRNA and rRNA fragments, and signs of fungal RNA phasing in the barley-powdery mildew interaction

The establishment of host-microbe interactions requires molecular communication between both partners, which may involve the mutual transfer of noncoding small RNAs. Previous evidence suggests that this is also true for powdery mildew disease in barley, which is caused by the fungal pathogen Blumeria hordei. However, previous studies lacked spatial resolution regarding the accumulation of small RNAs upon host infection by B. hordei. Here, we analysed site-specific small RNA repertoires in the context of the barley-B. hordei interaction. To this end, we dissected infected leaves into separate fractions representing different sites that are key to the pathogenic process: epiphytic fungal mycelium, infected plant epidermis, isolated haustoria, a vesicle-enriched fraction from infected epidermis, and extracellular vesicles. Unexpectedly, we discovered enrichment of specific 31-33-base 5'-terminal fragments of barley 5.8S ribosomal RNA in extracellular vesicles and infected epidermis, as well as particular B. hordei transfer RNA fragments in haustoria. We describe canonical small RNAs from both the plant host and the fungal pathogen that may confer cross-kingdom RNA interference activity. Interestingly, we found first evidence of phased small interfering RNAs in B. hordei, a feature usually attributed to plants, which may be associated with the posttranscriptional control of fungal coding genes, pseudogenes, and transposable elements. Our data suggest a key and possibly site-specific role for cross-kingdom RNA interference and noncoding RNA fragments in the host-pathogen communication between B. hordei and its host barley.

The Biological Roles of Puccinia striiformis f. sp. tritici Effectors during Infection of Wheat

Puccinia striiformis f. sp. tritici (Pst) is the causative agent of wheat stripe rust, which can lead to a significant loss in annual wheat yields. Therefore, there is an urgent need for a deeper comprehension of the basic mechanisms underlying Pst infection. Effectors are known as the agents that plant pathogens deliver into host tissues to promote infection, typically by interfering with plant physiology and biochemistry. Insights into effector activity can significantly aid the development of future strategies to generate disease-resistant crops. However, the functional analysis of Pst effectors is still in its infancy, which hinders our understanding of the molecular mechanisms of the interaction between Pst and wheat. In this review, we summarize the potential roles of validated and proposed Pst effectors during wheat infection, including proteinaceous effectors, non-coding RNAs (sRNA effectors), and secondary metabolites (SMs effectors). Further, we suggest specific countermeasures against Pst pathogenesis and future research directions, which may promote our understanding of Pst effector functions during wheat immunity attempts.

Two pathogen loci determine Blumeria graminis f. sp. tritici virulence to wheat resistance gene Pm1a

Blumeria graminis f. sp. tritici (Bgt) is a globally important fungal pathogen of wheat that can rapidly evolve to defeat wheat powdery mildew (Pm) resistance genes. Despite periodic regional deployment of the Pm1a resistance gene in US wheat production, Bgt strains that overcome Pm1a have been notably nonpersistent in the United States, while on other continents, they are more widely established. A genome-wide association study (GWAS) was conducted to map sequence variants associated with Pm1a virulence in 216 Bgt isolates from six countries, including the United States. A virulence variant apparently unique to Bgt isolates from the United States was detected in the previously mapped gene AvrPm1a (BgtE-5612) on Bgt chromosome 6; an in vitro growth assay suggested no fitness reduction associated with this variant. A gene on Bgt chromosome 8, Bgt-51526, was shown to function as a second determinant of Pm1a virulence, and despite < 30% amino acid identity, BGT-51526 and BGTE-5612 were predicted to share > 85% of their secondary structure. A co-expression study in Nicotiana benthamiana showed that BGTE-5612 and BGT-51526 each produce a PM1A-dependent hypersensitive response. More than one member of a B. graminis effector family can be recognized by a single wheat immune receptor, and a two-gene model is necessary to explain virulence to Pm1a.

Haplotype-phased and chromosome-level genome assembly of Puccinia polysora, a giga-scale fungal pathogen causing southern corn rust

Rust fungi are characterized by large genomes with high repeat content and have two haploid nuclei in most life stages, which makes achieving high-quality genome assemblies challenging. Here, we described a pipeline using HiFi reads and Hi-C data to assemble a gigabase-sized fungal pathogen, Puccinia polysora f.sp. zeae, to haplotype-phased and chromosome-scale. The final assembled genome is 1.71 Gbp, with ~850 Mbp and 18 chromosomes in each haplotype, being currently one of the two giga-scale fungi assembled to chromosome level. Transcript-based annotation identified 47,512 genes for the dikaryotic genome with a similar number for each haplotype. A high level of interhaplotype variation was found with 10% haplotype-specific BUSCO genes, 5.8 SNPs/kbp, and structural variation accounting for 3% of the genome size. The P. polysora genome displayed over 85% repeat contents, with genome-size expansion and copy number increasing of species-specific orthogroups. Interestingly, these features did not affect overall synteny with other Puccinia species having smaller genomes. Fine-time-point transcriptomics revealed seven clusters of coexpressed secreted proteins that are conserved between two haplotypes. The fact that candidate effectors interspersed with all genes indicated the absence of a "two-speed genome" evolution in P. polysora. Genome resequencing of 79 additional isolates revealed a clonal population structure of P. polysora in China with low geographic differentiation. Nevertheless, a minor population differentiated from the major population by having mutations on secreted proteins including AvrRppC, indicating the ongoing virulence to evade recognition by RppC, a major resistance gene in Chinese corn cultivars. The high-quality assembly provides valuable genomic resources for future studies on disease management and the evolution of P. polysora.

Genome-Wide Characterization of Effector Protein-Encoding Genes in Sclerospora graminicola and Its Validation in Response to Pearl Millet Downy Mildew Disease Stress

Pearl millet [Pennisetum glaucum (L.) R. Br.] is the essential food crop for over ninety million people living in drier parts of India and South Africa. Pearl millet crop production is harshly hindered by numerous biotic stresses. Sclerospora graminicola causes downy mildew disease in pearl millet. Effectors are the proteins secreted by several fungi and bacteria that manipulate the host cell structure and function. This current study aims to identify genes encoding effector proteins from the S. graminicola genome and validate them through molecular techniques. In silico analyses were employed for candidate effector prediction. A total of 845 secretory transmembrane proteins were predicted, out of which 35 proteins carrying LxLFLAK (Leucine–any amino acid–Phenylalanine–Leucine–Alanine–Lysine) motif were crinkler, 52 RxLR (Arginine, any amino acid, Leucine, Arginine), and 17 RxLR-dEER putative effector proteins. Gene validation analysis of 17 RxLR-dEER effector protein-producing genes was carried out, of which 5genes were amplified on the gel. These novel gene sequences were submitted to NCBI. This study is the first report on the identification and characterization of effector genes in Sclerospora graminicola. This dataset will aid in the integration of effector classes that act independently, paving the way to investigate how pearl millet responds to effector protein interactions. These results will assist in identifying functional effector proteins involving the omic approach using newer bioinformatics tools to protect pearl millet plants against downy mildew stress. Considered together, the identified effector protein-encoding functional genes can be utilized in screening oomycetes downy mildew diseases in other crops across the globe.

Sclerotinia sclerotiorum (Lib.) de Bary: Insights into the Pathogenomic Features of a Global Pathogen

Sclerotinia sclerotiorum (Lib.) de Bary is a broad host-range fungus that infects an inclusive array of plant species and afflicts significant yield losses globally. Despite being a notorious pathogen, it has an uncomplicated life cycle consisting of either basal infection from myceliogenically germinated sclerotia or aerial infection from ascospores of carpogenically germinated sclerotia. This fungus is unique among necrotrophic pathogens in that it inevitably colonizes aging tissues to initiate an infection, where a saprophytic stage follows the pathogenic phase. The release of cell wall-degrading enzymes, oxalic acid, and effector proteins are considered critical virulence factors necessary for the effective pathogenesis of S. sclerotiorum. Nevertheless, the molecular basis of S. sclerotiorum pathogenesis is still imprecise and remains a topic of continuing research. Previous comprehensive sequencing of the S. sclerotiorum genome has revealed new insights into its genome organization and provided a deeper comprehension of the sophisticated processes involved in its growth, development, and virulence. This review focuses on the genetic and genomic aspects of fungal biology and molecular pathogenicity to summarize current knowledge of the processes utilized by S. sclerotiorum to parasitize its hosts. Understanding the molecular mechanisms regulating the infection process of S. sclerotiorum will contribute to devising strategies for preventing infections caused by this destructive pathogen.

The nuclear effector ArPEC25 from the necrotrophic fungus Ascochyta rabiei targets the chickpea transcription factor CaβLIM1a and negatively modulates lignin biosynthesis, increasing host susceptibility

Fungal pathogens deploy a barrage of secreted effectors to subvert host immunity, often by evading, disrupting, or altering key components of transcription, defense signaling, and metabolic pathways. However, the underlying mechanisms of effectors and their host targets are largely unexplored in necrotrophic fungal pathogens. Here, we describe the effector protein Ascochyta rabiei PEXEL-like Effector Candidate 25 (ArPEC25), which is secreted by the necrotroph A. rabiei, the causal agent of Ascochyta blight disease in chickpea (Cicer arietinum), and is indispensable for virulence. After entering host cells, ArPEC25 localizes to the nucleus and targets the host LIM transcription factor CaβLIM1a. CaβLIM1a is a transcriptional regulator of CaPAL1, which encodes phenylalanine ammonia lyase (PAL), the regulatory, gatekeeping enzyme of the phenylpropanoid pathway. ArPEC25 inhibits the transactivation of CaβLIM1a by interfering with its DNA-binding ability, resulting in negative regulation of the phenylpropanoid pathway and decreased levels of intermediates of lignin biosynthesis, thereby suppressing lignin production. Our findings illustrate the role of fungal effectors in enhancing virulence by targeting a key defense pathway that leads to the biosynthesis of various secondary metabolites and antifungal compounds. This study provides a template for the study of less explored necrotrophic effectors and their host target functions.

Mitochondrial genomes in the iconic reindeer lichens: Architecture, variation, and synteny across multiple evolutionary scales

Variation in mitochondrial genome composition across intraspecific, interspecific, and higher taxonomic scales has been little studied in lichen obligate symbioses. Cladonia is one of the most diverse and ecologically important lichen genera, with over 500 species representing an array of unique morphologies and chemical profiles. Here, we assess mitochondrial genome diversity and variation in this flagship genus, with focused sampling of two clades of the "true" reindeer lichens, Cladonia subgenus Cladina, and additional genomes from nine outgroup taxa. We describe composition and architecture at the gene and the genome scale, examining patterns in organellar genome size in larger taxonomic groups in Ascomycota. Mitochondrial genomes of Cladonia, Pilophorus, and Stereocaulon were consistently larger than those of Lepraria and contained more introns, suggesting a selective pressure in asexual morphology in Lepraria driving it toward genomic simplification. Collectively, lichen mitochondrial genomes were larger than most other fungal life strategies, reaffirming the notion that coevolutionary streamlining does not correlate to genome size reductions. Genomes from Cladonia ravenelii and Stereocaulon pileatum exhibited ATP9 duplication, bearing paralogs that may still be functional. Homing endonuclease genes (HEGs), though scarce in Lepraria, were diverse and abundant in Cladonia, exhibiting variable evolutionary histories that were sometimes independent of the mitochondrial evolutionary history. Intraspecific HEG diversity was also high, with C. rangiferina especially bearing a range of HEGs with one unique to the species. This study reveals a rich history of events that have transformed mitochondrial genomes of Cladonia and related genera, allowing future study alongside a wealth of assembled genomes.

Barley RIC157, a potential RACB scaffold protein, is involved in susceptibility to powdery mildew

CRIB motif-containing barley RIC157 is a novel ROP scaffold protein that interacts directly with barley RACB, promotes susceptibility to fungal penetration, and colocalizes with RACB at the haustorial neck. Successful obligate pathogens benefit from host cellular processes. For the biotrophic ascomycete fungus Blumeria hordei (Bh) it has been shown that barley RACB, a small monomeric G-protein (ROP, Rho of plants), is required for full susceptibility to fungal penetration. The susceptibility function of RACB probably lies in its role in cell polarity, which may be co-opted by the pathogen for invasive ingrowth of its haustorium. However, how RACB supports fungal penetration success and which other host proteins coordinate this process is incompletely understood. RIC (ROP-Interactive and CRIB-(Cdc42/Rac Interactive Binding) motif-containing) proteins are considered scaffold proteins which can interact directly with ROPs via a conserved CRIB motif. Here we describe a previously uncharacterized barley RIC protein, RIC157, which can interact directly with RACB in planta. We show that, in the presence of constitutively activated RACB, RIC157 shows a localization at the cell periphery/plasma membrane, whereas it otherwise localizes to the cytoplasm. RIC157 appears to mutually stabilize the plasma membrane localization of the activated ROP. During fungal infection, RIC157 and RACB colocalize at the penetration site, particularly at the haustorial neck. Additionally, transiently overexpressed RIC157 renders barley epidermal cells more susceptible to fungal penetration. We discuss that RIC157 may promote fungal penetration into barley epidermal cells by operating probably downstream of activated RACB.

Comparative genomics of five Valsa species gives insights on their pathogenicity evolution

Valsa is a genus of ascomycetes within the Valsaceae family. This family includes many wood destructive pathogens such as the well known Valsa mali and Valsa pyri which cause canker diseases in fruit trees and threaten the global fruit production. Lack of genomic information of this family is impeding our understandings about their evolution and genetic basis of their pathogenicity divergence. Here, we report genome assemblies of Valsa malicola, Valsa persoonii, and Valsa sordida which represent close relatives of Valsa mali and Valsa pyri with different host preferences. Comparative genomics analysis revealed that segmental rearrangements, inversions, and translocations frequently occurred among Valsa spp. genomes. Gene families that exhibited gene copy expansions tended to be associated with secondary metabolism, transmembrane transport, and pyrophosphatase activities. Orthologous genes in regions lost synteny exhibited significantly higher rate of synonymous substitution (KS) than those in regions retained synteny. Moreover, among these genes, membrane transporter families associated with antidrug (MFS, DHA) activities and nutrient transportation (SP and APCs) activities were significantly over-represented. Lineage specific synonymous substitution (KS) and nonsynonymous substitution (KA) analysis based on the phylogeny constructed from 11 fungal species identified a set of genes with selection signatures in Valsa clade and these genes were significantly enriched in functions associated with fatty acid beta-oxidation, DNA helicase activity, and ATPase activity. Furthermore, unique genes that possessed or retained by each of the five Valsa species are more likely part of the secondary metabolic (SM) gene clusters. SM gene clusters conserved across five Valsa species showed various degrees of diversification in both identity and completeness. All 11 syntenically conserved SM clusters showed differential expression during the infection of apple branch with Valsa mali suggesting involvements of secondary metabolism in the pathogenicity of Valsa species.

A Genomic Resource for the Strawberry Powdery Mildew Pathogen Podosphaera aphanis

Powdery mildew is one of the most economically destructive diseases in protected strawberry production. Here we present the first genome assembly for Podosphaera aphanis, the causal agent of powdery mildew on strawberry. This obligate-biotrophic fungal pathogen was sampled from a naturally occurring outbreak on Fragaria × ananassa 'Malling Centenary' plants grown under cover in the United Kingdom. Assembled reads resolved a 55.6 Mb genome, composed of 12,357 contigs whose annotation led to prediction of 17,239 genes encoding 17,328 proteins. The genome is highly-complete, with 97.5% of conserved single-copy Ascomycete genes shown to be present. This annotated P. aphanis genome provides a molecular resource for further investigation into host-pathogen interactions in the strawberry powdery mildew pathosystem.

Warmer and drier ecosystems select for smaller bacterial genomes in global soils

Bacterial genome size reflects bacterial evolutionary processes and metabolic lifestyles, with implications for microbial community assembly and ecosystem functions. However, to understand the extent of genome-mediated microbial responses to environmental selections, we require studies that observe genome size distributions along environmental gradients representing different conditions that soil bacteria normally encounter. In this study, we used surface soils collected from 237 sites across the globe and analyzed how environmental conditions (e.g., soil carbon and nutrients, aridity, pH, and temperature) affect soil bacterial occurrences and genome size at the community level using bacterial community profiling. We used a joint species distribution model to quantify the effects of environments on species occurrences and found that aridity was a major regulator of genome size with warmer and drier environments selecting bacteria with smaller genomes. Drought-induced physiological constraints on bacterial growth (e.g., water scarcity for cell metabolisms) may have led to these correlations. This finding suggests that increasing cover by warmer and drier ecosystems may result in bacterial genome simplifications by a reduction of genome size.

Biotrophic Fungal Pathogens: a Critical Overview

Biotrophic fungi are one group of heterogeneous organisms and these fungi differ in their traits like mode of nutrition, types of reproduction, and dispersal systems. Generally, based on the nutritional mode, fungi are classified into three broad categories, viz. biotrophs, necrotrophs, and hemi-biotrophs. Biotrophs derive their nutrients and energy from living plant cells and survive within the interstitial space of the cells. Biotrophic fungi cause serious crop diseases but are highly challenging to investigate and develop a treatment strategy. Blumeria (Erysiphe) graminis, Uromyces fabae, Ustilago maydis, Cladosporium fulvum, Puccinia graminis, and Phytophthora infestans are some of the significant biotrophic fungi that affect mainly plants. One among the biotrophic fungus, Pneumocystis jirovecii (Taphrinomycotina subphylum of the Ascomycota) exclusively a human pathogen, can cause lung diseases such as "pneumocystis." Biotrophic fungus widely parasitizing Solanaceae family crops (Tomato and potato) has done massive damage to the crops and has led to economic impact worldwide. During infection and for nutrient absorption, biotrophs develops external appendages such as appressoria or haustoria. The hyphae or appressorium adheres to the plant cell wall and collapses the layers for their nutrient absorption. The pathogen also secretes effector molecules to escape from the plant defense mechanism. Later, plants activate their primary and secondary defense mechanisms; however, the pathogen induces virulence genes to escape the host immune responses. Obligate biotrophic fungi pathogenicity has not been fully understood at the molecular level because of the complex interaction, recognition, and signaling with the host. This review summarizes the mechanism of infection in the host, and immune response to emphasize the understanding of the biotrophic fungal biology and pathogenesis in crops. Thus, the detailed review will pave the way to design methods to overcome the resistance of biotrophic fungi and develop disease-free crops.

Barley endosomal MONENSIN SENSITIVITY1 is a target of the powdery mildew effector CSEP0162 and plays a role in plant immunity

Encasements formed around haustoria and biotrophic hyphae as well as hypersensitive reaction (HR) cell death are essential plant immune responses to filamentous pathogens. In this study we examine the components that may contribute to the absence of these responses in susceptible barley attacked by the powdery mildew fungus. We find that the effector CSEP0162 from this pathogen targets plant MONENSIN SENSITIVITY1 (MON1), which is important for the fusion of multivesicular bodies to their target membranes. Overexpression of CSEP0162 and silencing of barley MON1 both inhibit encasement formation. We find that the Arabidopsis ecotype No-0 has resistance to powdery mildew, and that this is partially dependent on MON1. Surprisingly, we find the MON1-dependent resistance in No-0 not only includes an encasement response, but also an effective HR. Similarly, silencing of MON1 in barley also blocks Mla3-mediated HR-based powdery mildew resistance. Our results indicate that MON1 is a vital plant immunity component, and we speculate that the barley powdery mildew fungus introduces the effector CSEP0162 to target MON1 and hence reduce encasement formation and HR.

The endosome as an effector target to mediate plant immunity?

This article comments on:

Liao W, Nielsen ME, Pedersen C, Xie W, Thordal-Christensen H. 2023. Barley endosomal MONENSIN SENSITIVITY1 is a target of the powdery mildew effector CSEP0162 and plays a role in plant immunity. Journal of Experimental Botany 74, 118–129.

Fungal Secondary Metabolites and Small RNAs Enhance Pathogenicity during Plant-Fungal Pathogen Interactions

Fungal plant pathogens use proteinaceous effectors as well as newly identified secondary metabolites (SMs) and small non-coding RNA (sRNA) effectors to manipulate the host plant's defense system via diverse plant cell compartments, distinct organelles, and many host genes. However, most molecular studies of plant-fungal interactions have focused on secreted effector proteins without exploring the possibly equivalent functions performed by fungal (SMs) and sRNAs, which are collectively known as "non-proteinaceous effectors". Fungal SMs have been shown to be generated throughout the plant colonization process, particularly in the early biotrophic stages of infection. The fungal repertoire of non-proteinaceous effectors has been broadened by the discovery of fungal sRNAs that specifically target plant genes involved in resistance and defense responses. Many RNAs, particularly sRNAs involved in gene silencing, have been shown to transmit bidirectionally between fungal pathogens and their hosts. However, there are no clear functional approaches to study the role of these SM and sRNA effectors. Undoubtedly, fungal SM and sRNA effectors are now a treasured land to seek. Therefore, understanding the role of fungal SM and sRNA effectors may provide insights into the infection process and identification of the interacting host genes that are targeted by these effectors. This review discusses the role of fungal SMs and sRNAs during plant-fungal interactions. It will also focus on the translocation of sRNA effectors across kingdoms, the application of cross-kingdom RNA interference in managing plant diseases and the tools that can be used to predict and study these non-proteinaceous effectors.

Comparative genomic analysis reveals contraction of gene families with putative roles in pathogenesis in the fungal boxwood pathogens Calonectria henricotiae and C. pseudonaviculata

Background

Boxwood blight disease caused by Calonectria henricotiae and C. pseudonaviculata is of ecological and economic significance in cultivated and native ecosystems worldwide. Prior research has focused on understanding the population genetic and genomic diversity of C. henricotiae and C. pseudonaviculata, but gene family evolution in the context of host adaptation, plant pathogenesis, and trophic lifestyle is poorly understood. This study applied bioinformatic and phylogenetic methods to examine gene family evolution in C. henricotiae, C. pseudonaviculata and 22 related fungi in the Nectriaceae that vary in pathogenic and saprobic (apathogenic) lifestyles.

Results

A total of 19,750 gene families were identified in the 24 genomes, of which 422 were rapidly evolving. Among the six Calonectria species, C. henricotiae and C. pseudonaviculata were the only species to experience high levels of rapid contraction of pathogenesis-related gene families (89% and 78%, respectively). In contrast, saprobic species Calonectria multiphialidica and C. naviculata, two of the closest known relatives of C. henricotiae and C. pseudonaviculata, showed rapid expansion of pathogenesis-related gene families.

Conclusions

Our results provide novel insight into gene family evolution within C. henricotiae and C. pseudonaviculata and suggest gene family contraction may have contributed to limited host-range expansion of these pathogens within the plant family Buxaceae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02035-4.

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.

The Venturia inaequalis effector repertoire is dominated by expanded families with predicted structural similarity, but unrelated sequence, to avirulence proteins from other plant-pathogenic fungi

BACKGROUND

Scab, caused by the biotrophic fungus Venturia inaequalis, is the most economically important disease of apples worldwide. During infection, V. inaequalis occupies the subcuticular environment, where it secretes virulence factors, termed effectors, to promote host colonization. Consistent with other plant-pathogenic fungi, many of these effectors are expected to be non-enzymatic proteins, some of which can be recognized by corresponding host resistance proteins to activate plant defences, thus acting as avirulence determinants. To develop durable control strategies against scab, a better understanding of the roles that these effector proteins play in promoting subcuticular growth by V. inaequalis, as well as in activating, suppressing, or circumventing resistance protein-mediated defences in apple, is required.

RESULTS

We generated the first comprehensive RNA-seq transcriptome of V. inaequalis during colonization of apple. Analysis of this transcriptome revealed five temporal waves of gene expression that peaked during early, mid, or mid-late infection. While the number of genes encoding secreted, non-enzymatic proteinaceous effector candidates (ECs) varied in each wave, most belonged to waves that peaked in expression during mid-late infection. Spectral clustering based on sequence similarity determined that the majority of ECs belonged to expanded protein families. To gain insights into function, the tertiary structures of ECs were predicted using AlphaFold2. Strikingly, despite an absence of sequence similarity, many ECs were predicted to have structural similarity to avirulence proteins from other plant-pathogenic fungi, including members of the MAX, LARS, ToxA and FOLD effector families. In addition, several other ECs, including an EC family with sequence similarity to the AvrLm6 avirulence effector from Leptosphaeria maculans, were predicted to adopt a KP6-like fold. Thus, proteins with a KP6-like fold represent another structural family of effectors shared among plant-pathogenic fungi.

CONCLUSIONS

Our study reveals the transcriptomic profile underpinning subcuticular growth by V. inaequalis and provides an enriched list of ECs that can be investigated for roles in virulence and avirulence. Furthermore, our study supports the idea that numerous sequence-unrelated effectors across plant-pathogenic fungi share common structural folds. In doing so, our study gives weight to the hypothesis that many fungal effectors evolved from ancestral genes through duplication, followed by sequence diversification, to produce sequence-unrelated but structurally similar proteins.

Evolution of pathogenicity in obligate fungal pathogens and allied genera

Obligate fungal pathogens (ascomycetes and basidiomycetes) and oomycetes are known to cause diseases in cereal crop plants. They feed on living cells and most of them have learned to bypass the host immune machinery. This paper discusses some of the factors that are associated with pathogenicity drawing examples from ascomycetes, basidiomycetes and oomycetes, with respect to their manifestation in crop plants. The comparisons have revealed a striking similarity in the three groups suggesting convergent pathways that have arisen from three lineages independently leading to an obligate lifestyle. This review has been written with the intent, that new information on adaptation strategies of biotrophs, modifications in pathogenicity strategies and population dynamics will improve current strategies for breeding with stable resistance.

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.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga's genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host-pathogen interactions.