Evolution of the Toxb Gene in Pyrenophora tritici-repentis and Related Species.

Pyrenophora tritici-repentis (tan spot) is a destructive foliar pathogen of wheat with global impact. This ascomycete fungus possesses a highly plastic open pangenome shaped by the gain and loss of effector genes. This study investigated the allelic variations in the chlorosis-encoding gene ToxB across 422 isolates representing all identified pathotypes and worldwide origins. To gain better insights into ToxB evolution, we examined its presence and variability in other Pyrenophora spp. A ToxB haplotype network was constructed, revealing the evolutionary relationships of this gene (20 haplotypes) across four Pyrenophora species. Notably, toxb, the homolog of ToxB, was detected for the first time in the barley pathogen Pyrenophora teres. The ToxB/toxb genes display evidence of selection that is characterized by loss of function, duplication, and diverse mutations. Within the ToxB/toxb open reading frame, 72 mutations were identified, including 14 synonymous, 55 nonsynonymous, and 3 indel mutations. Remarkably, a, ∼5.6-kb Copia-like retrotransposon, named Copia-1_Ptr, was found inserted in the toxb gene of a race 3 isolate. This insert disrupted the ToxB gene's function, a first case of effector gene disruption by a transposable element in P. tritici-repentis. Additionally, a microsatellite with 25 nucleotide repeats (0 to 10) in the upstream region of ToxB suggested a potential mechanism influencing ToxB expression and regulation. Exploring ToxB-like protein distribution in other ascomycetes revealed the presence of ToxB-like proteins in 19 additional species, including the Leotiomycetes class for the first time. The presence/absence pattern of ToxB-like proteins defied species relatedness compared with a phylogenetic tree, suggesting a past horizontal gene transfer event during the evolution of the ToxB gene. [Formula: see text] Copyright © 2024 His Majesty the King in Right of Canada, as represented by the Minister of Agriculture and Agri-Food. 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.

Spatiotemporal patterns of wheat response to Pyrenophora tritici-repentis in asymptomatic regions revealed by transcriptomic and X-ray fluorescence microscopy analyses.

Pathogen attacks elicit dynamic and widespread molecular responses in plants. While our understanding of plant responses has advanced considerably, little is known of the molecular responses in the asymptomatic 'green' regions adjoining lesions. Here, we explore gene expression data and high-resolution elemental imaging to report the spatiotemporal changes in the asymptomatic green region of susceptible and moderately resistant wheat cultivars infected with a necrotrophic fungal pathogen, Pyrenophora tritici-repentis. We show, with improved spatiotemporal resolution, that calcium oscillations are modified in the susceptible cultivar, resulting in 'frozen' host defence signals at the mature disease stage, and silencing of the host's recognition and defence mechanisms that would otherwise protect it from further attacks. In contrast, calcium accumulation and a heightened defence response were observed in the moderately resistant cultivar in the later stage of disease development. Furthermore, in the susceptible interaction, the asymptomatic green region was unable to recover after disease disruption. Our targeted sampling technique also enabled detection of eight previously predicted proteinaceous effectors in addition to the known ToxA effector. Collectively, our results highlight the benefits of spatially resolved molecular analysis and nutrient mapping to provide high-resolution spatiotemporal snapshots of host-pathogen interactions, paving the way for disentangling complex disease interactions in plants.

A unique resistance mechanism is associated with RBgh2 barley powdery mildew adult plant resistance.

KEY MESSAGE: Gene expression at the RBgh2 locus indicates involvement in cAMP/G-protein-coupled signalling and innate immunity in barley powdery mildew adult plant resistance. Barley powdery mildew is a globally significant disease, responsible for reduced grain yield and quality. A major effect adult plant resistance gene, RBgh2, was previously found in a landrace from Azerbaijan. The atypical phenotype suggested different underlying genetic factors compared to conventional resistance genes and to investigate this, genome-wide gene expression was compared between sets of heterogeneous doubled haploids. RBgh2 resistance is recessive and induces both temporary genome-wide gene expression changes during powdery mildew infection together with constitutive changes, principally at the RBgh2 locus. Defence-related genes significantly induced included homologues of genes associated with innate immunity and pathogen recognition. Intriguingly, RBgh2 resistance does not appear to be dependent on salicylic acid signalling, a key pathway in plant resistance to biotrophs. Constitutive co-expression of resistance gene homologues was evident at the 7HS RBgh2 locus, while no expression was evident for a 6-transmembrane gene, predicted in silico to contain both G-protein- and calmodulin-binding domains. The gene was disrupted at the 5' end, and G-protein-binding activity was suppressed. RBgh2 appears to operate through a unique mechanism that co-opts elements of innate immunity.

A Revised Nomenclature for ToxA Haplotypes Across Multiple Fungal Species.

ToxA is one of the most studied proteinaceous necrotrophic effectors produced by plant pathogens. It has been identified in four pathogens (Pyrenophora tritici-repentis, Parastagonospora nodorum, Parastagonospora pseudonodorum [formerly Parastagonospora avenaria f. sp. tritici], and Bipolaris sorokiniana) causing leaf spot diseases on cereals worldwide. To date, 24 different ToxA haplotypes have been identified. Some P. tritici-repentis and related species also express ToxB, another small protein necrotrophic effector. We present here a revised and standardized nomenclature for these effectors, which could be extended to other poly-haplotypic genes found across multiple species.

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.

A Conserved Hypothetical Gene Is Required but Not Sufficient for Ptr ToxC Production in Pyrenophora tritici-repentis.

The fungus Pyrenophora tritici-repentis causes tan spot, an important foliar disease of wheat worldwide. The fungal pathogen produces three necrotrophic effectors, namely Ptr ToxA, Ptr ToxB, and Ptr ToxC to induce necrosis or chlorosis in wheat. Both Ptr ToxA and Ptr ToxB are proteins, and their encoding genes have been cloned. Ptr ToxC was characterized as a low-molecular weight molecule 20 years ago but the one or more genes controlling its production in P. tritici-repentis are unknown. Here, we report the genetic mapping, molecular cloning, and functional analysis of a fungal gene that is required for Ptr ToxC production. The genetic locus controlling the production of Ptr ToxC, termed ToxC, was mapped to a subtelomeric region using segregating biparental populations, genome sequencing, and association analysis. Additional marker analysis further delimited ToxC to a 173-kb region. The predicted genes in the region were examined for presence/absence polymorphism in different races and isolates leading to the identification of a single candidate gene. Functional validation showed that this gene was required but not sufficient for Ptr ToxC production, thus it is designated as ToxC1. ToxC1 encoded a conserved hypothetical protein likely located on the vacuole membrane. The gene was highly expressed during infection, and only one haplotype was identified among 120 isolates sequenced. Our work suggests that Ptr ToxC is not a protein and is likely produced through a cascade of biosynthetic pathway. The identification of ToxC1 is a major step toward revealing the Ptr ToxC biosynthetic pathway and studying its molecular interactions with host factors.[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.

Immunomic Investigation of Holocyclotoxins to Produce the First Protective Anti-Venom Vaccine Against the Australian Paralysis Tick, Ixodes holocyclus.

Venom producing animals are ubiquitously disseminated among vertebrates and invertebrates such as fish, snakes, scorpions, spiders, and ticks. Of the ~890 tick species worldwide, 27 have been confirmed to cause paralysis in mammalian hosts. The Australian paralysis tick (Ixodes holocyclus) is the most potent paralyzing tick species known. It is an indigenous three host tick species that secretes potent neurotoxins known as holocyclotoxins (HTs). Holocyclotoxins cause a severe and harmful toxicosis leading to a rapid flaccid paralysis which can result in death of susceptible hosts such as dogs. Antivenins are generally polyclonal antibody treatments developed in sheep, horses or camels to administer following bites from venomous creatures. Currently, the methods to prevent or treat tick paralysis relies upon chemical acaricide preventative treatments or prompt removal of all ticks attached to the host followed by the administration of a commercial tick-antiserum (TAS) respectively. However, these methods have several drawbacks such as poor efficacies, non-standardized dosages, adverse effects and are expensive to administer. Recently the I. holocyclus tick transcriptome from salivary glands and viscera reported a large family of 19 holocyclotoxins at 38-99% peptide sequence identities. A pilot trial demonstrated that correct folding of holocyclotoxins is needed to induce protection from paralysis. The immunogenicity of the holocyclotoxins were measured using commercial tick antiserum selecting HT2, HT4, HT8 and HT11 for inclusion into the novel cocktail vaccine. A further 4 HTs (HT1, HT12, HT14 and HT17) were added to the cocktail vaccine to ensure that the sequence variation among the HT protein family was encompassed in the formulation. A second trial comparing the cocktail of 8 HTs to a placebo group demonstrated complete protection from tick challenge. Here we report the first successful anti-venom vaccine protecting dogs from tick paralysis.

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

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

Remote homology clustering identifies lowly conserved families of effector proteins in plant-pathogenic fungi.

Plant diseases caused by fungal pathogens are typically initiated by molecular interactions between 'effector' molecules released by a pathogen and receptor molecules on or within the plant host cell. In many cases these effector-receptor interactions directly determine host resistance or susceptibility. The search for fungal effector proteins is a developing area in fungal-plant pathology, with more than 165 distinct confirmed fungal effector proteins in the public domain. For a small number of these, novel effectors can be rapidly discovered across multiple fungal species through the identification of known effector homologues. However, many have no detectable homology by standard sequence-based search methods. This study employs a novel comparison method (RemEff) that is capable of identifying protein families with greater sensitivity than traditional homology-inference methods, leveraging a growing pool of confirmed fungal effector data to enable the prediction of novel fungal effector candidates by protein family association. Resources relating to the RemEff method and data used in this study are available from https://figshare.com/projects/Effector_protein_remote_homology/87965.

The first genome assembly of fungal pathogen Pyrenophora tritici-repentis race 1 isolate using Oxford Nanopore MinION sequencing.

OBJECTIVES: The assembly of fungal genomes using short-reads is challenged by long repetitive and low GC regions. However, long-read sequencing technologies, such as PacBio and Oxford Nanopore, are able to overcome many problematic regions, thereby providing an opportunity to improve fragmented genome assemblies derived from short reads only. Here, a necrotrophic fungal pathogen Pyrenophora tritici-repentis (Ptr) isolate 134 (Ptr134), which causes tan spot disease on wheat, was sequenced on a MinION using Oxford Nanopore Technologies (ONT), to improve on a previous Illumina short-read genome assembly and provide a more complete genome resource for pan-genomic analyses of Ptr. RESULTS: The genome of Ptr134 sequenced on a MinION using ONT was assembled into 28 contiguous sequences with a total length of 40.79 Mb and GC content of 50.81%. The long-read assembly provided 6.79 Mb of new sequence and 2846 extra annotated protein coding genes as compared to the previous short-read assembly. This improved genome sequence represents near complete chromosomes, an important resource for large scale and pan genomic comparative analyses.

A genome-wide genetic linkage map and reference quality genome sequence for a new race in the wheat pathogen Pyrenophora tritici-repentis.

Pyrenophora tritici-repentis is an ascomycete fungus that causes tan spot of wheat. The disease has a worldwide distribution and can cause significant yield and quality losses in wheat production. The fungal pathogen is homothallic in nature, which means it can undergo sexual reproduction by selfing to produce pseudothecia on wheat stubble for seasonal survival. Since homothallism precludes the development of bi-parental fungal populations, no genetic linkage map has been developed for P. tritici-repentis for mapping and map-based cloning of fungal virulence genes. In this work, we created two heterothallic strains by deleting one of the mating type genes in each of two parental isolates 86-124 (race 2) and AR CrossB10 (a new race) and developed a bi-parental fungal population between them. The draft genome sequences of the two parental isolates were aligned to the Pt-1C-BFP reference sequence to mine single nucleotide polymorphisms (SNPs). A total of 225 SNP markers were developed for genotyping the entire population. Additionally, 75 simple sequence repeat, and two gene markers were also developed and used in the genotyping. The resulting linkage map consisted of 13 linkage groups spanning 5,075.83 cM in genetic distance. Because the parental isolate AR CrossB10 is a new race and produces Ptr ToxC, it was sequenced using long-read sequencing platforms and de novo assembled into contigs. The majority of the contigs were further anchored into chromosomes with the aid of the linkage maps. The whole genome comparison of AR CrossB10 to the reference genome of M4 revealed a few chromosomal rearrangements. The genetic linkage map and the new AR CrossB10 genome sequence are valuable tools for gene cloning in P. tritici-repentis.

Potentiators of Disease During Barley Infection by Pyrenophora teres f. teres in a Susceptible Interaction.

Pyrenophora teres f. teres is a necrotrophic fungal pathogen and causal agent of net form net blotch (NFNB), a significant disease of barley. RNA-seq data encompassing asymptomatic and subsequent necrotrophic phases of the pathogen was obtained for P. teres f. teres isolate W1-1 in NFNB-sensitive cultivar Baudin. Host genes notably regulated during infection included concerted induction of over half the repertoire of disease resistance genes, together with genes involved in oxidation-reduction processes, characteristic of a hypersensitive response. Several systemic acquired resistance response genes were suppressed and there was a complete absence of defense-related thionin gene expression. In P. teres f. teres, genes involved in hydrolase activities and cell-wall catabolic processes were induced during infection, while nitrate assimilation and response to oxidative stress processes were suppressed. Timecourse data allowed a number of predicted P. teres f. teres effector genes with differing expression profiles to be identified that may underlie barley sensitivity to NFNB. Candidate genes involved in the host-pathogen interaction provide a basis for functional characterization and control strategies based on fungicide or mutation targets, which will facilitate further research aimed at controlling NFNB disease.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Physiological Changes in Barley mlo-11 Powdery Mildew Resistance Conditioned by Tandem Repeat Copy Number.

Wild barley accessions have evolved broad-spectrum defence against barley powdery mildew through recessive mlo mutations. However, the mlo defence response is associated with deleterious phenotypes with a cost to yield and fertility, with implications for natural fitness and agricultural productivity. This research elucidates the mechanism behind a novel mlo allele, designated mlo-11(cnv2), which has a milder phenotype compared to standard mlo-11. Bisulphite sequencing and histone ChIP-seq analyses using near-isogenic lines showed pronounced repression of the Mlo promoter in standard mlo-11 compared to mlo-11(cnv2), with repression governed by 24 nt heterochromatic small interfering RNAs. The mlo-11(cnv2) allele appears to largely reduce the physiological effects of mlo while still endorsing a high level of powdery mildew resistance. RNA sequencing showed that this is achieved through only partly restricted expression of Mlo, allowing adequate temporal induction of defence genes during infection and expression close to wild-type Mlo levels in the absence of infection. The two mlo-11 alleles showed copy number proportionate oxidase and peroxidase expression levels during infection, but lower amino acid and aromatic compound biosynthesis compared to the null allele mlo-5. Examination of highly expressed genes revealed a common WRKY W-box binding motif (consensus ACCCGGGACTAAAGG) and a transcription factor more highly expressed in mlo-11 resistance. In conclusion, mlo-11(cnv2) appears to significantly mitigate the trade-off between mlo defence and normal gene expression.

PacBio genome sequencing reveals new insights into the genomic organisation of the multi-copy ToxB gene of the wheat fungal pathogen Pyrenophora tritici-repentis.

BACKGROUND: Necrotrophic effector proteins secreted by fungal pathogens are important virulence factors that mediate the development of disease in wheat. Pyrenophora tritici-repentis (Ptr), the causal agent of wheat tan spot, has a race structure dependent on the combination of effectors. In Ptr, ToxA and ToxB are known proteinaceous effectors responsible for necrosis and chlorosis respectively. While Ptr ToxA is encoded by the single gene ToxA, ToxB has multiple loci in the Ptr genome, which is postulated to be directly related to the level of ToxB production and leaf chlorosis. Although previous analysis has indicated that the majority of the ToxB loci lie on a single chromosome, the exact number and chromosomal locations for all the ToxB loci have not been fully identified. RESULTS: In this study, we have sequenced the genome of a race 5 ToxB-producing isolate (DW5), using PacBio long read technology, and found that ToxB duplications are nested in the complex subtelomeric chromosomal regions. A total of ten identical ToxB gene copies were identified and based on flanking sequence identity, nine loci appeared associated with chromosome 10 and a single copy with chromosome 5. Chromosome 10 multiple ToxB gene loci were separated by large sequence regions between 31 and 66 kb within larger segmental duplications in an alternating pattern related to loci strand, and flanked by transposable elements. CONCLUSION: This work provides for the first time the full accompaniment of ToxB loci and surrounding regions, and identifies the organization and distribution of ten ToxB loci to subtelomeric regions. To our knowledge, this is the first report of an interwoven strand-related duplication pattern event. This study further highlights the importance of resolving the highly complex distal chromosomal regions, that remain difficult to assemble, and can harbour important effectors and virulence factors.

Expansion and Conservation of Biosynthetic Gene Clusters in Pathogenic Pyrenophora spp.

Pyrenophora is a fungal genus responsible for a number of major cereal diseases. Although fungi produce many specialised or secondary metabolites for defence and interacting with the surrounding environment, the repertoire of specialised metabolites (SM) within Pyrenophora pathogenic species remains mostly uncharted. In this study, an in-depth comparative analysis of the P. teres f. teres, P teres f. maculata and P. tritici-repentis potential to produce SMs, based on in silico predicted biosynthetic gene clusters (BGCs), was conducted using genome assemblies from PacBio DNA reads. Conservation of BGCs between the Pyrenophora species included type I polyketide synthases, terpene synthases and the first reporting of a type III polyketide synthase in P teres f. maculata. P. teres isolates exhibited substantial expansion of non-ribosomal peptide synthases relative to P. tritici-repentis, hallmarked by the presence of tailoring cis-acting nitrogen methyltransferase domains. P. teres isolates also possessed unique non-ribosomal peptide synthase (NRPS)-indole and indole BGCs, while a P. tritici-repentis phytotoxin BGC for triticone production was absent in P. teres. These differences highlight diversification between the pathogens that reflects their different evolutionary histories, host adaption and lifestyles.

Hierarchical clustering of MS/MS spectra from the firefly metabolome identifies new lucibufagin compounds.

Metabolite identification is the greatest challenge when analysing metabolomics data, as only a small proportion of metabolite reference standards exist. Clustering MS/MS spectra is a common method to identify similar compounds, however interrogation of underlying signature fragmentation patterns within clusters can be problematic. Previously published high-resolution LC-MS/MS data from the bioluminescent beetle (Photinus pyralis) provided an opportunity to mine new specialized metabolites in the lucibufagin class, compounds important for defense against predation. We aimed to 1) provide a workflow for hierarchically clustering MS/MS spectra for metabolomics data enabling users to cluster, visualise and easily interrogate the identification of underlying cluster ion profiles, and 2) use the workflow to identify key fragmentation patterns for lucibufagins in the hemolymph of P. pyralis. Features were aligned to their respective MS/MS spectra, then product ions were dynamically binned and resulting spectra were hierarchically clustered and grouped based on a cutoff distance threshold. Using the simplified visualization and the interrogation of cluster ion tables the number of lucibufagins was expanded from 17 to a total of 29.

Genomic Regions Associated with Virulence in Pyrenophora teres f. teres Identified by Genome-Wide Association Analysis and Biparental Mapping.

Net form net blotch (NFNB), caused by the fungal pathogen Pyrenophora teres f. teres, is an important foliar disease present in all barley-producing regions of the world. This fungus is a hemibiotrophic and heterothallic ascomycete, where sexual recombination can lead to changes in disease expression in the host. Knowledge of the genetic architecture and genes involved in virulence is vital to increase the durability of NFNB resistance in barley cultivars. We used a genome-wide association mapping approach to characterize P. teres f. teres genomic regions associated with virulence in Australian barley cultivars. One hundred eighty-eight P. teres f. teres isolates collected across five Australian states were genotyped using Diversity Arrays Technology sequence markers and phenotyped across 20 different barley genotypes. Association mapping identified 14 different genomic regions associated with virulence, with the majority located on P. teres f. teres chromosomes 3 and 5 and one each present on chromosomes 1, 6, and 9. Four of the regions identified were confirmed by quantitative trait loci (QTL) mapping. The QTL regions are discussed in the context of their genomic architecture together with examination of their gene contents, which identified 20 predicted effectors. The number of QTL shown in this study at the population level clearly illustrates a complex genetic basis of P. teres f. teres virulence compared with pure necrotrophs, such as the wheat pathogens Parastagonospora nodorum and Parastagonospora tritici-repentis.