Genomic distribution of a novel Pyrenophora tritici-repentis ToxA insertion element

A New ToxA Haplotype in the Wheat Fungal Pathogen Bipolaris sorokiniana

The necrotrophic effector ToxA is a well-studied virulence factor produced by several fungal necrotrophs. Initially cloned from the wheat tan spot pathogen Pyrenophora tritici-repentis in 1996, ToxA was found almost a decade later in another fungal pathogen, Parastagonospora nodorum, and its sister species, Parastagonospora pseudonodorum. In 2018, ToxA was detected in a third wheat fungal pathogenic species, Bipolaris sorokiniana, which causes spot blotch disease. However, unlike the case with P. tritici-repentis and P. nodorum, the ToxA in B. sorokiniana has only been investigated in recent years. In this report, five Australian B. sorokiniana isolates were assessed for the presence of ToxA. Four isolates were found to contain ToxA. While one isolate harbored the previously reported ToxA haplotype sequence (ToxA19), three isolates contain a different haplotype, designated herein as ToxA25, which has a nonsynonymous mutation resulting in an amino acid change of glycine to arginine at position 168. Both B. sorokiniana ToxA isoforms, when heterologously expressed in Escherichia coli, exhibited the classic ToxA necrosis-inducing activity on ToxA-sensitive Tsn1 cultivars. Preliminary analysis of the B. sorokiniana isolates in Australian wheat cultivars showed that isolates with ToxA19, ToxA25, or ToxA-deficient displayed various degrees of virulence, with the most aggressive isolates observed for those producing ToxA. Differences in spot blotch disease severity between Tsn1 and tsn1 cultivars were observed; however, this was not limited to the ToxA-producing isolates. The overall results suggest that the virulence of the Australian B. sorokiniana isolates is diverse, with the significance of ToxA-Tsn1 interactions depending on individual isolates.

Population-level transposable element expression dynamics influence trait evolution in a fungal crop pathogen

The rapid adaptive evolution of microbes is driven by strong selection pressure acting on genetic variation. How adaptive genetic variation is generated within species and how such variation influences phenotypic trait expression is often not well understood though. We focused on the recent activity of transposable elements (TEs) using deep population genomics and transcriptomics analyses of a fungal plant pathogen with a highly active content of TEs in the genome. Zymoseptoria tritici causes one of the most damaging diseases on wheat, with recent adaptation to the host and environment being facilitated by TE-associated mutations. We obtained genomic and RNA-sequencing data from 146 isolates collected from a single wheat field. We established a genome-wide map of TE insertion polymorphisms in the population by analyzing recent TE insertions among individuals. We quantified the locus-specific transcription of individual TE copies and found considerable population variation at individual TE loci in the population. About 20% of all TE copies show transcription in the genome suggesting that genomic defenses such as repressive epigenetic marks and repeat-induced polymorphisms are at least partially ineffective at preventing the proliferation of TEs in the genome. A quarter of recent TE insertions are associated with expression variation of neighboring genes providing broad potential to influence trait expression. We indeed found that TE insertions are likely responsible for variation in virulence on the host and potentially diverse components of secondary metabolite production. Our large-scale transcriptomics study emphasizes how TE-derived polymorphisms segregate even in individual microbial populations and can broadly underpin trait variation in pathogens.IMPORTANCEPathogens can rapidly adapt to new hosts, antimicrobials, or changes in the environment. Adaptation arises often from mutations in the genome; however, how such variation is generated remains poorly understood. We investigated the most dynamic regions of the genome of Zymoseptoria tritici, a major fungal pathogen of wheat. We focused on the transcription of transposable elements. A large proportion of the transposable elements not only show signatures of potential activity but are also variable within a single population of the pathogen. We find that this variation in activity is likely influencing many important traits of the pathogen. Hence, our work provides insights into how a microbial species can adapt over the shortest time periods based on the activity of transposable elements.

A Race Profile of Tan Spot in Australia Reveals Race 2 Isolates Harboring ToxC1

Tan spot disease is caused by Pyrenophora tritici-repentis (Ptr), one of the major necrotrophic fungal pathogens that affects wheat crops globally. Extensive research has shown that the necrotrophic fungal effectors ToxA, ToxB, and ToxC underlie the genetic interactions of Ptr race classification. ToxA and ToxB are both small proteins secreted during infection; however, the structure of ToxC remains unknown. In line with the recent discovery of the ToxC1 gene that is involved in ToxC production, a subset of 68 isolates collected from the Australian wheat cropping regions were assessed for the presence of all three effectors by pathotyping against four tan spot wheat differential lines and PCR amplification of ToxA, ToxB, and ToxC1. Based on the disease phenotypes, the 68 isolates were grouped into two races with 63 classified as race 1 and five as race 2. A representative selection of each race was tested against eight Australian commercial wheat cultivars and showed no distinction between the virulence levels. Sequencing of ToxA showed that both races had identical gene sequences of haplotype PtrA1. All the race 1 isolates possessed ToxC1 but three race 2 isolates also contained ToxC1 despite being unable to induce a spreading chlorotic symptom on the ToxC differential line. Quantitative trait loci mapping confirmed the absence of the ToxC-Tsc1 association in disease response caused by the ToxC1-containing race 2 isolate; however, ToxC1 expression was detected during plant infection. Altogether, these results suggest that there is a complex regulatory process involved in the production of ToxC within the Australian race 2 isolates.

Fungal Effectoromics: A World in Constant Evolution

Effectors are small, secreted molecules that mediate the establishment of interactions in nature. While some concepts of effector biology have stood the test of time, this area of study is ever-evolving as new effectors and associated characteristics are being revealed. In the present review, the different characteristics that underly effector classifications are discussed, contrasting past and present knowledge regarding these molecules to foster a more comprehensive understanding of effectors for the reader. Research gaps in effector identification and perspectives for effector application in plant disease management are also presented, with a focus on fungal effectors in the plant-microbe interaction and interactions beyond the plant host. In summary, the review provides an amenable yet thorough introduction to fungal effector biology, presenting noteworthy examples of effectors and effector studies that have shaped our present understanding of the field.

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.

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.

Whole-genome sequence analysis of Bipolaris sorokiniana infecting wheat in India and characterization of ToxA gene in different isolates as pathogenicity determinants

Spot blotch disease of wheat caused by Bipolaris sorokiniana Boerma (Sacc.) is an emerging problem in South Asian countries. Whole genome of a highly virulent isolate of B. sorokiniana BS112 (BHU, Uttar Pradesh; Accession no. GCA_004329375.1) was sequenced using a hybrid assembly approach. Secreted proteins, virulence gene(s), pathogenicity-related gene(s) were identified and the role of ToxA gene present in this genome, in the development of disease was recognized. ToxA gene (535 bp) was analyzed and identified in the genome of B. sorokiniana (BS112) which revealed 100% homology with the ToxA gene of Pyrenophora tritici repentis (Accession no. MH017419). Furthermore, ToxA gene was amplified, sequenced and validated in 39 isolates of B. sorokiniana which confirmed the presence of ToxA gene in all the isolates taken for this study. All ToxA sequences were submitted in NCBI database (MN601358-MN601396). As ToxA gene interacts with Tsn1 gene of host, 13 wheat genotypes were evaluated out of which 5 genotypes (38.4%) were found to be Tsn1 positive with more severe necrotic lesions compared to Tsn1-negative wheat genotypes. In vitro expression analysis of ToxA gene in the pathogen B. sorokiniana using qPCR revealed maximum upregulation (14.67 fold) at 1st day after inoculation (DAI) in the medium. Furthermore, in planta expression analysis of ToxA gene in Tsn1-positive and Tsn1-negative genotypes, revealed maximum expression (7.89-fold) in Tsn1-positive genotype, Agra local at 5th DAI compared to Tsn1-negative genotype Chiriya 7 showing minimum expression (0.048-fold) at 5th DAI. In planta ToxA-Tsn1 interaction studies suggested that spot blotch disease is more severe in Tsn1-positive genotypes, which will be helpful in better understanding and management of spot blotch disease of wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03213-3.

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.

Categorization of Orthologous Gene Clusters in 92 Ascomycota Genomes Reveals Functions Important for Phytopathogenicity

Phytopathogenic Ascomycota are responsible for substantial economic losses each year, destroying valuable crops. The present study aims to provide new insights into phytopathogenicity in Ascomycota from a comparative genomic perspective. This has been achieved by categorizing orthologous gene groups (orthogroups) from 68 phytopathogenic and 24 non-phytopathogenic Ascomycota genomes into three classes: Core, (pathogen or non-pathogen) group-specific, and genome-specific accessory orthogroups. We found that (i) ~20% orthogroups are group-specific and accessory in the 92 Ascomycota genomes, (ii) phytopathogenicity is not phylogenetically determined, (iii) group-specific orthogroups have more enriched functional terms than accessory orthogroups and this trend is particularly evident in phytopathogenic fungi, (iv) secreted proteins with signal peptides and horizontal gene transfers (HGTs) are the two functional terms that show the highest occurrence and significance in group-specific orthogroups, (v) a number of other functional terms are also identified to have higher significance and occurrence in group-specific orthogroups. Overall, our comparative genomics analysis determined positive enrichment existing between orthogroup classes and revealed a prediction of what genomic characteristics make an Ascomycete phytopathogenic. We conclude that genes shared by multiple phytopathogenic genomes are more important for phytopathogenicity than those that are unique in each genome.

Research advances in the Pyrenophora teres-barley interaction

Pyrenophora teres f. teres and P. teres f. maculata are significant pathogens that cause net blotch of barley. An increased number of loci involved in P. teres resistance or susceptibility responses of barley as well as interacting P. teres virulence effector loci have recently been identified through biparental and association mapping studies of both the pathogen and host. Characterization of the resistance/susceptibility loci in the host and the interacting effector loci in the pathogen will provide a path for targeted gene validation for better-informed release of resistant barley cultivars. This review assembles concise consensus maps for all loci published for both the host and pathogen, providing a useful resource for the community to be used in pathogen characterization and barley breeding for resistance to both forms of P. teres.

Transposon-Mediated Horizontal Transfer of the Host-Specific Virulence Protein ToxA between Three Fungal Wheat Pathogens

This work dissects the tripartite horizontal transfer of ToxA, a gene that has a direct negative impact on global wheat yields. Defining the extent of horizontally transferred DNA is important because it can provide clues to the mechanisms that facilitate HGT. Our analysis of ToxA and its surrounding 14 kb suggests that this gene was horizontally transferred in two independent events, with one event likely facilitated by a type II DNA transposon. These horizontal transfer events are now in various processes of decay in each species due to the repeated insertion of new transposons and subsequent rounds of targeted mutation by a fungal genome defense mechanism known as repeat induced point mutation. This work highlights the role that HGT plays in the evolution of host adaptation in eukaryotic pathogens. It also increases the growing body of evidence indicating that transposons facilitate adaptive HGT events between fungi present in similar environments and hosts.

Most known examples of horizontal gene transfer (HGT) between eukaryotes are ancient. These events are identified primarily using phylogenetic methods on coding regions alone. Only rarely are there examples of HGT where noncoding DNA is also reported. The gene encoding the wheat virulence protein ToxA and the surrounding 14 kb is one of these rare examples. ToxA has been horizontally transferred between three fungal wheat pathogens (Parastagonospora nodorum, Pyrenophora tritici-repentis, and Bipolaris sorokiniana) as part of a conserved ∼14 kb element which contains coding and noncoding regions. Here we used long-read sequencing to define the extent of HGT between these three fungal species. Construction of near-chromosomal-level assemblies enabled identification of terminal inverted repeats on either end of the 14 kb region, typical of a type II DNA transposon. This is the first description of ToxA with complete transposon features, which we call ToxhAT. In all three species, ToxhAT resides in a large (140-to-250 kb) transposon-rich genomic island which is absent in isolates that do not carry the gene (annotated here as toxa⁻). We demonstrate that the horizontal transfer of ToxhAT between P. tritici-repentis and P. nodorum occurred as part of a large (∼80 kb) HGT which is now undergoing extensive decay. In B. sorokiniana, in contrast, ToxhAT and its resident genomic island are mobile within the genome. Together, these data provide insight into the noncoding regions that facilitate HGT between eukaryotes and into the genomic processes which mask the extent of HGT between these species.