Genome evolution in filamentous plant pathogens: why bigger can be better

PcAvh87, a virulence essential RxLR effector of Phytophthora cinnamomi suppresses host defense and induces cell death in plant nucleus

Plants have developed intricate immune mechanisms to impede Phytophthora colonization. In response, Phytophthora secretes RxLR effector proteins that disrupt plant defense and promote infection. The specific molecular interactions through which Phytophthora RxLR effectors undermine plant immunity, however, remain inadequately defined. In this study, we delineate the role of the nuclear-localized RxLR effector PcAvh87, which is pivotal for the full virulence of Phytophthora cinnamomi. Gene expression analysis indicates that PcAvh87 expression is significantly upregulated during the initial infection stages, interacting with the immune responses triggered by the elicitin protein INF1 and pro-apoptotic protein BAX. Utilizing PEG/CaCl2-mediated protoplast transformation and CRISPR/Cas9-mediated gene editing, we generated PcAvh87 knockout mutants, which demonstrated compromised hyphal growth, sporangium development, and zoospore release, along with a marked reduction in pathogenicity. This underscores PcAvh87's crucial role as a virulence determinant. Notably, PcAvh87, conserved across the Phytophthora genus, was found to modulate the activity of plant immune protein 113, thereby attenuating plant immune responses. This implies that the PcAvh87-mediated regulatory mechanism could be a common strategy in Phytophthora species to manipulate plant immunity. Our findings highlight the multifaceted roles of PcAvh87 in promoting P. cinnamomi infection, including its involvement in sporangia production, mycelial growth, and the targeting of plant immune proteins to enhance pathogen virulence.

Melanized root-associated fungi: key players in plant-soil systems

Melanized root-associated fungi are a group of fungi that produce melanized structures and form root associations, including different mycorrhizal and endophytic symbioses with plants. They are pervasive across terrestrial ecosystems and play an important role in the prevailing soil carbon (C) and nutrient cycling syndromes through direct and indirect mechanisms, where they may strongly modulate plant-microbe interactions and structure root and soil microbiomes. Furthermore, melanized root-associated fungi can confer on plants an enhanced ability to tolerate abiotic and biotic stressors such as drought, extreme temperatures, heavy metals, and pathogen attacks. We propose that melanized root-associated fungi are a cohesive and ecologically relevant grouping that can be an indicator of plant-soil system functioning, and considering them will advance research on plant-soil interactions.

Histone (de)acetylation in epigenetic regulation of Phytophthora pathobiology

Phytophthora species are oomycetes that have evolved a broad spectrum of biological processes and improved strategies to cope with host and environmental challenges. A growing body of evidence indicates that the high pathogen plasticity is based on epigenetic regulation of gene expression linked to Phytophthora's rapid adjustment to endogenous cues and various stresses. As 5mC DNA methylation has not yet been identified in Phytophthora, the reversible processes of acetylation/deacetylation of histone proteins seem to play a pivotal role in the epigenetic control of gene expression in oomycetes. To explore this issue, we review the structure, diversity, and phylogeny of histone acetyltransferases (HATs) and histone deacetylases (HDACs) in six plant-damaging Phytophthora species: P. capsici, P. cinnamomi, P. infestans, P. parasitica, P. ramorum, and P. sojae. To further integrate and improve our understanding of the phylogenetic classification, evolutionary relationship, and functional characteristics, we supplement this review with a comprehensive view of HATs and HDACs using recent genome- and proteome-level databases. Finally, the potential functional role of transcriptional reprogramming mediated by epigenetic changes during Phytophthora species saprophytic and parasitic phases under nitro-oxidative stress is also briefly discussed.

Oomycete Metabolism Is Highly Dynamic and Reflects Lifestyle Adaptations

The selective pressure of pathogen-host symbiosis drives adaptations. How these interactions shape the metabolism of pathogens is largely unknown. Here, we use comparative genomics to systematically analyze the metabolic networks of oomycetes, a diverse group of eukaryotes that includes saprotrophs as well as animal and plant pathogens, with the latter causing devastating diseases with significant economic and/or ecological impacts. In our analyses of 44 oomycete species, we uncover considerable variation in metabolism that can be linked to lifestyle differences. Comparisons of metabolic gene content reveal that plant pathogenic oomycetes have a bipartite metabolism consisting of a conserved core and an accessory set. The accessory set can be associated with the degradation of defense compounds produced by plants when challenged by pathogens. Obligate biotrophic oomycetes have smaller metabolic networks, and taxonomically distantly related biotrophic lineages display convergent evolution by repeated gene losses in both the conserved as well as the accessory set of metabolisms. When investigating to what extent the metabolic networks in obligate biotrophs differ from those in hemibiotrophic plant pathogens, we observe that the losses of metabolic enzymes in obligate biotrophs are not random and that gene losses predominantly influence the terminal branches of the metabolic networks. Our analyses represent the first metabolism-focused comparison of oomycetes at this scale and will contribute to a better understanding of the evolution of oomycete metabolism in relation to lifestyle adaptation. Numerous oomycete species are devastating plant pathogens that cause major damage in crops and natural ecosystems. Their interactions with hosts are shaped by strong selection, but how selection affects adaptation of the primary metabolism to a pathogenic lifestyle is not yet well established. By pan-genome and metabolic network analyses of distantly related oomycete pathogens and their nonpathogenic relatives, we reveal considerable lifestyle- and lineage-specific adaptations. This study contributes to a better understanding of metabolic adaptations in pathogenic oomycetes in relation to lifestyle, host, and environment, and the findings will help in pinpointing potential targets for disease control. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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 Comparative Genomics of Botryosphaeriaceae Suggests Gene Families of Botryosphaeria dothidea Related to Pathogenicity on Chinese Hickory Tree

Trunk canker poses a major threat to the production of Chinese hickory tree (Carya cathayensis Sarg.), which is primarily determined by Botryosphaeriaceae. In our previous work, we identified Botryosphaeria dothidea as the predominant pathogen of this disease. However, it is still unclear about corresponding gene families and mechanisms associated with B. dothidea's pathogenicity on Chinese hickory tree. Here, we present a comparative analysis of high-quality genome assemblies of Botryosphaeria dothidea and other isolated pathogens, showing highly syntenic relationships between B. dothidea and its closely related species and the conservative evolution of the Botryosphaeriaceae family. Higher GC contents were found in the genomes of B. dothidea and three other isolated pathogens (Botryshaeria cortices, Botryshaeria fabicerciana, and Botryshaeria qingyuanensis) compared to Macrophomina phaseolina, Neofusicoccum parvum, Diplodia corticola, and Lasiodiplodia theobromae. An investigation of genes specific to or expanded in B. dothidea revealed that one secreted glucanase, one orsellinic acid biosynthesis enzyme, and two MFS transporters positively regulated B. dothidea's pathogenicity. We also observed an overrepresentation of viral integrase like gene and heterokaryon incompatibility proteins in the B. dothidea's genome. In addition, we observed one LRR-domain-containing protein and two Sec-domain-containing proteins (Sec_1 and Sec_7) that underwent positive selection. This study will help to understand B. dothidea's pathogenicity and potential influence on the infection of Chinese hickory, which will help in the development of disease control and ensure the security of Chinese hickory production.

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.

Segmental duplications drive the evolution of accessory regions in a major crop pathogen

Many pathogens evolved compartmentalized genomes with conserved core and variable accessory regions (ARs) that carry effector genes mediating virulence. The fungal plant pathogen Fusarium oxysporum has such ARs, often spanning entire chromosomes. The presence of specific ARs influences the host range, and horizontal transfer of ARs can modify the pathogenicity of the receiving strain. However, how these ARs evolve in strains that infect the same host remains largely unknown. We defined the pan-genome of 69 diverse F. oxysporum strains that cause Fusarium wilt of banana, a significant constraint to global banana production, and analyzed the diversity and evolution of the ARs. Accessory regions in F. oxysporum strains infecting the same banana cultivar are highly diverse, and we could not identify any shared genomic regions and in planta-induced effectors. We demonstrate that segmental duplications drive the evolution of ARs. Furthermore, we show that recent segmental duplications specifically in accessory chromosomes cause the expansion of ARs in F. oxysporum. Taken together, we conclude that extensive recent duplications drive the evolution of ARs in F. oxysporum, which contribute to the evolution of virulence.

Convenient synthesis and delivery of a megabase-scale designer accessory chromosome empower biosynthetic capacity

Synthetic biology confers new functions to hosts by introducing exogenous genetic elements, yet rebuilding complex traits that are based on large-scale genetic information remains challenging. Here, we developed a CRISPR/Cas9-mediated haploidization method that bypasses the natural process of meiosis. Based on the programmed haploidization in yeast, we further developed an easy-to-use method designated HAnDy (Haploidization-based DNA Assembly and Delivery in yeast) that enables efficient assembly and delivery of large DNA, with no need for any fussy in vitro manipulations. Using HAnDy, a de novo designed 1.024 Mb synthetic accessory chromosome (synAC) encoding 542 exogenous genes was parallelly assembled and then directly transferred to six phylogenetically diverse yeasts. The synAC significantly promotes hosts' adaptations and increases the scope of the metabolic network, which allows the emergence of valuable compounds. Our approach should facilitate the assembly and delivery of large-scale DNA for expanding and deciphering complex biological functions.

Quantitative pathogenicity and host adaptation in a fungal plant pathogen revealed by whole-genome sequencing

Knowledge of genetic determinism and evolutionary dynamics mediating host-pathogen interactions is essential to manage fungal plant diseases. Studies on the genetic architecture of fungal pathogenicity often focus on large-effect effector genes triggering strong, qualitative resistance. It is not clear how this translates to predominately quantitative interactions. Here, we use the Zymoseptoria tritici-wheat model to elucidate the genetic architecture of quantitative pathogenicity and mechanisms mediating host adaptation. With a multi-host genome-wide association study, we identify 19 high-confidence candidate genes associated with quantitative pathogenicity. Analysis of genetic diversity reveals that sequence polymorphism is the main evolutionary process mediating differences in quantitative pathogenicity, a process that is likely facilitated by genetic recombination and transposable element dynamics. Finally, we use functional approaches to confirm the role of an effector-like gene and a methyltransferase in phenotypic variation. This study highlights the complex genetic architecture of quantitative pathogenicity, extensive diversifying selection and plausible mechanisms facilitating pathogen adaptation.

Implications of the three-dimensional chromatin organization for genome evolution in a fungal plant pathogen

The spatial organization of eukaryotic genomes is linked to their biological functions, although it is not clear how this impacts the overall evolution of a genome. Here, we uncover the three-dimensional (3D) genome organization of the phytopathogen Verticillium dahliae, known to possess distinct genomic regions, designated adaptive genomic regions (AGRs), enriched in transposable elements and genes that mediate host infection. Short-range DNA interactions form clear topologically associating domains (TADs) with gene-rich boundaries that show reduced levels of gene expression and reduced genomic variation. Intriguingly, TADs are less clearly insulated in AGRs than in the core genome. At a global scale, the genome contains bipartite long-range interactions, particularly enriched for AGRs and more generally containing segmental duplications. Notably, the patterns observed for V. dahliae are also present in other Verticillium species. Thus, our analysis links 3D genome organization to evolutionary features conserved throughout the Verticillium genus.

Two-speed genomes of Epichloe fungal pathogens show contrasting signatures of selection between species and across populations

Antagonistic selection between pathogens and their hosts can drive rapid evolutionary change and leave distinct molecular footprints of past and ongoing selection in the genomes of the interacting species. Despite an increasing availability of tools able to identify signatures of selection, the genetic mechanisms underlying coevolutionary interactions and the specific genes involved are still poorly understood, especially in heterogeneous natural environments. We searched the genomes of two species of Epichloe plant pathogen for evidence of recent selection. The Epichloe genus includes highly host-specific species that can sterilize their grass hosts. We performed selection scans using genome-wide SNP data from seven natural populations of two co-occurring Epichloe sibling species specialized on different hosts. We found evidence of recent (and ongoing) selective sweeps across the genome in both species. However, selective sweeps were more abundant in the species with a larger effective population size. Sweep regions often overlapped with highly polymorphic AT-rich regions supporting the role of these genome compartments in adaptive evolution. Although most loci under selection were specific to individual populations, we could also identify several candidate genes targeted by selection in sweep regions shared among populations. The genes encoded small secreted proteins typical of fungal effectors and cell wall-degrading enzymes. By investigating the genomic signatures of selection across multiple populations and species, this study contributes to our understanding of complex adaptive processes in natural plant pathogen systems.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

A systematic screen for co-option of transposable elements across the fungal kingdom

How novel protein functions are acquired is a central question in molecular biology. Key paths to novelty include gene duplications, recombination or horizontal acquisition. Transposable elements (TEs) are increasingly recognized as a major source of novel domain-encoding sequences. However, the impact of TE coding sequences on the evolution of the proteome remains understudied. Here, we analyzed 1237 genomes spanning the phylogenetic breadth of the fungal kingdom. We scanned proteomes for evidence of co-occurrence of TE-derived domains along with other conventional protein functional domains. We detected more than 13,000 predicted proteins containing potentially TE-derived domain, of which 825 were identified in more than five genomes, indicating that many host-TE fusions may have persisted over long evolutionary time scales. We used the phylogenetic context to identify the origin and retention of individual TE-derived domains. The most common TE-derived domains are helicases derived from Academ, Kolobok or Helitron. We found putative TE co-options at a higher rate in genomes of the Saccharomycotina, providing an unexpected source of protein novelty in these generally TE depleted genomes. We investigated in detail a candidate host-TE fusion with a heterochromatic transcriptional silencing function that may play a role in TE and gene regulation in ascomycetes. The affected gene underwent multiple full or partial losses within the phylum. Overall, our work establishes a kingdom-wide view of putative host-TE fusions and facilitates systematic investigations of candidate fusion proteins.

Evolutionary dynamics of the LTR-retrotransposon crapaud in the Podospora anserina species complex and the interaction with repeat-induced point mutations

BACKGROUND

The genome of the filamentous ascomycete Podospora anserina shows a relatively high abundance of retrotransposons compared to other interspersed repeats. The LTR-retrotransposon family crapaud is particularly abundant in the genome, and consists of multiple diverged sequence variations specifically localized in the 5' half of both long terminal repeats (LTRs). P. anserina is part of a recently diverged species-complex, which makes the system ideal to classify the crapaud family based on the observed LTR variation and to study the evolutionary dynamics, such as the diversification and bursts of the elements over recent evolutionary time.

RESULTS

We developed a sequence similarity network approach to classify the crapaud repeats of seven genomes representing the P. anserina species complex into 14 subfamilies. This method does not utilize a consensus sequence, but instead it connects any copies that share enough sequence similarity over a set sequence coverage. Based on phylogenetic analyses, we found that the crapaud repeats likely diversified in the ancestor of the complex and have had activity at different time points for different subfamilies. Furthermore, while we hypothesized that the evolution into multiple subfamilies could have been a direct effect of escaping the genome defense system of repeat induced point mutations, we found this not to be the case.

CONCLUSIONS

Our study contributes to the development of methods to classify transposable elements in fungi, and also highlights the intricate patterns of retrotransposon evolution over short timescales and under high mutational load caused by nucleotide-altering genome defense.

Genome-Wide Identification and Characterization of Effector Candidates with Conserved Motif in Falciphora oryzae

Microbes employ effectors to disrupt immune responses and promote host colonization. Conserved motifs including RXLR, LFLAK-HVLVxxP (CRN), Y/F/WxC, CFEM, LysM, Chitin-bind, DPBB_1 (PNPi), and Cutinase have been discovered to play crucial roles in the functioning of effectors in filamentous fungi. Nevertheless, little is known about effectors with conserved motifs in endophytes. This research aims to discover the effector genes with conserved motifs in the genome of rice endophyte Falciphora oryzae. SignalP identified a total of 622 secreted proteins, out of which 227 were predicted as effector candidates by EffectorP. By utilizing HMM features, we discovered a total of 169 effector candidates with conserved motifs and three novel motifs. Effector candidates containing LysM, CFEM, DPBB_1, Cutinase, and Chitin_bind domains were conserved across species. In the transient expression assay, it was observed that one CFEM and one LysM activated cell death in tobacco leaves. Moreover, two CFEM and one Chitin_bind inhibited cell death induced by Bax protein. At various points during the infection, the genes' expression levels were increased. These results will help to identify functional effector proteins involving omics methods using new bioinformatics tools, thus providing a basis for the study of symbiosis mechanisms.

Pervasive associations between dark septate endophytic fungi with tree root and soil microbiomes across Europe

Trees interact with a multitude of microbes through their roots and root symbionts such as mycorrhizal fungi and root endophytes. Here, we explore the role of fungal root symbionts as predictors of the soil and root-associated microbiomes of widespread broad-leaved trees across a European latitudinal gradient. Our results suggest that, alongside factors such as climate, soil, and vegetation properties, root colonization by ectomycorrhizal, arbuscular mycorrhizal, and dark septate endophytic fungi also shapes tree-associated microbiomes. Notably, the structure of root and soil microbiomes across our sites is more strongly and consistently associated with dark septate endophyte colonization than with mycorrhizal colonization and many abiotic factors. Root colonization by dark septate endophytes also has a consistent negative association with the relative abundance and diversity of nutrient cycling genes. Our study not only indicates that root-symbiotic interactions are an important factor structuring soil communities and functions in forest ecosystems, but also that the hitherto less studied dark septate endophytes are likely to be central players in these interactions.

Target degradation specificity of phytoplasma effector phyllogen is regulated by the recruitment of host proteasome shuttle protein

Phytoplasmas infect a wide variety of plants and can cause distinctive symptoms including the conversion of floral organs into leaf-like organs, known as phyllody. Phyllody is induced by an effector protein family called phyllogens, which interact with floral MADS-box transcription factors (MTFs) responsible for determining the identity of floral organs. The MTF/phyllogen complex then interacts with the proteasomal shuttle protein RADIATION SENSITIVE23 (RAD23), which facilitates delivery of the MTF/phyllogen complex to the host proteasome for MTF degradation. Previous studies have indicated that the MTF degradation specificity of phyllogens is determined by their ability to bind to MTFs. However, in the present study, we discovered a novel mechanism determining the degradation specificity through detailed functional analyses of a phyllogen homologue of rice yellow dwarf phytoplasma (PHYLRYD ). PHYLRYD degraded a narrower range of floral MTFs than other phyllody-inducing phyllogens, resulting in compromised phyllody phenotypes in plants. Interestingly, PHYLRYD was able to bind to some floral MTFs that PHYLRYD was unable to efficiently degrade. However, the complex of PHYLRYD and the non-degradable MTF could not interact with RAD23. These results indicate that the MTF degradation specificity of PHYLRYD is correlated with the ability to form the MTF/PHYLRYD /RAD23 ternary complex, rather than the ability to bind to MTF. This study elucidated that phyllogen target specificity is regulated by both the MTF-binding ability and RAD23 recruitment ability of the MTF/phyllogen complex.

The soil-borne white root rot pathogen Rosellinia necatrix expresses antimicrobial proteins during host colonization

Rosellinia necatrix is a prevalent soil-borne plant-pathogenic fungus that is the causal agent of white root rot disease in a broad range of host plants. The limited availability of genomic resources for R. necatrix has complicated a thorough understanding of its infection biology. Here, we sequenced nine R. necatrix strains with Oxford Nanopore sequencing technology, and with DNA proximity ligation we generated a gapless assembly of one of the genomes into ten chromosomes. Whereas many filamentous pathogens display a so-called two-speed genome with more dynamic and more conserved compartments, the R. necatrix genome does not display such genome compartmentalization. It has recently been proposed that fungal plant pathogens may employ effectors with antimicrobial activity to manipulate the host microbiota to promote infection. In the predicted secretome of R. necatrix, 26 putative antimicrobial effector proteins were identified, nine of which are expressed during plant colonization. Two of the candidates were tested, both of which were found to possess selective antimicrobial activity. Intriguingly, some of the inhibited bacteria are antagonists of R. necatrix growth in vitro and can alleviate R. necatrix infection on cotton plants. Collectively, our data show that R. necatrix encodes antimicrobials that are expressed during host colonization and that may contribute to modulation of host-associated microbiota to stimulate disease development.

Novel function of the PsDMAP1 protein in regulating the growth and pathogenicity of Phytophthora sojae

The DNA methyltransferase 1-associated protein (DMAP1) was initially identified as an activator of DNA methyltransferase 1 (DNMT1), a conserved eukaryotic enzyme involved in diverse molecular processes, including histone acetylation and chromatin remodeling. However, the roles and regulatory mechanisms of DMAP1 in filamentous pathogens are still largely unknown. Here, employing bioinformatic analysis, we identified PsDMAP1 in P. sojae, which features a canonical histone tail-binding domain, as the ortholog of the human DMAP1. A phylogenetic analysis of DMAP1 protein sequences across diverse eukaryotic organisms revealed the remarkable conservation and distinctiveness of oomycete DMAP1 orthologs. Homozygous knockout of PsDMAP1 resulted in the mortality of P. sojae. Furthermore, silencing of PsDMAP1 caused a pronounced reduction in mycelial growth, production of sporangia and zoospore, cystospore germination, and virulence. PsDMAP1 also played a crucial role in the response of P. sojae to reactive oxygen species (ROS) and osmotic stresses. Moreover, PsDMAP1 interacted with DNA N6-methyladenine (6 mA) methyltransferase PsDAMT1, thereby enhancing its catalytic activity and effectively regulating 6 mA abundance in P. sojae. Our findings reveal the functional importance of PsDAMP1 in the development and infection of P. sojae, and this marks the initial exploration of the novel 6 mA regulator PsDMAP1 in plant pathogens.

The MetaInvert soil invertebrate genome resource provides insights into below-ground biodiversity and evolution

Soil invertebrates are among the least understood metazoans on Earth. Thus far, the lack of taxonomically broad and dense genomic resources has made it hard to thoroughly investigate their evolution and ecology. With MetaInvert we provide draft genome assemblies for 232 soil invertebrate species, representing 14 common groups and 94 families. We show that this data substantially extends the taxonomic scope of DNA- or RNA-based taxonomic identification. Moreover, we confirm that theories of genome evolution cannot be generalised across evolutionarily distinct invertebrate groups. The soil invertebrate genomes presented here will support the management of soil biodiversity through molecular monitoring of community composition and function, and the discovery of evolutionary adaptations to the challenges of soil conditions.

Neuroprotective α-pyrones from Nigrospora oryzae, an endophytic fungus residing in Taxus chinensis var. mairei

Endophytes coevolve with plant hosts and thus are more probable to acquire the character (in favor) of producing undescribed bioactive metabolites. Consequently, the topic has been intensely investigated for over two decades, but endophytic metabolites with neuroprotective effect remain scarce. The study presents the discovery of eight undescribed (named solanapyrones U-Z and prosolanapyrones A and B) and six known pyrones (solanapyrones A-C and E-G) from the culture of Nigrospora oryzae, an endophytic fungus associated with Taxus chinensis var. mairei. The structures and absolute configurations of undescribed pyrones were elucidated by extensive spectroscopic analysis, modified Mosher's method, and induced circular dichroism (ICD) spectrum. Solanapyrones A and B and an undescribed pyrone (solanapyrone U) were demonstrated to be more neuroprotective than clenbuterol in inducing bone marrow mesenchymal stem cells (bMSCs) to secret nerve growth factor (NGF). The work updates the pyrone chemodiversity in nature and extends the biofunction repertoire of solanapyrone-related polyketides.

Genome Sequence Resource of Cladosporium velox Strain C4 Causing Cotton Boll Disease in Xinjiang, China

Cladosporium spp., as one of the largest and most heterogeneous genera of hyphomycetes, are widely distributed worldwide. This genus is usually adaptable to a wide variety of extreme environments. However, only 11 genomes of Cladosporium genus have been publicly released. From 2017, we found for the first time that Cladosporium velox could cause cotton boll disease and lead to stiffness and cracking boll in Xinjiang, China. Herein, we provide a high-quality reference genome for the C. velox strain C4 isolated from cotton boll in Xinjiang, China. The genome size and encoding gene number of the C. velox strain C4 and C. cucumerinum strain CCNX2, which was recently released and caused the cucumber scab, showed minor differences. This resource will contribute to future research that aims to elucidate the genetic basis of C. velox pathogenicity and could expand our knowledge of Cladosporium spp. genomic characteristics that will be valuable for the development of Cladosporium disease control measures.

Intraspecific Variation of Transposable Elements Reveals Differences in the Evolutionary History of Fungal Phytopathogen Pathotypes

Transposable elements (TEs) contribute to intraspecific variation and play important roles in the evolution of fungal genomes. However, our understanding of the processes that shape TE landscapes is limited, as is our understanding of the relationship between TE content, population structure, and evolutionary history of fungal species. Fungal plant pathogens, which often have host-specific populations, are useful systems in which to study intraspecific TE content diversity. Here, we describe TE dynamics in five lineages of Magnaporthe oryzae, the fungus that causes blast disease of rice, wheat, and many other grasses. We identified differences in TE content across these lineages and showed that recent lineage-specific expansions of certain TEs have contributed to overall greater TE content in rice-infecting and Setaria-infecting lineages. We reconstructed the evolutionary histories of long terminal repeat-retrotransposon expansions and found that in some cases they were caused by complex proliferation dynamics of one element and in others by multiple elements from an older population of TEs multiplying in parallel. Additionally, we found evidence suggesting the recent transfer of a DNA transposon between rice- and wheat-infecting M. oryzae lineages and a region showing evidence of homologous recombination between those lineages, which could have facilitated such a transfer. By investigating intraspecific TE content variation, we uncovered key differences in the proliferation dynamics of TEs in various pathotypes of a fungal plant pathogen, giving us a better understanding of the evolutionary history of the pathogen itself.

Application of CRISPR/Cas9-based genome editing in ecotoxicology

Chemicals are widely used and released into the environment, and their degradation, accumulation, migration, and transformation processes in the environment can pose a threat to the ecosystem. The advancement in analytical methods with high-throughput screening of biomolecules has revolutionized the way toxicologists used to explore the effects of chemicals on organisms. CRISPR/Cas is a newly developed tool, widely used in the exploration of basic science and biologically engineered products given its high efficiency and low cost. For example, it can edit target genes efficiently, and save loss of the crop yield caused by environmental pollution as well as gain a better understanding of the toxicity mechanisms from various chemicals. This review briefly introduces the development history of CRISPR/Cas and summarizes the current application of CRISPR/Cas in ecotoxicology, including its application on improving crop yield and drug resistance towards agricultural pollution, antibiotic pollution and other threats. The benefits by applying the CRISPR/Cas9 system in conventional toxicity mechanism studies are fully demonstrated here together with its foreseeable expansions in other area of ecotoxicology. Finally, the prospects and disadvantages of CRISPR/Cas system in the field of ecotoxicology are also discussed.

Annotation of transposable elements in the transcriptome of the Neotropical brown stink bug Euschistus heros and its chromosomal distribution

Transposable elements (TEs) are DNA sequences capable of moving within the genome. Their distribution is very dynamic among organisms, and despite advances, there are still gaps in the understanding of the diversity and evolution of TEs in many insect species. In the case of Euschistus heros, considered the main stink bug in the soybean crop in Brazil, little is known about the participation of these elements. Therefore, the objective of the current work was to identify the different groups of transposable elements present in the E. heros transcriptome, evidencing their chromosomal distribution. Through RNA-Seq and de novo assembly, 60,009 transcripts were obtained, which were annotated locally via Blastn against specific databases. Of the 367 transcripts identified as TEs, 202 belong to Class II, with emphasis on the TIR order. Among Class I elements or retrotransposons, most were characterized as LINE. Phylogenetic analyses were performed with the protein domains, evidencing differences between Tc1-mariner sequences, which may be related to possible horizontal transfer events. The transposable elements that stood out in the transcriptome were selected for fluorescent in situ hybridization. DNA transposon probes hAT, Helitron, and Tc1-mariner showed mostly scattered signals, with the presence of some blocks. Retrotransposon probes Copia, Gypsy, Jockey, and RTE showed a more pulverized hybridization pattern, with the presence of small interstitial and/or terminal blocks. Studies like this one, integrating functional genomics and molecular cytogenetic tools, are essential to expanding knowledge about transcriptionally active mobile elements, and their behavior in the chromosomes.

Chromosome-scale assembly uncovers genomic compartmentation of Fusarium oxysporum f. sp. albedinis, the causal agent of Bayoud disease in date palm

Date palm (Phoenixdactylifera) is the most significant crop across North Africa and the Middle East. However, the crop faces a severe threat from Bayoud disease caused by the fungal pathogen Fusarium oxysporum f. sp. albedinis (FOA). FOA is a soil-borne fungus that infects the roots and vascular system of date palms, leading to widespread destruction of date palm plantations in North Africa over the last century. This is considered the most devastating pathogen of oasis agriculture in North Africa and responsible for loss of 13 million trees in Algeria and Morocco alone. In this study, we present a chromosome-scale high-quality genome assembly of the virulent isolate Foa 44, which provides valuable insights into understanding the genetic basis of Bayoud disease. The genome assembly consists of 11 chromosomes and 40 unplaced contigs, totalling 65,971,825 base pairs in size. It exhibits a GC ratio of 47.77% and a TE (transposable element) content of 17.30%. Through prediction and annotation, we identified 20,416 protein-coding genes. By combining gene and repeat densities analysis with alignment to Fusarium oxysporum f. sp. lycopersici (FOL) 4287 isolate genome sequence, we determined the core and lineage-specific compartments in Foa 44, shedding light on the genome structure of this pathogen. Furthermore, a phylogenomic analysis based on the 3,292 BUSCOs core genome revealed a distinct clade of FOA isolates within the Fusarium oxysporum species complex (FOSC). Notably, the genealogies of the five identified Secreted In Xylem (SIX) genes (1, 6, 9, 11 and 14) in FOA displayed a polyphyletic pattern, suggesting a horizontal inheritance of these effectors. These findings provide a valuable genomics toolbox for further research aimed at combatting the serious biotic constraints posed by FOA to date palm. This will pave the way for a deeper understanding of Bayoud disease and facilitate the development of effective diagnostic tools and control measures.

Dynamics of transposable element accumulation in the non-recombining regions of mating-type chromosomes in anther-smut fungi

In the absence of recombination, the number of transposable elements (TEs) increases due to less efficient selection, but the dynamics of such TE accumulations are not well characterized. Leveraging a dataset of 21 independent events of recombination cessation of different ages in mating-type chromosomes of Microbotryum fungi, we show that TEs rapidly accumulated in regions lacking recombination, but that TE content reached a plateau at ca. 50% of occupied base pairs by 1.5 million years following recombination suppression. The same TE superfamilies have expanded in independently evolved non-recombining regions, in particular rolling-circle replication elements (Helitrons). Long-terminal repeat (LTR) retrotransposons of the Copia and Ty3 superfamilies also expanded, through transposition bursts (distinguished from gene conversion based on LTR divergence), with both non-recombining regions and autosomes affected, suggesting that non-recombining regions constitute TE reservoirs. This study improves our knowledge of genome evolution by showing that TEs can accumulate through bursts, following non-linear decelerating dynamics.

A practical approach to genome assembly and annotation of Basidiomycota using the example of Armillaria

Technological advancements in genome sequencing, assembly and annotation platforms and algorithms that resulted in several genomic studies have created an opportunity to further our understanding of the biology of phytopathogens, including Armillaria species. Most Armillaria species are facultative necrotrophs that cause root- and stem-rot, usually on woody plants, significantly impacting agriculture and forestry worldwide. Genome sequencing, assembly and annotation in terms of samples used and methods applied in Armillaria genome projects are evaluated in this review. Infographic guidelines and a database of resources to facilitate future Armillaria genome projects were developed. Knowledge gained from genomic studies of Armillaria species is summarized and prospects for further research are provided. This guide can be applied to other diploid and dikaryotic fungal genomics.

The evolution of mini-chromosomes in the fungal genus Colletotrichum

Anthracnose diseases caused by Colletotrichum species are among the most common fungal diseases. These symptoms typically manifest as dark, sunken lesions on leaves, stems, and fruit. In China, mango anthracnose seriously affects fruit yield and quality. Genome sequencing of several species shows the presence of mini-chromosomes. These are thought to contribute to virulence, but their formation and activity remain to be fully elucidated. Here, we assembled 17 Colletotrichum genomes (16 isolated from mango plus one from persimmon) through PacBio long-read sequencing. Half of the assembled scaffolds had telomeric repeats at both ends indicating full-length chromosomes. Based on comparative genomics analysis at interspecies and intraspecies levels, we identified extensive chromosomal rearrangements events. We analyzed mini-chromosomes of Colletotrichum spp. and found large variation among close relatives. In C. fructicola, homology between core chromosomes and mini-chromosomes suggested that some mini-chromosomes were generated by recombination of core chromosomes. In C. musae GZ23-3, we found 26 horizontally transferred genes arranged in clusters on mini-chromosomes. In C. asianum FJ11-1, several potential pathogenesis-related genes on mini-chromosomes were upregulated, especially in strains with highly pathogenic phenotypes. Mutants of these upregulated genes showed obvious defects in virulence. Our findings provide insights into the evolution and potential relationships to virulence associated with mini-chromosomes. IMPORTANCE Colletotrichum is a cosmopolitan fungal genus that seriously affects fruit yield and quality of many plant species. Mini-chromosomes have been found to be related to virulence in Colletotrichum. Further examination of mini-chromosomes can help us elucidate some pathogenic mechanisms of Colletotrichum. In this study, we generated novel assemblies of several Colletotrichum strains. Comparative genomic analyses within and between Colletotrichum species were conducted. We then identified mini-chromosomes in our sequenced strains systematically. The characteristics and generation of mini-chromosomes were investigated. Transcriptome analysis and gene knockout revealed pathogenesis-related genes located on mini-chromosomes of C. asianum FJ11-1. This study represents the most comprehensive investigation of chromosome evolution and potential pathogenicity of mini-chromosomes in the Colletotrichum genus.

A chromosome-scale genome assembly of the grape powdery mildew pathogen Erysiphe necator reveals its genomic architecture and previously unknown features of its biology

Erysiphe necator is an obligate fungal pathogen that causes grape powdery mildew, globally the most important disease on grapevines. Previous attempts to obtain a quality genome assembly for this pathogen were hindered by its high repetitive DNA content. Here, chromatin conformation capture (Hi-C) with long-read PacBio sequencing was combined to obtain a chromosome-scale assembly and a high-quality annotation for E. necator isolate EnFRAME01. The resulting 81.1 Mb genome assembly is 98% complete and consists of 34 scaffolds, 11 of which represent complete chromosomes. All chromosomes contain large centromeric-like regions and lack synteny to the 11 chromosomes of the cereal PM pathogen Blumeria graminis. Further analysis of their composition showed that repeats and transposable elements (TEs) occupy 62.7% of their content. TEs were almost evenly interspersed outside centromeric and telomeric regions and massively overlapped with regions of annotated genes, suggesting that they could have a significant functional impact. Abundant gene duplicates were observed as well, particularly in genes encoding candidate secreted effector proteins. Moreover, younger in age gene duplicates exhibited more relaxed selection pressure and were more likely to be located physically close in the genome than older duplicates. A total of 122 genes with copy number variations among six isolates of E. necator were also identified and were enriched in genes that were duplicated in EnFRAME01, indicating they may reflect an adaptive variation. Taken together, our study illuminates higher-order genomic architectural features of E. necator and provides a valuable resource for studying genomic structural variations in this pathogen. IMPORTANCE Grape powdery mildew caused by the ascomycete fungus Erysiphe necator is economically the most important and recurrent disease in vineyards across the world. The obligate biotrophic nature of E. necator hinders the use of typical genetic methods to elucidate its pathogenicity and adaptation to adverse conditions, and thus comparative genomics has been a major method to study its genome biology. However, the current reference genome of E. necator isolate C-strain is highly fragmented with many non-coding regions left unassembled. This incompleteness prohibits in-depth comparative genomic analyses and the study of genomic structural variations (SVs) that are known to affect several aspects of microbial life, including fitness, virulence, and host adaptation. By obtaining a chromosome-scale genome assembly and a high-quality gene annotation for E. necator, we reveal the organization of its chromosomal content, unearth previously unknown features of its biology, and provide a reference for studying genomic SVs in this pathogen.

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

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

Epigenetic regulation of nuclear processes in fungal plant pathogens

Through the association of protein complexes to DNA, the eukaryotic nuclear genome is broadly organized into open euchromatin that is accessible for enzymes acting on DNA and condensed heterochromatin that is inaccessible. Chemical and physical alterations to chromatin may impact its organization and functionality and are therefore important regulators of nuclear processes. Studies in various fungal plant pathogens have uncovered an association between chromatin organization and expression of in planta-induced genes that are important for pathogenicity. This review discusses chromatin-based regulation mechanisms as determined in the fungal plant pathogen Verticillium dahliae and relates the importance of epigenetic transcriptional regulation and other nuclear processes more broadly in fungal plant pathogens.

Genomic and Transcriptomic Approaches Provide a Predictive Framework for Sesquiterpenes Biosynthesis in Desarmillaria tabescens CPCC 401429

Terpenoids constitute a structurally diverse class of secondary metabolites with wide applications in the pharmaceutical, fragrance and flavor industries. Desarmillaria tabescens CPCC 401429 is a basidiomycetous mushroom that could produce anti-tumor melleolides. To date, no studies have been conducted to thoroughly investigate the sesquiterpenes biosynthetic potential in Desarmillaria or related genus. This study aims to unravel the phylogeny, terpenome, and functional characterization of unique sesquiterpene biosynthetic genes of the strain CPCC 401429. Herein, we report the genome of the fungus containing 15,145 protein-encoding genes. MLST-based phylogeny and comparative genomic analyses shed light on the precise reclassification of D. tabescens suggesting that it belongs to the genus Desarmillaria. Gene ontology enrichment and pathway analyses uncover the hidden capacity for producing polyketides and terpenoids. Genome mining directed predictive framework reveals a diverse network of sesquiterpene synthases (STSs). Among twelve putative STSs encoded in the genome, six ones are belonging to the novel minor group: diverse Clade IV. In addition, RNA-sequencing based transcriptomic profiling revealed differentially expressed genes (DEGs) of the fungus CPCC 401429 in three different fermentation conditions, that of which enable us to identify noteworthy genes exemplified as STSs coding genes. Among the ten sesquiterpene biosynthetic DEGs, two genes including DtSTS9 and DtSTS10 were selected for functional characterization. Yeast cells expressing DtSTS9 and DtSTS10 could produce diverse sesquiterpene compounds, reinforced that STSs in the group Clade IV might be highly promiscuous producers. This highlights the potential of Desarmillaria in generating novel terpenoids. To summarize, our analyses will facilitate our understanding of phylogeny, STSs diversity and functional significance of Desarmillaria species. These results will encourage the scientific community for further research on uncharacterized STSs of Basidiomycota phylum, biological functions, and potential application of this vast source of secondary metabolites.

Recent transposable element bursts are associated with the proximity to genes in a fungal plant pathogen

The activity of transposable elements (TEs) contributes significantly to pathogen genome evolution. TEs often destabilize genome integrity but may also confer adaptive variation in pathogenicity or resistance traits. De-repression of epigenetically silenced TEs often initiates bursts of transposition activity that may be counteracted by purifying selection and genome defenses. However, how these forces interact to determine the expansion routes of TEs within a pathogen species remains largely unknown. Here, we analyzed a set of 19 telomere-to-telomere genomes of the fungal wheat pathogen Zymoseptoria tritici. Phylogenetic reconstruction and ancestral state estimates of individual TE families revealed that TEs have undergone distinct activation and repression periods resulting in highly uneven copy numbers between genomes of the same species. Most TEs are clustered in gene poor niches, indicating strong purifying selection against insertions near coding sequences, or as a consequence of insertion site preferences. TE families with high copy numbers have low sequence divergence and strong signatures of defense mechanisms (i.e., RIP). In contrast, small non-autonomous TEs (i.e., MITEs) are less impacted by defense mechanisms and are often located in close proximity to genes. Individual TE families have experienced multiple distinct burst events that generated many nearly identical copies. We found that a Copia element burst was initiated from recent copies inserted substantially closer to genes compared to older copies. Overall, TE bursts tended to initiate from copies in GC-rich niches that escaped inactivation by genomic defenses. Our work shows how specific genomic environments features provide triggers for TE proliferation in pathogen genomes.

A highly polymorphic effector protein promotes fungal virulence through suppression of plant-associated Actinobacteria

Plant pathogens secrete effector proteins to support host colonization through a wide range of molecular mechanisms, while plant immune systems evolved receptors to recognize effectors or their activities to mount immune responses to halt pathogens. Importantly, plants do not act as single organisms, but rather as holobionts that actively shape their microbiota as a determinant of health. The soil-borne fungal pathogen Verticillium dahliae was recently demonstrated to exploit the VdAve1 effector to manipulate the host microbiota to promote vascular wilt disease in the absence of the corresponding immune receptor Ve1. We identify a multiallelic V. dahliae gene displaying c. 65% sequence similarity to VdAve1, named VdAve1-like (VdAve1L), which shows extreme sequence variation, including alleles that encode dysfunctional proteins, indicative of selection pressure to overcome host recognition. We show that the orphan cell surface receptor Ve2, encoded at the Ve locus, does not recognize VdAve1L. Additionally, we demonstrate that the full-length variant VdAve1L2 possesses antimicrobial activity, like VdAve1, yet with a divergent activity spectrum, that is exploited by V. dahliae to mediate tomato colonization through the direct suppression of antagonistic Actinobacteria in the host microbiota. Our findings open up strategies for more targeted biocontrol against microbial plant pathogens.

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.

Engineering plant immune circuit: walking to the bright future with a novel toolbox

Plant pathogens destroy crops and cause severe yield losses, leading to an insufficient food supply to sustain the human population. Apart from relying on natural plant immune systems to combat biological agents or waiting for the appropriate evolutionary steps to occur over time, researchers are currently seeking new breakthrough methods to boost disease resistance in plants through genetic engineering. Here, we summarize the past two decades of research in disease resistance engineering against an assortment of pathogens through modifying the plant immune components (internal and external) with several biotechnological techniques. We also discuss potential strategies and provide perspectives on engineering plant immune systems for enhanced pathogen resistance and plant fitness.

Chromosome-level analysis of the Colletotrichum graminicola genome reveals the unique characteristics of core and minichromosomes

The fungal pathogen Colletotrichum graminicola causes the anthracnose of maize (Zea mays) and is responsible for significant yield losses worldwide. The genome of C. graminicola was sequenced in 2012 using Sanger sequencing, 454 pyrosequencing, and an optical map to obtain an assembly of 13 pseudochromosomes. We re-sequenced the genome using a combination of short-read (Illumina) and long-read (PacBio) technologies to obtain a chromosome-level assembly. The new version of the genome sequence has 13 chromosomes with a total length of 57.43 Mb. We detected 66 (23.62 Mb) structural rearrangements in the new assembly with respect to the previous version, consisting of 61 (21.98 Mb) translocations, 1 (1.41 Mb) inversion, and 4 (221 Kb) duplications. We annotated the genome and obtained 15,118 predicted genes and 3,614 new gene models compared to the previous version of the assembly. We show that 25.88% of the new assembly is composed of repetitive DNA elements (13.68% more than the previous assembly version), which are mostly found in gene-sparse regions. We describe genomic compartmentalization consisting of repeat-rich and gene-poor regions vs. repeat-poor and gene-rich regions. A total of 1,140 secreted proteins were found mainly in repeat-rich regions. We also found that ~75% of the three smallest chromosomes (minichromosomes, between 730 and 551 Kb) are strongly affected by repeat-induced point mutation (RIP) compared with 28% of the larger chromosomes. The gene content of the minichromosomes (MCs) comprises 121 genes, of which 83.6% are hypothetical proteins with no predicted function, while the mean percentage of Chr1-Chr10 is 36.5%. No predicted secreted proteins are present in the MCs. Interestingly, only 2% of the genes in Chr11 have homologs in other strains of C. graminicola, while Chr12 and 13 have 58 and 57%, respectively, raising the question as to whether Chrs12 and 13 are dispensable. The core chromosomes (Chr1-Chr10) are very different with respect to the MCs (Chr11-Chr13) in terms of the content and sequence features. We hypothesize that the higher density of repetitive elements and RIPs in the MCs may be linked to the adaptation and/or host co-evolution of this pathogenic fungus.

A Bioinformatic Guide to Identify Protein Effectors from Phytopathogens

Phytopathogenic fungi are a diverse and widespread group that has a significant detrimental impact on crops with an estimated annual average loss of 15% worldwide. Understanding the interaction between host plants and pathogenic fungi is critical to delineate underlying mechanisms of plant defense to mitigate agricultural losses. Fungal pathogens utilize suites of secreted molecules, called effectors, to modulate plant metabolism and immune response to overcome host defenses and promote colonization. Effectors come in many flavors including proteinaceous products, small RNAs, and metabolites such as mycotoxins. This review will focus on methods for identifying protein effectors from fungi. Excellent reviews have been published to identify secondary metabolites and small RNAs from fungi and therefore will not be part of this review.

Pathogenicity Variation in Two Genomes of Cercospora Species Causing Gray Leaf Spot in Maize

The gray leaf spots caused by Cercospora spp. severely affect the yield and quality of maize. However, the evolutionary relation and pathogenicity variation between species of the Cercospora genus is largely unknown. In this study, we constructed high-quality reference genomes by nanopore sequencing two Cercospora species, namely, C. zeae-maydis and C. zeina, with differing pathogenicity, collected from northeast (Liaoning [LN]) and southeast (Yunnan [YN]) China, respectively. The genome size of C. zeae-maydis-LN is 45.08 Mb, containing 10,839 annotated genes, whereas that of Cercospora zeina-YN is 42.18 Mb, containing 10,867 annotated genes, of which approximately 86.58% are common in the two species. The difference in their genome size is largely attributed to increased long terminal repeat retrotransposons of 3.8 Mb in total length in C. zeae-maydis-LN. There are 41 and 30 carbohydrate-binding gene subfamilies identified in C. zeae-maydis-LN and C. zeina-YN, respectively. A higher number of carbohydrate-binding families found in C. zeae-maydis-LN, and its unique CBM4, CBM37, and CBM66, in particular, may contribute to variation in pathogenicity between the two species, as the carbohydrate-binding genes are known to encode cell wall-degrading enzymes. Moreover, there are 114 and 107 effectors predicted, with 47 and 46 having unique potential pathogenicity in C. zeae-maydis-LN and C. zeina-YN, respectively. Of eight effectors randomly selected for pathogenic testing, five were found to inhibit cell apoptosis induced by Bcl-2-associated X. Taken together, our results provide genomic insights into variation in pathogenicity between C. zeae-maydis and C. zeina. [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.

Comparative Genomic Analysis of 31 Phytophthora Genomes Reveals Genome Plasticity and Horizontal Gene Transfer

Phytophthora species are oomycete plant pathogens that cause great economic and ecological impacts. The Phytophthora genus includes over 180 known species, infecting a wide range of plant hosts, including crops, trees, and ornamentals. We sequenced the genomes of 31 individual Phytophthora species and 24 individual transcriptomes to study genetic relationships across the genus. De novo genome assemblies revealed variation in genome sizes, numbers of predicted genes, and in repetitive element content across the Phytophthora genus. A genus-wide comparison evaluated orthologous groups of genes. Predicted effector gene counts varied across Phytophthora species by effector family, genome size, and plant host range. Predicted numbers of apoplastic effectors increased as the host range of Phytophthora species increased. Predicted numbers of cytoplasmic effectors also increased with host range but leveled off or decreased in Phytophthora species that have enormous host ranges. With extensive sequencing across the Phytophthora genus, we now have the genomic resources to evaluate horizontal gene transfer events across the oomycetes. Using a machine-learning approach to identify horizontally transferred genes with bacterial or fungal origin, we identified 44 candidates over 36 Phytophthora species genomes. Phylogenetic reconstruction indicates that the transfers of most of these 44 candidates happened in parallel to major advances in the evolution of the oomycetes and Phytophthora spp. We conclude that the 31 genomes presented here are essential for investigating genus-wide genomic associations in genus Phytophthora. [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.

Folding Features and Dynamics of 3D Genome Architecture in Plant Fungal Pathogens

The folding and dynamics of three-dimensional (3D) genome organization are fundamental for eukaryotes executing genome functions but have been largely unexplored in nonmodel fungi. Using high-throughput sequencing coupled with chromosome conformation capture (Hi-C) data, we generated two chromosome-level assemblies for Puccinia striiformis f. sp. tritici, a fungus causing stripe rust disease on wheat, for studying 3D genome architectures of plant pathogenic fungi. The chromatin organization of the fungus followed a combination of the fractal globule model and the equilibrium globule model. Surprisingly, chromosome compartmentalization was not detected. Dynamics of 3D genome organization during two developmental stages of P. striiformis f. sp. tritici indicated that regulation of gene activities might be independent of the changes of genome organization. In addition, chromatin conformation conservation was found to be independent of genome sequence synteny conservation among different fungi. These results highlighted the distinct folding principles of fungal 3D genomes. Our findings should be an important step toward a holistic understanding of the principles and functions of genome architecture across different eukaryotic kingdoms. IMPORTANCE Previously, our understanding of 3D genome architecture has mainly come from model mammals, insects, and plants. However, the organization and regulatory functions of 3D genomes in fungi are largely unknown. In this study, we comprehensively investigated P. striiformis f. sp. tritici, a plant fungal pathogen, and revealed distinct features of the 3D genome, comparing it with the universal folding feature of 3D genomes in higher eukaryotic organisms. We further suggested that there might be distinct regulatory mechanisms of gene expression that are independent of chromatin organization changes during the developmental stages of this rust fungus. Moreover, we showed that the evolutionary pattern of 3D genomes in this fungus is also different from the cases in mammalian genomes. In addition, the genome assembly pipeline and the generated two chromosome-level genomes will be valuable resources. These results highlighted the unexplored distinct features of 3D genome organization in fungi. Therefore, our study provided complementary knowledge to holistically understand the organization and functions of 3D genomes across different eukaryotes.

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.

Genetic co-option into plant-filamentous pathogen interactions

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

High Nucleotide Substitution Rates Associated with Retrotransposon Proliferation Drive Dynamic Secretome Evolution in Smut Pathogens

Transposable elements (TEs) play a pivotal role in shaping diversity in eukaryotic genomes. The covered smut pathogen on barley, Ustilago hordei, encountered a recent genome expansion. Using long reads, we assembled genomes of 6 U. hordei strains and 3 sister species, to study this genome expansion. We found that larger genome sizes can mainly be attributed to a higher genome fraction of long terminal repeat retrotransposons (LTR-RTs). In the studied smut genomes, LTR-RTs fractions are the largest in U. hordei and are positively correlated with the mating-type locus sizes, which is up to ~560 kb in U. hordei. Furthermore, LTR-RTs were found to be associated with higher nucleotide substitution levels, as these occur in specific genome regions of smut species with a recent LTR-RT proliferation. Moreover, genes in genome regions with higher nucleotide substitution levels generally reside closer to LTR-RTs than other genome regions. Genome regions with many nucleotide substitutions encountered an especially high fraction of CG substitutions, which is not observed for LTR-RT sequences. The high nucleotide substitution levels particularly accelerate the evolution of secretome genes, as their more accessory nature results in substitutions that often lead to amino acid alterations. IMPORTANCE Genomic alteration can be generated through various means, in which transposable elements (TEs) can play a pivotal role. Their mobility causes mutagenesis in itself and can disrupt the function of the sequences they insert into. They also impact genome evolution as their repetitive nature facilitates nonhomologous recombination. Furthermore, TEs have been linked to specific epigenetic genome organizations. We report a recent TE proliferation in the genome of the barley covered smut fungus, Ustilago hordei. This proliferation is associated with a distinct nucleotide substitution regime that has a higher rate and a higher fraction of CG substitutions. This different regime shapes the evolution of genes in subjected genome regions. We hypothesize that TEs may influence the error-rate of DNA polymerase in a hitherto unknown fashion.

Population genomics of Puccinia graminis f.sp. tritici highlights the role of admixture in the origin of virulent wheat rust races

Puccinia graminis f.sp. tritici (Pgt) causes stem rust disease in wheat that can result in severe yield losses. The factors driving the evolution of its virulence and adaptation remain poorly characterized. We utilize long-read sequencing to develop a haplotype-resolved genome assembly of a U.S. isolate of Pgt. Using Pgt haplotypes as a reference, we characterize the structural variants (SVs) and single nucleotide polymorphisms in a diverse panel of isolates. SVs impact the repertoire of predicted effectors, secreted proteins involved in host-pathogen interaction, and show evidence of purifying selection. By analyzing global and local genomic ancestry we demonstrate that the origin of 8 out of 12 Pgt clades is linked with either somatic hybridization or sexual recombination between the diverged donor populations. Our study shows that SVs and admixture events appear to play an important role in broadening Pgt virulence and the origin of highly virulent races, creating a resource for studying the evolution of Pgt virulence and preventing future epidemic outbreaks.

Genetic variation along an altitudinal gradient in the Phytophthora infestans effector gene Pi02860

Effector genes, together with climatic and other environmental factors, play multifaceted roles in the development of plant diseases. Understanding the role of environmental factors, particularly climate conditions affecting the evolution of effector genes, is important for predicting the long-term value of the genes in controlling agricultural diseases. Here, we collected Phytophthora infestans populations from five locations along a mountainous hill in China and sequenced the effector gene Pi02860 from >300 isolates. To minimize the influence of other ecological factors, isolates were sampled from the same potato cultivar on the same day. We also expressed the gene to visualise its cellular location, assayed its pathogenicity and evaluated its response to experimental temperatures. We found that Pi02860 exhibited moderate genetic variation at the nucleotide level which was mainly generated by point mutation. The mutations did not change the cellular location of the effector gene but significantly modified the fitness of P. infestans. Genetic variation and pathogenicity of the effector gene were positively associated with the altitude of sample sites, possibly due to increased mutation rate induced by the vertical distribution of environmental factors such as UV radiation and temperature. We further found that Pi02860 expression was regulated by experimental temperature with reduced expression as experimental temperature increased. Together, these results indicate that UV radiation and temperature are important environmental factors regulating the evolution of effector genes and provide us with considerable insight as to their future sustainable action under climate and other environmental change.

Current Status and Future Perspectives of Genomics Research in the Rust Fungi

Rust fungi in Pucciniales have caused destructive plant epidemics, have become more aggressive with new virulence, rapidly adapt to new environments, and continually threaten global agriculture. With the rapid advancement of genome sequencing technologies and data analysis tools, genomics research on many of the devastating rust fungi has generated unprecedented insights into various aspects of rust biology. In this review, we first present a summary of the main findings in the genomics of rust fungi related to variations in genome size and gene composition between and within species. Then we show how the genomics of rust fungi has promoted our understanding of the pathogen virulence and population dynamics. Even with great progress, many questions still need to be answered. Therefore, we introduce important perspectives with emphasis on the genome evolution and host adaptation of rust fungi. We believe that the comparative genomics and population genomics of rust fungi will provide a further understanding of the rapid evolution of virulence and will contribute to monitoring the population dynamics for disease management.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga’s genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host–pathogen interactions.