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

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.

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.

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.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01442-9.

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.

The unsung roles of microbial secondary metabolite effectors in the plant disease cacophony

Plants counter disease with an array of responses to styme pathogen ingress. In contrast to this cacophony, plant pathogens orchestrate a finely tuned repertoire of virulence mechanisms in their attempt to cause disease. One such example is the production of secondary metabolite effectors (SMEs). Despite many attempts to functionally categorize SMEs, their many roles in plant disease have proven they march to the beat of their producer's drum. Some lesser studied features of SMEs in plant disease include self-resistance (SR) and manipulation of the microbiome to enhance pathogen virulence. SR can be accomplished in three general compositions, with the first being the transport of the SME to a benign location; the second being modification of the SME so it cannot harm the producer; and the third being metabolic regulation of the SME or the producer homolog of the SME target. SMEs may also play an interlude prior to disease by shaping the plant microbial community, allowing producers to better establish themselves. Taken together, SMEs are integral players in the phytopathology canon.

Genome-wide comparison deciphers lifestyle adaptation and glass biodeterioration property of Curvularia eragrostidis C52

Glass biodeterioration by fungi has caused irreversible damage to valuable glass materials such as cultural heritages and optical devices. To date, knowledge about metabolic potential and genomic profile of biodeteriorative fungi is still scarce. Here, we report for the first time the whole genome sequence of Curvularia eragrostidis C52 that strongly degraded silica-based glasses coated with fluorine and hafnium, as expressed by the hyphal surface coverage of 46.16 ± 3.3% and reduced light transmission of 50.93 ± 1.45%. The genome of C. eragrostidis C52 is 36.9 Mb long with a GC content of 52.1% and contains 14,913 protein-coding genes, which is the largest genome ever recorded in the genus Curvularia. Phylogenomic analysis revealed C. eragrostidis C52 formed a distinct cluster with Curvularia sp. IFB-Z10 and was not evolved from compared genomes. Genome-wide comparison showed that strain C52 harbored significantly higher proportion of proteins involved in carbohydrate-active enzymes, peptidases, secreted proteins, and transcriptional factors, which may be potentially attributed to a lifestyle adaptation. Furthermore, 72 genes involved in the biosynthesis of 6 different organic acids were identified and expected to be crucial for the fungal survival in the glass environment. To form biofilm against stress, the fungal strain utilized 32 genes responsible for exopolysaccharide production. These findings will foster a better understanding of the biology of C. eragrostidis and the mechanisms behind fungal biodeterioration in the future.

Comparative genomics reveals a dynamic genome evolution in the ectomycorrhizal milk-cap (Lactarius) mushrooms

Ectomycorrhizal fungi play a key role in forests by establishing mutualistic symbioses with woody plants. Genome analyses have identified conserved symbiosis-related traits among ectomycorrhizal fungal species, but the molecular mechanisms underlying host specificity remain poorly known. We sequenced and compared the genomes of seven species of milk-cap fungi (Lactarius, Russulales) with contrasting host specificity. We also compared these genomes with those of symbiotic and saprotrophic Russulales species, aiming to identify genes involved in their ecology and host specificity. The size of Lactarius genomes is significantly larger than other Russulales species, owing to a massive accumulation of transposable elements and duplication of dispensable genes. As expected, their repertoire of genes coding for plant cell wall-degrading enzymes is restricted, but they retained a substantial set of genes involved in microbial cell wall degradation. Notably, Lactarius species showed a striking expansion of genes encoding proteases, such as secreted ectomycorrhiza-induced sedolisins. A high copy number of genes coding for small secreted LysM proteins and Lactarius-specific lectins were detected, which may be linked to host specificity. This study revealed a large diversity in the genome landscapes and gene repertoires within Russulaceae. The known host specificity of Lactarius symbionts may be related to mycorrhiza-induced species-specific genes, including secreted sedolisins.

Organic acids and glucose prime late-stage fungal biotrophy in maize

Many plant-associated fungi are obligate biotrophs that depend on living hosts to proliferate. However, little is known about the molecular basis of the biotrophic lifestyle, despite the impact of fungi on the environment and food security. In this work, we show that combinations of organic acids and glucose trigger phenotypes that are associated with the late stage of biotrophy for the maize pathogen Ustilago maydis. These phenotypes include the expression of a set of effectors normally observed only during biotrophic development, as well as the formation of melanin associated with sporulation in plant tumors. U. maydis and other hemibiotrophic fungi also respond to a combination of carbon sources with enhanced proliferation. Thus, the response to combinations of nutrients from the host may be a conserved feature of fungal biotrophy.

Beyond Nuclear Ribosomal DNA Sequences: Evolution, Taxonomy, and Closest Known Saprobic Relatives of Powdery Mildew Fungi (Erysiphaceae) Inferred From Their First Comprehensive Genome-Scale Phylogenetic Analyses

Powdery mildew fungi (Erysiphaceae), common obligate biotrophic pathogens of many plants, including important agricultural and horticultural crops, represent a monophyletic lineage within the Ascomycota. Within the Erysiphaceae, molecular phylogenetic relationships and DNA-based species and genera delimitations were up to now mostly based on nuclear ribosomal DNA (nrDNA) phylogenies. This is the first comprehensive genome-scale phylogenetic analysis of this group using 751 single-copy orthologous sequences extracted from 24 selected powdery mildew genomes and 14 additional genomes from Helotiales, the fungal order that includes the Erysiphaceae. Representative genomes of all powdery mildew species with publicly available whole-genome sequencing (WGS) data that were of sufficient quality were included in the analyses. The 24 powdery mildew genomes included in the analysis represented 17 species belonging to eight out of 19 genera recognized within the Erysiphaceae. The epiphytic genera, all but one represented by multiple genomes, belonged each to distinct, well-supported lineages. Three hemiendophytic genera, each represented by a single genome, together formed the hemiendophytic lineage. Out of the 14 other taxa from the Helotiales, Arachnopeziza araneosa, a saprobic species, was the only taxon that grouped together with the 24 genome-sequenced powdery mildew fungi in a monophyletic clade. The close phylogenetic relationship between the Erysiphaceae and Arachnopeziza was revealed earlier by a phylogenomic study of the Leotiomycetes. Further analyses of powdery mildew and Arachnopeziza genomes may discover signatures of the evolutionary processes that have led to obligate biotrophy from a saprobic way of life. A separate phylogeny was produced using the 18S, 5.8S, and 28S nrDNA sequences of the same set of powdery mildew specimens and compared to the genome-scale phylogeny. The nrDNA phylogeny was largely congruent to the phylogeny produced using 751 orthologs. This part of the study has revealed multiple contamination and other quality issues in some powdery mildew genomes. We recommend that the presence of 28S, internal transcribed spacer (ITS), and 18S nrDNA sequences in powdery mildew WGS datasets that are identical to those determined by Sanger sequencing should be used to assess the quality of assemblies, in addition to the commonly used Benchmarking Universal Single-Copy Orthologs (BUSCO) values.

Infection Strategies and Pathogenicity of Biotrophic Plant Fungal Pathogens

Biotrophic plant pathogenic fungi are widely distributed and are among the most damaging pathogenic organisms of agriculturally important crops responsible for significant losses in quality and yield. However, the pathogenesis of obligate parasitic pathogenic microorganisms is still under investigation because they cannot reproduce and complete their life cycle on an artificial medium. The successful lifestyle of biotrophic fungal pathogens depends on their ability to secrete effector proteins to manipulate or evade plant defense response. By integrating genomics, transcriptomics, and effectoromics, insights into how the adaptation of biotrophic plant fungal pathogens adapt to their host populations can be gained. Efficient tools to decipher the precise molecular mechanisms of rust–plant interactions, and standardized routines in genomics and functional pipelines have been established and will pave the way for comparative studies. Deciphering fungal pathogenesis not only allows us to better understand how fungal pathogens infect host plants but also provides valuable information for plant diseases control, including new strategies to prevent, delay, or inhibit fungal development. Our review provides a comprehensive overview of the efforts that have been made to decipher the effector proteins of biotrophic fungal pathogens and demonstrates how rapidly research in the field of obligate biotrophy has progressed.

An interolog-based barley interactome as an integration framework for immune signaling

The barley MLA nucleotide-binding leucine-rich-repeat (NLR) receptor and its orthologs confer recognition specificity to many fungal diseases, including powdery mildew, stem-, and stripe rust. We used interolog inference to construct a barley protein interactome (Hordeum vulgare predicted interactome, HvInt) comprising 66,133 edges and 7,181 nodes, as a foundation to explore signaling networks associated with MLA. HvInt was compared with the experimentally validated Arabidopsis interactome of 11,253 proteins and 73,960 interactions, verifying that the 2 networks share scale-free properties, including a power-law distribution and small-world network. Then, by successive layering of defense-specific "omics" datasets, HvInt was customized to model cellular response to powdery mildew infection. Integration of HvInt with expression quantitative trait loci (eQTL) enabled us to infer disease modules and responses associated with fungal penetration and haustorial development. Next, using HvInt and infection-time-course RNA sequencing of immune signaling mutants, we assembled resistant and susceptible subnetworks. The resulting differentially coexpressed (resistant - susceptible) interactome is essential to barley immunity, facilitates the flow of signaling pathways and is linked to mildew resistance locus a (Mla) through trans eQTL associations. Lastly, we anchored HvInt with new and previously identified interactors of the MLA coiled coli + nucleotide-binding domains and extended these to additional MLA alleles, orthologs, and NLR outgroups to predict receptor localization and conservation of signaling response. These results link genomic, transcriptomic, and physical interactions during MLA-specified immunity.

Virulence Structure and Genetic Diversity of Blumeria graminis f. sp. avenae from Different Regions of Europe

The structure and dynamics of changes in pathogen populations can be analysed by assessing the level of virulence and genetic diversity. The aim of the present study was to determine the diversity of Blumeria graminis f. sp. avenae populations. Diversity and virulence of B. graminis f. sp. avenae was assessed based on 80 single-spore isolates collected in different European countries such as Poland (40 isolates), Germany (10), Finland (10), Czech Republic (10) and Ireland (10) using ISSR (Inter-Simple Sequence Repeats) and SCoT (Start Codon Targeted) markers. This work demonstrated differences in virulence of B. graminis f. sp. avenae isolates sampled from different countries. Molecular analysis showed that both systems were useful for assessing genetic diversity, but ISSR markers were superior and generated more polymorphic products, as well as higher PIC and RP values. UPMGA and PCoA divided the isolates into groups corresponding with their geographical origin. In conclusion, the low level of genetic differentiation of the analysed isolates has suggested that the evolution of B. graminis f. sp. Avenae population is slow, and thus the evolutionary potential of the pathogen is low. This work paves the way for future studies on B. graminis f. sp. Avenae population structure and dynamics based on genetic variability.

Population genomics of the grapevine pathogen Eutypa lata reveals evidence for population expansion and intraspecific differences in secondary metabolite gene clusters

Eutypa dieback of grapevine is an important disease caused by the generalist Ascomycete fungus Eutypa lata. Despite the relevance of this species to the global wine industry, its genomic diversity remains unknown, with only a single publicly available genome assembly. Whole-genome sequencing and comparative genomics was performed on forty Australian E. lata isolates to understand the genome evolution, adaptation, population size and structure of these isolates. Phylogenetic and linkage disequilibrium decay analyses provided evidence of extensive gene flow through sexual recombination between isolates obtained from different geographic locations and hosts. Investigation of the genetic diversity of these isolates suggested rapid population expansion, likely as a consequence of the recent growth of the Australian wine industry. Genomic regions affected by selective sweeps were shown to be enriched for genes associated with secondary metabolite clusters and included genes encoding proteins with a role in nutrient acquisition, degradation of host cell wall and metal and drug resistance, suggesting recent adaptation to both abiotic factors and potentially host genotypes. Genome synteny analysis using long-read genome assemblies showed significant intraspecific genomic plasticity with extensive chromosomal rearrangements impacting the secondary metabolite production potential of this species. Finally, k-mer based GWAS analysis identified a potential locus associated with mycelia recovery in canes of Vitis vinifera that will require further investigations.

Eutypa dieback of grapevine, caused by the Ascomycete fungus Eutypa lata, is responsible for significant economic losses to the wine industry. Despite the worldwide prevalence of this pathogen, its genomic diversity remains unknown, with only a single publicly available genome assembly. This knowledge gap was addressed by performing whole-genome sequencing of 40 E. lata isolates sourced from different hosts and geographical locations around Australia. Investigation of the genetic diversity of this population showed a high degree of gene-flow and sexual recombination as well as demographic expansion. Through the inspection of signatures of selective sweeps, repeat-mediated chromosomal rearrangements, and pan-genomic elements, it was shown that this species has a highly dynamic secondary metabolite production potential that could have important implications for its pathogenicity and lifestyle. In addition, application of a k-mer based GWAS methodology, identified a potential locus associated with the growth of this species within canes of Vitis vinifera.

First Draft Genome Assemblies of Pleochaeta shiraiana and Phyllactinia moricola, Two Tree-Parasitic Powdery Mildew Fungi with Hemiendophytic Mycelia

Powdery mildew fungi (Erysiphaceae) are widespread obligate biotrophic plant pathogens. Thus, applying genetic and omics approaches to study these fungi remains a major challenge, particularly for species with hemiendophytic mycelium. These belong to a distinct phylogenetic lineage within the family Erysiphaceae. To date, only a single draft genome assembly is available for this clade, obtained for Leveillula taurica. Here, we generated the first draft genome assemblies of Pleochaeta shiraiana and Phyllactinia moricola, two tree-parasitic powdery mildew species with hemiendophytic mycelium, representing two genera that have not yet been investigated with genomics tools. The Pleochaeta shiraiana assembly was 96,769,103 bp in length and consisted of 14,447 scaffolds, and the Phyllactinia moricola assembly was 180,382,532 bp in length on 45,569 scaffolds. Together with the draft genome of L. taurica, these resources will be pivotal for understanding the molecular basis of the lifestyle of these fungi, which is unique within the family Erysiphaceae.

A myxozoan genome reveals mosaic evolution in a parasitic cnidarian

Background

Parasite evolution has been conceptualized as a process of genetic loss and simplification. Contrary to this model, there is evidence of expansion and conservation of gene families related to essential functions of parasitism in some parasite genomes, reminiscent of widespread mosaic evolution—where subregions of a genome have different rates of evolutionary change. We found evidence of mosaic genome evolution in the cnidarian Myxobolus honghuensis, a myxozoan parasite of fish, with extremely simple morphology.

Results

We compared M. honghuensis with other myxozoans and free-living cnidarians, and determined that it has a relatively larger myxozoan genome (206 Mb), which is less reduced and less compact due to gene retention, large introns, transposon insertion, but not polyploidy. Relative to other metazoans, the M. honghuensis genome is depleted of neural genes and has only the simplest animal immune components. Conversely, it has relatively more genes involved in stress resistance, tissue invasion, energy metabolism, and cellular processes compared to other myxozoans and free-living cnidarians. We postulate that the expansion of these gene families is the result of evolutionary adaptations to endoparasitism. M. honghuensis retains genes found in free-living Cnidaria, including a reduced nervous system, myogenic components, ANTP class Homeobox genes, and components of the Wnt and Hedgehog pathways.

Conclusions

Our analyses suggest that the M. honghuensis genome evolved as a mosaic of conservative, divergent, depleted, and enhanced genes and pathways. These findings illustrate that myxozoans are not as genetically simple as previously regarded, and the evolution of some myxozoans is driven by both genomic streamlining and expansion.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01249-8.

In Silico Predictions of Ecological Plasticity Mediated by Protein Family Expansions in Early-Diverging Fungi

Early-diverging fungi (EDF) are ubiquitous and versatile. Their diversity is reflected in their genome sizes and complexity. For instance, multiple protein families have been reported to expand or disappear either in particular genomes or even whole lineages. The most commonly mentioned are CAZymes (carbohydrate-active enzymes), peptidases and transporters that serve multiple biological roles connected to, e.g., metabolism and nutrients intake. In order to study the link between ecology and its genomic underpinnings in a more comprehensive manner, we carried out a systematic in silico survey of protein family expansions and losses among EDF with diverse lifestyles. We found that 86 protein families are represented differently according to EDF ecological features (assessed by median count differences). Among these there are 19 families of proteases, 43 CAZymes and 24 transporters. Some of these protein families have been recognized before as serine and metallopeptidases, cellulases and other nutrition-related enzymes. Other clearly pronounced differences refer to cell wall remodelling and glycosylation. We hypothesize that these protein families altogether define the preliminary fungal adaptasome. However, our findings need experimental validation. Many of the protein families have never been characterized in fungi and are discussed in the light of fungal ecology for the first time.

Comparative genome analysis of plant ascomycete fungal pathogens with different lifestyles reveals distinctive virulence strategies

Background

Pathogens have evolved diverse lifestyles and adopted pivotal new roles in both natural ecosystems and human environments. However, the molecular mechanisms underlying their adaptation to new lifestyles are obscure. Comparative genomics was adopted to determine distinct strategies of plant ascomycete fungal pathogens with different lifestyles and to elucidate their distinctive virulence strategies.

Results

We found that plant ascomycete biotrophs exhibited lower gene gain and loss events and loss of CAZyme-encoding genes involved in plant cell wall degradation and biosynthesis gene clusters for the production of secondary metabolites in the genome. Comparison with the candidate effectome detected distinctive variations between plant biotrophic pathogens and other groups (including human, necrotrophic and hemibiotrophic pathogens). The results revealed the biotroph-specific and lifestyle-conserved candidate effector families. These data have been configured in web-based genome browser applications for public display (http://lab.malab.cn/soft/PFPG). This resource allows researchers to profile the genome, proteome, secretome and effectome of plant fungal pathogens.

Conclusions

Our findings demonstrated different genome evolution strategies of plant fungal pathogens with different lifestyles and explored their lifestyle-conserved and specific candidate effectors. It will provide a new basis for discovering the novel effectors and their pathogenic mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08165-1.