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

Genome comparison of long-circulating field CnmeGV isolates from the same region

Cnaphalocrocis medinalis granulovirus (CnmeGV), belonging to Betabaculovirus cnamedinalis, can infect the rice pest, the rice leaf roller. In 1979, a CnmeGV isolate, CnmeGV-EP, was collected from Enping County, China. In 2014, we collected another CnmeGV isolate, CnmeGV-EPDH3, at the same location and obtained the complete virus genome sequence using Illumina and ONT sequencing technologies. By combining these two virus isolates, we updated the genome annotation of CnmeGV and conducted an in-depth analysis of its genome features. CnmeGV genome contains abundant tandem repeat sequences, and the repeating units in the homologous regions (hrs) exhibit overlapping and nested patterns. The genetic variations within EPDH3 population show the high stability of CnmeGV genome, and tandem repeats are the only region of high genetic variation in CnmeGV genome replication. Some defective viral genomes formed by recombination were found within the population. Comparison analysis of the two virus isolates collected from Enping showed that the proteins encoded by the CnmeGV-specific genes were less conserved relative to the baculovirus core genes. At the genomic level, there are a large number of SNPs and InDels between the two virus isolates, especially in and around the bro genes and hrs. Additionally, we discovered that CnmeGV acquired a segment of non-ORF sequence from its host, which does not provide any new proteins but rather serves as redundant genetic material integrated into the viral genome. Furthermore, we observed that the host's transposon piggyBac has inserted into some virus genes. Together, dsDNA viruses could acquire non-coding genetic material from their hosts to expand the size of their genomes. These findings provide new insights into the evolution of dsDNA viruses.

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.

Systematic identification of cargo-mobilizing genetic elements reveals new dimensions of eukaryotic diversity

Cargo-mobilizing mobile elements (CMEs) are genetic entities that faithfully transpose diverse protein coding sequences. Although common in bacteria, we know little about eukaryotic CMEs because no appropriate tools exist for their annotation. For example, Starships are giant fungal CMEs whose functions are largely unknown because they require time-intensive manual curation. To address this knowledge gap, we developed starfish, a computational workflow for high-throughput eukaryotic CME annotation. We applied starfish to 2 899 genomes of 1 649 fungal species and found that starfish recovers known Starships with 95% combined precision and recall while expanding the number of annotated elements ten-fold. Extant Starship diversity is partitioned into 11 families that differ in their enrichment patterns across fungal classes. Starship cargo changes rapidly such that elements from the same family differ substantially in their functional repertoires, which are predicted to contribute to diverse biological processes such as metabolism. Many elements have convergently evolved to insert into 5S rDNA and AT-rich sequence while others integrate into random locations, revealing both specialist and generalist strategies for persistence. Our work establishes a framework for advancing mobile element biology and provides the means to investigate an emerging dimension of eukaryotic genetic diversity, that of genomes within genomes.

Genome comparisons reveal accessory genes crucial for the evolution of apple Glomerella leaf spot pathogenicity in Colletotrichum fungi

Apple Glomerella leaf spot (GLS) is an emerging fungal disease caused by Colletotrichum fructicola and other Colletotrichum species. These species are polyphyletic and it is currently unknown how these pathogens convergently evolved to infect apple. We generated chromosome-level genome assemblies of a GLS-adapted isolate and a non-adapted isolate in C. fructicola using long-read sequencing. Additionally, we resequenced 17 C. fructicola and C. aenigma isolates varying in GLS pathogenicity using short-read sequencing. Genome comparisons revealed a conserved bipartite genome architecture involving minichromosomes (accessory chromosomes) shared by C. fructicola and other closely related species within the C. gloeosporioides species complex. Moreover, two repeat-rich genomic regions (1.61 Mb in total) were specifically conserved among GLS-pathogenic isolates in C. fructicola and C. aenigma. Single-gene deletion of 10 accessory genes within the GLS-specific regions of C. fructicola identified three that were essential for GLS pathogenicity. These genes encoded a putative non-ribosomal peptide synthetase, a flavin-binding monooxygenase and a small protein with unknown function. These results highlight the crucial role accessory genes play in the evolution of Colletotrichum pathogenicity and imply the significance of an unidentified secondary metabolite in GLS pathogenesis.

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.

Frequent horizontal chromosome transfer between asexual fungal insect pathogens

Entire chromosomes are typically only transmitted vertically from one generation to the next. The horizontal transfer of such chromosomes has long been considered improbable, yet gained recent support in several pathogenic fungi where it may affect the fitness or host specificity. To date, it is unknown how these transfers occur, how common they are, and whether they can occur between different species. In this study, we show multiple independent instances of horizontal transfers of the same accessory chromosome between two distinct strains of the asexual entomopathogenic fungus Metarhizium robertsii during experimental co-infection of its insect host, the Argentine ant. Notably, only the one chromosome-but no other-was transferred from the donor to the recipient strain. The recipient strain, now harboring the accessory chromosome, exhibited a competitive advantage under certain host conditions. By phylogenetic analysis, we further demonstrate that the same accessory chromosome was horizontally transferred in a natural environment between M. robertsii and another congeneric insect pathogen, Metarhizium guizhouense. Hence, horizontal chromosome transfer is not limited to the observed frequent events within species during experimental infections but also occurs naturally across species. The accessory chromosome that was transferred contains genes that may be involved in its preferential horizontal transfer or support its establishment. These genes encode putative histones and histone-modifying enzymes, as well as putative virulence factors. Our study reveals that both intra- and interspecies horizontal transfer of entire chromosomes is more frequent than previously assumed, likely representing a not uncommon mechanism for gene exchange.

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.

Unraveling the dynamics of wheat leaf blight complex: isolation, characterization, and insights into pathogen population under Indian conditions

Wheat, a staple food crop for 35% of the global population, faces a threat from Helminthosporium leaf blight (HLB), a complex of spot blotch (Bipolaris sorokiniana) and tan spot (Pyrenophora-tritici-repentis) diseases under warm and humid conditions. However, in Indian conditions, the knowledge of existing pathogen populations associated with the HLB complex is limited and largely dominated by only B. sorokiniana (spot blotch). To address this, diseased samples were collected from all six wheat growing zones during 2020-2022. The pathogenic species were identified through in-depth morphological characterization, supplemented with ITS-rDNA and GAPDH sequence analysis, a diagnostic SCAR marker, and pathogenicity studies on two wheat varieties: Sonalika and HD2733. The 32 isolates collected from 10 different states consist of B. spicifera (12.5% of all isolates), Exserohilum rostratum (9.3%), Bipolaris oryzae (3.1%), and B. sorokiniana (75%). B. sorokiniana exhibited the highest disease severity on both varieties. Other lesser-known pathogenic species also produced comparable disease severity as B. sorokiniana isolates and, therefore are economically important. Unraveling pathogen composition and biology aids in disease control and resistance breeding. Our study highlights economically impactful and lesser-known pathogenic species causing wheat leaf blight/spot blotch in India, guiding both current management and future resistance breeding strategies in plant pathology.

Whole-genome sequencing and evolutionary analysis of the wild edible mushroom, Morchella eohespera

Morels (Morchella, Ascomycota) are an extremely desired group of edible mushrooms with worldwide distribution. Morchella eohespera is a typical black morel species, belonging to the Elata clade of Morchella species. The biological and genetic studies of this mushroom are rare, largely hindering the studies of molecular breeding and evolutionary aspects. In this study, we performed de novo sequencing and assembly of the M. eohespera strain m200 genome using the third-generation nanopore sequencing platform. The whole-genome size of M. eohespera was 53.81 Mb with a contig N50 of 1.93 Mb, and the GC content was 47.70%. A total of 9,189 protein-coding genes were annotated. Molecular dating showed that M. eohespera differentiated from its relative M. conica at ~19.03 Mya (million years ago) in Burdigalian. Evolutionary analysis showed that 657 gene families were contracted and 244 gene families expanded in M. eohespera versus the related morel species. The non-coding RNA prediction results showed that there were 336 tRNAs, 76 rRNAs, and 45 snRNAs in the M. eohespera genome. Interestingly, there was a high degree of repetition (20.93%) in the M. eohespera genome, and the sizes of long interspersed nuclear elements, short interspersed nuclear elements, and long terminal repeats were 0.83 Mb, 0.009 Mb, and 4.56 Mb, respectively. Additionally, selection pressure analysis identified that a total of 492 genes in the M. eohespera genome have undergone signatures of positive selection. The results of this study provide new insights into the genome evolution of M. eohespera and lay the foundation for in-depth research into the molecular biology of the genus Morchella in the future.

Analysis of five near-complete genome assemblies of the tomato pathogen Cladosporium fulvum uncovers additional accessory chromosomes and structural variations induced by transposable elements effecting the loss of avirulence genes

BACKGROUND

Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogen Cladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. fulvum isolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.

RESULTS

Here, we obtained near-complete genome assemblies of four additional C. fulvum isolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genes Avr4E, Avr5, and Avr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome of C. fulvum.

CONCLUSIONS

Our results reveal new evolutionary aspects of the C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

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.

Surface frustration re-patterning underlies the structural landscape and evolvability of fungal orphan candidate effectors

Pathogens secrete effector proteins to subvert host physiology and cause disease. Effectors are engaged in a molecular arms race with the host resulting in conflicting evolutionary constraints to manipulate host cells without triggering immune responses. The molecular mechanisms allowing effectors to be at the same time robust and evolvable remain largely enigmatic. Here, we show that 62 conserved structure-related families encompass the majority of fungal orphan effector candidates in the Pezizomycotina subphylum. These effectors diversified through changes in patterns of thermodynamic frustration at surface residues. The underlying mutations tended to increase the robustness of the overall effector protein structure while switching potential binding interfaces. This mechanism could explain how conserved effector families maintained biological activity over long evolutionary timespans in different host environments and provides a model for the emergence of sequence-unrelated effector families with conserved structures.

Increased Conidia Production and Germination In Vitro Correlate with Virulence Enhancement in Fusarium oxysporum f. sp. cucumerinum

Cucumber plants commonly suffer from Fusarium wilt disease, which is caused by Fusarium oxysporum f. sp. cucumerinum (Foc). Although resistant cultivars assist with Fusarium wilt disease control, enhancement of the virulence of Foc has been identified after monoculture of wilt-resistant cultivars. To investigate the biological characteristics that contribute to the virulence evolution of Foc, a wildtype strain foc-3b (WT) and its virulence-enhanced variant Ra-4 (InVir) were compared in terms of their growth, reproduction, stress tolerance, and colonization in cucumber plants. The InVir strain showed similar culture characteristics on PDA media to the WT strain but produced significantly more conidia (>two fold), with a distinctly higher germination rate (>four fold) than the WT strain. The colony diameter of the InVir strain increased faster than the WT strain on PDA plates; however, the mycelia dry weight of the InVir was significantly lower (<70%) than that of the WT harvested from PDB. The InVir strain exhibited a significant increase in tolerance to osmolality (1 M NaCl, 1 M KCl, etc.). The GFP-labeled InVir strain propagated in the cucumber vascular faster than the WT strain. These results suggest that increased conidia production and germination in vitro may correlate with virulence enhancement in Fusarium oxysporum f. sp. cucumerinum. This study will provide an insight into its virulence evolution and help us understand the mechanisms underlying the evolutionary biology of F. oxysporum.

Prevalence and Importance of the Necrotrophic Effector Gene ToxA in Bipolaris sorokiniana Populations Collected from Spring Wheat and Barley

Bipolaris sorokiniana is a necrotrophic fungal pathogen that causes foliar and root diseases on wheat and barley. These diseases are common in all wheat- and barley-growing regions, with more severe outbreaks occurring under warm and humid conditions. B. sorokiniana can also infect a wide range of grass species in the family Poaceae and secrete ToxA, an important necrotrophic effector also identified other wheat leaf spotting pathogens. In this study, the prevalence and virulence role of ToxA were investigated in a collection of 278 B. sorokiniana isolates collected from spring wheat and barley in the Upper Midwest of the United States or other places, including 169 from wheat leaves, 75 from wheat roots, 30 from barley leaves, and 4 from wild quack grass leaves. ToxA was present in the isolates from wheat leaves, wheat roots, and wild grass leaves but was absent from isolates collected from barley leaves. Prevalence of ToxA in wheat leaf isolates (34.3%) was much higher than that in wheat root isolates (16%). Sequencing analysis revealed the presence of two haplotypes, with the majority being BsH2. All ToxA+ isolates produced the functional effector in liquid cultures. Pathogenicity assays revealed that ToxA+ isolates caused significantly more disease on spring wheat lines harboring Tsn1 than their tsn1 mutants, suggesting that the ToxA-Tsn1 interaction plays an important role in spot blotch development. This work confirms the importance of ToxA in B. sorokiniana populations infecting wheat and, thus, the need to eliminate Tsn1 from spring wheat cultivars to reduce susceptibility to spot blotch.

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.

The distribution and characteristic of two transposable elements in the genome of Cydia pomonella granulovirus and codling moth

Baculoviruses are capable to acquire insect host transposable elements (TEs) in their genomes and are hypothesized as possible vectors of insect transposons between Lepidopteran species. Here, we investigated the host origin of two TEs, namely the Tc1/mariner-like element TCp3.2 and a 0.7 kbp insertion sequence (IS07), found in the genome of different isolates of Cydia pomonella granulovirus (CpGV), a member of the Betabaculovirus genus. The sequences of both TEs were searched for in the full genome sequence database of codling moth (CM, Cydia pomonella L.). A total of eleven TCp3.2 TE copies and 76 copies of the IS07 fragments were identified in the CM genome. These TEs were distributed over the 22 autosomes and the Z chromosome (chr1) of CM, except chr6, chr12, chr16, chr23, chr27 and the W chromosome (chr29). TCp3.2 copies with two transposase genes in opposite direction, representing a novel feature, were identified on chr10 and chr18. The TCp3.2 transposase was characterized by DD41D motif of classic Tc1/mariner transposons, consisting of DNA-binding domain, catalytic domain and nuclear localization signal (NLS). Transcription analyses of uninfected and CpGV-infected CM larvae suggested a doubling of the TCp3.2 transposase transcription rate in virus infected larvae. Furthermore, IS07 insertion into the CpGV genome apparently added new transcription initiation sites to the viral genome. The global analysis of the distribution of two TEs in the genome of CM addressed the influx of mobile TEs from CM to CpGV, a genetic process that contributes to the population diversity of baculoviruses.

A small cysteine-rich fungal effector, BsCE66 is essential for the virulence of Bipolaris sorokiniana on wheat plants

The Spot Blotch (SB) caused by hemibiotrophic fungal pathogen Bipolaris sorokiniana is one of the most devastating wheat diseases leading to 15-100% crop loss. However, the biology of Triticum-Bipolaris interactions and host immunity modulation by secreted effector proteins remain underexplored. Here, we identified a total of 692 secretory proteins including 186 predicted effectors encoded by B. sorokiniana genome. Gene Ontology categorization showed that these proteins belong to cellular, metabolic and signaling processes, and exhibit catalytic and binding activities. Further, we functionally characterized a cysteine-rich, B. sorokiniana Candidate Effector 66 (BsCE66) that was induced at 24-96 hpi during host colonization. The Δbsce66 mutant did not show vegetative growth defects or stress sensitivity compared to wild-type, but developed drastically reduced necrotic lesions upon infection in wheat plants. The loss-of-virulence phenotype was rescued upon complementing the Δbsce66 mutant with BsCE66 gene. Moreover, BsCE66 does not form homodimer and conserved cysteine residues form intra-molecular disulphide bonds. BsCE66 localizes to the host nucleus and cytosol, and triggers a strong oxidative burst and cell death in Nicotiana benthamiana. Overall, our findings demonstrate that BsCE66 is a key virulence factor that is necessary for host immunity modulation and SB disease progression. These findings would significantly improve our understanding of Triticum-Bipolaris interactions and assist in the development of SB resistant wheat varieties.

Starships are active eukaryotic transposable elements mobilized by a new family of tyrosine recombinases

Transposable elements in eukaryotic organisms have historically been considered "selfish," at best conferring indirect benefits to their host organisms. The Starships are a recently discovered feature in fungal genomes that are, in some cases, predicted to confer beneficial traits to their hosts and also have hallmarks of being transposable elements. Here, we provide experimental evidence that Starships are indeed autonomous transposons, using the model Paecilomyces variotii, and identify the HhpA "Captain" tyrosine recombinase as essential for their mobilization into genomic sites with a specific target site consensus sequence. Furthermore, we identify multiple recent horizontal gene transfers of Starships, implying that they jump between species. Fungal genomes have mechanisms to defend against mobile elements, which are frequently detrimental to the host. We discover that Starships are also vulnerable to repeat-induced point mutation defense, thereby having implications on the evolutionary stability of such elements.

How has horizontal gene transfer shaped the evolution of insect genomes?

As the most diverse group of animals on Earth, insects are key organisms in ecosystems. Horizontal gene transfer (HGT) refers to the transfer of genetic material between species by non-reproductive means. HGT is a major evolutionary force in prokaryotic genome evolution, but its importance in different eukaryotic groups, such as insects, has only recently begun to be understood. Genomic data from hundreds of insect species have enabled the detection of large numbers of HGT events and the elucidation of the functions of some of these foreign genes. Although quantification of the extent of HGT in insects broadens our understanding of its role in insect evolution, the scope of its influence and underlying mechanism(s) of its occurrence remain open questions for the field.

Genome-level analyses resolve an ancient lineage of symbiotic ascomycetes

Ascomycota account for about two-thirds of named fungal species.1 Over 98% of known Ascomycota belong to the Pezizomycotina, including many economically important species as well as diverse pathogens, decomposers, and mutualistic symbionts.2 Our understanding of Pezizomycotina evolution has until now been based on sampling traditionally well-defined taxonomic classes.3,4,5 However, considerable diversity exists in undersampled and uncultured, putatively early-diverging lineages, and the effect of these on evolutionary models has seldom been tested. We obtained genomes from 30 putative early-diverging lineages not included in recent phylogenomic analyses and analyzed these together with 451 genomes covering all available ascomycete genera. We show that 22 of these lineages, collectively representing over 600 species, trace back to a single origin that diverged from the common ancestor of Eurotiomycetes and Lecanoromycetes over 300 million years BP. The new clade, which we recognize as a more broadly defined Lichinomycetes, includes lichen and insect symbionts, endophytes, and putative mycorrhizae and encompasses a range of morphologies so disparate that they have recently been placed in six different taxonomic classes. To test for shared hidden features within this group, we analyzed genome content and compared gene repertoires to related groups in Ascomycota. Regardless of their lifestyle, Lichinomycetes have smaller genomes than most filamentous Ascomycota, with reduced arsenals of carbohydrate-degrading enzymes and secondary metabolite gene clusters. Our expanded genome sample resolves the relationships of numerous "orphan" ascomycetes and establishes the independent evolutionary origins of multiple mutualistic lifestyles within a single, morphologically hyperdiverse clade of fungi.

Minos and Restless transposon insertion mutagenesis of psychrotrophic fungus for red pigment synthesis adaptive to normal temperature

The polar psychrotrophic fungus Geomyces sp. WNF-15A can produce high-quality natural red pigment for the potential use as edible pigment. However, it shows low-temperature-dependent synthesis of red pigment, which limits its large-scale industrial applications due to the difficult and high-cost bioprocess control. This study aims to develop transposon-mediated mutagenesis methods to generate mutants that are able to synthesize red pigment at normal temperature. Four transposable systems, including single and dual transposable systems, were established in this fungus based on the Minos from Drosophila hydei and the Restless from Tolypocladium inflatum. A total of 23 production-dominant mutants and 12 growth-dominant mutants were thus obtained by constructed transposable systems. At 14 °C and 20 °C, the MPS1 mutant strain achieved the highest level of red pigment (OD520 of 43.3 and 29.7, respectively), which was increased by 78.4% and 128.7% compared to the wild-type, respectively. Of note, 4 mutants (MPS1, MPS3, MPS4 and MPD1) successfully synthesized red pigment (OD520 of 5.0, 5.3, 4.7 and 4.9, respectively) at 25 °C, which broke the limit of the wild-type production under normal temperature. Generally, the dual transposable systems of Minos and Restless were more efficient than their single transposable systems for mutagenesis in this fungus. However, the positive mutation ratios were similar between the dual and single transposable systems for either Minos or Restless. This study provides alternative tools for genetic mutagenesis breeding of fungi from extreme environments.

Genetic co-option into plant-filamentous pathogen interactions

Plants are engaged in a coevolutionary arms race with their pathogens that drives rapid diversification and specialization of genes involved in resistance and virulence. However, some major innovations in plant-pathogen interactions, such as molecular decoys, trans-kingdom RNA interference, two-speed genomes, and receptor networks, evolved through the expansion of the functional landscape of genes. This is a typical outcome of genetic co-option, the evolutionary process by which available genes are recruited into new biological functions. Co-option into plant-pathogen interactions emerges generally from (i) cis-regulatory variation, (ii) horizontal gene transfer (HGT), (iii) mutations altering molecular promiscuity, and (iv) rewiring of gene networks and protein complexes. Understanding these molecular mechanisms is key for the functional and predictive biology of plant-pathogen interactions.

The pangenome of the wheat pathogen Pyrenophora tritici-repentis reveals novel transposons associated with necrotrophic effectors ToxA and ToxB

BACKGROUND

In fungal plant pathogens, genome rearrangements followed by selection pressure for adaptive traits have facilitated the co-evolutionary arms race between hosts and their pathogens. Pyrenophora tritici-repentis (Ptr) has emerged recently as a foliar pathogen of wheat worldwide and its populations consist of isolates that vary in their ability to produce combinations of different necrotrophic effectors. These effectors play vital roles in disease development. Here, we sequenced the genomes of a global collection (40 isolates) of Ptr to gain insights into its gene content and genome rearrangements.

RESULTS

A comparative genome analysis revealed an open pangenome, with an abundance of accessory genes (~ 57%) reflecting Ptr's adaptability. A clear distinction between pathogenic and non-pathogenic genomes was observed in size, gene content, and phylogenetic relatedness. Chromosomal rearrangements and structural organization, specifically around effector coding genes, were detailed using long-read assemblies (PacBio RS II) generated in this work in addition to previously assembled genomes. We also discovered the involvement of large mobile elements associated with Ptr's effectors: ToxA, the gene encoding for the necrosis effector, was found as a single copy within a 143-kb 'Starship' transposon (dubbed 'Horizon') with a clearly defined target site and target site duplications. 'Horizon' was located on different chromosomes in different isolates, indicating mobility, and the previously described ToxhAT transposon (responsible for horizontal transfer of ToxA) was nested within this newly identified Starship. Additionally, ToxB, the gene encoding the chlorosis effector, was clustered as three copies on a 294-kb element, which is likely a different putative 'Starship' (dubbed 'Icarus') in a ToxB-producing isolate. ToxB and its putative transposon were missing from the ToxB non-coding reference isolate, but the homolog toxb and 'Icarus' were both present in a different non-coding isolate. This suggests that ToxB may have been mobile at some point during the evolution of the Ptr genome which is contradictory to the current assumption of ToxB vertical inheritance. Finally, the genome architecture of Ptr was defined as 'one-compartment' based on calculated gene distances and evolutionary rates.

CONCLUSIONS

These findings together reflect on the highly plastic nature of the Ptr genome which has likely helped to drive its worldwide adaptation and has illuminated the involvement of giant transposons in facilitating the evolution of virulence in Ptr.

Horizontal gene transfer from genetically modified plants - Regulatory considerations

Gene technology regulators receive applications seeking permission for the environmental release of genetically modified (GM) plants, many of which possess beneficial traits such as improved production, enhanced nutrition and resistance to drought, pests and diseases. The regulators must assess the risks to human and animal health and to the environment from releasing these GM plants. One such consideration, of many, is the likelihood and potential consequence of the introduced or modified DNA being transferred to other organisms, including people. While such gene transfer is most likely to occur to sexually compatible relatives (vertical gene transfer), horizontal gene transfer (HGT), which is the acquisition of genetic material that has not been inherited from a parent, is also a possibility considered during these assessments. Advances in HGT detection, aided by next generation sequencing, have demonstrated that HGT occurrence may have been previously underestimated. In this review, we provide updated evidence on the likelihood, factors and the barriers for the introduced or modified DNA in GM plants to be horizontally transferred into a variety of recipients. We present the legislation and frameworks the Australian Gene Technology Regulator adheres to with respect to the consideration of risks posed by HGT. Such a perspective may generally be applicable to regulators in other jurisdictions as well as to commercial and research organisations who develop GM plants.

Identification of a Novel ToxA Haplotype of Pyrenophora tritici-repentis from Japan

Pyrenophora tritici-repentis was described first as a pathogen of wheat (tan spot) in Japan in the 1920s, but, since then, no reports on P. tritici-repentis race structure or its effectors in Japan have been published. In this study, 10 single-spore isolates of P. tritici-repentis were collected from bread wheat in Japan. These isolates were evaluated for virulence on four differential wheat genotypes and tested for the presence/absence of the effector-encoding genes, ToxA and ToxB, in multiplex PCR assays. These isolates were identified as ToxA producers, of which eight were designated as race 2 (ToxA producers) and two were classified as race 1 (ToxA and ToxC producers) based on their virulence patterns. Sequence analysis of the ToxA amplicons from these 10 isolates indicated the presence of a novel ToxA haplotype (denoted PtrA2). A comparative sequence analysis and resequencing of ToxA from reference P. tritici-repentis isolates showed that all previously published ToxA haplotypes in P. tritici-repentis were identical, and are hence denoted PtrA1 in this study. A total of 163 PtrToxA sequences from global origins were already deposited in GenBank and were confirmed identical to PtrA1. Sequence variation in PtrA1 and PtrA2 open reading frames were found at three positions: one synonymous mutation at position 412 (C/G) and two nonsynonymous mutations at positions 342 and 362 that alter amino acid sequence. These mutations did not seem to affect the necrosis development on a ToxA-sensitive wheat genotype when rated for symptoms 5 to 7 days after inoculation. This is the first report correctly confirming the presence of an additional novel ToxA haplotype in P. tritici-repentis for which we have predicted its isoform and updated the ToxA haplotype evolutionary network.

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.

Giant Starship Elements Mobilize Accessory Genes in Fungal Genomes

Accessory genes are variably present among members of a species and are a reservoir of adaptive functions. In bacteria, differences in gene distributions among individuals largely result from mobile elements that acquire and disperse accessory genes as cargo. In contrast, the impact of cargo-carrying elements on eukaryotic evolution remains largely unknown. Here, we show that variation in genome content within multiple fungal species is facilitated by Starships, a newly discovered group of massive mobile elements that are 110 kb long on average, share conserved components, and carry diverse arrays of accessory genes. We identified hundreds of Starship-like regions across every major class of filamentous Ascomycetes, including 28 distinct Starships that range from 27 to 393 kb and last shared a common ancestor ca. 400 Ma. Using new long-read assemblies of the plant pathogen Macrophomina phaseolina, we characterize four additional Starships whose activities contribute to standing variation in genome structure and content. One of these elements, Voyager, inserts into 5S rDNA and contains a candidate virulence factor whose increasing copy number has contrasting associations with pathogenic and saprophytic growth, suggesting Voyager's activity underlies an ecological trade-off. We propose that Starships are eukaryotic analogs of bacterial integrative and conjugative elements based on parallels between their conserved components and may therefore represent the first dedicated agents of active gene transfer in eukaryotes. Our results suggest that Starships have shaped the content and structure of fungal genomes for millions of years and reveal a new concerted route for evolution throughout an entire eukaryotic phylum.

A chromosome-scale genome assembly of the tomato pathogen Cladosporium fulvum reveals a compartmentalized genome architecture and the presence of a dispensable chromosome

Cladosporium fulvum is a fungal pathogen that causes leaf mould of tomato. The reference genome of this pathogen was released in 2012 but its high repetitive DNA content prevented a contiguous assembly and further prohibited the analysis of its genome architecture. In this study, we combined third generation sequencing technology with the Hi-C chromatin conformation capture technique, to produce a high-quality and near complete genome assembly and gene annotation of a Race 5 isolate of C. fulvum. The resulting genome assembly contained 67.17 Mb organized into 14 chromosomes (Chr1-to-Chr14), all of which were assembled telomere-to-telomere. The smallest of the chromosomes, Chr14, is only 460 kb in size and contains 25 genes that all encode hypothetical proteins. Notably, PCR assays revealed that Chr14 was absent in 19 out of 24 isolates of a world-wide collection of C. fulvum, indicating that Chr14 is dispensable. Thus, C. fulvum is currently the second species of Capnodiales shown to harbour dispensable chromosomes. The genome of C. fulvum Race 5 is 49.7 % repetitive and contains 14 690 predicted genes with an estimated completeness of 98.9%, currently one of the highest among the Capnodiales. Genome structure analysis revealed a compartmentalized architecture composed of gene-dense and repeat-poor regions interspersed with gene-sparse and repeat-rich regions. Nearly 39.2 % of the C. fulvum Race 5 genome is affected by Repeat-Induced Point (RIP) mutations and evidence of RIP leakage toward non-repetitive regions was observed in all chromosomes, indicating the RIP plays an important role in the evolution of this pathogen. Finally, 345 genes encoding candidate effectors were identified in C. fulvum Race 5, with a significant enrichment of their location in gene-sparse regions, in accordance with the 'two-speed genome' model of evolution. Overall, the new reference genome of C. fulvum presents several notable features and is a valuable resource for studies in plant pathogens.

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.

Using structural variants to understand the ecological and evolutionary dynamics of fungal plant pathogens

Deletions, duplications, insertions, inversions and translocations are commonly referred to as structural variants (SVs). Fungal plant pathogens have compact genomes, facilitating the generation of accurate maps of SVs for these species in recent studies. Structural variants have been found to constitute a significant proportion of the standing genetic variation in fungal plant pathogen populations, potentially leading to the generation of accessory genes, regions or chromosomes enriched in pathogenicity factors. Structural variants are involved in the rapid adaptation and ecological traits of pathogens, including host specialization and mating. Long-read sequencing techniques coupled with theoretical and experimental approaches have considerable potential for elucidating the phenotypic effects of SVs and deciphering the evolutionary and genomic mechanisms underlying the formation of SVs in fungal plant pathogens.

A large transposable element mediates metal resistance in the fungus Paecilomyces variotii

The horizontal transfer of large gene clusters by mobile elements is a key driver of prokaryotic adaptation in response to environmental stresses. Eukaryotic microbes face similar stresses; however, a parallel role for mobile elements has not been established. A stress faced by many microorganisms is toxic metal ions in their environment. In fungi, identified mechanisms for protection against metals generally rely on genes that are dispersed within an organism's genome. Here, we discover a large (∼85 kb) region that confers tolerance to five metal/metalloid ions (arsenate, cadmium, copper, lead, and zinc) in the genomes of some, but not all, strains of a fungus, Paecilomyces variotii. We name this region HEPHAESTUS (Hφ) and present evidence that it is mobile within the P. variotii genome with features characteristic of a transposable element. HEPHAESTUS contains the greatest complement of host-beneficial genes carried by a transposable element in eukaryotes, suggesting that eukaryotic transposable elements might play a role analogous to bacteria in the horizontal transfer of large regions of host-beneficial DNA. Genes within HEPHAESTUS responsible for individual metal tolerances include those encoding a P-type ATPase transporter-PcaA-required for cadmium and lead tolerance, a transporter-ZrcA-providing tolerance to zinc, and a multicopper oxidase-McoA-conferring tolerance to copper. In addition, a subregion of Hφ confers tolerance to arsenate. The genome sequences of other fungi in the Eurotiales contain further examples of HEPHAESTUS, suggesting that it is responsible for independently assembling tolerance to a diverse array of ions, including chromium, mercury, and sodium.

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.

Population genomics of transposable element activation in the highly repressive genome of an agricultural pathogen

The activity of transposable elements (TEs) can be an important driver of genetic diversity with TE-mediated mutations having a wide range of fitness consequences. To avoid deleterious effects of TE activity, some fungi have evolved highly sophisticated genomic defences to reduce TE proliferation across the genome. Repeat-induced point mutation (RIP) is a fungal-specific TE defence mechanism efficiently targeting duplicated sequences. The rapid accumulation of RIPs is expected to deactivate TEs over the course of a few generations. The evolutionary dynamics of TEs at the population level in a species with highly repressive genome defences is poorly understood. Here, we analyse 366 whole-genome sequences of Parastagonospora nodorum, a fungal pathogen of wheat with efficient RIP. A global population genomics analysis revealed high levels of genetic diversity and signs of frequent sexual recombination. Contrary to expectations for a species with RIP, we identified recent TE activity in multiple populations. The TE composition and copy numbers showed little divergence among global populations regardless of the demographic history. Miniature inverted-repeat transposable elements (MITEs) and terminal repeat retrotransposons in miniature (TRIMs) were largely underlying recent intra-species TE expansions. We inferred RIP footprints in individual TE families and found that recently active, high-copy TEs have possibly evaded genomic defences. We find no evidence that recent positive selection acted on TE-mediated mutations rather that purifying selection maintained new TE insertions at low insertion frequencies in populations. Our findings highlight the complex evolutionary equilibria established by the joint action of TE activity, selection and genomic repression.

Dynamics of transposable elements in recently diverged fungal pathogens: lineage-specific transposable element content and efficiency of genome defenses

Transposable elements (TEs) impact genome plasticity, architecture, and evolution in fungal plant pathogens. The wide range of TE content observed in fungal genomes reflects diverse efficacy of host-genome defense mechanisms that can counter-balance TE expansion and spread. Closely related species can harbor drastically different TE repertoires. The evolution of fungal effectors, which are crucial determinants of pathogenicity, has been linked to the activity of TEs in pathogen genomes. Here, we describe how TEs have shaped genome evolution of the fungal wheat pathogen Zymoseptoria tritici and four closely related species. We compared de novo TE annotations and repeat-induced point mutation signatures in 26 genomes from the Zymoseptoria species-complex. Then, we assessed the relative insertion ages of TEs using a comparative genomics approach. Finally, we explored the impact of TE insertions on genome architecture and plasticity. The 26 genomes of Zymoseptoria species reflect different TE dynamics with a majority of recent insertions. TEs associate with accessory genome compartments, with chromosomal rearrangements, with gene presence/absence variation, and with effectors in all Zymoseptoria species. We find that the extent of RIP-like signatures varies among Z. tritici genomes compared to genomes of the sister species. The detection of a reduction of RIP-like signatures and TE recent insertions in Z. tritici reflects ongoing but still moderate TE mobility.

The identification of a transposon affecting the asexual reproduction of the wheat pathogen Zymoseptoria tritici

Zymoseptoria tritici, the causal agent of Septoria tritici blotch, is a fungal wheat pathogen that causes significant global yield losses. Within Z. tritici populations, quantitative differences in virulence among different isolates are commonly observed; however, the genetic components that underpin these differences remain elusive. In this study, intraspecific comparative transcriptomic analysis was used to identify candidate genes that contribute to differences in virulence on the wheat cultivar WW2449. This led to the identification of a multicopy gene that was not expressed in the high-virulence isolate when compared to the medium- and low-virulence isolates. Further investigation suggested this gene resides in a 7.9-kb transposon. Subsequent long-read sequencing of the isolates used in the transcriptomic analysis confirmed that this gene did reside in an active Class II transposon, which is composed of four genes named REP9-1 to -4. Silencing and overexpression of REP9-1 in two distinct genetic backgrounds demonstrated that its expression alone reduces the number of pycnidia produced by Z. tritici during infection. The REP9-1 gene identified within a Class II transposon is the first discovery of a gene in a transposable element that influences the virulence of Z. tritici. This discovery adds further complexity to genetic loci that contribute to quantitative virulence in this important pathogen.

Historical genomics reveals the evolutionary mechanisms behind multiple outbreaks of the host-specific coffee wilt pathogen Fusarium xylarioides

BACKGROUND

Nearly 50% of crop yields are lost to pests and disease, with plants and pathogens locked in an amplified co-evolutionary process of disease outbreaks. Coffee wilt disease, caused by Fusarium xylarioides, decimated coffee production in west and central Africa following its initial outbreak in the 1920s. After successful management, it later re-emerged and by the 2000s comprised two separate epidemics on arabica coffee in Ethiopia and robusta coffee in east and central Africa.

RESULTS

Here, we use genome sequencing of six historical culture collection strains spanning 52 years to identify the evolutionary processes behind these repeated outbreaks. Phylogenomic reconstruction using 13,782 single copy orthologs shows that the robusta population arose from the initial outbreak, whilst the arabica population is a divergent sister clade to the other strains. A screen for putative effector genes involved in pathogenesis shows that the populations have diverged in gene content and sequence mainly by vertical processes within lineages. However, 15 putative effector genes show evidence of horizontal acquisition, with close homology to genes from F. oxysporum. Most occupy small regions of homology within wider scaffolds, whereas a cluster of four genes occupy a 20Kb scaffold with strong homology to a region on a mobile pathogenicity chromosome in F. oxysporum that houses known effector genes. Lacking a match to the whole mobile chromosome, we nonetheless found close associations with DNA transposons, especially the miniature impala type previously proposed to facilitate horizontal transfer of pathogenicity genes in F. oxysporum. These findings support a working hypothesis that the arabica and robusta populations partly acquired distinct effector genes via transposition-mediated horizontal transfer from F. oxysporum, which shares coffee as a host and lives on other plants intercropped with coffee.

CONCLUSION

Our results show how historical genomics can help reveal mechanisms that allow fungal pathogens to keep pace with our efforts to resist them. Our list of putative effector genes identifies possible future targets for fungal control. In turn, knowledge of horizontal transfer mechanisms and putative donor taxa might help to design future intercropping strategies that minimize the risk of transfer of effector genes between closely-related Fusarium taxa.

How Do Pathogens Evolve Novel Virulence Activities?

This article is part of the Top 10 Unanswered Questions in MPMI invited review series.We consider the state of knowledge on pathogen evolution of novel virulence activities, broadly defined as anything that increases pathogen fitness with the consequence of causing disease in either the qualitative or quantitative senses, including adaptation of pathogens to host immunity and physiology, host species, genotypes, or tissues, or the environment. The evolution of novel virulence activities as an adaptive trait is based on the selection exerted by hosts on variants that have been generated de novo or arrived from elsewhere. In addition, the biotic and abiotic environment a pathogen experiences beyond the host may influence pathogen virulence activities. We consider host-pathogen evolution, host range expansion, and external factors that can mediate pathogen evolution. We then discuss the mechanisms by which pathogens generate and recombine the genetic variation that leads to novel virulence activities, including DNA point mutation, transposable element activity, gene duplication and neofunctionalization, and genetic exchange. In summary, if there is an (epi)genetic mechanism that can create variation in the genome, it will be used by pathogens to evolve virulence factors. Our knowledge of virulence evolution has been biased by pathogen evolution in response to major gene resistance, leaving other virulence activities underexplored. Understanding the key driving forces that give rise to novel virulence activities and the integration of evolutionary concepts and methods with mechanistic research on plant-microbe interactions can help inform crop protection.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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

Background

The fungus Parastagonospora nodorum causes septoria nodorum blotch (SNB) of wheat (Triticum aestivum) and is a model species for necrotrophic plant pathogens. The genome assembly of reference isolate Sn15 was first reported in 2007. P. nodorum infection is promoted by its production of proteinaceous necrotrophic effectors, three of which are characterised – ToxA, Tox1 and Tox3.

Results

A chromosome-scale genome assembly of P. nodorum Australian reference isolate Sn15, which combined long read sequencing, optical mapping and manual curation, produced 23 chromosomes with 21 chromosomes possessing both telomeres. New transcriptome data were combined with fungal-specific gene prediction techniques and manual curation to produce a high-quality predicted gene annotation dataset, which comprises 13,869 high confidence genes, and an additional 2534 lower confidence genes retained to assist pathogenicity effector discovery. Comparison to a panel of 31 internationally-sourced isolates identified multiple hotspots within the Sn15 genome for mutation or presence-absence variation, which was used to enhance subsequent effector prediction. Effector prediction resulted in 257 candidates, of which 98 higher-ranked candidates were selected for in-depth analysis and revealed a wealth of functions related to pathogenicity. Additionally, 11 out of the 98 candidates also exhibited orthology conservation patterns that suggested lateral gene transfer with other cereal-pathogenic fungal species. Analysis of the pan-genome indicated the smallest chromosome of 0.4 Mbp length to be an accessory chromosome (AC23). AC23 was notably absent from an avirulent isolate and is predominated by mutation hotspots with an increase in non-synonymous mutations relative to other chromosomes. Surprisingly, AC23 was deficient in effector candidates, but contained several predicted genes with redundant pathogenicity-related functions.

Conclusions

We present an updated series of genomic resources for P. nodorum Sn15 – an important reference isolate and model necrotroph – with a comprehensive survey of its predicted pathogenicity content.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07699-8.

Comparison of mitochondrial genomes provides insights into intron dynamics and evolution in Botryosphaeria dothidea and B. kuwatsukai

Botryosphaeria dothidea is one of the most common fungal pathogens on a large number of hosts worldwide. Botryosphaeria dothidea and B. kuwatsukai are also the main causal agents of apple ring rot. In this study, we sequenced, assembled and annotated the circular mitogenomes of 12 diverse B. dothidea isolates (105.7-114.8 kb) infecting various plants including apple, and five diverse B. kuwatsukai isolates (118.0-124.6 kb) from apple. B. dothidea mitogenomes harboured a set of 29-31 introns and 48-52 ORFs. In contrast, B. kuwatsukai mitogenomes harboured more introns (32-34) and ORFs (51-54). The variation in mitogenome sizes was associated mainly with different numbers of introns and insertions of mobile genetic elements. Interestingly, B. dothidea and B. kuwatsukai displayed distinct intron distribution patterns, with three intron loci showing presence/absence dynamics in each species. Large numbers of introns (57% in B. dothidea and 49% in B. kuwatsukai) were most likely obtained through horizontal transfer from non-Dothideomycetes. The mitochondrial gene phylogeny supported the differentiation of the two species. Overall, this study sheds light into the mitochondrial evolution of the plant pathogens B. dothidea and B. kuwatsukai, and intron distribution patterns could be useful markers for studies on population diversity.

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.

The Enterprise, a massive transposon carrying Spok meiotic drive genes

The genomes of eukaryotes are full of parasitic sequences known as transposable elements (TEs). Here, we report the discovery of a putative giant tyrosine-recombinase-mobilized DNA transposon, Enterprise, from the model fungus Podospora anserina Previously, we described a large genomic feature called the Spok block which is notable due to the presence of meiotic drive genes of the Spok gene family. The Spok block ranges from 110 kb to 247 kb and can be present in at least four different genomic locations within P. anserina, despite what is an otherwise highly conserved genome structure. We propose that the reason for its varying positions is that the Spok block is not only capable of meiotic drive but is also capable of transposition. More precisely, the Spok block represents a unique case where the Enterprise has captured the Spoks, thereby parasitizing a resident genomic parasite to become a genomic hyperparasite. Furthermore, we demonstrate that Enterprise (without the Spoks) is found in other fungal lineages, where it can be as large as 70 kb. Lastly, we provide experimental evidence that the Spok block is deleterious, with detrimental effects on spore production in strains which carry it. This union of meiotic drivers and a transposon has created a selfish element of impressive size in Podospora, challenging our perception of how TEs influence genome evolution and broadening the horizons in terms of what the upper limit of transposition may be.

Evolution and Regulation of a Large Effector Family of Pyricularia oryzae

Plant pathogen effectors play important roles in parasitism, including countering plant immunity. However, investigations of the emergence and diversification of fungal effectors across host-adapted populations has been limited. We previously identified a gene encoding a suppressor of plant cell death in Pyricularia oryzae (syn. Magnaporthe oryzae). Here, we report the gene is one of a 21-member gene family and we characterize sequence diversity in different populations. Within the rice pathogen population, nucleotide diversity is low, however; the majority of gene family members display presence-absence polymorphism or other null alleles. Gene family allelic diversity is greater between host-adapted populations and, thus, we named them host-adapted genes (HAGs). Multiple copies of HAGs were found in some genome assemblies and sequence divergence between the alleles in two cases suggested they were the result of repeat-induced point mutagenesis. Transfer of family members between populations and novel HAG haplotypes resulting from apparent recombination were observed. HAG family transcripts were induced in planta and a subset of HAGs are dependent on a key regulator of pathogenesis, PMK1. We also found differential intron splicing for some HAGs that would prevent ex planta protein expression. For some genes, spliced transcript was expressed in antiphase with an overlapping antisense transcript. Characterization of HAG expression patterns and allelic diversity reveal novel mechanisms for HAG regulation and mechanisms generating sequence diversity and novel allele combinations. This evidence of strong in planta-specific expression and selection operating on the HAG family is suggestive of a role in parasitism.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Disarming the Host: Detoxification of Plant Defense Compounds During Fungal Necrotrophy

While fungal biotrophs are dependent on successfully suppressing/subverting host defenses during their interaction with live cells, necrotrophs, due to their lifestyle are often confronted with a suite of toxic metabolites. These include an assortment of plant defense compounds (PDCs) which can demonstrate broad antifungal activity. These PDCs can be either constitutively present in plant tissue or induced in response to infection, but are nevertheless an important obstacle which needs to be overcome for successful pathogenesis. Fungal necrotrophs have developed a number of strategies to achieve this goal, from the direct detoxification of these compounds through enzymatic catalysis and modification, to the active transport of various PDCs to achieve toxin sequestration and efflux. Studies have shown across multiple pathogens that the efficient detoxification of host PDCs is both critical for successful infection and often a determinant factor in pathogen host range. Here, we provide a broad and comparative overview of the various mechanisms for PDC detoxification which have been identified in both fungal necrotrophs and fungal pathogens which depend on detoxification during a necrotrophic phase of infection. Furthermore, the effect that these mechanisms have on fungal host range, metabolism, and disease control will be discussed.

Metarhizium: jack of all trades, master of many

The genus Metarhizium and Pochonia chlamydosporia comprise a monophyletic clade of highly abundant globally distributed fungi that can transition between long-term beneficial associations with plants to transitory pathogenic associations with frequently encountered protozoans, nematodes or insects. Some very common 'specialist generalist' species are adapted to particular soil and plant ecologies, but can overpower a wide spectrum of insects with numerous enzymes and toxins that result from extensive gene duplications made possible by loss of meiosis and associated genome defence mechanisms. These species use parasexuality instead of sex to combine beneficial mutations from separate clonal individuals into one genome (Vicar of Bray dynamics). More weakly endophytic species which kill a narrow range of insects retain sexuality to facilitate host-pathogen coevolution (Red Queen dynamics). Metarhizium species can fit into numerous environments because they are very flexible at the genetic, physiological and ecological levels, providing tractable models to address how new mechanisms for econutritional heterogeneity, host switching and virulence are acquired and relate to diverse sexual life histories and speciation. Many new molecules and functions have been discovered that underpin Metarhizium associations, and have furthered our understanding of the crucial ecology of these fungi in multiple habitats.

Dynamics in Secondary Metabolite Gene Clusters in Otherwise Highly Syntenic and Stable Genomes in the Fungal Genus Botrytis

Fungi of the genus Botrytis infect >1,400 plant species and cause losses in many crops. Besides the broad host range pathogen Botrytis cinerea, most other species are restricted to a single host. Long-read technology was used to sequence genomes of eight Botrytis species, mostly pathogenic on Allium species, and the related onion white rot fungus, Sclerotium cepivorum. Most assemblies contained <100 contigs, with the Botrytis aclada genome assembled in 16 gapless chromosomes. The core genome and pan-genome of 16 Botrytis species were defined and the secretome, effector, and secondary metabolite repertoires analyzed. Among those genes, none is shared among all Allium pathogens and absent from non-Allium pathogens. The genome of each of the Allium pathogens contains 8-39 predicted effector genes that are unique for that single species, none stood out as potential determinant for host specificity. Chromosome configurations of common ancestors of the genus Botrytis and family Sclerotiniaceae were reconstructed. The genomes of B. cinerea and B. aclada were highly syntenic with only 19 rearrangements between them. Genomes of Allium pathogens were compared with ten other Botrytis species (nonpathogenic on Allium) and with 25 Leotiomycetes for their repertoire of secondary metabolite gene clusters. The pattern was complex, with several clusters displaying patchy distribution. Two clusters involved in the synthesis of phytotoxic metabolites are at distinct genomic locations in different Botrytis species. We provide evidence that the clusters for botcinic acid production in B. cinerea and Botrytis sinoallii were acquired by horizontal transfer from taxa within the same genus.

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 11 and a single copy with chromosome 5. Chromosome 11 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.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12864-020-07029-4.

Plant genes hijacked by necrotrophic fungal pathogens

Plant fungal pathogens can be classified according to their lifestyles. Biotrophs feed on living tissue and constitute an economically significant group of pathogens historically. Necrotrophs, which feed on dead tissue, have become economically significant over recent decades, especially those of the Dothideomycetes, which produce necrotrophic effectors (NEs) to modulate the host response. Some of these pathogens interact with their hosts in an inverse gene-for-gene manner, where NEs are recognized by specific dominant genes in the host leading to host-mediated programmed cell death allowing the pathogen to cause disease. Whereas the NE genes tend to be unique, several of the plant 'susceptibility' genes belong to the nucleotide-binding leucine-rich repeat class of disease 'resistance' genes, and one is a wall-associated kinase. These susceptible interactions exhibit hallmarks of defense responses to biotrophic pathogens. Therefore, there is now accumulating evidence that many necrotrophic specialists hijack the resistance mechanisms that are effective against biotrophic pathogens.