Chromosome-level genome assembly and manually-curated proteome of model necrotroph Parastagonospora nodorum Sn15 reveals a genome-wide trove of candidate effector homologs, and redundancy of virulence-related functions within an accessory chromosome

Loss of the accessory chromosome converts a pathogenic tree-root fungus into a mutualistic endophyte

Some fungal accessory chromosomes (ACs) may contribute to virulence in plants. However, the mechanisms by which ACs determine specific traits associated with lifestyle transitions along a symbiotic continuum are not clear. Here we delineated the genetic divergence in two sympatric but considerably variable isolates (16B and 16W) of the poplar-associated fungus Stagonosporopsis rhizophilae. We identified a ∼0.6-Mb horizontally acquired AC in 16W that resulted in a mildly parasitic lifestyle in plants. Complete deletion of the AC (Δ16W) significantly altered the fungal phenotype. Specifically, Δ16W was morphologically more similar to 16B, showed enhanced melanization, and established beneficial interactions with poplar plants, thereby acting as a dark septate endophyte. RNA sequencing (RNA-seq) analysis showed that AC loss induced the upregulation of genes related to root colonization and biosynthesis of indole acetic acid and melanin. We observed that the AC maintained a more open status of chromatin across the genome, indicating an impressive remodeling of cis-regulatory elements upon AC loss, which potentially enhanced symbiotic effectiveness. We demonstrated that the symbiotic capacities were non-host-specific through comparable experiments on Triticum- and Arabidopsis-fungus associations. Furthermore, the three isolates generated symbiotic interactions with a nonvascular liverwort. In summary, our study suggests that the AC is a suppressor of symbiosis and provides insights into the underlying mechanisms of mutualism with vascular plants in the absence of traits encoded by the AC. We speculate that AC-situated effectors and other potential secreted molecules may have evolved to specifically target vascular plants and promote mild virulence.

The Necrotrophic Pathogen Parastagonospora nodorum Is a Master Manipulator of Wheat Defense

Parastagonospora nodorum is a necrotrophic pathogen of wheat that is particularly destructive in major wheat-growing regions of the United States, northern Europe, Australia, and South America. P. nodorum secretes necrotrophic effectors that target wheat susceptibility genes to induce programmed cell death (PCD), resulting in increased colonization of host tissue and, ultimately, sporulation to complete its pathogenic life cycle. Intensive research over the last two decades has led to the functional characterization of five proteinaceous necrotrophic effectors, SnTox1, SnToxA, SnTox267, SnTox3, and SnTox5, and three wheat susceptibility genes, Tsn1, Snn1, and Snn3D-1. Functional characterization has revealed that these effectors, in addition to inducing PCD, have additional roles in pathogenesis, including chitin binding that results in protection from wheat chitinases, blocking defense response signaling, and facilitating plant colonization. There are still large gaps in our understanding of how this necrotrophic pathogen is successfully manipulating wheat defense to complete its life cycle. This review summarizes our current knowledge, identifies knowledge gaps, and provides a summary of well-developed tools and resources currently available to study the P. nodorum-wheat interaction, which has become a model for necrotrophic specialist interactions. Further functional characterization of the effectors involved in this interaction and work toward a complete understanding of how P. nodorum manipulates wheat defense will provide fundamental knowledge about this and other necrotrophic interactions. Additionally, a broader understanding of this interaction will contribute to the successful management of Septoria nodorum blotch disease on wheat. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

CoreDetector: a flexible and efficient program for core-genome alignment of evolutionary diverse genomes

MOTIVATION

Whole genome alignment of eukaryote species remains an important method for the determination of sequence and structural variations and can also be used to ascertain the representative non-redundant core-genome sequence of a population. Many whole genome alignment tools were first developed for the more mature analysis of prokaryote species with few current tools containing the functionality to process larger genomes of eukaryotes as well as genomes of more divergent species. In addition, the functionality of these tools becomes computationally prohibitive due to the significant compute resources needed to handle larger genomes.

RESULTS

In this research, we present CoreDetector, an easy-to-use general-purpose program that can align the core-genome sequences for a range of genome sizes and divergence levels. To illustrate the flexibility of CoreDetector, we conducted alignments of a large set of closely related fungal pathogen and hexaploid wheat cultivar genomes as well as more divergent fly and rodent species genomes. In all cases, compared to existing multiple genome alignment tools, CoreDetector exhibited improved flexibility, efficiency, and competitive accuracy in tested cases.

AVAILABILITY AND IMPLEMENTATION

CoreDetector was developed in the cross platform, and easily deployable, Java language. A packaged pipeline is readily executable in a bash terminal without any external need for Perl or Python environments. Installation, example data, and usage instructions for CoreDetector are freely available from https://github.com/mfruzan/CoreDetector.

The First Telomere-to-Telomere Chromosome-Level Genome Assembly of Stagonospora tainanensis Causing Sugarcane Leaf Blight

The sexual morph Leptosphaeria taiwanensis Yen and Chi and its asexual morph Stagonospora tainanensis W. H. Hsieh is an important necrotrophic fungal phytopathogen, which causes sugarcane leaf blight, resulting in loss of cane tonnage and sucrose in susceptible sugarcane varieties. Decoding the genome and understanding of the basis of virulence is vitally important for devising effective disease control strategies. Here, we present a 38.25-Mb high-quality genome assembly of S. tainanensis strain StFZ01, denovo assembled with 10.19 Gb Nanopore sequencing long reads (~267×) and 3.82 Gb Illumina short reads (~100×). The genome assembly consists of 12 contigs with N50 of 2.86 Mb of which 5 belong to the telomere to telomere (T2T) chromosome. It contains 13.20% repeat sequences, 12,543 proteins, and 12,206 protein-coding genes with the BUSCO completeness 99.18% at fungi (n = 758) and 99.87% at ascomycota (n = 1706), indicating the high accuracy and completeness of our gene annotations. The virulence analysis in silico revealed the presence of 2379 PHIs, 599 CAZys, 248 membrane transport proteins, 191 cytochrome P450 enzymes, 609 putative secreted proteins, and 333 effectors in the StFZ01 genome. The genomic resources presented here will not only be helpful for development of specific molecular marker and diagnosis technique, population genetics, molecular taxonomy, and disease managements, it can also provide a significant precise genomic reference for investigating the ascomycetous genome, the necrotrophic lifestyle, and pathogenicity in the future.

A global pangenome for the wheat fungal pathogen Pyrenophora tritici-repentis and prediction of effector protein structural homology

The adaptive potential of plant fungal pathogens is largely governed by the gene content of a species, consisting of core and accessory genes across the pathogen isolate repertoire. To approximate the complete gene repertoire of a globally significant crop fungal pathogen, a pan genomic analysis was undertaken for Pyrenophora tritici-repentis (Ptr), the causal agent of tan (or yellow) spot disease in wheat. In this study, 15 new Ptr genomes were sequenced, assembled and annotated, including isolates from three races not previously sequenced. Together with 11 previously published Ptr genomes, a pangenome for 26 Ptr isolates from Australia, Europe, North Africa and America, representing nearly all known races, revealed a conserved core-gene content of 57 % and presents a new Ptr resource for searching natural homologues (orthologues not acquired by horizontal transfer from another species) using remote protein structural homology. Here, we identify for the first time a non-synonymous mutation in the Ptr necrotrophic effector gene ToxB, multiple copies of the inactive toxb within an isolate, a distant natural Pyrenophora homologue of a known Parastagonopora nodorum necrotrophic effector (SnTox3), and clear genomic break points for the ToxA effector horizontal transfer region. This comprehensive genomic analysis of Ptr races includes nine isolates sequenced via long read technologies. Accordingly, these resources provide a more complete representation of the species, and serve as a resource to monitor variations potentially involved in pathogenicity.

Association mapping reveals a reciprocal virulence/avirulence locus within diverse US Pyrenophora teres f. maculata isolates

Background

Spot form net blotch (SFNB) caused by the necrotrophic fungal pathogen Pyrenophora teres f. maculata (Ptm) is an economically important disease of barley that also infects wheat. Using genetic analysis to characterize loci in Ptm genomes associated with virulence or avirulence is an important step to identify pathogen effectors that determine compatible (virulent) or incompatible (avirulent) interactions with cereal hosts. Association mapping (AM) is a powerful tool for detecting virulence loci utilizing phenotyping and genotyping data generated for natural populations of plant pathogenic fungi.

Results

Restriction-site associated DNA genotyping-by-sequencing (RAD-GBS) was used to generate 4,836 single nucleotide polymorphism (SNP) markers for a natural population of 103 Ptm isolates collected from Idaho, Montana and North Dakota. Association mapping analyses were performed utilizing the genotyping and infection type data generated for each isolate when challenged on barley seedlings of thirty SFNB differential barley lines. A total of 39 marker trait associations (MTAs) were detected across the 20 barley lines corresponding to 30 quantitative trait loci (QTL); 26 novel QTL and four that were previously mapped in Ptm biparental populations. These results using diverse US isolates and barley lines showed numerous barley-Ptm genetic interactions with seven of the 30 Ptm virulence/avirulence loci falling on chromosome 3, suggesting that it is a reservoir of diverse virulence effectors. One of the loci exhibited reciprocal virulence/avirulence with one haplotype predominantly present in isolates collected from Idaho increasing virulence on barley line MXB468 and the alternative haplotype predominantly present in isolates collected from North Dakota and Montana increasing virulence on barley line CI9819.

Conclusions

Association mapping provided novel insight into the host pathogen genetic interactions occurring in the barley-Ptm pathosystem. The analysis suggests that chromosome 3 of Ptm serves as an effector reservoir in concordance with previous reports for Pyrenophora teres f. teres, the causal agent of the closely related disease net form net blotch. Additionally, these analyses identified the first reported case of a reciprocal pathogen virulence locus. However, further investigation of the pathosystem is required to determine if multiple genes or alleles of the same gene are responsible for this genetic phenomenon.

Genome Analysis of the Broad Host Range Necrotroph Nalanthamala psidii Highlights Genes Associated With Virulence

Guava wilt disease is caused by the fungus Nalanthamala psidii. The wilt disease results in large-scale destruction of orchards in South Africa, Taiwan, and several Southeast Asian countries. De novo assembly, annotation, and in-depth analysis of the N. psidii genome were carried out to facilitate the identification of characteristics associated with pathogenicity and pathogen evolution. The predicted secretome revealed a range of CAZymes, proteases, lipases and peroxidases associated with plant cell wall degradation, nutrient acquisition, and disease development. Further analysis of the N. psidii carbohydrate-active enzyme profile exposed the broad-spectrum necrotrophic lifestyle of the pathogen, which was corroborated by the identification of putative effectors and secondary metabolites with the potential to induce tissue necrosis and cell surface-dependent immune responses. Putative regulatory proteins including transcription factors and kinases were identified in addition to transporters potentially involved in the secretion of secondary metabolites. Transporters identified included important ABC and MFS transporters involved in the efflux of fungicides. Analysis of the repetitive landscape and the detection of mechanisms linked to reproduction such as het and mating genes rendered insights into the biological complexity and evolutionary potential of N. psidii as guava pathogen. Hence, the assembly and annotation of the N. psidii genome provided a valuable platform to explore the pathogenic potential and necrotrophic lifestyle of the guava wilt pathogen.

EffectorP 3.0: Prediction of Apoplastic and Cytoplasmic Effectors in Fungi and Oomycetes

Many fungi and oomycete species are devasting plant pathogens. These eukaryotic filamentous pathogens secrete effector proteins to facilitate plant infection. Fungi and oomycete pathogens have diverse infection strategies and their effectors generally do not share sequence homology. However, they occupy similar host environments, either the plant apoplast or plant cytoplasm, and, therefore, may share some unifying properties based on the requirements of these host compartments. Here, we exploit these biological signals and present the first classifier (EffectorP 3.0) that uses two machine-learning models: one trained on apoplastic effectors and one trained on cytoplasmic effectors. EffectorP 3.0 accurately predicts known apoplastic and cytoplasmic effectors in fungal and oomycete secretomes with low estimated false-positive rates of 3 and 8%, respectively. Cytoplasmic effectors have a higher proportion of positively charged amino acids, whereas apoplastic effectors are enriched for cysteine residues. The combination of fungal and oomycete effectors in training leads to a higher number of predicted cytoplasmic effectors in biotrophic fungi. EffectorP 3.0 expands predicted effector repertoires beyond small, cysteine-rich secreted proteins in fungi and RxLR-motif containing secreted proteins in oomycetes. We show that signal peptide prediction is essential for accurate effector prediction, because EffectorP 3.0 recognizes a cytoplasmic signal also in intracellular, nonsecreted proteins.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Biology and molecular interactions of Parastagonospora nodorum blotch of wheat

Parastagonospora nodorum is one of the important necrotrophic pathogens of wheat which causes severe economical loss to crop yield. So far, a number of effectors of Parastagonospora nodorum origin and their target interacting genes on the host plant have been characterized. Since targeting effector-sensitive gene carefully can be helpful in breeding for resistance. Therefore, constant efforts are required to further characterize the effectors, their interacting genes, and underlying biochemical pathways. Furthermore, to develop effective counter-strategies against emerging diseases, continuous efforts are required to determine the qualitative resistance that demands to screen of diverse genotypes for host resistance. Stagonospora nodorum blotch also refers to as Stagonospora glume blotch and leaf is caused by Parastagonospora nodorum. The pathogen deploys necrotrophic effectors for the establishment and development on wheat plants. The necrotrophic effectors and their interaction with host receptors lead to the establishment of infection on leaves and extensive lesions formation which either results in host cell death or suppression/activation of host defence mechanisms. The wheat Stagonospora nodorum interaction involves a set of nine host gene-necrotrophic effector interactions. Out of these, Snn1-SnTox1, Tsn1-SnToxA and Snn-SnTox3 are one of the most studied interaction, due to its role in the suppression of reactive oxygen species production, regulating the cytokinin content through ethylene-dependent wayduring initial infection stage. Further, although the molecular basis is not fully unveiled, these effectors regulate the redox state and influence the ethylene biosynthesis in infected wheat plants. Here, we have discussed the biology of the wheat pathogen Parastagonospora nodorum, role of its necrotrophic effectors and their interacting sensitivity genes on the redox state, how they hijack the resistance mechanisms, hormonal regulated immunity and other signalling pathways in susceptible wheat plants. The information generated from effectors and their corresponding sensitivity genes and other biological processes could be utilized effectively for disease management strategies.