Emergence of a hybrid PKS-NRPS secondary metabolite cluster in a clonal population of the rice blast fungus Magnaporthe oryzae

Discovery of fungal onoceroid triterpenoids through domainless enzyme-targeted global genome mining

Genomics-guided methodologies have revolutionized the discovery of natural products. However, a major challenge in the field of genome mining is determining how to selectively extract biosynthetic gene clusters (BGCs) for untapped natural products from numerous available genome sequences. In this study, we developed a fungal genome mining tool that extracts BGCs encoding enzymes that lack a detectable protein domain (i.e., domainless enzymes) and are not recognized as biosynthetic proteins by existing bioinformatic tools. We searched for BGCs encoding a homologue of Pyr4-family terpene cyclases, which are representative examples of apparently domainless enzymes, in approximately 2000 fungal genomes and discovered several BGCs with unique features. The subsequent characterization of selected BGCs led to the discovery of fungal onoceroid triterpenoids and unprecedented onoceroid synthases. Furthermore, in addition to the onoceroids, a previously unreported sesquiterpene hydroquinone, of which the biosynthesis involves a Pyr4-family terpene cyclase, was obtained. Our genome mining tool has broad applicability in fungal genome mining and can serve as a beneficial platform for accessing diverse, unexploited natural products.

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.

Genomic Comparative Analysis of Cordyceps pseudotenuipes with Other Species from Cordyceps

The whole genome of Cordyceps pseudotenuipes was sequenced, annotated, and compared with three related species to characterize the genome. The antibiotics and Secondary Metabolites Analysis Shell (antiSMASH) and local BLAST analysis were used to explore the secondary metabolites (SMs) and biosynthesis gene clusters (BGCs) of the genus Cordyceps. The genome-wide basic characteristics of C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris revealed unequal genome size, with C. cicadae as the largest (34.11 Mb), followed by C. militaris (32.27 Mb). However, the total gene lengths of C. pseudotenuipes and C. tenuipes were similar (30.1 Mb and 30.06 Mb). The GC contents of C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris genomes differed slightly (51.40% to 54.11%). AntiSMASH and local BLAST analysis showed that C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris had 31, 28, 31, and 29 putative SM BGCs, respectively. The SM BGCs contained different quantities of polyketide synthetase (PKS), nonribosomal peptide synthetase (NRPS), terpene, hybrid PKS + NRPS, and hybrid NRPS + Other. Moreover, C. pseudotenuipes, C. tenuipes, C. cicadae, and C. militaris had BGCs for the synthesis of dimethylcoprogen. C. pseudotenuipes, C. tenuipes, and C. cicadae had BGCs for the synthesis of leucinostatin A/B, neosartorin, dimethylcoprogen, wortmanamide A/B, and beauvericin. In addition, the SM BGCs unique to C. pseudotenuipes were clavaric acid, communesin, and deoxynivalenol. Synteny analysis indicated that the scaffolds where the SM BGC was located were divided into more than 70 collinear blocks, and there might be rearrangements. Altogether, these findings improved our understanding of the molecular biology of the genus Cordyceps and will facilitate the discovery of new biologically active SMs from the genus Cordyceps using heterologous expression and gene knockdown methods.

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

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

Genetic Variation Bias toward Noncoding Regions and Secreted Proteins in the Rice Blast Fungus Magnaporthe oryzae

The genomes of plant pathogens are highly variable and plastic. Pathogen gene repertoires change quickly with the plant environment, which results in a rapid loss of plant resistance shortly after the pathogen emerges in the field. Extensive studies have evaluated natural pathogen populations to understand their evolutionary effects; however, the number of studies that have examined the dynamic processes of the mutation and adaptation of plant pathogens to host plants remains limited. Here, we applied experimental evolution and high-throughput pool sequencing to Magnaporthe oryzae, a fungal pathogen that causes massive losses in rice production, to observe the evolution of genome variation. We found that mutations, including single-nucleotide variants (SNVs), insertions and deletions (indels), and transposable element (TE) insertions, accumulated very rapidly throughout the genome of M. oryzae during sequential plant inoculation and preferentially in noncoding regions, while such mutations were not frequently found in coding regions. However, we also observed that new TE insertions accumulated with time and preferentially accumulated at the proximal region of secreted protein (SP) coding genes in M. oryzae populations. Taken together, these results revealed a bias in genetic variation toward noncoding regions and SP genes in M. oryzae and may contribute to the rapid adaptive evolution of the blast fungal effectors under host selection.IMPORTANCE Plants "lose" resistance toward pathogens shortly after their widespread emergence in the field because plant pathogens mutate and adapt rapidly under resistance selection. Thus, the rapid evolution of pathogens is a serious threat to plant health. Extensive studies have evaluated natural pathogen populations to understand their evolutionary effects; however, the study of the dynamic processes of the mutation and adaptation of plant pathogens to host plants remains limited. Here, by performing an experimental evolution study, we found a bias in genetic variation toward noncoding regions and SPs in the rice blast fungus Magnaporthe oryzae, which explains the ability of the rice blast fungus to maintain high virulence variation to overcome rice resistance in the field.