Blast Fungal Genomes Show Frequent Chromosomal Changes, Gene Gains and Losses, and Effector Gene Turnover

Transposable elements impact the population divergence of rice blast fungus Magnaporthe oryzae

Dynamic transposition of transposable elements (TEs) in fungal pathogens has significant impact on genome stability, gene expression, and virulence to the host. In Magnaporthe oryzae, genome plasticity resulting from TE insertion is a major driving force leading to the rapid evolution and diversification of this fungus. Despite their importance in M. oryzae population evolution and divergence, our understanding of TEs in this context remains limited. Here, we conducted a genome-wide analysis of TE transposition dynamics in the 11 most abundant TE families in M. oryzae populations. Our results show that these TEs have specifically expanded in recently isolated M. oryzae rice populations, with the presence/absence polymorphism of TE insertions highly concordant with population divergence on Geng/Japonica and Xian/Indica rice cultivars. Notably, the genes targeted by clade-specific TEs showed clade-specific expression patterns and are involved in the pathogenic process, suggesting a transcriptional regulation of TEs on targeted genes. Our study provides a comprehensive analysis of TEs in M. oryzae populations and demonstrates a crucial role of recent TE bursts in adaptive evolution and diversification of the M. oryzae rice-infecting lineage.

IMPORTANCE

Magnaporthe oryzae is the causal agent of the destructive blast disease, which caused massive loss of yield annually worldwide. The fungus diverged into distinct clades during adaptation toward the two rice subspecies, Xian/Indica and Geng/Japonica. Although the role of TEs in the adaptive evolution was well established, mechanisms underlying how TEs promote the population divergence of M. oryzae remain largely unknown. In this study, we reported that TEs shape the population divergence of M. oryzae by differentially regulating gene expression between Xian/Indica-infecting and Geng/Japonica-infecting populations. Our results revealed a TE insertion-mediated gene expression adaption that led to the divergence of M. oryzae population infecting different rice subspecies.

The massive 340 megabase genome of Anisogramma anomala, a biotrophic ascomycete that causes eastern filbert blight of hazelnut

BACKGROUND

The ascomycete fungus Anisogramma anomala causes Eastern Filbert Blight (EFB) on hazelnut (Corylus spp.) trees. It is a minor disease on its native host, the American hazelnut (C. americana), but is highly destructive on the commercially important European hazelnut (C. avellana). In North America, EFB has historically limited commercial production of hazelnut to west of the Rocky Mountains. A. anomala is an obligately biotrophic fungus that has not been grown in continuous culture, rendering its study challenging. There is a 15-month latency before symptoms appear on infected hazelnut trees, and only a sexual reproductive stage has been observed. Here we report the sequencing, annotation, and characterization of its genome.

RESULTS

The genome of A. anomala was assembled into 108 scaffolds totaling 342,498,352 nt with a GC content of 34.46%. Scaffold N50 was 33.3 Mb and L50 was 5. Nineteen scaffolds with lengths over 1 Mb constituted 99% of the assembly. Telomere sequences were identified on both ends of two scaffolds and on one end of another 10 scaffolds. Flow cytometry estimated the genome size of A. anomala at 370 Mb. The genome exhibits two-speed evolution, with 93% of the assembly as AT-rich regions (32.9% GC) and the other 7% as GC-rich (57.1% GC). The AT-rich regions consist predominantly of repeats with low gene content, while 90% of predicted protein coding genes were identified in GC-rich regions. Copia-like retrotransposons accounted for more than half of the genome. Evidence of repeat-induced point mutation (RIP) was identified throughout the AT-rich regions, and two copies of the rid gene and one of dim-2, the key genes in the RIP mutation pathway, were identified in the genome. Consistent with its homothallic sexual reproduction cycle, both MAT1-1 and MAT1-2 idiomorphs were found. We identified a large suite of genes likely involved in pathogenicity, including 614 carbohydrate active enzymes, 762 secreted proteins and 165 effectors.

CONCLUSIONS

This study reveals the genomic structure, composition, and putative gene function of the important pathogen A. anomala. It provides insight into the molecular basis of the pathogen's life cycle and a solid foundation for studying EFB.

A fungal core effector exploits the OsPUX8B.2-OsCDC48-6 module to suppress plant immunity

Proteins containing a ubiquitin regulatory X (UBX) domain are cofactors of Cell Division Cycle 48 (CDC48) and function in protein quality control. However, whether and how UBX-containing proteins participate in host-microbe interactions remain unclear. Here we show that MoNLE1, an effector from the fungal pathogen Magnaporthe oryzae, is a core virulence factor that suppresses rice immunity by specifically interfering with OsPUX8B.2. The UBX domain of OsPUX8B.2 is required for its binding to OsATG8 and OsCDC48-6 and controls its 26 S proteasome-dependent stability. OsPUX8B.2 and OsCDC48-6 positively regulate plant immunity against blast fungus, while the high-temperature tolerance heat-shock protein OsBHT, a putative cytoplasmic substrate of OsPUX8B.2-OsCDC48-6, negatively regulates defense against blast infection. MoNLE1 promotes the nuclear migration and degradation of OsPUX8B.2 and disturbs its association with OsBHT. Given the high conservation of MoNLE1 among fungal isolates, plants with broad and durable blast resistance might be generated by engineering intracellular proteins resistant to MoNLE1.

Recombination and repeat-induced point mutation landscapes reveal trade-offs between the sexual and asexual cycles of Magnaporthe oryzae

The fungal disease caused by Magnaporthe oryzae is one of the most devastating diseases that endanger many crops worldwide. Evidence shows that sexual reproduction can be advantageous for fungal diseases as hybridization facilitates host-jumping. However, the pervasive clonal lineages of M. oryzae observed in natural fields contradict this expectation. A better understanding of the roles of recombination and the fungi-specific repeat-induced point mutation (RIP) in shaping its evolutionary trajectory is essential to bridge this knowledge gap. Here we systematically investigate the RIP and recombination landscapes in M. oryzae using a whole genome sequencing data from 252 population samples and 92 cross progenies. Our data reveal that the RIP can robustly capture the population history of M. oryzae, and we provide accurate estimations of the recombination and RIP rates across different M. oryzae clades. Significantly, our results highlight a parent-of-origin bias in both recombination and RIP rates, tightly associating with their sexual potential and variations of effector proteins. This bias suggests a critical trade-off between generating novel allelic combinations in the sexual cycle to facilitate host-jumping and stimulating transposon-associated diversification of effectors in the asexual cycle to facilitate host coevolution. These findings provide unique insights into understanding the evolution of blast fungus.

Loss of PWT7, Located on a Supernumerary Chromosome, Is Associated with Parasitic Specialization of Pyricularia oryzae on Wheat

Pyricularia oryzae, a blast fungus of gramineous plants, is composed of various host genus-specific pathotypes. The avirulence of an Avena isolate on wheat is conditioned by PWT3 and PWT4. We isolated the third avirulence gene from the Avena isolate and designated it as PWT7. PWT7 was effective as an avirulence gene only at the seedling stage or on leaves. PWT7 homologs were widely distributed in a subpopulation of the Eleusine pathotype and the Lolium pathotype but completely absent in the Triticum pathotype (the wheat blast fungus). The PWT7 homolog found in the Eleusine pathotype was one of the five genes involved in its avirulence on wheat. A comparative analysis of distribution of PWT7 and the other two genes previously identified in the Eleusine pathotype suggested that, in the course of parasitic specialization toward the wheat blast fungus, a common ancestor of the Eleusine, Lolium, Avena, and Triticum pathotypes first lost PWT6, secondly PWT7, and, finally, the function of PWT3. PWT7 or its homologs were located on core chromosomes in Setaria and Eleusine isolates but on supernumerary chromosomes in Lolium and Avena isolates. This is an example of interchromosomal translocations of effector genes between core and supernumerary chromosomes. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Pyricularia's Capability of Infecting Different Grasses in Two Regions of Mexico

The genus Pyricularia includes species that are phytopathogenic fungi, which infect different species of Poaceae, such as rice and sorghum. However, few isolates have been genetically characterized in North America. The current study addresses this lack of information by characterizing an additional 57 strains of three grasses (Stenotaphrum secundatum, Cenchrus ciliaris and Digitaria ciliaris) from two distant regions of Mexico. A Pyricularia dataset with ITS sequences retrieved from GenBank and the studied sequences were used to build a haplotype network that allowed us to identify a few redundant haplotypes highly related to P. oryzae species. An analysis considering only the Mexican sequences allowed us to identify non-redundant haplotypes in the isolates of C. ciliaris and D. ciliaris, with a high identity with P. pennisetigena. The Pot2-TIR genomic fingerprinting technique resulted in high variability and allowed for the isolates to be grouped according to their host grass, whilst the ERIC-PCR technique was able to separate the isolates according to their host grass and their region of collection. Representative isolates from different host grasses were chosen to explore the pathogenic potential of these isolates. The selected isolates showed a differential pathogenic profile. Cross-infection with representative isolates from S. secundatum and C. ciliaris showed that these were unable to infect D. ciliaris grass and that the DY1 isolate from D. ciliaris was only able to infect its host grass. The results support the identification of pathogenic strains of Pyricularia isolates and their cross-infection potential in different grasses surrounding important crops in Mexico.

Comparison of Fungal Genera Isolated from Cucumber Plants and Rhizosphere Soil by Using Various Cultural Media

Plant endophytic fungi and rhizosphere soil fungi are often reported as biocontrol agents against plant pathogens or with plant growth promotion potential. Four treatments were performed in field and greenhouse experiments where cucumber plants were inoculated with Trichoderma harzianum and Fusarium oxysporum in 2022. The roots, stems and leaves of cucumber plants and their rhizosphere soil were collected twice individually from the field and greenhouse for isolation of cucumber endophytic and rhizosphere soil fungi. All fungal strains were identified through sequence similarity of the ITS1-5.8s-ITS2 rDNA region. The potato dextrose agar (PDA) media yielded the highest number of genera isolated from cucumber plants, rhizosphere soil and both compared to other media. There were no significant differences among the four media for the isolation of all cucumber endophytic fungi. However, in the roots, the number of endophytic fungi isolated by MRBA was significantly higher than that isolated on malt extract agar (MEA), while in the stems, the number of fungi isolated with PDA was significantly higher than that isolated with Martin's rose bengal agar medium (MRBA). PDA had significantly higher isolation efficiency for the rhizosphere soil fungi than MRBA. The 28 fungal genera had high isolation efficiency, and the endophytic Trichoderma strains were significantly more isolated by MEA than those of MRBA. It is suggested that PDA can be used as a basic medium, and different cultural media can be considered for specific fungal genera.

Adaptive evolution in virulence effectors of the rice blast fungus Pyricularia oryzae

Plant pathogens secrete proteins called effectors that target host cellular processes to promote disease. Recently, structural genomics has identified several families of fungal effectors that share a similar three-dimensional structure despite remarkably variable amino-acid sequences and surface properties. To explore the selective forces that underlie the sequence variability of structurally-analogous effectors, we focused on MAX effectors, a structural family of effectors that are major determinants of virulence in the rice blast fungus Pyricularia oryzae. Using structure-informed gene annotation, we identified 58 to 78 MAX effector genes per genome in a set of 120 isolates representing seven host-associated lineages. The expression of MAX effector genes was primarily restricted to the early biotrophic phase of infection and strongly influenced by the host plant. Pangenome analyses of MAX effectors demonstrated extensive presence/absence polymorphism and identified gene loss events possibly involved in host range adaptation. However, gene knock-in experiments did not reveal a strong effect on virulence phenotypes suggesting that other evolutionary mechanisms are the main drivers of MAX effector losses. MAX effectors displayed high levels of standing variation and high rates of non-synonymous substitutions, pointing to widespread positive selection shaping the molecular diversity of MAX effectors. The combination of these analyses with structural data revealed that positive selection acts mostly on residues located in particular structural elements and at specific positions. By providing a comprehensive catalog of amino acid polymorphism, and by identifying the structural determinants of the sequence diversity, our work will inform future studies aimed at elucidating the function and mode of action of MAX effectors.

The completed genome sequence of Pestalotiopsis versicolor, a pathogenic ascomycete fungus with implications for bayberry production

The pathogenic fungus Pestalotiopsis versicolor is a major etiological agent of fungal twig blight disease affecting bayberry trees. However, the lack of complete genome sequence information for this crucial pathogenic fungus hinders the molecular and genetic investigation of its pathogenic mechanism. To address this knowledge gap, we have generated the complete genome sequence of P. versicolor strain XJ27, employing a combination of Illumina, PacBio, and Hi-C sequencing technologies. This comprehensive genome sequence, comprising 7 chromosomes with an N50 contig size of 7,275,017 bp, a GC content ratio of 50.16%, and a total size of 50.80 Mb, encompasses 13,971 predicted coding genes. By performing comparative genomic analysis between P. versicolor and the genomes of eleven plant-pathogenic fungi, as well as three closely related fungi within the same group, we have gained initial insights into its evolutionary trajectory, particularly through gene family analysis. These findings shed light on the distinctive characteristics and evolutionary history of P. versicolor. Importantly, the availability of this high-quality genetic resource will serve as a foundational tool for investigating the biology, molecular pathogenesis, and virulence of P. versicolor. Furthermore, it will facilitate the development of more potent antifungal medications by uncovering potential vulnerabilities in its genetic makeup.

Horizontally Transferred DNA in the Genome of the Fungus Pyricularia oryzae is Associated With Repressive Histone Modifications

Horizontal gene transfer (HGT) is a means of exchanging genetic material asexually. The process by which horizontally transferred genes are domesticated by the host genome is of great interest but is not well understood. In this study, we determined the telomere-to-telomere genome sequence of the wheat-infecting Pyricularia oryzae strain Br48. SNP analysis indicated that the Br48 strain is a hybrid of wheat- and Brachiaria-infecting strains by a sexual or parasexual cross. Comparative genomic analysis identified several megabase-scale "insertions" in the Br48 genome, some of which were possibly gained by HGT-related events from related species, such as P. pennisetigena or P. grisea. Notably, the mega-insertions often contained genes whose phylogeny is not congruent with the species phylogeny. Moreover, some of the genes have a close homolog even in distantly related organisms, such as basidiomycetes or prokaryotes, implying the involvement of multiple HGT events. Interestingly, the levels of the silent epigenetic marks H3K9me3 and H3K27me3 in a genomic region tended to be negatively correlated with the phylogenetic concordance of genes in the same region, suggesting that horizontally transferred DNA is preferentially targeted for epigenetic silencing. Indeed, the putative HGT-derived genes were activated when MoKmt6, the gene responsible for H3K27me3 modification, was deleted. Notably, these genes also tended to be up-regulated during infection, suggesting that they are now under host control and have contributed to establishing a fungal niche. In conclusion, this study suggests that epigenetic modifications have played an important role in the domestication of HGT-derived genes in the P. oryzae genome.

Host Specificity Controlled by PWL1 and PWL2 Effector Genes in the Finger Millet Blast Pathogen Magnaporthe oryzae in Eastern Africa

Magnaporthe oryzae, a devastating pathogen of finger millet (Eleusine coracana), secretes effector molecules during infection to manipulate host immunity. This study determined the presence of avirulence effector genes PWL1 and PWL2 in 221 Eleusine blast isolates from eastern Africa. Most Ethiopian isolates carried both PWL1 and PWL2. Kenyan and Ugandan isolates largely lacked both genes, and Tanzanian isolates carried either PWL1 or lacked both. The roles of PWL1 and PWL2 towards pathogenicity on alternative chloridoid hosts, including weeping lovegrass (Eragrostis curvula), were also investigated. PWL1 and PWL2 were cloned from Ethiopian isolate E22 and were transformed separately into Ugandan isolate U34, which lacked both genes. Resulting transformants harboring either gene gained varying degrees of avirulence on Eragrostis curvula but remained virulent on finger millet. Strains carrying one or both PWL1 and PWL2 infected the chloridoid species Sporobolus phyllotrichus and Eleusine tristachya, indicating the absence of cognate resistance (R) genes for PWL1 and PWL2 in these species. Other chloridoid grasses, however, were fully resistant, regardless of the presence of one or both PWL1 and PWL2, suggesting the presence of effective R genes against PWL and other effectors. Partial resistance in some Eragrostis curvula accessions to some blast isolates lacking PWL1 and PWL2 also indicated the presence of other interactions between fungal avirulence (AVR) genes and host resistance (R) genes. Related chloridoid species thus harbor resistance genes that could be useful to improve finger millet for blast resistance. Conversely, loss of AVR genes in the fungus could expand its host range, as demonstrated by the susceptibility of Eragrostis curvula to finger millet blast isolates that had lost PWL1 and PWL2. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

A high-quality reference genome for the fission yeast Schizosaccharomyces osmophilus

Fission yeasts are an ancient group of fungal species that diverged from each other from tens to hundreds of million years ago. Among them is the preeminent model organism Schizosaccharomyces pombe, which has significantly contributed to our understandings of molecular mechanisms underlying fundamental cellular processes. The availability of the genomes of S. pombe and 3 other fission yeast species S. japonicus, S. octosporus, and S. cryophilus has enabled cross-species comparisons that provide insights into the evolution of genes, pathways, and genomes. Here, we performed genome sequencing on the type strain of the recently identified fission yeast species S. osmophilus and obtained a complete mitochondrial genome and a nuclear genome assembly with gaps only at rRNA gene arrays. A total of 5,098 protein-coding nuclear genes were annotated and orthologs for more than 95% of them were identified. Genome-based phylogenetic analysis showed that S. osmophilus is most closely related to S. octosporus and these 2 species diverged around 16 million years ago. To demonstrate the utility of this S. osmophilus reference genome, we conducted cross-species comparative analyses of centromeres, telomeres, transposons, the mating-type region, Cbp1 family proteins, and mitochondrial genomes. These analyses revealed conservation of repeat arrangements and sequence motifs in centromere cores, identified telomeric sequences composed of 2 types of repeats, delineated relationships among Tf1/sushi group retrotransposons, characterized the evolutionary origins and trajectories of Cbp1 family domesticated transposases, and discovered signs of interspecific transfer of 2 types of mitochondrial selfish elements.

Disentangling the complex gene interaction networks between rice and the blast fungus identifies a new pathogen effector

Studies focused solely on single organisms can fail to identify the networks underlying host-pathogen gene-for-gene interactions. Here, we integrate genetic analyses of rice (Oryza sativa, host) and rice blast fungus (Magnaporthe oryzae, pathogen) and uncover a new pathogen recognition specificity of the rice nucleotide-binding domain and leucine-rich repeat protein (NLR) immune receptor Pik, which mediates resistance to M. oryzae expressing the avirulence effector gene AVR-Pik. Rice Piks-1, encoded by an allele of Pik-1, recognizes a previously unidentified effector encoded by the M. oryzae avirulence gene AVR-Mgk1, which is found on a mini-chromosome. AVR-Mgk1 has no sequence similarity to known AVR-Pik effectors and is prone to deletion from the mini-chromosome mediated by repeated Inago2 retrotransposon sequences. AVR-Mgk1 is detected by Piks-1 and by other Pik-1 alleles known to recognize AVR-Pik effectors; recognition is mediated by AVR-Mgk1 binding to the integrated heavy metal-associated (HMA) domain of Piks-1 and other Pik-1 alleles. Our findings highlight how complex gene-for-gene interaction networks can be disentangled by applying forward genetics approaches simultaneously to the host and pathogen. We demonstrate dynamic coevolution between an NLR integrated domain and multiple families of effector proteins.

Fungal Effectoromics: A World in Constant Evolution

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

HAG Effector Evolution in Pyricularia Species and Plant Cell Death Suppression by HAG4

Seventy host-adapted gene (HAG) effector family members from Pyricularia species are found in P. oryzae and three closely related species (isolates LS and 18-2 from an unknown Pyricularia sp., P. grisea, and P. pennisetigena) that share at least eight orthologous HAG family members with P. oryzae. The genome sequence of a more distantly related species, P. penniseti, lacks HAG genes, suggesting a time frame for the origin of the gene family in the genus. In P. oryzae, HAG4 is uniquely found in the genetic lineage that contains populations adapted to Setaria and Oryza hosts. We find a nearly identical HAG4 allele in a P. grisea isolate, suggesting transfer of HAG4 from P. grisea to P. oryzae. HAG4 encodes a suppressor of plant cell death. Yeast two-hybrid screens with several HAG genes independently identify common interacting clones from a rice complementary DNA library, suggesting conservation of protein surface motifs between HAG homologs with as little as 40% protein sequence identity. HAG family orthologs have diverged rapidly and HAG15 orthologs display unusually high rates of sequence divergence compared with adjacent genes suggesting gene-specific accelerated divergence. The sequence diversity of the HAG homologs in Pyricularia species provides a resource for examining mechanisms of gene family evolution and the relationship to structural and functional evolution of HAG effector family activity. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Identification of Differentially Expressed Genes Reveal Conserved Mechanisms in the Rice-Magnaporthe oryzae Interaction

Magnaporthe oryzae causes rice blast disease and is responsible for major losses in rice production worldwide. Although numerous studies have focused on the interactions between Oryza sativa and M. oryzae, to date, the conserved mechanisms remain in part unclear. In this study, a comparative analysis of transcriptomes of O. sativa L. ssp. japonica cv. 'Nipponbare' interacting with three M. oryzae strains (248, 235, and 163) were performed to explore the conserved molecular mechanisms. Differentially expressed genes with similar expression patterns in the interactions between cultivar 'Nipponbare' and three M. oryzae strains were defined as Conserved Differentially Expressed Genes (CDEGs). These included 3,647 O. sativa CDEGs and 3,655 M. oryzae CDEGs. Four rice CDEGs (LOC_Os03g19270, LOC_Os07g36600, LOC_Os05g28740, and LOC_Os01g32780) encoding universal stress protein (USP) were induced within 24 h post-inoculation (hpi) by three M. oryzae strains. Meanwhile, overexpression of LOC_Os07g36600 resulted in enhanced rice resistance against M. oryzae. Furthermore, four rice genes coding light-harvesting chlorophyll a/b-binding (LHC) protein (LOC_Os02g52650, LOC_Os09g12540, LOC_Os11g13850, LOC_Os05g22730) were also identified as CDEGs and were induced at 48 hpi, which might contribute to blast resistance through reactive oxygen species (ROS) accumulation. MoCDIP4 is M. oryzae effector inducing rice cell death and were verified that include AA9 CAZy domain (namely GH61 domain). In this study, we found seven MoCDIP4-homologous genes coding proteins with signal peptides and AA9 CAZy domains, which were continuously up-regulated across all infection stages relative to uninoculated control. This study uncovered that genes are required for conserved mechanisms of rice-M. oryzae interaction, which includes rice genes encoding USP proteins and LHC proteins, as well as M. oryzae genes encoding AA9 proteins. This study will help us to understand how O. sativa responds to M. oryzae infections and the molecular mechanisms of M. oryzae pathogenicity.

Leaf blast, a new disease on Pennisetum sinese caused by Pyricularia pennisetigena in China

Pennisetum sinese Roxb, named king grass or 'sugarcane grass', belonging to the Poaceae family, is widely cultivated in China because of its use for livestock feed, bioenergy, water-soil conservation and phytoremediation. In September 2018, leaf blast on P. sinese was observed in a pasture in Zhanjiang, Guangdong province (21.15°N, 110.30°E ). The lesions initially were round, water-soaked, oval or spindle-shaped, then later turned to typical eye-shaped, whitish in the center with brown-black necrotic borders surrounded by a yellow halo (Fig. S1). Almost all of plants were infected and most of the spots appeared on the blades in the lower part of the plants. The leaves gradually turned yellow and withered, which affected the growth and quality of the herbage. Diseased leaves were surface-sterilized with 75% ethanol for 30s, followed by 3% hydrogen peroxide solution for 3 mins, and rinsed three times with sterile water, then placed in a hermetic container with moistened filter paper, and kept at 21-24℃ to stimulate conidiogenesis. After conidiation, conidial suspensions were prepared with sterile water and evenly spread onto PDA, and incubated at 25 ℃ in the dark to collect single-spore isolates. Colonies on PDA were fuscous black with grey centers, broad white and flat. Conidiophores were solitary, erect, straight or curved, unbranched, medium brown, smooth, 2 septate, and sometimes up to 5 septate. Conidia with a protruding and truncating hilum on basal-cells were pyriform to obclavate, 2 septate, 22.2-28.8 × 7.1-9.6 μm ( 24.9±1.7 × 8.5±0.6 μm) (Fig S2). The internal transcribed spacer region (ITS), nrRNA gene large subunit (LSU), actin gene (ACT) and RNA polymerase II largest subunit gene (RPB2) were PCR amplified, and the amplicons were sequenced by Sangon Biotech (Shanghai) Co., Ltd. (ITS: OM294657, LSU: OM294656, ACT: OM304642, RPB2: OM304643). BLASTn result showed above 99% nucleotide identity with the ITS (KM484935 465/465 bp), LSU (MH412641 862/863 bp), ACT (XM_029891619 418/421 bp) and RPB2 (XM_029892603 894/898 bp) of Pyricularia pennisetigena strains. Maximum likelihood (ML) analysis based on ITS using MEGA revealed that the fungal isolate clustered in P. pennisetigena clade with 99% bootstrap support (Fig S3). Based on morpho-molecular criteria, the fungus was identified as P. pennisetigena (Klaubauf et al., 2014; Gómez Luciano et al, 2019; Wang et al., 2019). Four healthy leaves of a 2-month-old plant were surface disinfected, and 2 leaves were wounded with a sterile needle (3 wound sites each leaf), then one leaf was inoculated with 10 μl 105 conidia/ml suspension on each site, and on the other leaf, sterile water was used as the control inoculation. The similar inoculations were conducted on the other 2 unwounded leaves. Pathogenicity tests were repeated three times. Seven days later, all inoculated leaves (wounded and unwounded) showed typical symptoms, while no symptoms were observed with the controls on both the wounded and unwounded leaves. The same fungus was recovered and confirmed by morphology and the ITS and LSU sequences. P. pennisetigena is pathogenic on a wide range of Poaceae plants in Brazil, Japan, Mali, Philippines, United States, China (Klaubauf et al., 2014; Reges et al., 2016, Wang et al., 2019). Han et al (2020) reported bacterial leaf blast on king grass while this is first report of P. pennisetigena causing leaf blast. This finding is a warning to prevent and control of this disease and its infectivity to Poaceae plants.

EffectorP 3.0: Prediction of Apoplastic and Cytoplasmic Effectors in Fungi and Oomycetes

Many fungi and oomycete species are devasting plant pathogens. These eukaryotic filamentous pathogens secrete effector proteins to facilitate plant infection. Fungi and oomycete pathogens have diverse infection strategies and their effectors generally do not share sequence homology. However, they occupy similar host environments, either the plant apoplast or plant cytoplasm, and, therefore, may share some unifying properties based on the requirements of these host compartments. Here, we exploit these biological signals and present the first classifier (EffectorP 3.0) that uses two machine-learning models: one trained on apoplastic effectors and one trained on cytoplasmic effectors. EffectorP 3.0 accurately predicts known apoplastic and cytoplasmic effectors in fungal and oomycete secretomes with low estimated false-positive rates of 3 and 8%, respectively. Cytoplasmic effectors have a higher proportion of positively charged amino acids, whereas apoplastic effectors are enriched for cysteine residues. The combination of fungal and oomycete effectors in training leads to a higher number of predicted cytoplasmic effectors in biotrophic fungi. EffectorP 3.0 expands predicted effector repertoires beyond small, cysteine-rich secreted proteins in fungi and RxLR-motif containing secreted proteins in oomycetes. We show that signal peptide prediction is essential for accurate effector prediction, because EffectorP 3.0 recognizes a cytoplasmic signal also in intracellular, nonsecreted proteins.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A single amino acid polymorphism in a conserved effector of the multihost blast fungus pathogen expands host-target binding spectrum

Accelerated gene evolution is a hallmark of pathogen adaptation and specialization following host-jumps. However, the molecular processes associated with adaptive evolution between host-specific lineages of a multihost plant pathogen remain poorly understood. In the blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae), host specialization on different grass hosts is generally associated with dynamic patterns of gain and loss of virulence effector genes that tend to define the distinct genetic lineages of this pathogen. Here, we unravelled the biochemical and structural basis of adaptive evolution of APikL2, an exceptionally conserved paralog of the well-studied rice-lineage specific effector AVR-Pik. Whereas AVR-Pik and other members of the six-gene AVR-Pik family show specific patterns of presence/absence polymorphisms between grass-specific lineages of M. oryzae, APikL2 stands out by being ubiquitously present in all blast fungus lineages from 13 different host species. Using biochemical, biophysical and structural biology methods, we show that a single aspartate to asparagine polymorphism expands the binding spectrum of APikL2 to host proteins of the heavy-metal associated (HMA) domain family. This mutation maps to one of the APikL2-HMA binding interfaces and contributes to an altered hydrogen-bonding network. By combining phylogenetic ancestral reconstruction with an analysis of the structural consequences of allelic diversification, we revealed a common mechanism of effector specialization in the AVR-Pik/APikL2 family that involves two major HMA-binding interfaces. Together, our findings provide a detailed molecular evolution and structural biology framework for diversification and adaptation of a fungal pathogen effector family following host-jumps.

Transcriptional Landscapes of Long Non-coding RNAs and Alternative Splicing in Pyricularia oryzae Revealed by RNA-Seq

Pyricularia oryzae causes the rice blast, which is one of the most devastating crop diseases worldwide, and is a model fungal pathogen widely used for dissecting the molecular mechanisms underlying fungal virulence/pathogenicity. Although the whole genome sequence of P. oryzae is publicly available, its current transcriptomes remain incomplete, lacking the information on non-protein coding genes and alternative splicing. Here, we performed and analyzed RNA-Seq of conidia and hyphae, resulting in the identification of 3,374 novel genes. Interestingly, the vast majority of these novel genes likely transcribed long non-coding RNAs (lncRNAs), and most of them were localized in the intergenic regions. Notably, their expressions were concomitant with the transcription of neighboring genes thereof in conidia and hyphae. In addition, 2,358 genes were found to undergo alternative splicing events. Furthermore, we exemplified that a lncRNA was important for hyphal growth likely by regulating the neighboring protein-coding gene and that alternative splicing of the transcription factor gene CON7 was required for appressorium formation. In summary, results from this study indicate that lncRNA transcripts and alternative splicing events are two important mechanisms for regulating the expression of genes important for conidiation, hyphal growth, and pathogenesis, and provide new insights into transcriptomes and gene regulation in the rice blast fungus.

Genome-Wide Analysis of Nutrient Signaling Pathways Conserved in Arbuscular Mycorrhizal Fungi

Arbuscular mycorrhizal (AM) fungi form a mutualistic symbiosis with a majority of terrestrial vascular plants. To achieve an efficient nutrient trade with their hosts, AM fungi sense external and internal nutrients, and integrate different hierarchic regulations to optimize nutrient acquisition and homeostasis during mycorrhization. However, the underlying molecular networks in AM fungi orchestrating the nutrient sensing and signaling remain elusive. Based on homology search, we here found that at least 72 gene components involved in four nutrient sensing and signaling pathways, including cAMP-dependent protein kinase A (cAMP-PKA), sucrose non-fermenting 1 (SNF1) protein kinase, target of rapamycin kinase (TOR) and phosphate (PHO) signaling cascades, are well conserved in AM fungi. Based on the knowledge known in model yeast and filamentous fungi, we outlined the possible gene networks functioning in AM fungi. These pathways may regulate the expression of downstream genes involved in nutrient transport, lipid metabolism, trehalase activity, stress resistance and autophagy. The RNA-seq analysis and qRT-PCR results of some core genes further indicate that these pathways may play important roles in spore germination, appressorium formation, arbuscule longevity and sporulation of AM fungi. We hope to inspire further studies on the roles of these candidate genes involved in these nutrient sensing and signaling pathways in AM fungi and AM symbiosis.

Correlation of Genomic Compartments with Contrastive Modes of Functional Losses of Host Specificity Determinants During Pathotype Differentiation in Pyricularia oryzae

The specificity between pathotypes of Pyricularia oryzae and genera of gramineous plants is governed by gene-for-gene interactions. Here, we show that avirulence genes involved in this host specificity have undergone different modes of functional losses dependent on or affected by genomic compartments harboring them. The avirulence of an Eleusine pathotype on wheat is controlled by five genes, including PWT3, which played a key role in the evolution of the Triticum pathotype (the wheat blast fungus). We cloned another gene using an association of its presence or absence with pathotypes and designated it as PWT6. PWT6 was widely distributed in a lineage composed of Eleusine and Eragrostis isolates but was completely absent in a lineage composed of Lolium and Triticum isolates. On the other hand, PWT3 homologs were present in all isolates, and their loss of function in Triticum isolates was caused by insertions of transposable elements or nucleotide substitutions. Analyses of whole-genome sequences of representative isolates revealed that these two genes were located in different genomic compartments; PWT6 was located in a repeat-rich region, while PWT3 was located in a repeat-poor region. These results suggest that the course of differentiation of the pathotypes in P. oryzae appears to be illustrated as processes of functional losses of avirulence genes but that modes of the losses are affected by genomic compartments in which they reside.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Mining oomycete proteomes for metalloproteases leads to identification of candidate virulence factors in Phytophthora infestans

Pathogens deploy a wide range of pathogenicity factors, including a plethora of proteases, to modify host tissue or manipulate host defences. Metalloproteases (MPs) have been implicated in virulence in several animal and plant pathogens. Here we investigated the repertoire of MPs in 46 stramenopile species including 37 oomycetes, 5 diatoms, and 4 brown algae. Screening their complete proteomes using hidden Markov models (HMMs) trained for MP detection resulted in over 4,000 MPs, with most species having between 65 and 100 putative MPs. Classification in clans and families according to the MEROPS database showed a highly diverse MP repertoire in each species. Analyses of domain composition, orthologous groups, distribution, and abundance within the stramenopile lineage revealed a few oomycete-specific MPs and MPs potentially related to lifestyle. In-depth analyses of MPs in the plant pathogen Phytophthora infestans revealed 91 MPs, divided over 21 protein families, including 25 MPs with a predicted signal peptide or signal anchor. Expression profiling showed different patterns of MP gene expression during pre-infection and infection stages. When expressed in leaves of Nicotiana benthamiana, 12 MPs changed the sizes of lesions caused by inoculation with P. infestans; with 9 MPs the lesions were larger, suggesting a positive effect on the virulence of P. infestans, while 3 MPs had a negative effect, resulting in smaller lesions. To the best of our knowledge, this is the first systematic inventory of MPs in oomycetes and the first study pinpointing MPs as potential pathogenicity factors in Phytophthora.

Evolution and Regulation of a Large Effector Family of Pyricularia oryzae

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

Genomic rearrangements generate hypervariable mini-chromosomes in host-specific isolates of the blast fungus

Supernumerary mini-chromosomes–a unique type of genomic structural variation–have been implicated in the emergence of virulence traits in plant pathogenic fungi. However, the mechanisms that facilitate the emergence and maintenance of mini-chromosomes across fungi remain poorly understood. In the blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae), mini-chromosomes have been first described in the early 1990s but, until very recently, have been overlooked in genomic studies. Here we investigated structural variation in four isolates of the blast fungus M. oryzae from different grass hosts and analyzed the sequences of mini-chromosomes in the rice, foxtail millet and goosegrass isolates. The mini-chromosomes of these isolates turned out to be highly diverse with distinct sequence composition. They are enriched in repetitive elements and have lower gene density than core-chromosomes. We identified several virulence-related genes in the mini-chromosome of the rice isolate, including the virulence-related polyketide synthase Ace1 and two variants of the effector gene AVR-Pik. Macrosynteny analyses around these loci revealed structural rearrangements, including inter-chromosomal translocations between core- and mini-chromosomes. Our findings provide evidence that mini-chromosomes emerge from structural rearrangements and segmental duplication of core-chromosomes and might contribute to adaptive evolution of the blast fungus.

The genomes of plant pathogens often exhibit an architecture that facilitates high rates of dynamic rearrangements and genetic diversification in virulence associated regions. These regions, which tend to be gene sparse and repeat rich, are thought to serve as a cradle for adaptive evolution. Supernumerary chromosomes, i.e. chromosomes that are only present in some but not all individuals of a species, are a special type of structural variation that have been observed in plants, animals, and fungi. Here we identified and studied supernumerary mini-chromosomes in the blast fungus Magnaporthe oryzae, a pathogen that causes some of the most destructive plant diseases. We found that rice, foxtail millet and goosegrass isolates of this pathogen contain mini-chromosomes with distinct sequence composition. All mini-chromosomes are rich in repetitive genetic elements and have lower gene densities than core-chromosomes. Further, we identified virulence-related genes on the mini-chromosome of the rice isolate. We observed large-scale genomic rearrangements around these loci, indicative of a role of mini-chromosomes in facilitating genome dynamics. Taken together, our results indicate that mini-chromosomes contribute to genome rearrangements and possibly adaptive evolution of the blast fungus.

The Impact of Blast Disease: Past, Present, and Future

Rice blast disease is both the most explosive and potentially damaging disease of the world's rice (Oryza sativa) crop and a model system for research on the molecular mechanisms that fungi use to cause plant disease. The blast fungus, Magnaporthe oryzae, is highly evolved to sense when it is on a leaf surface; to develop a pressurized cell, the appressorium, to punch through the leaf cuticle; and then to hijack living rice cells to assist it in causing disease. Host specificity, determining which plants particular fungal strains can infect, is also an important topic for research. The blast fungus is a moving target, quickly overcoming rice resistance genes we deploy to control it, and recently emerging to cause devastating disease on an entirely new cereal crop, wheat. M. oryzae is highly adaptable, with multiple examples of genetic instability at certain gene loci and in certain genomic regions. Understanding the biology of the fungus in the field, and its potential for genetic and genome variability, is key to keep it from adapting to life in the research laboratory and losing relevance to the significant impact it has on global food security.

Suppression of wheat blast resistance by an effector of Pyricularia oryzae is counteracted by a host specificity resistance gene in wheat

Wheat blast caused by the Triticum pathotype of Pyricularia oryzae poses a serious threat to wheat production in South America and Asia and is now becoming a pandemic disease. Here, we show that Rmg8, a promising wheat gene for resistance breeding, is suppressed by PWT4, an effector gene of P. oryzae, and in turn that the suppression is counteracted by Rwt4, a wheat gene recognizing PWT4. When PWT4 was introduced into a wheat blast isolate carrying AVR-Rmg8 (an avirulence gene corresponding to Rmg8), PWT4 suppressed wheat resistance conferred by Rmg8. PWT4 did not alter the expression of AVR-Rmg8, but higher expression of PWT4 led to more efficient suppression. This suppression was observed in rwt4 carriers, but not in Rwt4 carriers, indicating that it is counteracted by Rwt4. PWT4 was assumed to have been horizontally transferred from a weed-associated cryptic species, P. pennisetigena, to an Avena isolate of P. oryzae in Brazil. This implies a potential risk of the acquisition of PWT4 by the wheat blast fungus and the 'breakdown' of Rmg8. We suggest that Rmg8 should be introduced together with Rwt4 into a wheat cultivar when it is used for resistance breeding.

Genome Sequence Resource of Magnaporthe oryzae Laboratory Strain 2539

Magnaporthe oryzae causes blast disease on more than 50 species of monocot plants, including important crops such as rice, millet, and most recently wheat. Additionally, it is an important model system for studying host-pathogen interaction. Here, we report a high-quality genome assembly and annotation of a laboratory strain 2539 of M. oryzae, which is a widely used progeny of a rice-infecting isolate and a grass-infecting isolate. The genome sequence of strain 2539 will be useful for studying the evolution, host adaption, and pathogenicity of M. oryzae, which will be beneficial for a better understanding of the mechanisms of host-pathogen interaction.