Regulation of proteinaceous effector expression in phytopathogenic fungi

Uncovering the Mechanisms: The Role of Biotrophic Fungi in Activating or Suppressing Plant Defense Responses

This paper discusses the mechanisms by which fungi manipulate plant physiology and suppress plant defense responses by producing effectors that can target various host proteins. Effector-triggered immunity and effector-triggered susceptibility are pivotal elements in the complex molecular dialogue underlying plant-pathogen interactions. Pathogen-produced effector molecules possess the ability to mimic pathogen-associated molecular patterns or hinder the binding of pattern recognition receptors. Effectors can directly target nucleotide-binding domain, leucine-rich repeat receptors, or manipulate downstream signaling components to suppress plant defense. Interactions between these effectors and receptor-like kinases in host plants are critical in this process. Biotrophic fungi adeptly exploit the signaling networks of key plant hormones, including salicylic acid, jasmonic acid, abscisic acid, and ethylene, to establish a compatible interaction with their plant hosts. Overall, the paper highlights the importance of understanding the complex interplay between plant defense mechanisms and fungal effectors to develop effective strategies for plant disease management.

The Velvet transcription factor PnVeA regulates necrotrophic effectors and secondary metabolism in the wheat pathogen Parastagonospora nodorum

The fungus Parastagonospora nodorum causes septoria nodorum blotch on wheat. The role of the fungal Velvet-family transcription factor VeA in P. nodorum development and virulence was investigated here. Deletion of the P. nodorum VeA ortholog, PnVeA, resulted in growth abnormalities including pigmentation, abolished asexual sporulation and highly reduced virulence on wheat. Comparative RNA-Seq and RT-PCR analyses revealed that the deletion of PnVeA also decoupled the expression of major necrotrophic effector genes. In addition, the deletion of PnVeA resulted in an up-regulation of four predicted secondary metabolite (SM) gene clusters. Using liquid-chromatography mass-spectrometry, it was observed that one of the SM gene clusters led to an accumulation of the mycotoxin alternariol. PnVeA is essential for asexual sporulation, full virulence, secondary metabolism and necrotrophic effector regulation.

HacA, a key transcription factor for the unfolded protein response, is required for fungal development, aflatoxin biosynthesis and pathogenicity of Aspergillus flavus

Aspergillus flavus is a fungus notorious for contaminating food and feed with aflatoxins. As a saprophytic fungus, it secretes large amounts of enzymes to access nutrients, making endoplasmic reticulum (ER) homeostasis important for protein folding and secretion. The role of HacA, a key transcription factor in the unfolded protein response pathway, remains poorly understood in A. flavus. In this study, the hacA gene in A. flavus was knockout. Results showed that the absence of hacA led to a decreased pathogenicity of the strain, as it failed to colonize intact maize kernels. This may be due to retarded vegetable growth, especially the abnormal development of swollen tips and shorter hyphal septa. Deletion of hacA also hindered conidiogenesis and sclerotial development. Notably, the mutant strain failed to produce aflatoxin B1. Moreover, compared to the wild type, the mutant strain showed increased sensitivity to ER stress inducer such as Dithiothreitol (DTT), and heat stress. It also displayed heightened sensitivity to other environmental stresses, including cell wall, osmotic, and pH stresses. Further transcriptomic analysis revealed the involvement of the hacA in numerous biological processes, including filamentous growth, asexual reproduction, mycotoxin biosynthetic process, signal transduction, budding cell apical bud growth, invasive filamentous growth, response to stimulus, and so on. Taken together, HacA plays a vital role in fungal development, pathogenicity and aflatoxins biosynthesis. This highlights the potential of targeting hacA as a novel approach for early prevention of A. flavus contamination.

A MYB-related transcription factor regulates effector gene expression in an oomycete pathogen

Phytophthora pathogens possess hundreds of effector genes that exhibit diverse expression patterns during infection, yet how the expression of effector genes is precisely regulated remains largely elusive. Previous studies have identified a few potential conserved transcription factor binding sites (TFBSs) in the promoters of Phytophthora effector genes. Here, we report a MYB-related protein, PsMyb37, in Phytophthora sojae, the major causal agent of root and stem rot in soybean. Yeast one-hybrid and electrophoretic mobility shift assays showed that PsMyb37 binds to the TACATGTA motif, the most prevalent TFBS in effector gene promoters. The knockout mutant of PsMyb37 exhibited significantly reduced virulence on soybean and was more sensitive to oxidative stress. Consistently, transcriptome analysis showed that numerous effector genes associated with suppressing plant immunity or scavenging reactive oxygen species were down-regulated in the PsMyb37 knockout mutant during infection compared to the wild-type P. sojae. Several promoters of effector genes were confirmed to drive the expression of luciferase in a reporter assay. These results demonstrate that a MYB-related transcription factor contributes to the expression of effector genes in P. sojae.

Regulation of effector gene expression as concerted waves in Leptosphaeria maculans: a two-player game

Effector genes, encoding molecules involved in disease establishment, are concertedly expressed throughout the lifecycle of plant-pathogenic fungi. However, little is known about how effector gene expression is regulated. Since many effector genes are located in repeat-rich regions, the role of chromatin remodeling in their regulation was recently investigated, notably establishing that the repressive histone modification H3K9me3, deposited by KMT1, was involved in several fungal species including Leptosphaeria maculans. Nevertheless, previous data suggest that a second regulatory layer, probably involving a specific transcription factor (TF), might be required. In L. maculans, a Dothideomycete causing stem canker of oilseed rape, we identified the ortholog of Pf2, a TF belonging to the Zn2Cys6 fungal-specific family, and described as essential for pathogenicity and effector gene expression. We investigated its role together with KMT1, by inactivating and over-expressing LmPf2 in a wild-type strain and a ∆kmt1 mutant. Functional analyses of the corresponding transformants highlighted an essential role of LmPf2 in the establishment of pathogenesis and we found a major effect of LmPf2 on the induction of effector gene expression once KMT1 repression is lifted. Our results show, for the first time, a dual control of effector gene expression.

Distribution and conservation of simple sequence repeats in plant pathogenic species of Zymoseptoria and development of genomic resources for its orphaned species

To better understand the structure and evolution of the genomes of four plant pathogenic species of Zymoseptoria, we analyzed the occurrence, relative abundance (RA), and density (RD) of simple sequence repeats (SSRs) in their whole genome and transcriptome sequences. In this study, SSRs are defined as repeats of more than 12 bases in length. The genome and transcriptome sequences of Zymoseptoria ardabiliae show the highest RA (201.1 and 129.9) and RD (3229.4 and 1928.2) of SSRs, while those of Zymoseptoria pseudotritici show the lowest RA (167.2 and 118.5) and RD (2482.2 and 1687.0). The majority of SSRs in the genomic and transcriptome sequences of species were trinucleotide SSRs, while dinucleotide SSRs were the least common. The most common trinucleotide motifs in the transcriptomic sequences across all species were those that encoded the amino acid arginine. As per our motif conservation study, Zymoseptoria tritici (12.4%) possessed the most unique motifs, while Z. pseudotritici (3.9%) had the fewest. Overall, only 38.1% of the motifs were found to be conserved among the species. Gene enrichment studies reveal that three of the species, Z. ardabiliae, Zymoseptoria brevis, and Z. pseudotritici, have SSRs in their genes related to cellular metabolism, while the remaining Z. tritici harbors SSRs in genes related to DNA synthesis and gene expression. In an effort to improve the genetic resources for the orphan species of pathogenic Zymoseptoria, a total of 73,134 primers were created. The genomic resources developed in this study could help with analyses of genetic relatedness within the population and the development of species-specific markers.

The biological relevance of the FspTF transcription factor, homologous of Bqt4, in Fusarium sp. associated with the ambrosia beetle Xylosandrus morigerus

Transcription factors in phytopathogenic fungi are key players due to their gene expression regulation leading to fungal growth and pathogenicity. The KilA-N family encompasses transcription factors unique to fungi, and the Bqt4 subfamily is included in it and is poorly understood in filamentous fungi. In this study, we evaluated the role in growth and pathogenesis of the homologous of Bqt4, FspTF, in Fusarium sp. isolated from the ambrosia beetle Xylosandrus morigerus through the characterization of a CRISPR/Cas9 edited strain in Fsptf. The phenotypic analysis revealed that TF65-6, the edited strain, modified its mycelia growth and conidia production, exhibited affectation in mycelia and culture pigmentation, and in the response to certain stress conditions. In addition, the plant infection process was compromised. Untargeted metabolomic and transcriptomic analysis, clearly showed that FspTF may regulate secondary metabolism, transmembrane transport, virulence, and diverse metabolic pathways such as lipid metabolism, and signal transduction. These data highlight for the first time the biological relevance of an orthologue of Bqt4 in Fusarium sp. associated with an ambrosia beetle.

The SET domain protein PsKMT3 regulates histone H3K36 trimethylation and modulates effector gene expression in the soybean pathogen Phytophthora sojae

Plant pathogens secrete effector proteins to overcome host immunity and promote colonization. In oomycete plant pathogens, the expression of many effector genes is altered upon infection; however, the regulatory mechanisms are unclear. In this study, we identified a su(var)3-9, enhancer of zeste, and trithorax (SET) domain protein-encoding gene, PsKMT3, that was highly induced at early infection stages in Phytophthora sojae. Deletion of PsKMT3 led to asexual development and pathogenicity defects. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) and western blot analyses demonstrated that histone H3K36 trimethylation (H3K36me3) was significantly reduced genome-wide in mutants. RNA-seq analysis identified 374 genes encoding secreted proteins that were differentially expressed in pskmt3 at the mycelium stage. The significantly altered genes encompassed the RxLR (Arg-x-Lys-Arg) effector gene family, including the essential effector genes Avh23, Avh181, Avh240, and Avh241. Transcriptome analysis at early infection stages showed misregulation of effector gene expression waves in pskmt3. H3K36me3 was directly and indirectly associated with RxLR effector gene activation. Our results reveal a role of a SET domain protein in regulating effector gene expression and modulating histone methylation in P. sojae.

Transcription factor lineages in plant-pathogenic fungi, connecting diversity with fungal virulence

Plant-pathogenic fungi span diverse taxonomic lineages. Their host-infection strategies are often specialised and require the coordinated regulation of molecular virulence factors. Transcription factors (TFs) are fundamental regulators of gene expression, yet relatively few virulence-specific regulators are characterised in detail and their evolutionary trajectories are not well understood. Hence, this study compared the full range of TFs across taxonomically-diverse fungal proteomes and classified their lineages through an orthology analysis. The primary aims were to characterise differences in the range and profile of TF lineages broadly linked to plant-host association or pathogenic lifestyles, and to better characterise the evolutionary origin and trajectory of experimentally-validated virulence regulators. We observed significantly fewer TFs among obligate, host-associated pathogens, largely attributed to contractions in several Zn2Cys6 TF-orthogroup lineages. We also present novel insight into the key virulence-regulating TFs Ste12, Pf2 and EBR1, providing evidence for their ancestral origins, expansion and/or loss. Ultimately, the analysis presented here provides both primary evidence for TF evolution in fungal phytopathogenicity, as well as a practical phylogenetic resource to guide further detailed investigation on the regulation of virulence within key pathogen lineages.

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

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

Variability in an effector gene promoter of a necrotrophic fungal pathogen dictates epistasis and effector-triggered susceptibility in wheat

The fungus Parastagonospora nodorum uses proteinaceous necrotrophic effectors (NEs) to induce tissue necrosis on wheat leaves during infection, leading to the symptoms of septoria nodorum blotch (SNB). The NEs Tox1 and Tox3 induce necrosis on wheat possessing the dominant susceptibility genes Snn1 and Snn3B1/Snn3D1, respectively. We previously observed that Tox1 is epistatic to the expression of Tox3 and a quantitative trait locus (QTL) on chromosome 2A that contributes to SNB resistance/susceptibility. The expression of Tox1 is significantly higher in the Australian strain SN15 compared to the American strain SN4. Inspection of the Tox1 promoter region revealed a 401 bp promoter genetic element in SN4 positioned 267 bp upstream of the start codon that is absent in SN15, called PE401. Analysis of the world-wide P. nodorum population revealed that a high proportion of Northern Hemisphere isolates possess PE401 whereas the opposite was observed in representative P. nodorum isolates from Australia and South Africa. The presence of PE401 removed the epistatic effect of Tox1 on the contribution of the SNB 2A QTL but not Tox3. PE401 was introduced into the Tox1 promoter regulatory region in SN15 to test for direct regulatory roles. Tox1 expression was markedly reduced in the presence of PE401. This suggests a repressor molecule(s) binds PE401 and inhibits Tox1 transcription. Infection assays also demonstrated that P. nodorum which lacks PE401 is more pathogenic on Snn1 wheat varieties than P. nodorum carrying PE401. An infection competition assay between P. nodorum isogenic strains with and without PE401 indicated that the higher Tox1-expressing strain rescued the reduced virulence of the lower Tox1-expressing strain on Snn1 wheat. Our study demonstrated that Tox1 exhibits both 'selfish' and 'altruistic' characteristics. This offers an insight into a complex NE-NE interaction that is occurring within the P. nodorum population. The importance of PE401 in breeding for SNB resistance in wheat is discussed.

A Novel DCL2-Dependent Micro-Like RNA Vm-PC-3p-92107_6 Affects Pathogenicity by Regulating the Expression of Vm-VPS10 in Valsa mali

Dicer proteins are mainly responsible for generating small RNAs (sRNAs), which are involved in gene silencing in most eukaryotes. In previous research, two DCL proteins in Valsa mali, the pathogenic fungus causing apple tree Valsa canker, were found associated with both the pathogenicity and generation of sRNAs. In this study, the differential expression of small interfering RNAs (siRNAs) and miRNA-like RNAs (milRNAs) was analyzed based on the deep sequencing of the wild type and Vm-DCL2 mutant, respectively. Overall, the generation of 40 siRNAs and 18 milRNAs was evidently associated with Vm-DCL2. The target genes of milRNAs were then identified using degradome sequencing; according to the prediction results, most candidate targets are related to pathogenicity. Further, expression of Vm-PC-3p-92107_6 was confirmed in the wild type but not in the Vm-DCL2 mutant. Moreover, the pathogenicity of Vm-PC-3p-92107_6 deletion mutants (ΔVm-PC-3p-92107_6) and the over-expression transformants (Vm-PC-3p-92107_6-OE) was significantly increased and decreased, respectively. Based on those degradome results, vacuolar protein sorting 10 (Vm-VPS10) was identified as the target of Vm-PC-3p-92107_6. Co-expression analysis in tobacco leaves further confirmed that Vm-PC-3p-92107_6 could suppress the expression of Vm-VPS10. Meanwhile, the expression levels of Vm-PC-3p-92107_6 and Vm-VPS10 displayed divergent trends in ΔVm-PC-3p-92107_6 and Vm-PC-3p-92107_6-OE, respectively. Perhaps most importantly, ΔVm-VPS10 featured a significant reduction in pathogenicity. Taken together, our results indicate that a DCL2-dependent milRNA Vm-PC-3p-92107_6 plays roles in pathogenicity by regulating the expression of Vm-VPS10. This study lays a foundation for the comprehensive analysis of pathogenic mechanisms of V. mali and deepens our understanding of the generation and function of fungal sRNA.

SsPEP1, an Effector with Essential Cellular Functions in Sugarcane Smut Fungus

Biotrophic fungi have to infect their host to obtain nutrients and must establish an interaction with the host to complete their life cycle. In this process, effectors play important roles in manipulating the host's immune system to avoid being attacked. Sporisorium scitamineum is the causative agent of sugarcane smut, the most important disease in sugarcane-producing regions worldwide. In this work, we functionally characterized the conserved effector PEP1 in S. scitamineum. The mating process and the expression of genes in the MAPK signaling pathway and the a and b loci were adversely affected in Sspep1-null mutants. The requirement for SsPEP1 in pathogenicity and symptom development was allele dosage-dependent, i.e., deleting one Sspep1 allele in the mating pair turned a normal black whip with abundant teliospores into a white whip with few teliospores; however, deleting both alleles almost abolished infectivity and whip development. ΔSspep1 mutants produced significantly less mycelium mass within infected plants. Additionally, SsPEP1 was identified as a potent inhibitor of sugarcane POD-1a peroxidase activity, implying that SsPEP1 may function to relieve reactive oxygen species-related stress within the host plant. Taken together, our work demonstrated that SsPEP1 is a multifaceted effector essential for S. scitamineum growth, development, and pathogenicity.

Reactive Oxygen Species in Host Plant Are Required for an Early Defense Response against Attack of Stagonospora nodorum Berk. Necrotrophic Effectors SnTox

Reactive oxygen species (ROS) play a central role in plant immune responses. The most important virulence factors of the Stagonospora nodorum Berk. are multiple fungal necrotrophic effectors (NEs) (SnTox) that affect the redox-status and cause necrosis and/or chlorosis in wheat lines possessing dominant susceptibility genes (Snn). However, the effect of NEs on ROS generation at the early stages of infection has not been studied. We studied the early stage of infection of various wheat genotypes with S nodorum isolates -Sn4VD, SnB, and Sn9MN, carrying a different set of NE genes. Our results indicate that all three NEs of SnToxA, SnTox1, SnTox3 significantly contributed to cause disease, and the virulence of the isolates depended on their differential expression in plants (Triticum aestivum L.). The Tsn1–SnToxA, Snn1–SnTox1and Snn3–SnTox3 interactions played an important role in inhibition ROS production at the initial stage of infection. The Snn3–SnTox3 inhibited ROS production in wheat by affecting NADPH-oxidases, peroxidases, superoxide dismutase and catalase. The Tsn1–SnToxA inhibited ROS production in wheat by affecting peroxidases and catalase. The Snn1–SnTox1 inhibited the production of ROS in wheat by mainly affecting a peroxidase. Collectively, these results show that the inverse gene-for gene interactions between effector of pathogen and product of host sensitivity gene suppress the host’s own PAMP-triggered immunity pathway, resulting in NE-triggered susceptibility (NETS). These results are fundamentally changing our understanding of the development of this economical important wheat disease.

Genome-wide mapping of histone modifications during axenic growth in two species of Leptosphaeria maculans showing contrasting genomic organization

Leptosphaeria maculans 'brassicae' (Lmb) and Leptosphaeria maculans 'lepidii' (Lml) are closely related phytopathogenic species that exhibit a large macrosynteny but contrasting genome structure. Lmb has more than 30% of repeats clustered in large repeat-rich regions, while the Lml genome has only a small amount of evenly distributed repeats. Repeat-rich regions of Lmb are enriched in effector genes, expressed during plant infection. The distinct genome structures of Lmb and Lml provide an excellent model for comparing the organization of pathogenicity genes in relation to the chromatin landscape in two closely related phytopathogenic fungi. Here, we performed chromatin immunoprecipitation (ChIP) during axenic culture, targeting histone modifications typical for heterochromatin or euchromatin, combined with transcriptomic analysis to analyze the influence of chromatin organization on gene expression. In both species, we found that facultative heterochromatin is enriched with genes lacking functional annotation, including numerous effector and species-specific genes. Notably, orthologous genes located in H3K27me3 domains are enriched with effector genes. Compared to other fungal species, including Lml, Lmb is distinct in having large H3K9me3 domains associated with repeat-rich regions that contain numerous species-specific effector genes. Discovery of these two distinctive heterochromatin landscapes now raises questions about their involvement in the regulation of pathogenicity, the dynamics of these domains during plant infection and the selective advantage to the fungus to host effector genes in H3K9me3 or H3K27me3 domains.

Research Progress on Phytopathogenic Fungi and Their Role as Biocontrol Agents

Phytopathogenic fungi decrease crop yield and quality and cause huge losses in agricultural production. To prevent the occurrence of crop diseases and insect pests, farmers have to use many synthetic chemical pesticides. The extensive use of these pesticides has resulted in a series of environmental and ecological problems, such as the increase in resistant weed populations, soil compaction, and water pollution, which seriously affect the sustainable development of agriculture. This review discusses the main advances in research on plant-pathogenic fungi in terms of their pathogenic factors such as cell wall-degrading enzymes, toxins, growth regulators, effector proteins, and fungal viruses, as well as their application as biocontrol agents for plant pests, diseases, and weeds. Finally, further studies on plant-pathogenic fungal resources with better biocontrol effects can help find new beneficial microbial resources that can control diseases.

GWAS analysis reveals distinct pathogenicity profiles of Australian Parastagonospora nodorum isolates and identification of marker-trait-associations to septoria nodorum blotch

The fungus Parastagonospora nodorum is the causal agent of septoria nodorum leaf blotch (SNB) and glume blotch which are common in many wheat growing regions in the world. The disease is complex and could be explained by multiple interactions between necrotrophic effectors secreted by the pathogen and matching susceptibility genes in wheat. An Australian P. nodorum population was clustered into five groups with contrasting properties. This study was set to identify their pathogenicity profiles using a diverse wheat panel of 134 accessions which are insensitive to SnToxA and SnTox1 in both in vitro and in vivo conditions. SNB seedling resistance/susceptibility to five representative isolates from the five clusters, responses to crude culture-filtrates (CFs) of three isolates and sensitivity to SnTox3 semi-purified effector together with 11,455 SNP markers have been used for linkage disequilibrium (LD) and association analyses. While quantitative trait loci (QTL) on 1D, 2A, 2B, 4B, 5B, 6A, 6B, 7A, 7D chromosomes were consistently detected across isolates and conditions, distinct patterns and isolate specific QTL were also observed among these isolates. In this study, SnTox3–Snn3-B1 interaction for the first time in Australia and SnTox3–Snn3-D1 interaction for the first time in bread wheat were found active using wild-type isolates. These findings could be due to new SnTox3 haplotype/isoform and exotic CIMMYT/ICARDA and Vavilov germplasm used, respectively. This study could provide useful information for dissecting novel and different SNB disease components, helping to prioritise research targets and contributing valuable information on genetic loci/markers for marker-assisted selection in SNB resistance wheat breeding programme.

Tandem DNA repeats contain cis-regulatory sequences that activate biotrophy-specific expression of Magnaporthe effector gene PWL2

During plant infection, fungi secrete effector proteins in coordination with distinct infection stages. Thus, the success of plant infection is determined by precise control of effector gene expression. We analysed the PWL2 effector gene of the rice blast fungus Magnaporthe oryzae to understand how effector genes are activated specifically during the early biotrophic stages of rice infection. Here, we used confocal live-cell imaging of M. oryzae transformants with various PWL2 promoter fragments fused to sensitive green fluorescent protein reporter genes to determine the expression patterns of PWL2 at the cellular level, together with quantitative reverse transcription PCR analyses at the tissue level. We found PWL2 expression was coupled with sequential biotrophic invasion of rice cells. PWL2 expression was induced in the appressorium upon penetration into a living rice cell but greatly declined in the highly branched hyphae when the first-invaded rice cell was dead. PWL2 expression then increased again as the hyphae penetrate into living adjacent cells. The expression of PWL2 required fungal penetration into living plant cells of either host rice or nonhost onion. Deletion and mutagenesis experiments further revealed that the tandem repeats in the PWL2 promoter contain 12-base pair sequences required for expression. We conclude that PWL2 expression is (a) activated by an unknown signal commonly present in living plant cells, (b) specific to biotrophic stages of fungal infection, and (c) requires 12-base pair cis-regulatory sequences in the promoter.

Every cloud has a silver lining: how abiotic stresses affect gene expression in plant-pathogen interactions

Current environmental and climate changes are having a pronounced influence on the outcome of plant-pathogen interactions, further highlighting the fact that abiotic stresses strongly affect biotic interactions at various levels. For instance, physiological parameters such as plant architecture and tissue organization together with primary and specialized metabolism are affected by environmental constraints, and these combine to make an individual plant either a more or less suitable host for a given pathogen. In addition, abiotic stresses can affect the timely expression of plant defense and pathogen virulence. Indeed, several studies have shown that variations in temperature, and in water and mineral nutrient availability affect the expression of plant defense genes. The expression of virulence genes, known to be crucial for disease outbreak, is also affected by environmental conditions, potentially modifying existing pathosystems and paving the way for emerging pathogens. In this review, we summarize our current knowledge on the impact of abiotic stress on biotic interactions at the transcriptional level in both the plant and the pathogen side of the interaction. We also perform a metadata analysis of four different combinations of abiotic and biotic stresses, which identifies 197 common modulated genes with strong enrichment in Gene Ontology terms related to defense . We also describe the multistress-specific responses of selected defense-related genes.

Vm-milR37 contributes to pathogenicity by regulating glutathione peroxidase gene VmGP in Valsa mali

MicroRNAs play important roles in various biological processes by regulating their corresponding target genes. However, the function and regulatory mechanism of fungal microRNA-like RNAs (milRNAs) are still largely unknown. In this study, a milRNA (Vm-milR37) was isolated and identified from Valsa mali, which causes the most serious disease on the trunk of apple trees in China. Based on the results of deep sequencing and quantitative reverse transcription PCR, Vm-milR37 was found to be expressed in the mycelium, while it was not expressed during the V. mali infection process. Overexpression of Vm-milR37 did not affect vegetative growth, but significantly decreased pathogenicity. Based on degradome sequencing, the target of Vm-milR37 was identified as VmGP, a glutathione peroxidase. The expression of Vm-milR37 and VmGP showed a divergent trend in V. mali-apple interaction samples and Vm-milR37 overexpression transformants. The expression of VmGP could be suppressed significantly by Vm-milR37 when coexpressed in tobacco leaves. Deletion of VmGP showed significantly reduced pathogenicity compared with the wild type. VmGP deletion mutants showed more sensitivity to hydrogen peroxide. Apple leaves inoculated with Vm-milR37 overexpression transformants and VmGP deletion mutant displayed increased accumulation of reactive oxygen species compared with the wild type. Thus, Vm-milR37 plays a critical role in pathogenicity by regulating VmGP, which contributes to the oxidative stress response during V. mali infection. These results provide important evidence to define the roles of milRNAs and their corresponding target genes in pathogenicity.

Leaping into the Unknown World of Sporisorium scitamineum Candidate Effectors

Sporisorium scitamineum is a biotrophic fungus causing sugarcane smut disease. In this study, we set up a pipeline and used genomic and dual transcriptomic data previously obtained by our group to identify candidate effectors of S. scitamineum and their expression profiles in infected smut-resistant and susceptible sugarcane plants. The expression profile of different genes after infection in contrasting sugarcane genotypes assessed by RT-qPCR depended on the plant genotypes and disease progression. Three candidate effector genes expressed earlier only in resistant plants, four expressed in both genotypes, and three later in susceptible plants. Ten genes were cloned and transiently expressed in N. benthamiana leaves to determine their subcellular location, while four localized in more than one compartment. Two candidates, g3890 having a nucleoplasmic and mitochondrial location and g5159 targeting the plant cell wall, were selected to obtain their possible corresponding host targets using co-immunoprecipitation (CoIP) experiments and mass spectrometry. Various potential interactors were identified, including subunits of the protein phosphatase 2A and an endochitinase. We investigated the presence of orthologs in sugarcane and using transcriptome data present their expression profiles. Orthologs of sugarcane shared around 70% similarity. Identifying a set of putative fungal effectors and their plant targets provides a valuable resource for functional characterization of the molecular events leading to smut resistance in sugarcane plants and uncovers further opportunities for investigation.

Verticillium longisporum Elicits Media-Dependent Secretome Responses With Capacity to Distinguish Between Plant-Related Environments

Verticillia cause a vascular wilt disease affecting a broad range of economically valuable crops. The fungus enters its host plants through the roots and colonizes the vascular system. It requires extracellular proteins for a successful plant colonization. The exoproteomes of the allodiploid Verticillium longisporum upon cultivation in different media or xylem sap extracted from its host plant Brassica napus were compared. Secreted fungal proteins were identified by label free liquid chromatography-tandem mass spectrometry screening. V. longisporum induced two main secretion patterns. One response pattern was elicited in various non-plant related environments. The second pattern includes the exoprotein responses to the plant-related media, pectin-rich simulated xylem medium and pure xylem sap, which exhibited similar but additional distinct features. These exoproteomes include a shared core set of 221 secreted and similarly enriched fungal proteins. The pectin-rich medium significantly induced the secretion of 143 proteins including a number of pectin degrading enzymes, whereas xylem sap triggered a smaller but unique fungal exoproteome pattern with 32 enriched proteins. The latter pattern included proteins with domains of known pathogenicity factors, metallopeptidases and carbohydrate-active enzymes. The most abundant proteins of these different groups are the necrosis and ethylene inducing-like proteins Nlp2 and Nlp3, the cerato-platanin proteins Cp1 and Cp2, the metallopeptidases Mep1 and Mep2 and the carbohydrate-active enzymes Gla1, Amy1 and Cbd1. Their pathogenicity contribution was analyzed in the haploid parental strain V. dahliae. Deletion of the majority of the corresponding genes caused no phenotypic changes during ex planta growth or invasion and colonization of tomato plants. However, we discovered that the MEP1, NLP2, and NLP3 deletion strains were compromised in plant infections. Overall, our exoproteome approach revealed that the fungus induces specific secretion responses in different environments. The fungus has a general response to non-plant related media whereas it is able to fine-tune its exoproteome in the presence of plant material. Importantly, the xylem sap-specific exoproteome pinpointed Nlp2 and Nlp3 as single effectors required for successful V. dahliae colonization.

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.

Long-Read-Based de novo Genome Assembly and Comparative Genomics of the Wheat Leaf Rust Pathogen Puccinia triticina Identifies Candidates for Three Avirulence Genes

Leaf rust, caused by Puccinia triticina (Pt), is one of the most devastating diseases of wheat, affecting production in nearly all wheat-growing regions worldwide. Despite its economic importance, genomic resources for Pt are very limited. In the present study, we have used long-read sequencing (LRS) and the pipeline of FALCON and FALCON-Unzip (v4.1.0) to carry out the first LRS-based de novo genome assembly for Pt. Using 22.4-Gb data with an average read length of 11.6 kb and average coverage of 150-fold, we generated a genome assembly for Pt104 [strain 104-2,3,(6),(7),11; isolate S423], considered to be the founding isolate of a clonal lineage of Pt in Australia. The Pt104 genome contains 162 contigs with a total length of 140.5 Mb and N50 of 2 Mb, with the associated haplotigs providing haplotype information for 91% of the genome. This represents the best quality of Pt genome assembly to date, which reduces the contig number by 91-fold and improves the N50 by 4-fold as compared to the previous Pt race1 assembly. An annotation pipeline that combined multiple lines of evidence including the transcriptome assemblies derived from RNA-Seq, previously identified expressed sequence tags and Pt race 1 protein sequences predicted 29,043 genes for Pt104 genome. Based on the presence of a signal peptide, no transmembrane segment, and no target location to mitochondria, 2,178 genes were identified as secreted proteins (SPs). Whole-genome sequencing (Illumina paired-end) was performed for Pt104 and six additional strains with differential virulence profile on the wheat leaf rust resistance genes Lr26, Lr2a, and Lr3ka. To identify candidates for the corresponding avirulence genes AvrLr26, AvrLr2a, and AvrLr3ka, genetic variation within each strain was first identified by mapping to the Pt104 genome. Variants within predicted SP genes between the strains were then correlated to the virulence profiles, identifying 38, 31, and 37 candidates for AvrLr26, AvrLr2a, and AvrLr3ka, respectively. The identification of these candidate genes lays a good foundation for future studies on isolating these avirulence genes, investigating the molecular mechanisms underlying host-pathogen interactions, and the development of new diagnostic tools for pathogen monitoring.

Functional characterization of an inducible bidirectional promoter from Fusarium oxysporum f. sp. cubense

Bidirectional promoters (BDPs) are regulatory DNA sequences (~1000 bp long) intervening two genes arranged on opposite strands with their 5' ends in close proximity. These genes are mostly co-expressed; but, instances of anti-correlation and independent transcription have been observed. In fungal systems, BDPs have shown to provide an improved genetic circuit by assembling and regulating transcription of different genes of a common metabolic pathway. We have identified an intergenic region (1063 bp) from the genome of Fusarium oxysporum f. sp. cubense (Foc), a banana root pathogen. This intergenic region regulates the expression of a gene pair required for the breakdown of hemicellulose. For characterization, it was cloned into pCSN44 vector backbone between two reporter genes, namely β-glucuronidase (GUS) and enhanced green fluorescent protein (EGFP). The newly formed vector was transformed into Foc and tested for its bidirectional expression activity. Using histochemical staining and fluorescence microscopy, the kinetics for both, GUS and EGFP expression were tested under different growth conditions respectively. The activity was differentially regulated by inducers such as xylan, arabinogalactan and pectin. This is the first report on the isolation of the intergenic region with inducible bidirectional promoter activity from Fusarium. Characterization of such BDPs will find applications in genetic engineering, metabolic engineering and synthetic biology using fungal systems.

Low Amplitude Boom-and-Bust Cycles Define the Septoria Nodorum Blotch Interaction

INTRODUCTION

Septoria nodorum blotch (SNB) is a complex fungal disease of wheat caused by the Dothideomycete fungal pathogen Parastagonospora nodorum. The fungus infects through the use of necrotrophic effectors (NEs) that cause necrosis on hosts carrying matching dominant susceptibility genes. The Western Australia (WA) wheatbelt is a SNB "hot spot" and experiences significant under favorable conditions. Consequently, SNB has been a major target for breeders in WA for many years.

MATERIALS AND METHODS

In this study, we assembled a panel of 155 WA P. nodorum isolates collected over a 44-year period and compared them to 23 isolates from France and the USA using 28 SSR loci.

RESULTS

The WA P. nodorum population was clustered into five groups with contrasting properties. 80% of the studied isolates were assigned to two core groups found throughout the collection location and time. The other three non-core groups that encompassed transient and emergent populations were found in restricted locations and time. Changes in group genotypes occurred during periods that coincided with the mass adoption of a single or a small group of widely planted wheat cultivars. When introduced, these cultivars had high scores for SNB resistance. However, the field resistance of these new cultivars often declined over subsequent seasons prompting their replacement with new, more resistant varieties. Pathogenicity assays showed that newly emerged isolates non-core are more pathogenic than old isolates. It is likely that the non-core groups were repeatedly selected for increased virulence on the contemporary popular cultivars.

DISCUSSION

The low level of genetic diversity within the non-core groups, difference in virulence, low abundance, and restriction to limited locations suggest that these populations more vulnerable to a population crash when the cultivar was replaced by one that was genetically different and more resistant. We characterize the observed pattern as a low-amplitude boom-and-bust cycle in contrast with the classical high amplitude boom-and-bust cycles seen for biotrophic pathogens where the contrast between resistance and susceptibility is typically much greater. Implications of the results are discussed relating to breeding strategies for more sustainable SNB resistance and more generally for pathogens with NEs.

Effect of the host-specific toxin SnTOX3 from Stagonospora nodorum on ethylene signaling pathway regulation and redox-state in common wheat

The fungus Stagonospora nodorum Berk. is the causative agent of Septoria nodorum blotch (SNB) of wheat. The most important factors of Stagonospora nodorum virulence include numerous fungal necrotrophic effectors (NEs) encoded by SnTox genes. They interact with the matching products of host susceptibility genes (Snn). SnTox-Snn interactions are mirror images of classical gene-for-gene interactions and lead to the development of disease. We have studied the SnTox3-Snn3 interaction, resulting in the development of infection on leaves and formation of extensive lesions. The mechanism of SnTox3 action is likely to be linked to the regulation of redox metabolism and the influence on ethylene synthesis in the wheat plants, although the molecular mechanisms are not fully unveiled. To characterize the SnTox3-Snn3 interaction, we used S. nodorum isolates differing in the expression of the NEs genes SnTox3 (SnB (Tox3+), Sn4VD (Tox3–)) and two soft spring wheat (Triticum aestivum L.) cultivars, contrasting in resistance to the SNB agent and differing in the allelic composition of the susceptibility locus Snn3-B1: Kazakhstanskaya 10 (susceptible) and Omskaya 35 (resistant). We carried out a comparative assessment of the transcriptional activity patterns of genes responsible for ethylene biosynthesis (TaACS1, TaACО) and signaling pathway (TaEIL1, TaERF1) by real-time PCR and estimated the redox state of wheat plants infected with different isolates of S. nodorum by spectrometry. The induction of ethylene biosynthesis and signaling has been shown to result from gene-for-gene interaction between Snn3-B1 and SnTox3. The results of plant redox status estimation showed that ethylene inhibited accumulation of hydrogen peroxide in SnTox3-sensitive genotypes by regulating the operation of various pro-/antioxidant enzymes at the transcriptional and posttranslational levels. Our results suggest that NE SnTox3 influences ethylene biosynthesis and signaling, thereby regulating redox metabolism in infected wheat plants as necessary for successful host colonization at the initial phases of infection, which ultimately leads to extensive lesions due to fast pathogen reproduction.

Genome-wide identification of long non-coding RNAs suggests a potential association with effector gene transcription in Phytophthora sojae

Numerous long non-coding RNAs (lncRNAs) identified and characterized in mammals, plants and fungi have been found to play critical regulatory roles in biological processes. However, little is known about the role of lncRNAs in oomycete plant pathogens, which cause devastating damage to the economy and ecosystems. We used strand-specific RNA sequencing (RNA-seq) to generate a computational pipeline to identify lncRNAs in Phytophthora sojae, a model oomycete plant pathogen. In total, 940 lncRNAs with 1010 isoforms were identified from RNA-seq data obtained from four representative stages of P. sojae. The lncRNAs had shorter transcript lengths, longer exon lengths, fewer numbers of exons, lower GC content and higher minimum free energy values compared with protein-coding genes. lncRNAs in P. sojae exhibited low sequence conservation amongst oomycetes and P. sojae isolates. Transcriptional data indicated that P. sojae lncRNAs tended to be transcribed in a stage-specific manner; representative lncRNAs were validated by semi-quantitative reverse transcription-polymerase chain reaction. Phytophthora sojae lncRNAs were concentrated in gene-sparse regions, and lncRNAs were associated with secreted protein and effector coding genes. The neighbouring genes of lncRNAs encoded various effector family members, and RNA-seq data revealed a correlation between the transcription level of lncRNAs and their neighbouring genes. Our results provide the first comprehensive identification of lncRNAs in oomycetes and suggest a potential association between lncRNAs and effector genes.

The Genome Biology of Effector Gene Evolution in Filamentous Plant Pathogens

Filamentous pathogens, including fungi and oomycetes, pose major threats to global food security. Crop pathogens cause damage by secreting effectors that manipulate the host to the pathogen's advantage. Genes encoding such effectors are among the most rapidly evolving genes in pathogen genomes. Here, we review how the major characteristics of the emergence, function, and regulation of effector genes are tightly linked to the genomic compartments where these genes are located in pathogen genomes. The presence of repetitive elements in these compartments is associated with elevated rates of point mutations and sequence rearrangements with a major impact on effector diversification. The expression of many effectors converges on an epigenetic control mediated by the presence of repetitive elements. Population genomics analyses showed that rapidly evolving pathogens show high rates of turnover at effector loci and display a mosaic in effector presence-absence polymorphism among strains. We conclude that effective pathogen containment strategies require a thorough understanding of the effector genome biology and the pathogen's potential for rapid adaptation.

The genome assembly of the fungal pathogen Pyrenochaeta lycopersici from Single-Molecule Real-Time sequencing sheds new light on its biological complexity

The first draft genome sequencing of the non-model fungal pathogen Pyrenochaeta lycopersici showed an expansion of gene families associated with heterokaryon incompatibility and lacking of mating-type genes, providing insights into the genetic basis of this “imperfect” fungus which lost the ability to produce the sexual stage. However, due to the Illumina short-read technology, the draft genome was too fragmented to allow a comprehensive characterization of the genome, especially of the repetitive sequence fraction. In this work, the sequencing of another P. lycopersici isolate using long-read Single Molecule Real-Time sequencing technology was performed with the aim of obtaining a gapless genome. Indeed, a gapless genome assembly of 62.7 Mb was obtained, with a fraction of repetitive sequences representing 30% of the total bases. The gene content of the two P. lycopersici isolates was very similar, and the large difference in genome size (about 8 Mb) might be attributable to the high fraction of repetitive sequences detected for the new sequenced isolate. The role of repetitive elements, including transposable elements, in modulating virulence effectors is well established in fungal plant pathogens. Moreover, transposable elements are of fundamental importance in creating and re-modelling genes, especially in imperfect fungi. Their abundance in P. lycopersici, together with the large expansion of heterokaryon incompatibility genes in both sequenced isolates, suggest the presence of possible mechanisms alternative to gene re-assorting mediated by sexual recombination. A quite large fraction (~9%) of repetitive elements in P. lycopersici, has no homology with known classes, strengthening this hypothesis. The availability of a gapless genome of P. lycopersici allowed the in-depth analysis of its genome content, by annotating functional genes and TEs. This goal will be an important resource for shedding light on the evolution of the reproductive and pathogenic behaviour of this soilborne pathogen and the onset of a possible speciation within this species.

The global regulator of pathogenesis PnCon7 positively regulates Tox3 effector gene expression through direct interaction in the wheat pathogen Parastagonospora nodorum

To investigate effector gene regulation in the wheat pathogenic fungus Parastagonospora nodorum, the promoter and expression of Tox3 was characterised through a series of complementary approaches. Promoter deletion and DNase I footprinting experiments identified a 25 bp region in the Tox3 promoter as being required for transcription. Subsequent yeast one-hybrid analysis using the DNA sequence as bait identified that interacting partner as the C2H2 zinc finger transcription factor PnCon7, a putative master regulator of pathogenesis. Silencing of PnCon7 resulted in the down-regulation of Tox3 demonstrating that the transcription factor has a positive regulatory role on gene expression. Analysis of Tox3 expression in the PnCon7 silenced strains revealed a strong correlation with PnCon7 transcript levels, supportive of a direct regulatory role. Subsequent pathogenicity assays using PnCon7-silenced isolates revealed that the transcription factor was required for Tox3-mediated disease. The expression of two other necrotrophic effectors (ToxA and Tox1) was also affected but in a non-dose dependent manner suggesting that the regulatory role of PnCon7 on these genes was indirect. Collectively, these data have advanced our fundamental understanding of the Con7 master regulator of pathogenesis by demonstrating its positive regulatory role on the Tox3 effector in P. nodorum through direct interaction.

De Novo Assembly and Phasing of Dikaryotic Genomes from Two Isolates of Puccinia coronata f. sp. avenae, the Causal Agent of Oat Crown Rust

Oat crown rust, caused by the fungus Pucinnia coronata f. sp. avenae, is a devastating disease that impacts worldwide oat production. For much of its life cycle, P. coronata f. sp. avenae is dikaryotic, with two separate haploid nuclei that may vary in virulence genotype, highlighting the importance of understanding haplotype diversity in this species. We generated highly contiguous de novo genome assemblies of two P. coronata f. sp. avenae isolates, 12SD80 and 12NC29, from long-read sequences. In total, we assembled 603 primary contigs for 12SD80, for a total assembly length of 99.16 Mbp, and 777 primary contigs for 12NC29, for a total length of 105.25 Mbp; approximately 52% of each genome was assembled into alternate haplotypes. This revealed structural variation between haplotypes in each isolate equivalent to more than 2% of the genome size, in addition to about 260,000 and 380,000 heterozygous single-nucleotide polymorphisms in 12SD80 and 12NC29, respectively. Transcript-based annotation identified 26,796 and 28,801 coding sequences for isolates 12SD80 and 12NC29, respectively, including about 7,000 allele pairs in haplotype-phased regions. Furthermore, expression profiling revealed clusters of coexpressed secreted effector candidates, and the majority of orthologous effectors between isolates showed conservation of expression patterns. However, a small subset of orthologs showed divergence in expression, which may contribute to differences in virulence between 12SD80 and 12NC29. This study provides the first haplotype-phased reference genome for a dikaryotic rust fungus as a foundation for future studies into virulence mechanisms in P. coronata f. sp. avenaeIMPORTANCE Disease management strategies for oat crown rust are challenged by the rapid evolution of Puccinia coronata f. sp. avenae, which renders resistance genes in oat varieties ineffective. Despite the economic importance of understanding P. coronata f. sp. avenae, resources to study the molecular mechanisms underpinning pathogenicity and the emergence of new virulence traits are lacking. Such limitations are partly due to the obligate biotrophic lifestyle of P. coronata f. sp. avenae as well as the dikaryotic nature of the genome, features that are also shared with other important rust pathogens. This study reports the first release of a haplotype-phased genome assembly for a dikaryotic fungal species and demonstrates the amenability of using emerging technologies to investigate genetic diversity in populations of P. coronata f. sp. avenae.

Transcription Factor Repertoire of Necrotrophic Fungal Phytopathogen Ascochyta rabiei: Predominance of MYB Transcription Factors As Potential Regulators of Secretome

Transcription factors (TFs) are the key players in gene expression and their study is highly significant for shedding light on the molecular mechanisms and evolutionary history of organisms. During host-pathogen interaction, extensive reprogramming of gene expression facilitated by TFs is likely to occur in both host and pathogen. To date, the knowledge about TF repertoire in filamentous fungi is in infancy. The necrotrophic fungus Ascochyta rabiei, that causes destructive Ascochyta blight (AB) disease of chickpea (Cicer arietinum), demands more comprehensive study for better understanding of Ascochyta-legume pathosystem. In the present study, we performed the genome-wide identification and analysis of TFs in A. rabiei. Taking advantage of A. rabiei genome sequence, we used a bioinformatic approach to predict the TF repertoire of A. rabiei. For identification and classification of A. rabiei TFs, we designed a comprehensive pipeline using a combination of BLAST and InterProScan software. A total of 381 A. rabiei TFs were predicted and divided into 32 fungal specific families of TFs. The gene structure, domain organization and phylogenetic analysis of abundant families of A. rabiei TFs were also carried out. Comparative study of A. rabiei TFs with that of other necrotrophic, biotrophic, hemibiotrophic, symbiotic, and saprotrophic fungi was performed. It suggested presence of both conserved as well as unique features among them. Moreover, cis-acting elements on promoter sequences of earlier predicted A. rabiei secretome were also identified. With the help of published A. rabiei transcriptome data, the differential expression of TF and secretory protein coding genes was analyzed. Furthermore, comprehensive expression analysis of few selected A. rabiei TFs using quantitative real-time polymerase chain reaction revealed variety of expression patterns during host colonization. These genes were expressed in at least one of the time points tested post infection. Overall, this study illustrates the first genome-wide identification and analysis of TF repertoire of A. rabiei. This work would provide the basis for further studies to dissect role of TFs in the molecular mechanisms during A. rabiei-chickpea interactions.