EffectorP: predicting fungal effector proteins from secretomes using machine learning

Large-scale transcriptomics to dissect 2 years of the life of a fungal phytopathogen interacting with its host plant

BACKGROUND

The fungus Leptosphaeria maculans has an exceptionally long and complex relationship with its host plant, Brassica napus, during which it switches between different lifestyles, including asymptomatic, biotrophic, necrotrophic, and saprotrophic stages. The fungus is also exemplary of "two-speed" genome organisms in the genome of which gene-rich and repeat-rich regions alternate. Except for a few stages of plant infection under controlled conditions, nothing is known about the genes mobilized by the fungus throughout its life cycle, which may last several years in the field.

RESULTS

We performed RNA-seq on samples corresponding to all stages of the interaction of L. maculans with its host plant, either alive or dead (stem residues after harvest) in controlled conditions or in field experiments under natural inoculum pressure, over periods of time ranging from a few days to months or years. A total of 102 biological samples corresponding to 37 sets of conditions were analyzed. We show here that about 9% of the genes of this fungus are highly expressed during its interactions with its host plant. These genes are distributed into eight well-defined expression clusters, corresponding to specific infection lifestyles or to tissue-specific genes. All expression clusters are enriched in effector genes, and one cluster is specific to the saprophytic lifestyle on plant residues. One cluster, including genes known to be involved in the first phase of asymptomatic fungal growth in leaves, is re-used at each asymptomatic growth stage, regardless of the type of organ infected. The expression of the genes of this cluster is repeatedly turned on and off during infection. Whatever their expression profile, the genes of these clusters are enriched in heterochromatin regions associated with H3K9me3 or H3K27me3 repressive marks. These findings provide support for the hypothesis that part of the fungal genes involved in niche adaptation is located in heterochromatic regions of the genome, conferring an extreme plasticity of expression.

CONCLUSION

This work opens up new avenues for plant disease control, by identifying stage-specific effectors that could be used as targets for the identification of novel durable disease resistance genes, or for the in-depth analysis of chromatin remodeling during plant infection, which could be manipulated to interfere with the global expression of effector genes at crucial stages of plant infection.

Large-Scale Genome Scanning within Exonic Regions Revealed the Contributions of Selective Sweep Prone Genes to Host Divergence and Adaptation in Magnaporthe oryzae Species Complex

Magnaporthe oryzae, one of the most notorious plant pathogens in the agronomic ecosystem, causes a destructive rice blast disease around the world. The blast fungus infects wide arrays of cultivated and non-cultivated plants within the Poaceae. Studies have shown that host speciation exerts selection pressure that drives the evolution and divergence of the M. oryzae population. Population genetic relationship deducted by genome-wide single nucleotide polymorphisms showed that M. oryzae differentiation is highly consistent with the host speciation process. In particular, the rice-infecting population of M. oryzae is distinct from populations from other hosts. However, how genome regions prone to host-mediated selection pressures associated with speciation in M. oryzae, especially at a large-scale population level, has not been extensively characterized. Here, we detected strong evidence of sweep selection throughout the genomes of rice and non-rice pathotypes of M. oryzae population using integrated haplotype score (iHS), cross population extended haplotype homozygosity (XPEHH), and cross population composite likelihood ratio (XPCLR) tests. Functional annotation analyses of the genes associated with host-mediated selection pressure showed that 14 pathogenicity-related genes are under positive selection pressure. Additionally, we showed that 17 candidate effector proteins are under positive and divergent selection among the blast fungus population through sweep selection analysis. Specifically, we find that a divergent selective gene, MGG_13871, is experiencing host-directed mutation in two amino acid residues in rice and non-rice infecting populations. These results provide a crucial insight into the impact of selective sweeping on the differentiation of M. oryzae populations and the dynamic influences of genomic regions in promoting host adaptation and speciation among M. oryzae species.

Size-exclusion chromatography allows the isolation of EVs from the filamentous fungal plant pathogen Fusarium oxysporum f. sp. vasinfectum (Fov)

Extracellular vesicles (EVs) are nano-sized compartments involved in cell communication and macromolecule transport that are well characterized in mammalian organisms. Fungal EVs transport virulence-related cargo and modulate the host immune response, but most work has been focused on human yeast pathogens. Additionally, the study of EVs from filamentous fungi has been hindered by the lack of protein markers and efficient isolation methods. In this study we performed the isolation and proteomic characterization of EVs from the filamentous cotton pathogen Fusarium oxysporum f. sp. vasinfectum (Fov). EVs were recovered from two different growth media, Czapek Dox and Saboraud's dextrose broth, and purified by size-exclusion chromatography. Our results show that the EV proteome changes depending on the growth medium but EV production remains constant. EVs contained proteins involved in polyketide synthesis, cell wall modifications, proteases and potential effectors. These results support a role in modulation of host-pathogen interactions for Fov EVs.

Comparative analysis of extracellular proteomes reveals putative effectors of the boxwood blight pathogens, Calonectria henricotiae and C. pseudonaviculata

Calonectria henricotiae (Che) and C. pseudonaviculata (Cps) are destructive fungal pathogens causing boxwood blight, a persistent threat to horticultural production, landscape industries, established gardens, and native ecosystems. Although extracellular proteins including effectors produced by fungal pathogens are known to play a fundamental role in pathogenesis, the composition of Che and Cps extracellular proteins has not been examined. Using liquid chromatography-tandem mass spectrometry (LC-MS/MS) and bioinformatics prediction tools, 630 extracellular proteins and 251 cell membrane proteins of Che and Cps were identified in the classical secretion pathway in the present study. In the non-classical secretion pathway, 79 extracellular proteins were identified. The cohort of proteins belonged to 364 OrthoMCL clusters, with the majority (62%) present in both species, and a subset unique to Che (19%) and Cps (20%). These extracellular proteins were predicted to play important roles in cell structure, regulation, metabolism, and pathogenesis. A total of 124 proteins were identified as putative effectors. Many of them are orthologs of proteins with documented roles in suppressing host defense and facilitating infection processes in other pathosystems, such as SnodProt1-like proteins in the OrthoMCL cluster OG5_152723 and PhiA-like cell wall proteins in the cluster OG5_155754. This exploratory study provides a repository of secreted proteins and putative effectors that can provide insights into the virulence mechanisms of the boxwood blight pathogens.

The Cytospora chrysosperma Virulence Effector CcCAP1 Mainly Localizes to the Plant Nucleus To Suppress Plant Immune Responses

Canker disease is caused by the fungus Cytospora chrysosperma and damages a wide range of woody plants, causing major losses to crops and native plants. Plant pathogens secrete virulence-related effectors into host cells during infection to regulate plant immunity and promote colonization. However, the functions of C. chrysosperma effectors remain largely unknown. In this study, we used Agrobacterium tumefaciens-mediated transient expression system in Nicotiana benthamiana and confocal microscopy to investigate the immunoregulation roles and subcellular localization of CcCAP1, a virulence-related effector identified in C. chrysosperma CcCAP1 was significantly induced in the early stages of infection and contains cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 proteins (CAP) superfamily domain with four cysteines. CcCAP1 suppressed the programmed cell death triggered by Bcl-2-associated X protein (BAX) and the elicitin infestin1 (INF1) in transient expression assays with Nicotiana benthamiana The CAP superfamily domain was sufficient for its cell death-inhibiting activity and three of the four cysteines in the CAP superfamily domain were indispensable for its activity. Pathogen challenge assays in N. benthamiana demonstrated that transient expression of CcCAP1 promoted Botrytis cinerea infection and restricted reactive oxygen species accumulation, callose deposition, and defense-related gene expression. In addition, expression of green fluorescent protein-labeled CcCAP1 in N. benthamiana showed that it localized to both the plant nucleus and the cytoplasm, but the nuclear localization was essential for its full immune inhibiting activity. These results suggest that this virulence-related effector of C. chrysosperma modulates plant immunity and functions mainly via its nuclear localization and the CAP domain.IMPORTANCE The data presented in this study provide a key resource for understanding the biology and molecular basis of necrotrophic pathogen responses to Nicotiana benthamiana resistance utilizing effector proteins, and CcCAP1 may be used in future studies to understand effector-triggered susceptibility processes in the Cytospora chrysosperma-poplar interaction system.

Near-Chromosome-Level Genome Assembly of the Dark Septate Endophyte Laburnicola rhizohalophila: A Model for Investigating Root-Fungus Symbiosis

The novel DSE Laburnicola rhizohalophila (Pleosporales, Ascomycota) is frequently found in the halophytic seepweed (Suaeda salsa). In this article, we report a near-chromosome-level hybrid assembly of this fungus using a combination of short-read Illumina data to polish assemblies generated from long-read Nanopore data. The reference genome for L. rhizohalophila was assembled into 26 scaffolds with a total length of 64.0 Mb and a N50 length of 3.15 Mb. Of them, 17 scaffolds approached the length of intact chromosomes, and 5 had telomeres at one end only. A total of 10,891 gene models were predicted. Intriguingly, 27.5 Mb of repeat sequences that accounted for 42.97% of the genome was identified, and long terminal repeat retrotransposons were the most frequent known transposable elements, indicating that transposable element proliferation contributes to its increased genome size. BUSCO analyses using the Fungi_odb10 data set showed that 95.0% of genes were complete. In addition, 292 carbohydrate active enzymes, 33 secondary metabolite clusters, and 84 putative effectors were identified in silico. The resulting high-quality assembly and genome features are not only an important resource for further research on understanding the mechanism of root-fungi symbiotic interactions but will also contribute to comparative analyses of genome biology and evolution within Pleosporalean species.

A small secreted protein from Zymoseptoria tritici interacts with a wheat E3 ubiquitin ligase to promote disease

Septoria tritici blotch (STB), caused by the ascomycete fungus Zymoseptoria tritici, is a major threat to wheat production worldwide. The Z. tritici genome encodes many small secreted proteins (ZtSSPs) that are likely to play a key role in the successful colonization of host tissues. However, few of these ZtSSPs have been functionally characterized for their role during infection. In this study, we identified and characterized a small, conserved cysteine-rich secreted effector from Z. tritici which has homologues in other plant pathogens in the Dothideomycetes. ZtSSP2 was expressed throughout Z. tritici infection in wheat, with the highest levels observed early during infection. A yeast two-hybrid assay revealed an interaction between ZtSSP2 and wheat E3 ubiquitin ligase (TaE3UBQ) in yeast, and this was further confirmed in planta using bimolecular fluorescence complementation and co-immunoprecipitation. Down-regulation of this wheat E3 ligase using virus-induced gene silencing increased the susceptibility of wheat to STB. Together, these results suggest that TaE3UBQ is likely to play a role in plant immunity to defend against Z. tritici.

The Epichloë festucae Antifungal Protein Efe-AfpA Is also a Possible Effector Protein Required for the Interaction of the Fungus with Its Host Grass Festuca rubra subsp. rubra

Strong creeping red fescue (Festuca rubra subsp. rubra) is a commercially important low-maintenance turfgrass and is often naturally infected with the fungal endophyte Epichloë festucae. Epichloë spp. are endophytes of several cool-season grass species, often conferring insect resistance to the grass hosts due to the production of toxic alkaloids. In addition to insect resistance, a unique feature of the strong creeping red fescue/E. festucae symbiosis is the endophyte-mediated disease resistance to the fungal pathogen Clarireedia jacksonii, the causal agent of dollar spot disease. Such disease resistance is not a general feature of other grass/ Epichloë interactions. E. festucae isolates infecting red fescue have an antifungal protein gene Efe-afpA, whereas most other Epichloë spp. do not have a similar gene. The uniqueness of this gene suggests it may, therefore, be a component of the unique disease resistance seen in endophyte-infected red fescue. Here, we report the generation of CRISPR-Cas9 Efe-afpA gene knockouts with the goal of determining if absence of the protein in endophyte-infected Festuca rubra leads to disease susceptibility. However, it was not possible to infect plants with the knockout isolates, although infection was possible with the wild type E. festucae and with complemented isolates. This raises the interesting possibility that, in addition to having antifungal activity, the protein is required for the symbiotic interaction. The antifungal protein is a small secreted protein with high expression in planta relative to its expression in culture, all characteristics consistent with effector proteins. If Efe-AfpA is an effector protein it must be specific to certain interactions, since most Epichloë spp. do not have such a gene in their genomes.

Regulation of host-infection ability in the grass-symbiotic fungus Epichloë festucae by histone H3K9 and H3K36 methyltransferases

Recent studies have identified key genes that control the symbiotic interaction between Epichloë festucae and Lolium perenne. Here we report on the identification of specific E. festucae genes that control host infection. Deletion of setB, which encodes a homologue of the H3K36 histone methyltransferase Set2/KMT3, reduced histone H3K36 trimethylation and led to severe defects in colony growth and hyphal development. The E. festucae ΔclrD mutant, which lacks the gene encoding the homologue of the H3K9 methyltransferase KMT1, displays similar developmental defects. Both mutants are completely defective in their ability to infect L. perenne. Alleles that complement the culture and plant phenotypes of both mutants also complement the histone methylation defects. Co-inoculation of either ΔsetB or ΔclrD with the wild-type strain enables these mutants to colonize the host. However, successful colonization by the mutants resulted in death or stunting of the host plant. Transcriptome analysis at the early infection stage identified four fungal candidate genes, three of which encode small-secreted proteins, that are differentially regulated in these mutants compared to wild type. Deletion of crbA, which encodes a putative carbohydrate binding protein, resulted in significantly reduced host infection rates by E. festucae.

In-depth secretome analysis of Puccinia striiformis f. sp. tritici in infected wheat uncovers effector functions

The importance of wheat yellow rust disease, caused by Puccinia striiformis f. sp. tritici (Pst), has increased substantially due to the emergence of aggressive new Pst races in the last couple of decades. In an era of escalating human populations and climate change, it is vital to understand the infection mechanism of Pst in order to develop better strategies to combat wheat yellow disease. The present study focuses on the identification of small secreted proteins (SSPs) and candidate-secreted effector proteins (CSEPs) that are used by the pathogen to support infection and control disease development. We generated de novo assembled transcriptomes of Pst collected from wheat fields in central Anatolia. We inoculated both susceptible and resistant seedlings with Pst and analyzed haustoria formation. At 10 days post-inoculation (dpi), we analyzed the transcriptomes and identified 10550 Differentially Expressed Unigenes (DEGs), of which 6072 were Pst-mapped. Among those Pst-related genes, 227 were predicted as PstSSPs. In silico characterization was performed using an approach combining the transcriptomic data and data mining results to provide a reliable list to narrow down the ever-expanding repertoire of predicted effectorome. The comprehensive analysis detected 14 Differentially Expressed Small-Secreted Proteins (DESSPs) that overlapped with the genes in available literature data to serve as the best CSEPs for experimental validation. One of the CSEPs was cloned and studied to test the reliability of the presented data. Biological assays show that the randomly selected CSEP, Unigene17495 (PSTG_10917), localizes in the chloroplast and is able to suppress cell death induced by INF1 in a Nicotiana benthamiana heterologous expression system.

First Draft Genome Resource for the Tomato Black Leaf Mold Pathogen Pseudocercospora fuligena

Pseudocercospora fuligena is a fungus that causes black leaf mold, an important disease of tomato in tropical and subtropical regions of the world. Despite its economic importance, genomic resources for this pathogen are scarce and no reference genome was available thus far. Here, we report a 50.6-Mb genome assembly for P. fuligena, consisting of 348 contigs with an N₅₀ value of 0.407 Mb. In total, 13,764 protein-coding genes were predicted with an estimated BUSCO completeness of 98%. Among the predicted genes there were 179 candidate effectors, 445 carbohydrate-active enzymes, and 30 secondary metabolite gene clusters. The resources presented in this study will allow genome-wide comparative analyses and population genomic studies of this pathogen, ultimately improving management strategies for black leaf mold of tomato.

Comparative Genome Analyses of 18 Verticillium dahliae Tomato Isolates Reveals Phylogenetic and Race Specific Signatures

Host resistance is one of the few strategies available to combat the soil borne pathogenic fungus Verticillium dahliae. Understanding pathogen diversity in populations is key to successfully deploying host resistance. In this study the genomes of 18 V. dahliae isolates of races 1 (n = 2), 2 (n = 4), and 3 (n = 12) from Japan, California, and North Carolina were sequenced and mapped to the reference genome of JR2 (from tomato). The genomes were analyzed for phylogenetic and pathogen specific signatures to classify specific strains or genes for future research. Four highly clonal lineages/groups were discovered, including a lineage unique to North Carolina isolates, which had the rare MAT1-1 mating type. No evidence for recombination between isolates of different mating types was observed, even in isolates of different mating types discovered in the same field. By mapping these 18 isolates genomes to the JR2 reference genome, 193 unique candidate effectors were found using SignalP and EffectorP. Within these effectors, 144 highly conserved effectors, 42 mutable effectors (truncated or present in some isolates but absent in others), and 7 effectors present in highly variable regions of the chromosomes were discovered. Of the 144 core effectors, 21 were highly conserved in V. alfalfae and V. longisporum, 7 of which have no known function. Within the non-core effectors 30 contained large numbers of non-synonymous mutations, while 15 of them contained indels, frameshift mutations, or were present on highly variable regions of the chromosome. Two of these highly variable region effectors (HVREs) were only present in race 2 isolates, but not in race 3 isolates. The race 1 effector Ave1 was also present in a highly variable region. These data may suggest that these highly variable regions are enriched in race determinant genes, consistent with the two-speed genome hypothesis.

Deciphering the Infectious Process of Colletotrichum lupini in Lupin through Transcriptomic and Proteomic Analysis

The fungal phytopathogen Colletotrichum lupini is responsible for lupin anthracnose, resulting in significant yield losses worldwide. The molecular mechanisms underlying this infectious process are yet to be elucidated. This study proposes to evaluate C. lupini gene expression and protein synthesis during lupin infection, using, respectively, an RNAseq-based transcriptomic approach and a mass spectrometry-based proteomic approach. Patterns of differentially-expressed genes in planta were evaluated from 24 to 84 hours post-inoculation, and compared to in vitro cultures. A total of 897 differentially-expressed genes were identified from C. lupini during interaction with white lupin, of which 520 genes were predicted to have a putative function, including carbohydrate active enzyme, effector, protease or transporter-encoding genes, commonly described as pathogenicity factors for other Colletotrichum species during plant infection, and 377 hypothetical proteins. Simultaneously, a total of 304 proteins produced during the interaction were identified and quantified by mass spectrometry. Taken together, the results highlight that the dynamics of symptoms, gene expression and protein synthesis shared similarities to those of hemibiotrophic pathogens. In addition, a few genes with unknown or poorly-described functions were found to be specifically associated with the early or late stages of infection, suggesting that they may be of importance for pathogenicity. This study, conducted for the first time on a species belonging to the Colletotrichum acutatum species complex, presents an opportunity to deepen functional analyses of the genes involved in the pathogenicity of Colletotrichum spp. during the onset of plant infection.

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.

Chromatin Dynamics Contribute to the Spatiotemporal Expression Pattern of Virulence Genes in a Fungal Plant Pathogen

Fungal plant pathogens possess a large repertoire of genes encoding putative effectors, which are crucial for infection. Many of these genes are expressed at low levels in the absence of the host but are strongly induced at specific stages of the infection. The mechanisms underlying this transcriptional reprogramming remain largely unknown. We investigated the role of the genomic environment and associated chromatin modifications of effector genes in controlling their expression pattern in the fungal wheat pathogen Zymoseptoria tritici. Depending on their genomic location, effector genes are epigenetically repressed in the absence of the host and during the initial stages of infection. Derepression of effector genes occurs mainly during and after penetration of plant leaves and is associated with changes in histone modifications. Our work demonstrates the role of chromatin in shaping the expression of virulence components and, thereby, the interaction between fungal pathogens and their plant hosts.

Dynamic changes in transcription profiles are key for the success of pathogens in colonizing their hosts. In many pathogens, genes associated with virulence, such as effector genes, are located in regions of the genome that are rich in transposable elements and heterochromatin. The contribution of chromatin modifications to gene expression in pathogens remains largely unknown. Using a combination of a reporter gene-based approach and chromatin immunoprecipitation, we show that the heterochromatic environment of effector genes in the fungal plant pathogen Zymoseptoria tritici is a key regulator of their specific spatiotemporal expression patterns. Enrichment in trimethylated lysine 27 of histone H3 dictates the repression of effector genes in the absence of the host. Chromatin decondensation during host colonization, featuring a reduction in this repressive modification, indicates a major role for epigenetics in effector gene induction. Our results illustrate that chromatin modifications triggered during host colonization determine the specific expression profile of effector genes at the cellular level and, hence, provide new insights into the regulation of virulence in fungal plant pathogens.

Machine learning in plant-pathogen interactions: empowering biological predictions from field scale to genome scale

Machine learning (ML) encompasses statistical methods that learn to identify patterns in complex datasets. Here, I review application areas in plant-pathogen interactions that have recently benefited from ML, such as disease monitoring, the discovery of gene regulatory networks, genomic selection for disease resistance and prediction of pathogen effectors. However, achieving robust performance from ML is not trivial and requires knowledge of both the methodology and the biology. I discuss common pitfalls and challenges in using ML approaches. Finally, I highlight future opportunities for ML as a tool for dissecting plant-pathogen interactions using high-throughput data, for example, through integration of diverse data sources and the analysis with higher resolution, such as from individual cells or on elaborate spatial and temporal scales.

Predicating the Effector Proteins Secreted by Puccinia triticina Through Transcriptomic Analysis and Multiple Prediction Approaches

Wheat leaf rust caused by Puccinia triticina is one of the most common and serious diseases in wheat production. The constantly changing pathogens overcome the plant resistance to P. triticina. Plant pathogens secrete effector proteins that alter the structure of the host cell, interfere plant defenses, or modify the physiology of plant cells. Therefore, the identification of effector proteins is critical to reveal the pathogenic mechanism. We used SignalP v4.1, TargetP v1.1, TMHMM v2.0, and EffectorP v2.0 to screen the candidate effector proteins in P. triticina isolates – KHTT, JHKT, and THSN. As a result, a total of 635 candidate effector proteins were obtained. Structural analysis showed that effector proteins were small in size (50AA to 422AA) and of diverse sequences, and the conserved sequential elements or clear common elements were not involved, regardless of their secretion from the pathogen to the host. There were 427 candidate effector proteins that contain more than or equal to 4 cysteine residues, and 339 candidate effector proteins contained the known motifs. Sixteen families, 9 domains, and 53 other known functional types were found in 186 candidate effector proteins using the Pfam search. Three novel motifs were found by MEME. Heterogeneous expression system was performed to verify the functions of 30 candidate effectors by inhibiting the programmed cell death (PCD) induced by BAX (the mouse-apoptotic gene elicitor) on Nicotiana benthamiana. Hypersensitive response (HR) can be induced by the six effectors in the wheat leaf rust resistance near isogenic lines, and this would be shown by the method of transient expression through Agrobacterium tumefaciens infiltration. The quantitative reverse transcription PCR (qRT-PCR) analysis of 14 candidate effector proteins secreted after P. triticina inoculation showed that the tested effectors displayed different expression patterns in different stages, suggesting that they may be involved in the wheat–P. triticina interaction. The results showed that the prediction of P. triticina effector proteins based on transcriptomic analysis and multiple bioinformatics software is effective and more accurate, laying the foundation of revealing the pathogenic mechanism of Pt and controlling disease.

Identification of effector candidate genes of Rhizoctonia solani AG-1 IA expressed during infection in Brachypodium distachyon

Rhizoctonia solani is a necrotrophic phytopathogen belonging to basidiomycetes. It causes rice sheath blight which inflicts serious damage in rice production. The infection strategy of this pathogen remains unclear. We previously demonstrated that salicylic acid-induced immunity could block R. solani AG-1 IA infection in both rice and Brachypodium distachyon. R. solani may undergo biotrophic process using effector proteins to suppress host immunity before necrotrophic stage. To identify pathogen genes expressed at the early infection process, here we developed an inoculation method using B. distachyon which enables to sample an increased amount of semi-synchronous infection hyphae. Sixty-one R. solani secretory effector-like protein genes (RsSEPGs) were identified using in silico approach with the publicly available gene annotation of R. solani AG-1 IA genome and our RNA-sequencing results obtained from hyphae grown on agar medium. Expression of RsSEPGs was analyzed at 6, 10, 16, 24, and 32 h after inoculation by a quantitative reverse transcription-polymerase chain reaction and 52 genes could be detected at least on a single time point tested. Their expressions showed phase-specific patterns which were classified into 6 clusters. The 23 RsSEPGs in the cluster 1–3 and 29 RsSEPGs in the cluster 4–6 are expected to be involved in biotrophic and necrotrophic interactions, respectively.

Effector Biology of Biotrophic Plant Fungal Pathogens: Current Advances and Future Prospects

The interaction of fungal pathogens with their host requires a novel invading mechanism and the presence of various virulence-associated components responsible for promoting the infection. The small secretory proteins, explicitly known as effector proteins, are one of the prime mechanisms of host manipulation utilized by the pathogen to disarm the host. Several effector proteins are known to translocate from fungus to the plant cell for host manipulation. Many fungal effectors have been identified using genomic, transcriptomic, and bioinformatics approaches. Most of the effector proteins are devoid of any conserved signatures, and their prediction based on sequence homology is very challenging, therefore by combining the sequence consensus based upon machine learning features, multiple tools have also been developed for predicting apoplastic and cytoplasmic effectors. Various post-genomics approaches like transcriptomics of virulent isolates have also been utilized for identifying active consortia of effectors. Significant progress has been made in understanding biotrophic effectors; however, most of it is underway due to their complex interaction with host and complicated recognition and signaling networks. This review discusses advances, and challenges in effector identification and highlighted various features of the potential effector proteins and approaches for understanding their genetics and strategies for regulation.

Apoplastic effector proteins of plant-associated fungi and oomycetes

The outcome of an interaction between a plant and a fungus or an oomycete, whether compatibility or incompatibility, is often determined in the hostile extracellular spaces and matrices of the apoplast. Indeed, for compatibility to occur, many plant-associated fungi and oomycetes must first neutralize the apoplast, which is both monitored by plant cell-surface immune receptors, and enriched in plant (and frequently, competitor)-derived antimicrobial compounds. Research is highlighting the diverse roles that fungal and oomycete effector proteins play in the apoplast to promote compatibility, with most recent progress made towards understanding the role of these proteins in evading chitin-triggered immunity. Research is also showcasing the ability of apoplastic effector proteins to bring about incompatibility upon recognition by diverse plant cell-surface immune receptors, and the use of effectoromics to rapidly identify apoplastic effector protein-cell-surface immune receptor interactions.

Reference Genome Assembly for Australian Ascochyta rabiei Isolate ArME14

Ascochyta rabiei is the causal organism of ascochyta blight of chickpea and is present in chickpea crops worldwide. Here we report the release of a high-quality PacBio genome assembly for the Australian A. rabiei isolate ArME14. We compare the ArME14 genome assembly with an Illumina assembly for Indian A. rabiei isolate, ArD2. The ArME14 assembly has gapless sequences for nine chromosomes with telomere sequences at both ends and 13 large contig sequences that extend to one telomere. The total length of the ArME14 assembly was 40,927,385 bp, which was 6.26 Mb longer than the ArD2 assembly. Division of the genome by OcculterCut into GC-balanced and AT-dominant segments reveals 21% of the genome contains gene-sparse, AT-rich isochores. Transposable elements and repetitive DNA sequences in the ArME14 assembly made up 15% of the genome. A total of 11,257 protein-coding genes were predicted compared with 10,596 for ArD2. Many of the predicted genes missing from the ArD2 assembly were in genomic regions adjacent to AT-rich sequence. We compared the complement of predicted transcription factors and secreted proteins for the two A. rabiei genome assemblies and found that the isolates contain almost the same set of proteins. The small number of differences could represent real differences in the gene complement between isolates or possibly result from the different sequencing methods used. Prediction pipelines were applied for carbohydrate-active enzymes, secondary metabolite clusters and putative protein effectors. We predict that ArME14 contains between 450 and 650 CAZymes, 39 putative protein effectors and 26 secondary metabolite clusters.

Genome-Wide Identification of Effector Candidates With Conserved Motifs From the Wheat Leaf Rust Fungus Puccinia triticina

Rust fungi secrete various specialized effectors into host cells to manipulate the plant defense response. Conserved motifs, including RXLR, LFLAK-HVLVxxP (CRN), Y/F/WxC, CFEM, LysM, EAR, [SG]-P-C-[KR]-P, DPBB_1 (PNPi), and ToxA, have been identified in various oomycete and fungal effectors and are reported to be crucial for effector translocation or function. However, little is known about potential effectors containing any of these conserved motifs in the wheat leaf rust fungus (Puccinia triticina, Pt). In this study, sequencing was performed on RNA samples collected from the germ tubes (GT) of uredospores of an epidemic Pt pathotype PHTT(P) and Pt-infected leaves of a susceptible wheat cultivar "Chinese Spring" at 4, 6, and 8 days post-inoculation (dpi). The assembled transcriptome data were compared to the reference genome of "Pt 1-1 BBBD Race 1." A total of 17,976 genes, including 2,284 "novel" transcripts, were annotated. Among all these genes, we identified 3,149 upregulated genes upon Pt infection at all time points compared to GT, whereas 1,613 genes were more highly expressed in GT. A total of 464 secreted proteins were encoded by those upregulated genes, with 79 of them also predicted as possible effectors by EffectorP. Using hmmsearch and Regex, we identified 719 RXLR-like, 19 PNPi-like, 19 CRN-like, 138 Y/F/WxC, and 9 CFEM effector candidates from the deduced protein database including data based on the "Pt 1-1 BBBD Race 1" genome and the transcriptome data collected here. Four of the PNPi-like effector candidates with DPBB_1 conserved domain showed physical interactions with wheat NPR1 protein in yeast two-hybrid assay. Nine Y/F/WxC and seven CFEM effector candidates were transiently expressed in Nicotiana benthamiana. None of these effector candidates showed induction or suppression of cell death triggered by BAX protein, but the expression of one CFEM effector candidate, PTTG_08198, accelerated the progress of cell death and promoted the accumulation of reactive oxygen species (ROS). In conclusion, we profiled genes associated with the infection process of the Pt pathotype PHTT(P). The identified effector candidates with conserved motifs will help guide the investigation of virulent mechanisms of leaf rust fungus.

Genome mining of the citrus pathogen Elsinoë fawcettii; prediction and prioritisation of candidate effectors, cell wall degrading enzymes and secondary metabolite gene clusters

Elsinoë fawcettii, a necrotrophic fungal pathogen, causes citrus scab on numerous citrus varieties around the world. Known pathotypes of E. fawcettii are based on host range; additionally, cryptic pathotypes have been reported and more novel pathotypes are thought to exist. E. fawcettii produces elsinochrome, a non-host selective toxin which contributes to virulence. However, the mechanisms involved in potential pathogen-host interactions occurring prior to the production of elsinochrome are unknown, yet the host-specificity observed among pathotypes suggests a reliance upon such mechanisms. In this study we have generated a whole genome sequencing project for E. fawcettii, producing an annotated draft assembly 26.01 Mb in size, with 10,080 predicted gene models and low (0.37%) coverage of transposable elements. A small proportion of the assembly showed evidence of AT-rich regions, potentially indicating genomic regions with increased plasticity. Using a variety of computational tools, we mined the E. fawcettii genome for potential virulence genes as candidates for future investigation. A total of 1,280 secreted proteins and 276 candidate effectors were predicted and compared to those of other necrotrophic (Botrytis cinerea, Parastagonospora nodorum, Pyrenophora tritici-repentis, Sclerotinia sclerotiorum and Zymoseptoria tritici), hemibiotrophic (Leptosphaeria maculans, Magnaporthe oryzae, Rhynchosporium commune and Verticillium dahliae) and biotrophic (Ustilago maydis) plant pathogens. Genomic and proteomic features of known fungal effectors were analysed and used to guide the prioritisation of 120 candidate effectors of E. fawcettii. Additionally, 378 carbohydrate-active enzymes were predicted and analysed for likely secretion and sequence similarity with known virulence genes. Furthermore, secondary metabolite prediction indicated nine additional genes potentially involved in the elsinochrome biosynthesis gene cluster than previously described. A further 21 secondary metabolite clusters were predicted, some with similarity to known toxin producing gene clusters. The candidate virulence genes predicted in this study provide a comprehensive resource for future experimental investigation into the pathogenesis of E. fawcettii.

AcDCXR Is a Cowpea Aphid Effector With Putative Roles in Altering Host Immunity and Physiology

Cowpea, Vigna unguiculata, is a crop that is essential to semiarid areas of the world like Sub-Sahara Africa. Cowpea is highly susceptible to cowpea aphid, Aphis craccivora, infestation that can lead to major yield losses. Aphids feed on their host plant by inserting their hypodermal needlelike flexible stylets into the plant to reach the phloem sap. During feeding, aphids secrete saliva, containing effector proteins, into the plant to disrupt plant immune responses and alter the physiology of the plant to their own advantage. Liquid chromatography tandem mass spectrometry (LC-MS/MS) was used to identify the salivary proteome of the cowpea aphid. About 150 candidate proteins were identified including diacetyl/L-xylulose reductase (DCXR), a novel enzyme previously unidentified in aphid saliva. DCXR is a member of short-chain dehydrogenases/reductases with dual enzymatic functions in carbohydrate and dicarbonyl metabolism. To assess whether cowpea aphid DCXR (AcDCXR) has similar functions, recombinant AcDCXR was purified and assayed enzymatically. For carbohydrate metabolism, the oxidation of xylitol to xylulose was tested. The dicarbonyl reaction involved the reduction of methylglyoxal, an α-β-dicarbonyl ketoaldehyde, known as an abiotic and biotic stress response molecule causing cytotoxicity at high concentrations. To assess whether cowpea aphids induce methylglyoxal in plants, we measured methylglyoxal levels in both cowpea and pea (Pisum sativum) plants and found them elevated transiently after aphid infestation. Agrobacterium-mediated transient overexpression of AcDCXR in pea resulted in an increase of cowpea aphid fecundity. Taken together, our results indicate that AcDCXR is an effector with a putative ability to generate additional sources of energy to the aphid and to alter plant defense responses. In addition, this work identified methylglyoxal as a potential novel aphid defense metabolite adding to the known repertoire of plant defenses against aphid pests.

Host interactors of effector proteins of the lettuce downy mildew Bremia lactucae obtained by yeast two-hybrid screening

Plant pathogenic bacteria, fungi and oomycetes secrete effector proteins to manipulate host cell processes to establish a successful infection. Over the last decade the genomes and transcriptomes of many agriculturally important plant pathogens have been sequenced and vast candidate effector repertoires were identified using bioinformatic analyses. Elucidating the contribution of individual effectors to pathogenicity is the next major hurdle. To advance our understanding of the molecular mechanisms underlying lettuce susceptibility to the downy mildew Bremia lactucae, we mapped physical interactions between B. lactucae effectors and lettuce candidate target proteins. Using a lettuce cDNA library-based yeast-two-hybrid system, 61 protein-protein interactions were identified, involving 21 B. lactucae effectors and 46 unique lettuce proteins. The top ten interactors based on the number of independent colonies identified in the Y2H and two interactors that belong to gene families involved in plant immunity, were further characterized. We determined the subcellular localization of the fluorescently tagged lettuce proteins and their interacting effectors. Importantly, relocalization of effectors or their interactors to the nucleus was observed for four protein-pairs upon their co-expression, supporting their interaction in planta.

Genome Sequence of Fusarium graminearum Strain CML3066, Isolated from a Wheat Spike in Southern Brazil

Fusarium graminearum is a global fungal pathogen of wheat and other small grains, causing Fusarium head blight (FHB) disease, also known as wheat scab. We report here the annotated genome of a deoxynivalenol/15-acetyl-deoxynivalenol-producing Brazilian strain called CML3066, isolated from FHB-symptomatic wheat spikes collected in 2009.

Fusarium graminearum is a global fungal pathogen of wheat and other small grains, causing Fusarium head blight (FHB) disease, also known as wheat scab. We report here the annotated genome of a deoxynivalenol/15-acetyl-deoxynivalenol-producing Brazilian strain called CML3066, isolated from FHB-symptomatic wheat spikes collected in 2009.

EffHunter: A Tool for Prediction of Effector Protein Candidates in Fungal Proteomic Databases

Pathogens are able to deliver small-secreted, cysteine-rich proteins into plant cells to enable infection. The computational prediction of effector proteins remains one of the most challenging areas in the study of plant fungi interactions. At present, there are several bioinformatic programs that can help in the identification of these proteins; however, in most cases, these programs are managed independently. Here, we present EffHunter, an easy and fast bioinformatics tool for the identification of effectors. This predictor was used to identify putative effectors in 88 proteomes using characteristics such as size, cysteine residue content, secretion signal and transmembrane domains.

Functional analyses of small secreted cysteine-rich proteins identified candidate effectors in Verticillium dahliae

Secreted small cysteine-rich proteins (SCPs) play a critical role in modulating host immunity in plant-pathogen interactions. Bioinformatic analyses showed that the fungal pathogen Verticillium dahliae encodes more than 100 VdSCPs, but their roles in host-pathogen interactions have not been fully characterized. Transient expression of 123 VdSCP-encoding genes in Nicotiana benthamiana identified three candidate genes involved in host-pathogen interactions. The expression of these three proteins, VdSCP27, VdSCP113, and VdSCP126, in N. benthamiana resulted in cell death accompanied by a reactive oxygen species burst, callose deposition, and induction of defence genes. The three VdSCPs mainly localized to the periphery of the cell. BAK1 and SOBIR1 (associated with receptor-like protein) were required for the immunity triggered by these three VdSCPs in N. benthamiana. Site-directed mutagenesis showed that cysteine residues that form disulphide bonds are essential for the functioning of VdSCP126, but not VdSCP27 and VdSCP113. VdSCP27, VdSCP113, and VdSCP126 individually are not essential for V. dahliae infection of N. benthamiana and Gossypium hirsutum, although there was a significant reduction of virulence on N. benthamiana and G. hirsutum when inoculated with the VdSCP27/VdSCP126 double deletion strain. These results illustrate that the SCPs play a critical role in the V. dahliae-plant interaction via an intrinsic virulence function and suppress immunity following infection.

Distinct Life Histories Impact Dikaryotic Genome Evolution in the Rust Fungus Puccinia striiformis Causing Stripe Rust in Wheat

Stripe rust of wheat, caused by the obligate biotrophic fungus Puccinia striiformis f.sp. tritici, is a major threat to wheat production worldwide with an estimated yearly loss of US $1 billion. The recent advances in long-read sequencing technologies and tailored-assembly algorithms enabled us to disentangle the two haploid genomes of Pst. This provides us with haplotype-specific information at a whole-genome level. Exploiting this novel information, we perform whole-genome comparative genomics of two P. striiformis f.sp. tritici isolates with contrasting life histories. We compare one isolate of the old European lineage (PstS0), which has been asexual for over 50 years, and a Warrior isolate (PstS7 lineage) from a novel incursion into Europe in 2011 from a sexual population in the Himalayan region. This comparison provides evidence that long-term asexual evolution leads to genome expansion, accumulation of transposable elements, and increased heterozygosity at the single nucleotide, structural, and allele levels. At the whole-genome level, candidate effectors are not compartmentalized and do not exhibit reduced levels of synteny. Yet we were able to identify two subsets of candidate effector populations. About 70% of candidate effectors are invariant between the two isolates, whereas 30% are hypervariable. The latter might be involved in host adaptation on wheat and explain the different phenotypes of the two isolates. Overall, this detailed comparative analysis of two haplotype-aware assemblies of P. striiformis f.sp. tritici is the first step in understanding the evolution of dikaryotic rust fungi at a whole-genome level.

Distinct Life Histories Impact Dikaryotic Genome Evolution in the Rust Fungus Puccinia striiformis Causing Stripe Rust in Wheat

Stripe rust of wheat, caused by the obligate biotrophic fungus Puccinia striiformis f.sp. tritici, is a major threat to wheat production worldwide with an estimated yearly loss of US $1 billion. The recent advances in long-read sequencing technologies and tailored-assembly algorithms enabled us to disentangle the two haploid genomes of Pst. This provides us with haplotype-specific information at a whole-genome level. Exploiting this novel information, we perform whole-genome comparative genomics of two P. striiformis f.sp. tritici isolates with contrasting life histories. We compare one isolate of the old European lineage (PstS0), which has been asexual for over 50 years, and a Warrior isolate (PstS7 lineage) from a novel incursion into Europe in 2011 from a sexual population in the Himalayan region. This comparison provides evidence that long-term asexual evolution leads to genome expansion, accumulation of transposable elements, and increased heterozygosity at the single nucleotide, structural, and allele levels. At the whole-genome level, candidate effectors are not compartmentalized and do not exhibit reduced levels of synteny. Yet we were able to identify two subsets of candidate effector populations. About 70% of candidate effectors are invariant between the two isolates, whereas 30% are hypervariable. The latter might be involved in host adaptation on wheat and explain the different phenotypes of the two isolates. Overall, this detailed comparative analysis of two haplotype-aware assemblies of P. striiformis f.sp. tritici is the first step in understanding the evolution of dikaryotic rust fungi at a whole-genome level.

Prediction of the secretomes of endophytic and nonendophytic fungi reveals similarities in host plant infection and colonization strategies

Endophytic fungi are microorganisms that inhabit internal plant tissues without causing apparent damage. During the infection process, both endophytic and phytopathogenic fungi secrete proteins to resist or supplant the plant's defense mechanisms. This study analyzed the predicted secretomes of six species of endophytic fungi and compared them with predicted secretomes of eight fungal species with different lifestyles: saprophytic, necrotrophic, hemibiotrophic, and biotrophic. The sizes of the predicted secretomes varied from 260 to 1640 proteins, and the predicted secretomes have a wide diversity of CAZymes, proteases, and conserved domains. Regarding the CAZymes in the secretomes of the analyzed fungi, the most abundant CAZyme families were glycosyl hydrolase and serine proteases. Several predicted proteins have characteristics similar to those found in small, secreted proteins with effector characteristics (SSPEC). The most abundant conserved domains, besides those found in the SSPEC, have oxidation activities, indicating that these proteins can protect the fungus against oxidative stress, against domains with protease activity, which may be involved in the mechanisms of nutrition, or against lytic enzymes secreted by the host plant. This study demonstrates that secretomes of endophytic and nonendophytic fungi share an arsenal of proteins important in the process of infection and colonization of host plants.

Label free proteomics and systematic analysis of secretome reveals effector candidates regulated by SGE1 and FTF1 in the plant pathogen Fusarium oxysporum f. sp. cubense tropical race 4

BACKGROUND

Phytopathogens secreted effectors during host colonization to suppress or trigger plant immunity. Identification of new effectors is one of the research focuses in recent years. There is only a limited knowledge about effectors of Fusarium oxysporum f. sp. Cubense tropical race 4 (Foc TR4), the causal agent of wilt disease in Cavendish banana.

RESULTS

Two transcription factors, SGE1 and FTF1, were constitutively over-expressed in Foc TR4 to partially mimic the in-planta state. Secreted proteins with high purity were prepared through a two-round extraction method. Then the secretome were analyzed via label free proteomics method. A total of 919 non-redundant proteins were detected, of which 74 proteins were predicted to be effector candidates. Among these candidates, 29 were up-regulated and 13 down-regulated in the strain over-expressing SGE1 and FTF1, 8 were up-regulated and 4 down-regulated in either SGE1 or FTF1 over expression strain.

CONCLUSIONS

Through label free proteomics analysis, a series of effector candidates were identified in secretome of Foc TR4. Our work put a foundation for functional research of these effectors.

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

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

Seeking the Roles for Fungal Small-Secreted Proteins in Affecting Saprophytic Lifestyles

Small secreted proteins (SSPs) comprise 40–60% of the total fungal secretome and are present in fungi of all phylogenetic groups, representing the entire spectrum of lifestyles. They are characteristically shorter than 300 amino acids in length and have a signal peptide. The majority of SSPs are coded by orphan genes, which lack known domains or similarities to known protein sequences. Effectors are a group of SSPs that have been investigated extensively in fungi that interact with living hosts, either pathogens or mutualistic systems. They are involved in suppressing the host defense response and altering its physiology. Here, we aim to delineate some of the potential roles of SSPs in saprotrophic fungi, that have been bioinformatically predicted as effectors, and termed in this mini-review as “effector-like” proteins. The effector-like Ssp1 from the white-rot fungus Pleurotus ostreatus is presented as a case study, and its potential role in regulating the ligninolytic system, secondary metabolism, development, and fruiting body initiation are discussed. We propose that deciphering the nature of effector-like SSPs will contribute to our understanding of development and communication in saprophytic fungi, as well as help, to elucidate the origin, regulation, and mechanisms of fungal-host, fungal-fungal, and fungal-bacterial interactions.

Time-Course RNAseq Reveals Exserohilum turcicum Effectors and Pathogenicity Determinants

Exserohilum turcicum (sexual stage Setosphaeria turcica) is the hemibiotrophic causal agent of northern leaf blight of maize and sorghum. This study aimed to identify the genes involved in host colonization during the biotrophic and necrotrophic phases of infection. It also aimed to identify race-specific differences in gene expression. RNAseq of maize seedlings inoculated with a race 13N or 23N E. turcicum isolate was conducted before inoculation and at 2, 5, 7, and 13 days post-inoculation (dpi). Biological replicates were pooled per time point for each race and sequenced. A bioinformatics pipeline was used to identify candidate effectors, and expression was validated for selected candidates. Fungal biomass was positively correlated with the percentages of E. turcicum reads mapped, which were low at early time points (2-7 dpi) with a significant increase at 13 dpi, indicating a lifestyle switch from biotrophy to necrotrophy between 7 and 13 dpi. AVRHt1 is the putative E. turcicum effector recognized by the maize resistance gene Ht1. Consistent with this, AVRHt1 was expressed in planta by race 23N, but transcripts were absent in race 13N. In addition, specific transposable elements were expressed in 23N only. Genes encoding the virulence-associated peptidases leupeptin-inhibiting protein 1 and fungalysin were expressed in planta. Transcriptional profiles of genes involved in secondary metabolite synthesis or cell wall degradation revealed the importance of these genes during late stages of infection (13 dpi). A total of 346 expressed candidate effectors were identified, including Ecp6 and proteins similar to the secreted in xylem (SIX) effectors common to formae speciales of Fusarium oxysporum, SIX13 and SIX5. Expression profiling of Ecp6 and SIX13-like indicated a peak in expression at 5 and 7 dpi compared to 2 and 13 dpi. Sequencing of SIX13-like from diverse isolates of E. turcicum revealed host-specific polymorphisms that were mostly non-synonymous, resulting in two groups of SIX13-like proteins that corresponded to the maize or sorghum origin of each isolate. This study suggests putative mechanisms whereby E. turcicum causes disease. Identification of the candidate effector SIX13-like is consistent with the infection mode of E. turcicum through the xylem of susceptible hosts.

Draft Genome Sequence of the Strawberry Anthracnose Pathogen Colletotrichum fructicola

Colletotrichum fructicola is a causal agent of strawberry anthracnose and a major economic pathogen of horticultural and ornamental crops worldwide. Here, we present an annotated draft genome sequence for a C. fructicola isolate previously used for transcriptomic analysis. The assembly totals 58.0 Mb in 477 contigs with 18,143 predicted genes.

Colletotrichum fructicola is a causal agent of strawberry anthracnose and a major economic pathogen of horticultural and ornamental crops worldwide. Here, we present an annotated draft genome sequence for a C. fructicola isolate previously used for transcriptomic analysis. The assembly totals 58.0 Mb in 477 contigs with 18,143 predicted genes.

Distinctive Gene Expression Profiles and Effectors Consistent With Host Specificity in Two Formae Speciales of Marssonina brunnea

The knowledge on the host specificity of a pathogen underlying an interaction is becoming an urgent necessity for global warming. In this study, the gene expression profiles and the roles of effectors in host specificity were integrally characterized with two formae speciales, multigermtubi and monogermtubi, of a hemibiotrophic pathogen Marssonina brunnea when they were infecting respective susceptible poplar hosts. With a functional genome comparison referring to a de novo transcriptome of M. brunnea and Pathogen-Host Interaction database functional annotations, the multigermtubi strain showed abundant and significant differentially expressed unigenes (DEGs) (more than 40%) in colonizing the initial invasion stage and in the necrotrophic stage. The monogermtubi strain induced less than 10% of DEGs in the initial invasion stage but which abruptly increased to more than 80% DEGs in the necrotrophic stage. Both strains induced the least DEGs in the biotrophic stage compared to the initial invasion and necrotrophic stages. The orthologs of the effector genes Ecp6, PemG1, XEG1, ACE1, and Mg3LysM were exclusively induced by one of the two formae speciales depending on the infection stages. Some unigenes homologous to carbohydrate lytic enzyme genes, especially pectate lyases, were notably induced with multigermtubi forma specialis infection but not expressed in the monogermtubi forma specialis at an earlier infection stage. The extraordinary differences in the functional genome level between the two formae speciales of M. brunnea could be fundamental to exploring their host specificity determinant and evolution. This study also firstly provided the fungal transcriptome of the monogermtubi forma specialis for M. brunnea.

Hybridization and introgression drive genome evolution of Dutch elm disease pathogens

Hybridization and the resulting introgression can drive the success of invasive species via the rapid acquisition of adaptive traits. The Dutch elm disease pandemics in the past 100 years were caused by three fungal lineages with permeable reproductive barriers: Ophiostoma ulmi, Ophiostoma novo-ulmi subspecies novo-ulmi and Ophiostoma novo-ulmi subspecies americana. Using whole-genome sequences and growth phenotyping of a worldwide collection of isolates, we show that introgression has been the main driver of genomic diversity and that it impacted fitness-related traits. Introgressions contain genes involved in host-pathogen interactions and reproduction. Introgressed isolates have enhanced growth rate at high temperature and produce different necrosis sizes on an in vivo model for pathogenicity. In addition, lineages diverge in many pathogenicity-associated genes and exhibit differential mycelial growth in the presence of a proxy of a host defence compound, implying an important role of host trees in the molecular and functional differentiation of these pathogens.

Different Phenotypes, Similar Genomes: Three Newly Sequenced Fusarium fujikuroi Strains Induce Different Symptoms in Rice Depending on Temperature

Bakanae, caused by the hemibiotrophic fungus Fusarium fujikuroi, is one of the most important diseases of rice and is attributed to up to 75% of losses, depending on the strain and environmental conditions. Some strains cause elongation and thin leaves, whereas others induce stunting and chlorotic seedlings. Differences in symptoms are attributed to genetic differences in the strains. F. fujikuroi strains Augusto2, CSV1, and I1.3 were sequenced with Illumina MiSeq, and pathogenicity trials were conducted on rice cultivar Galileo, which is susceptible to bakanae. By performing gene prediction, single nucleotide polymorphism (SNP) calling, and structural variant analysis with a reference genome, we show how an extremely limited number of polymorphisms in genes not commonly associated with bakanae disease can cause strong differences in phenotype. CSV1 and Augusto2 were particularly close, with only 21,887 SNPs between them, but they differed in virulence, reaction to temperature, induced symptoms, colony morphology and color, growth speed, fumonisin, and gibberellin production. Genes potentially involved in the shift in phenotype were identified. Furthermore, we show how temperature variation may result in different symptoms even in rice plants inoculated with the same F. fujikuroi strain. Moreover, all of the F. fujikuroi strains became more virulent at higher temperatures. Significant differences were likewise observed in gibberellic acid production and in the expression of both fungal and plant gibberellin biosynthetic genes.

Population genomics and evolution of a fungal pathogen after releasing exotic strains to control insect pests for 20 years

Entomopathogenic fungi are one of the key regulators of insect populations in nature. Some species such as Beauveria bassiana with a wide host range have been developed as promising alternatives to chemical insecticides for the biocontrol of insect pests. However, the long-term persistence of the released strains, the effect on non-target hosts and local fungal populations remains elusive, but they are considerable concerns with respect to environmental safety. Here we report the temporal features of the Beauveria population genomics and evolution over 20 years after releasing exotic strains to control pine caterpillar pests. We found that the isolates within the biocontrol site were mostly of clonal origins. The released strains could persist in the environment for a long time but with low recovery rates. Similar to the reoccurrence of host jumping by local isolates, the infection of non-target insects by the released strains was evident to endemically occur in association with host seasonality. No obvious dilution effect on local population structure was evident by the releases. However, the population was largely replaced by genetically divergent isolates once per decade but evolved with a pattern of balancing selection and towards expansion through adaptation, non-random outcrossing and isolate migration. This study not only unveils the real-time features of entomopathogenic fungal population genomics and evolution but also provides added values to alleviate the concerns of environmental safety regarding the biocontrol application of mycoinsecticides.

Characterization of Growth Morphology and Pathology, and Draft Genome Sequencing of Botrytis fabae, the Causal Organism of Chocolate Spot of Faba Bean (Vicia faba L.)

Chocolate spot is a major fungal disease of faba bean caused by the ascomycete fungus, Botrytis fabae. B. fabae is also implicated in botrytis gray mold disease in lentils, along with B. cinerea. Here we have isolated and characterized two B. fabae isolates from chocolate spot lesions on faba bean leaves. In plant disease assays on faba bean and lentil, B. fabae was more aggressive than B. cinerea and we observed variation in susceptibility among a small set of cultivars for both plant hosts. Using light microscopy, we observed a spreading, generalized necrosis response in faba bean toward B. fabae. In contrast, the plant response to B. cinerea was localized to epidermal cells underlying germinated spores and appressoria. In addition to the species characterization of B. fabae, we produced genome assemblies for both B. fabae isolates using Illumina sequencing. Genome sequencing coverage and assembly size for B. fabae isolates, were 27x and 45x, and 43.2 and 44.5 Mb, respectively. Following genome assembly and annotation, carbohydrate-active enzyme (CAZymes) and effector genes were predicted. There were no major differences in the numbers of each of the major classes of CAZymes. We predicted 29 effector genes for B. fabae, and using the same selection criteria for B. cinerea, we predicted 34 putative effector genes. For five of the predicted effector genes, the pairwise dN/dS ratio between orthologs from B. fabae and B. cinerea was greater than 1.0, suggesting positive selection and the potential evolution of molecular mechanisms for host specificity in B. fabae. Furthermore, a homology search of secondary metabolite clusters revealed the absence of the B. cinerea phytotoxin botrydial and several other uncharacterized secondary metabolite biosynthesis genes from B. fabae. Although there were no obvious differences in the number or proportional representation of different transposable element classes, the overall proportion of AT-rich DNA sequence in B. fabae was double that of B. cinerea.

A 19-isolate reference-quality global pangenome for the fungal wheat pathogen Zymoseptoria tritici

BACKGROUND

The gene content of a species largely governs its ecological interactions and adaptive potential. A species is therefore defined by both core genes shared between all individuals and accessory genes segregating presence-absence variation. There is growing evidence that eukaryotes, similar to bacteria, show intra-specific variability in gene content. However, it remains largely unknown how functionally relevant such a pangenome structure is for eukaryotes and what mechanisms underlie the emergence of highly polymorphic genome structures.

RESULTS

Here, we establish a reference-quality pangenome of a fungal pathogen of wheat based on 19 complete genomes from isolates sampled across six continents. Zymoseptoria tritici causes substantial worldwide losses to wheat production due to rapidly evolved tolerance to fungicides and evasion of host resistance. We performed transcriptome-assisted annotations of each genome to construct a global pangenome. Major chromosomal rearrangements are segregating within the species and underlie extensive gene presence-absence variation. Conserved orthogroups account for only ~ 60% of the species pangenome. Investigating gene functions, we find that the accessory genome is enriched for pathogenesis-related functions and encodes genes involved in metabolite production, host tissue degradation and manipulation of the immune system. De novo transposon annotation of the 19 complete genomes shows that the highly diverse chromosomal structure is tightly associated with transposable element content. Furthermore, transposable element expansions likely underlie recent genome expansions within the species.

CONCLUSIONS

Taken together, our work establishes a highly complex eukaryotic pangenome providing an unprecedented toolbox to study how pangenome structure impacts crop-pathogen interactions.

A 19-isolate reference-quality global pangenome for the fungal wheat pathogen Zymoseptoria tritici

BACKGROUND

The gene content of a species largely governs its ecological interactions and adaptive potential. A species is therefore defined by both core genes shared between all individuals and accessory genes segregating presence-absence variation. There is growing evidence that eukaryotes, similar to bacteria, show intra-specific variability in gene content. However, it remains largely unknown how functionally relevant such a pangenome structure is for eukaryotes and what mechanisms underlie the emergence of highly polymorphic genome structures.

RESULTS

Here, we establish a reference-quality pangenome of a fungal pathogen of wheat based on 19 complete genomes from isolates sampled across six continents. Zymoseptoria tritici causes substantial worldwide losses to wheat production due to rapidly evolved tolerance to fungicides and evasion of host resistance. We performed transcriptome-assisted annotations of each genome to construct a global pangenome. Major chromosomal rearrangements are segregating within the species and underlie extensive gene presence-absence variation. Conserved orthogroups account for only ~ 60% of the species pangenome. Investigating gene functions, we find that the accessory genome is enriched for pathogenesis-related functions and encodes genes involved in metabolite production, host tissue degradation and manipulation of the immune system. De novo transposon annotation of the 19 complete genomes shows that the highly diverse chromosomal structure is tightly associated with transposable element content. Furthermore, transposable element expansions likely underlie recent genome expansions within the species.

CONCLUSIONS

Taken together, our work establishes a highly complex eukaryotic pangenome providing an unprecedented toolbox to study how pangenome structure impacts crop-pathogen interactions.

Target the core: durable plant resistance against filamentous plant pathogens through effector recognition

Plant pathogens colonize their host through the secretion of effector proteins that modulate plant metabolism and immune responses to their benefit. Plants evolve towards effector recognition, leading to host immunity. Typically, pathogen effectors are targets for recognition through plant receptors that are encoded by resistance genes. Resistance gene mediated crop immunity puts a tremendous pressure on pathogens to adapt and alter their effector repertoire to overcome recognition. We argue that the type of effector that is recognized by the host may have considerable implications on the durability of resistance against filamentous plant pathogens. Effector genes that are conserved among pathogens and reside in core genome regions are most likely to hold indispensable virulence functions. Consequently, the cost for the pathogen to overcome recognition by the host is higher than for diversified, host-specific effectors with a quantitative impact on virulence. Consequently, resistance genes that directly target conserved effector proteins without the interception of other effector proteins are potentially excellent resistance resources. © 2019 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Population Genomic Analyses Reveal Connectivity via Human-Mediated Transport across Populus Plantations in North America and an Undescribed Subpopulation of Sphaerulina musiva

Domestication of plant species has affected the evolutionary dynamics of plant pathogens in agriculture and forestry. A model system for studying the consequences of plant domestication on the evolution of an emergent plant disease is the fungal pathogen Sphaerulina musiva. This ascomycete causes leaf spot and stem canker disease of Populus spp. and their hybrids. A population genomics approach was used to determine the degree of population structure and evidence for selection on the North American population of S. musiva. In total, 122 samples of the fungus were genotyped identifying 120,016 single-nucleotide polymorphisms after quality filtering. In North America, S. musiva has low to moderate degrees of differentiation among locations. Three main genetic clusters were detected: southeastern United States, midwestern United States and Canada, and a new British Columbia cluster (BC2). Population genomics suggest that BC2 is a novel genetic cluster from central British Columbia, clearly differentiated from previously reported S. musiva from coastal British Columbia, and the product of a single migration event. Phenotypic measurements from greenhouse experiments indicate lower aggressiveness of BC2 on Populus trichocarpa. In summary, S. musiva has geographic structure across broad regions indicative of gene flow among clusters. The interconnectedness of the North American S. musiva populations across large geographic distances further supports the hypothesis of anthropogenic-facilitated transport of the pathogen.

Chromosome-Level Reference Genome of Venturia effusa, Causative Agent of Pecan Scab

Pecan scab, caused by Venturia effusa, is a devastating disease of pecan (Carya illinoinensis), which results in economic losses on susceptible cultivars throughout the southeastern United States. To enhance our understanding of pathogenicity in V. effusa, we have generated a complete telomere-to-telomere reference genome of V. effusa isolate FRT5LL7-Albino. By combining Illumina MiSeq and Oxford Nanopore MinION data, we assembled a 45.2-Mb genome represented by 20 chromosomes and containing 10,820 putative genes, of which 7,619 have at least one functional annotation. The likely causative mutation of the albino phenotype was identified as a single base insertion and a resulting frameshift in the gene encoding the polyketide synthase ALM1. This genome represents the first full chromosome-level assembly of any Venturia sp.

A Comparative Genomic Analysis of the Barley Pathogen Pyrenophora teres f. teres Identifies Subtelomeric Regions as Drivers of Virulence

Pyrenophora teres f. teres causes net form net blotch of barley and is an economically important pathogen throughout the world. However, P. teres f. teres is lacking in the genomic resources necessary to characterize the mechanisms of virulence. Recently a high-quality reference genome was generated for P. teres f. teres isolate 0-1. Here, we present the reference quality sequence and annotation of four new isolates and we use the five available P. teres f. teres genomes for an in-depth comparison, resulting in the generation of hypotheses pertaining to the potential mechanisms and evolution of virulence. Comparative analyses were performed between all five P. teres f. teres genomes, examining genomic organization, structural variations, and core and accessory genomic content, specifically focusing on the genomic characterization of known virulence loci and the localization of genes predicted to encode secreted and effector proteins. We showed that 14 of 15 currently published virulence quantitative trait loci (QTL) span accessory genomic regions, consistent with these accessory regions being important drivers of host adaptation. Additionally, these accessory genomic regions were frequently found in subtelomeric regions of chromosomes, with 10 of the 14 accessory region QTL localizing to subtelomeric regions. Comparative analysis of the subtelomeric regions of P. teres f. teres chromosomes revealed translocation events in which homology was detected between nonhomologous chromosomes at a significantly higher rate than the rest of the genome. These results indicate that the subtelomeric accessory genomic compartments not only harbor most of the known virulence loci but, also, that these regions have the capacity to rapidly evolve.

Genomics Evolutionary History and Diagnostics of the Alternaria alternata Species Group Including Apple and Asian Pear Pathotypes

The Alternaria section alternaria (Alternaria alternata species group) represents a diverse group of saprotroph, human allergens, and plant pathogens. Alternaria taxonomy has benefited from recent phylogenetic revision but the basis of differentiation between major phylogenetic clades within the group is not yet understood. Furthermore, genomic resources have been limited for the study of host-specific pathotypes. We report near complete genomes of the apple and Asian pear pathotypes as well as draft assemblies for a further 10 isolates representing Alternaria tenuissima and Alternaria arborescens lineages. These assemblies provide the first insights into differentiation of these taxa as well as allowing the description of effector and non-effector profiles of apple and pear conditionally dispensable chromosomes (CDCs). We define the phylogenetic relationship between the isolates sequenced in this study and a further 23 Alternaria spp. based on available genomes. We determine which of these genomes represent MAT1-1-1 or MAT1-2-1 idiomorphs and designate host-specific pathotypes. We show for the first time that the apple pathotype is polyphyletic, present in both the A. arborescens and A. tenuissima lineages. Furthermore, we profile a wider set of 89 isolates for both mating type idiomorphs and toxin gene markers. Mating-type distribution indicated that gene flow has occurred since the formation of A. tenuissima and A. arborescens lineages. We also developed primers designed to AMT14, a gene from the apple pathotype toxin gene cluster with homologs in all tested pathotypes. These primers allow identification and differentiation of apple, pear, and strawberry pathotypes, providing new tools for pathogen diagnostics.

Comparative genomics of Alternaria species provides insights into the pathogenic lifestyle of Alternaria brassicae - a pathogen of the Brassicaceae family

BACKGROUND

Alternaria brassicae, a necrotrophic pathogen, causes Alternaria Leaf Spot, one of the economically important diseases of Brassica crops. Many other Alternaria spp. such as A. brassicicola and A. alternata are known to cause secondary infections in the A. brassicae-infected Brassicas. The genome architecture, pathogenicity factors, and determinants of host-specificity of A. brassicae are unknown. In this study, we annotated and characterised the recently announced genome assembly of A. brassicae and compared it with other Alternaria spp. to gain insights into its pathogenic lifestyle.

RESULTS

We also sequenced the genomes of two A. alternata isolates that were co-infecting B. juncea using Nanopore MinION sequencing for additional comparative analyses within the Alternaria genus. Genome alignments within the Alternaria spp. revealed high levels of synteny between most chromosomes with some intrachromosomal rearrangements. We show for the first time that the genome of A. brassicae, a large-spored Alternaria species, contains a dispensable chromosome. We identified 460 A. brassicae-specific genes, which included many secreted proteins and effectors. Furthermore, we have identified the gene clusters responsible for the production of Destruxin-B, a known pathogenicity factor of A. brassicae.

CONCLUSION

The study provides a perspective into the unique and shared repertoire of genes within the Alternaria genus and identifies genes that could be contributing to the pathogenic lifestyle of A. brassicae.

Chloroplasts and Plant Immunity: Where Are the Fungal Effectors?

Chloroplasts play a central role in plant immunity through the synthesis of secondary metabolites and defense compounds, as well as phytohormones, such as jasmonic acid and salicylic acid. Additionally, chloroplast metabolism results in the production of reactive oxygen species and nitric oxide as defense molecules. The impact of viral and bacterial infections on plastids and chloroplasts has been well documented. In particular, bacterial pathogens are known to introduce effectors specifically into chloroplasts, and many viral proteins interact with chloroplast proteins to influence viral replication and movement, and plant defense. By contrast, clear examples are just now emerging for chloroplast-targeted effectors from fungal and oomycete pathogens. In this review, we first present a brief overview of chloroplast contributions to plant defense and then discuss examples of connections between fungal interactions with plants and chloroplast function. We then briefly consider well-characterized bacterial effectors that target chloroplasts as a prelude to discussing the evidence for fungal effectors that impact chloroplast activities.