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

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.

Phenotypic and Genotypic Diversity of Ascochyta fabae Populations in Southern Australia

Ascochyta fabae Speg. is a serious foliar fungal disease of faba bean and a constraint to production worldwide. This study investigated the phenotypic and genotypic diversity of the A. fabae pathogen population in southern Australia and the pathogenic variability of the population was examined on a differential set of faba bean cultivars. The host set was inoculated with 154 A. fabae isolates collected from 2015 to 2018 and a range of disease reactions from high to low aggressiveness was observed. Eighty percent of isolates collected from 2015 to 2018 were categorized as pathogenicity group (PG) PG-2 (pathogenic on Farah) and were detected in every region in each year of collection. Four percent of isolates were non-pathogenic on Farah and designated as PG-1. A small group of isolates (16%) were pathogenic on the most resistant differential cultivars, PBA Samira or Nura, and these isolates were designated PG-3. Mating types of 311 isolates collected between 1991 and 2018 were determined and showed an equal ratio of MAT1-1 and MAT1-2 in the southern Australian population. The genetic diversity and population structure of 305 isolates were examined using DArTseq genotyping, and results suggest no association of genotype with any of the population descriptors viz.: collection year, region, host cultivar, mating type, or PG. A Genome-Wide Association Study (GWAS) was performed to assess genetic association with pathogenicity traits and a significant trait-associated genomic locus for disease in Farah AR and PBA Zahra, and PG was revealed. The high frequency of mating of A. fabae indicated by the wide distribution of the two mating types means changes to virulence genes would be quickly distributed to other genotypes. Continued monitoring of the A. fabae pathogen population through pathogenicity testing will be important to identify any increases in aggressiveness or emergence of novel PGs. GWAS and future genetic studies using biparental mating populations could be useful for identifying virulence genes responsible for the observed changes in pathogenicity.

Characterization of Host-Specific Genes from Pine- and Grass-Associated Species of the Fusarium fujikuroi Species Complex

The Fusarium fujikuroi species complex (FFSC) includes socioeconomically important pathogens that cause disease for numerous crops and synthesize a variety of secondary metabolites that can contaminate feedstocks and food. Here, we used comparative genomics to elucidate processes underlying the ability of pine-associated and grass-associated FFSC species to colonize tissues of their respective plant hosts. We characterized the identity, possible functions, evolutionary origins, and chromosomal positions of the host-range-associated genes encoded by the two groups of fungi. The 72 and 47 genes identified as unique to the respective genome groups were potentially involved in diverse processes, ranging from transcription, regulation, and substrate transport through to virulence/pathogenicity. Most genes arose early during the evolution of Fusarium/FFSC and were only subsequently retained in some lineages, while some had origins outside Fusarium. Although differences in the densities of these genes were especially noticeable on the conditionally dispensable chromosome of F. temperatum (representing the grass-associates) and F. circinatum (representing the pine-associates), the host-range-associated genes tended to be located towards the subtelomeric regions of chromosomes. Taken together, these results demonstrate that multiple mechanisms drive the emergence of genes in the grass- and pine-associated FFSC taxa examined. It also highlighted the diversity of the molecular processes potentially underlying niche-specificity in these and other Fusarium species.

Comparative Genome Analyses of Plant Rust Pathogen Genomes Reveal a Confluence of Pathogenicity Factors to Quell Host Plant Defense Responses

Switchgrass rust caused by Puccinia novopanici (P. novopanici) has the ability to significantly affect the biomass yield of switchgrass, an important biofuel crop in the United States. A comparative genome analysis of P. novopanici with rust pathogen genomes infecting monocot cereal crops wheat, barley, oats, maize and sorghum revealed the presence of larger structural variations contributing to their genome sizes. A comparative alignment of the rust pathogen genomes resulted in the identification of collinear and syntenic relationships between P. novopanici and P. sorghi; P. graminis tritici 21–0 (Pgt 21) and P. graminis tritici Ug99 (Pgt Ug99) and between Pgt 21 and P. triticina (Pt). Repeat element analysis indicated a strong presence of retro elements among different Puccinia genomes, contributing to the genome size variation between ~1 and 3%. A comparative look at the enriched protein families of Puccinia spp. revealed a predominant role of restriction of telomere capping proteins (RTC), disulfide isomerases, polysaccharide deacetylases, glycoside hydrolases, superoxide dismutases and multi-copper oxidases (MCOs). All the proteomes of Puccinia spp. share in common a repertoire of 75 secretory and 24 effector proteins, including glycoside hydrolases cellobiohydrolases, peptidyl-propyl isomerases, polysaccharide deacetylases and protein disulfide-isomerases, that remain central to their pathogenicity. Comparison of the predicted effector proteins from Puccinia spp. genomes to the validated proteins from the Pathogen–Host Interactions database (PHI-base) resulted in the identification of validated effector proteins PgtSR1 (PGTG_09586) from P. graminis and Mlp124478 from Melampsora laricis across all the rust pathogen genomes.

Phytophthora: an ancient, historic, biologically and structurally cohesive and evolutionarily successful generic concept in need of preservation

The considerable economic and social impact of the oomycete genus Phytophthora is well known. In response to evidence that all downy mildews (DMs) reside phylogenetically within Phytophthora, rendering Phytophthora paraphyletic, a proposal has been made to split the genus into multiple new genera. We have reviewed the status of the genus and its relationship to the DMs. Despite a substantial increase in the number of described species and improvements in molecular phylogeny the Phytophthora clade structure has remained stable since first demonstrated in 2000. Currently some 200 species are distributed across twelve major clades in a relatively tight monophyletic cluster. In our assessment of 196 species for twenty morphological and behavioural criteria the clades show good biological cohesion. Saprotrophy, necrotrophy and hemi-biotrophy of woody and non-woody roots, stems and foliage occurs across the clades. Phylogenetically less related clades often show strong phenotypic and behavioural similarities and no one clade or group of clades shows the synapomorphies that might justify a unique generic status. We propose the clades arose from the migration and worldwide radiation ~ 140 Mya (million years ago) of an ancestral Gondwanan Phytophthora population, resulting in geographic isolation and clade divergence through drift on the diverging continents combined with adaptation to local hosts, climatic zones and habitats. The extraordinary flexibility of the genus may account for its global 'success'. The 20 genera of the obligately biotrophic, angiosperm-foliage specialised DMs evolved from Phytophthora at least twice via convergent evolution, making the DMs as a group polyphyletic and Phytophthora paraphyletic in cladistic terms. The long phylogenetic branches of the DMs indicate this occurred rather rapidly, via paraphyletic evolutionary 'jumps'. Such paraphyly is common in successful organisms. The proposal to divide Phytophthora appears more a device to address the issue of the convergent evolution of the DMs than the structure of Phytophthora per se. We consider it non-Darwinian, putting the emphasis on the emergent groups (the DMs) rather than the progenitor (Phytophthora) and ignoring the evolutionary processes that gave rise to the divergence. Further, the generic concept currently applied to the DMs is narrower than that between some closely related Phytophthora species. Considering the biological and structural cohesion of Phytophthora, its historic and social impacts and its importance in scientific communication and biosecurity protocol, we recommend that the current broad generic concept is retained by the scientific community.

Using a Hybrid Mapping Population to Identify Genomic Regions of Pyrenophora teres Associated With Virulence

Net blotches caused by Pyrenophora teres are important foliar fungal diseases of barley and result in significant yield losses of up to 40%. The two types of net blotch, net-form net blotch and spot-form net blotch, are caused by P. teres f. teres (Ptt) and P. teres f. maculata (Ptm), respectively. This study is the first to use a cross between Ptt and Ptm to identify quantitative trait loci (QTL) associated with virulence and leaf symptoms. A genetic map consisting of 1,965 Diversity Arrays Technology (DArT) markers was constructed using 351 progenies of the Ptt/Ptm cross. Eight barley cultivars showing differential reactions to the parental isolates were used to phenotype the hybrid progeny isolates. Five QTL associated with virulence and four QTL associated with leaf symptoms were identified across five linkage groups. Phenotypic variation explained by these QTL ranged from 6 to 16%. Further phenotyping of selected progeny isolates on 12 more barley cultivars revealed that three progeny isolates are moderately to highly virulent across these cultivars. The results of this study suggest that accumulation of QTL in hybrid isolates can result in enhanced virulence.

Proteinaceous Effector Discovery and Characterization in Plant Pathogenic Colletotrichum Fungi

Anthracnose caused by plant pathogenic Colletotrichum fungi results in large economic losses in field crop production worldwide. To aid the establishment of plant host infection, Colletotrichum pathogens secrete numerous effector proteins either in apoplastic space or inside of host cells for effective colonization. Understanding these effector repertoires is critical for developing new strategies for resistance breeding and disease management. With the advance of genomics and bioinformatics tools, a large repertoire of putative effectors has been identified in Colletotrichum genomes, and the biological functions and molecular mechanisms of some studied effectors have been summarized. Here, we review recent advances in genomic identification, understanding of evolutional characteristics, transcriptional profiling, and functional characterization of Colletotrichum effectors. We also offer a perspective on future research.

Short Linear Motifs (SLiMs) in “Core” RxLR Effectors of Phytophthora parasitica var. nicotianae : a Case of PpRxLR1 Effector

Oomycetes of the genus Phytophthora encompass several of the most successful plant pathogens described to date. The success of infection by Phytophthora species is attributed to the pathogens’ ability to secrete effector proteins that alter the host’s physiological processes. Structural analyses of effector proteins mainly from bacterial and viral pathogens have revealed the presence of intrinsically disordered regions that host short linear motifs (SLiMs). These motifs play important biological roles by facilitating protein-protein interactions as well as protein translocation. Nonetheless, SLiMs in Phytophthora species RxLR effectors have not been investigated previously and their roles remain unknown. Using a bioinformatics pipeline, we identified 333 candidate RxLR effectors in the strain INRA 310 of Phytophthora parasitica. Of these, 71 (21%) were also found to be present in 10 other genomes of P. parasitica, and hence, these were designated core RxLR effectors (CREs). Within the CRE sequences, the N terminus exhibited enrichment in intrinsically disordered regions compared to the C terminus, suggesting a potential role of disorder in effector translocation. Although the disorder content was reduced in the C-terminal regions, it is important to mention that most SLiMs were in this terminus. PpRxLR1 is one of the 71 CREs identified in this study, and its genes encode a 6-amino acid (aa)-long SLiM at the C terminus. We showed that PpRxLR1 interacts with several host proteins that are implicated in defense. Structural analysis of this effector using homology modeling revealed the presence of potential ligand-binding sites. Among key residues that were predicted to be crucial for ligand binding, L¹⁰² and Y¹⁰⁶ were of interest since they form part of the 6-aa-long PpRxLR1 SLiM. In silico substitution of these two residues to alanine was predicted to have a significant effect on both the function and the structure of PpRxLR1 effector. Molecular docking simulations revealed possible interactions between PpRxLR1 effector and ubiquitin-associated proteins. The ubiquitin-like SLiM carried in this effector was shown to be a potential mediator of these interactions. Further studies are required to validate and elucidate the underlying molecular mechanism of action.

IMPORTANCE The continuous gain and loss of RxLR effectors makes the control of Phytophthora spp. difficult. Therefore, in this study, we endeavored to identify RxLR effectors that are highly conserved among species, also known as “core” RxLR effectors (CREs). We reason that these highly conserved effectors target conserved proteins or processes; thus, they can be harnessed in breeding for durable resistance in plants. To further understand the mechanisms of action of CREs, structural dissection of these proteins is crucial. Intrinsically disordered regions (IDRs) that do not adopt a fixed, three-dimensional fold carry short linear motifs (SLiMs) that mediate biological functions of proteins. The presence and potential role of these SLiMs in CREs of Phytophthora spp. have been overlooked. To our knowledge, we have effectively identified CREs as well as SLiMs with the potential of promoting effector virulence. Together, this work has advanced our comprehension of Phytophthora RxLR effector function and may facilitate the development of innovative and effective control strategies.

Convergence Analysis of Rust Fungi and Anther Smuts Reveals Their Common Molecular Adaptation to a Phytoparasitic Lifestyle

Convergent evolution between distantly related taxa often mirrors adaptation to similar environments. Rust fungi and anther smuts, which belong to different classes in Pucciniomycotina, have independently evolved a phytoparasitic lifestyle, representing an example of convergent evolution in the fungal kingdom. To investigate their adaptations and the genetic bases underlying their phytoparasitic lifestyles, we performed genome-wide convergence analysis of amino acid substitutions, evolutionary rates, and gene gains and losses. Convergent substitutions were detected in ATPeV0D and RP-S27Ae, two genes important for the generation of turgor pressure and ribosomal biosynthesis, respectively. A total of 51 positively selected genes were identified, including eight genes associated with translation and three genes related to the secretion pathway. In addition, rust fungi and anther smuts contained more proteins associated with oligopeptide transporters and vacuolar proteases than did other fungi. For rust fungi and anther smuts, these forms of convergence suggest four adaptive mechanisms for a phytoparasitic lifestyle: 1) reducing the metabolic demand for hyphal growth and penetration at the pre-penetration stage, 2) maintaining the efficiency of protein synthesis during colonization, 3) ensuring the normal secretion of rapidly evolving secreted proteins, and 4) improving the capacity for oligopeptide metabolism. Our results are the first to shed light on the genetic convergence mechanisms and molecular adaptation underlying phytoparasitic lifestyles in fungi.

A chromosome-scale genome assembly of the tomato pathogen Cladosporium fulvum reveals a compartmentalized genome architecture and the presence of a dispensable chromosome

Cladosporium fulvum is a fungal pathogen that causes leaf mould of tomato. The reference genome of this pathogen was released in 2012 but its high repetitive DNA content prevented a contiguous assembly and further prohibited the analysis of its genome architecture. In this study, we combined third generation sequencing technology with the Hi-C chromatin conformation capture technique, to produce a high-quality and near complete genome assembly and gene annotation of a Race 5 isolate of C. fulvum. The resulting genome assembly contained 67.17 Mb organized into 14 chromosomes (Chr1-to-Chr14), all of which were assembled telomere-to-telomere. The smallest of the chromosomes, Chr14, is only 460 kb in size and contains 25 genes that all encode hypothetical proteins. Notably, PCR assays revealed that Chr14 was absent in 19 out of 24 isolates of a world-wide collection of C. fulvum, indicating that Chr14 is dispensable. Thus, C. fulvum is currently the second species of Capnodiales shown to harbour dispensable chromosomes. The genome of C. fulvum Race 5 is 49.7 % repetitive and contains 14 690 predicted genes with an estimated completeness of 98.9%, currently one of the highest among the Capnodiales. Genome structure analysis revealed a compartmentalized architecture composed of gene-dense and repeat-poor regions interspersed with gene-sparse and repeat-rich regions. Nearly 39.2 % of the C. fulvum Race 5 genome is affected by Repeat-Induced Point (RIP) mutations and evidence of RIP leakage toward non-repetitive regions was observed in all chromosomes, indicating the RIP plays an important role in the evolution of this pathogen. Finally, 345 genes encoding candidate effectors were identified in C. fulvum Race 5, with a significant enrichment of their location in gene-sparse regions, in accordance with the 'two-speed genome' model of evolution. Overall, the new reference genome of C. fulvum presents several notable features and is a valuable resource for studies in plant pathogens.

Comparative Genome Analysis Across 128 Phytophthora Isolates Reveal Species-Specific Microsatellite Distribution and Localized Evolution of Compartmentalized Genomes

Phytophthora sp. are invasive groups of pathogens belonging to class Oomycetes. In order to contain and control them, a deep knowledge of their biology and infection strategy is imperative. With the availability of large-scale sequencing data, it has been possible to look directly into their genetic material and understand the strategies adopted by them for becoming successful pathogens. Here, we have studied the genomes of 128 Phytophthora species available publicly with reasonable quality. Our analysis reveals that the simple sequence repeats (SSRs) of all Phytophthora sp. follow distinct isolate specific patterns. We further show that TG/CA dinucleotide repeats are far more abundant in Phytophthora sp. than other classes of repeats. In case of tri- and tetranucleotide SSRs also, TG/CA-containing motifs always dominate over others. The GC content of the SSRs are stable without much variation across the isolates of Phytophthora. Telomeric repeats of Phytophthora follow a pattern of (TTTAGGG)n or (TTAGGGT)n rather than the canonical (TTAGGG)n. RxLR (arginine-any amino acid-leucine-arginine) motifs containing effectors diverge rapidly in Phytophthora and do not show any core common group. The RxLR effectors of some Phytophthora isolates have a tendency to form clusters with RxLRs from other species than within the same species. An analysis of the flanking intergenic distance clearly indicates a two-speed genome organization for all the Phytophthora isolates. Apart from effectors and the transposons, a large number of other virulence genes such as carbohydrate-active enzymes (CAZymes), transcriptional regulators, signal transduction genes, ATP-binding cassette transporters (ABC), and ubiquitins are also present in the repeat-rich compartments. This indicates a rapid co-evolution of this powerful arsenal for successful pathogenicity. Whole genome duplication studies indicate that the pattern followed is more specific to a geographic location. To conclude, the large-scale genomic studies of Phytophthora have thrown light on their adaptive evolution, which is largely guided by the localized host-mediated selection pressure.

Evasion of plant immunity by microbial pathogens

Plant pathogenic viruses, bacteria, fungi and oomycetes cause destructive diseases in natural habitats and agricultural settings, thereby threatening plant biodiversity and global food security. The capability of plants to sense and respond to microbial infection determines the outcome of plant-microorganism interactions. Host-adapted microbial pathogens exploit various infection strategies to evade or counter plant immunity and eventually establish a replicative niche. Evasion of plant immunity through dampening host recognition or the subsequent immune signalling and defence execution is a crucial infection strategy used by different microbial pathogens to cause diseases, underpinning a substantial obstacle for efficient deployment of host genetic resistance genes for sustainable disease control. In this Review, we discuss current knowledge of the varied strategies microbial pathogens use to evade the complicated network of plant immunity for successful infection. In addition, we discuss how to exploit this knowledge to engineer crop resistance.

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

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

Comparative genome analysis of plant ascomycete fungal pathogens with different lifestyles reveals distinctive virulence strategies

BACKGROUND

Pathogens have evolved diverse lifestyles and adopted pivotal new roles in both natural ecosystems and human environments. However, the molecular mechanisms underlying their adaptation to new lifestyles are obscure. Comparative genomics was adopted to determine distinct strategies of plant ascomycete fungal pathogens with different lifestyles and to elucidate their distinctive virulence strategies.

RESULTS

We found that plant ascomycete biotrophs exhibited lower gene gain and loss events and loss of CAZyme-encoding genes involved in plant cell wall degradation and biosynthesis gene clusters for the production of secondary metabolites in the genome. Comparison with the candidate effectome detected distinctive variations between plant biotrophic pathogens and other groups (including human, necrotrophic and hemibiotrophic pathogens). The results revealed the biotroph-specific and lifestyle-conserved candidate effector families. These data have been configured in web-based genome browser applications for public display ( http://lab.malab.cn/soft/PFPG ). This resource allows researchers to profile the genome, proteome, secretome and effectome of plant fungal pathogens.

CONCLUSIONS

Our findings demonstrated different genome evolution strategies of plant fungal pathogens with different lifestyles and explored their lifestyle-conserved and specific candidate effectors. It will provide a new basis for discovering the novel effectors and their pathogenic mechanisms.

Genetic Analysis of Plant Pathogens Natural Populations

Population genetics allow to address fundamental questions about the biology of plant pathogens. By testing specific hypotheses, population genetics provide insights into the population genetic variability of pathogens across different geographical areas, time, and associated plant hosts, as well as on the structure and differentiation of populations, and on the possibility that a population is introduced and from where it has originated. In this chapter, basic concepts of population genetics are introduced, as well as the five evolutionary factors affecting populations, that is, mutations, recombination, variation in population size, gene flow, and natural selection. A step-by-step workflow, from sampling to data analysis, on how to perform a genetic analysis of natural populations of plant pathogens is discussed. Increased knowledge of the population biology of pathogens is pivotal to improve management strategies of diseases in agricultural and forest ecosystems.

An Improved Assembly of the Albugo candida Ac2V Genome Reveals the Expansion of the "CCG" Class of Effectors

Albugo candida is an obligate oomycete pathogen that infects many plants in the Brassicaceae family. We resequenced the genome of isolate Ac2V using PacBio long reads and constructed an assembly augmented by Illumina reads. The Ac2VPB genome assembly is 10% larger and more contiguous compared with a previous version. Our annotation of the new assembly, aided by RNA-sequencing information, revealed a 175% expansion (40 to 110) in the CHxC effector class, which we redefined as "CCG" based on motif analysis. This class of effectors consist of arrays of phylogenetically related paralogs residing in gene sparse regions, and shows signatures of positive selection and presence/absence polymorphism. This work provides a resource that allows the dissection of the genomic components underlying A. candida adaptation and, particularly, the role of CCG effectors in virulence and avirulence on different hosts.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Transcriptional response to host chemical cues underpins the expansion of host range in a fungal plant pathogen lineage

The host range of parasites is an important factor in assessing the dynamics of disease epidemics. The evolution of pathogens to accommodate new hosts may lead to host range expansion, a process the molecular bases of which are largely enigmatic. The fungus Sclerotinia sclerotiorum has been reported to parasitize more than 400 plant species from diverse eudicot families while its close relative, S. trifoliorum, is restricted to plants from the Fabaceae family. We analyzed S. sclerotiorum global transcriptome reprogramming on hosts from six botanical families and reveal a flexible, host-specific transcriptional program. We generated a chromosome-level genome assembly for S. trifoliorum and found near-complete gene space conservation in two representative strains of broad and narrow host range Sclerotinia species. However, S. trifoliorum showed increased sensitivity to the Brassicaceae defense compound camalexin. Comparative analyses revealed a lack of transcriptional response to camalexin in the S. trifoliorum strain and suggest that regulatory variation in detoxification and effector genes at the population level may associate with the genetic accommodation of Brassicaceae in the Sclerotinia host range. Our work proposes transcriptional plasticity and the co-existence of signatures for generalist and polyspecialist adaptive strategies in the genome of a plant pathogen.

Functional Genomics and Comparative Lineage-Specific Region Analyses Reveal Novel Insights into Race Divergence in Verticillium dahliae

Verticillium dahliae is a widespread soilborne fungus that causes Verticillium wilt on numerous economically important plant species. In tomato, until now, three races have been characterized based on the response of differential cultivars to V. dahliae, but the genetic basis of race divergence in V. dahliae remains undetermined. To investigate the genetic basis of race divergence, we sequenced the genomes of two race 2 strains and four race 3 strains for comparative analyses with two known race 1 genomes. The genetic basis of race divergence was described by the pathogenicity-related genes among the three races, orthologue analyses, and genomic structural variations. Global comparative genomics showed that chromosomal rearrangements are not the only source of race divergence and that race 3 should be split into two genotypes based on orthologue clustering. Lineage-specific regions (LSRs), frequently observed between genomes of the three races, encode several predicted secreted proteins that potentially function as suppressors of immunity triggered by known effectors. These likely contribute to the virulence of the three races. Two genes in particular that can act as markers for race 2 and race 3 (VdR2e and VdR3e, respectively) contribute to virulence on tomato, and the latter acts as an avirulence factor of race 3. We elucidated the genetic basis of race divergence through global comparative genomics and identified secreted proteins in LSRs that could potentially play critical roles in the differential virulence among the races in V. dahliae.

IMPORTANCE Deciphering the gene-for-gene relationships during host-pathogen interactions is the basis of modern plant resistance breeding. In the Verticillium dahliae-tomato pathosystem, two races (races 1 and 2) and their corresponding avirulence (Avr) genes have been identified, but strains that lack these two Avr genes exist in nature. In this system, race 3 has been described, but the corresponding Avr gene has not been identified. We de novo-sequenced genomes of six strains and identified secreted proteins within the lineage-specific regions (LSRs) distributed among the genomes of the three races that could potentially function as manipulators of host immunity. One of the LSR genes, VdR3e, was confirmed as the Avr gene for race 3. The results indicate that differences in transcriptional regulation may contribute to race differentiation. This is the first study to describe these differences and elucidate roles of secreted proteins in LSRs that play roles in race differentiation.

Comparative analysis of Phytophthora genomes data

The data presented here are related to the article entitled “Comparative analysis of Phytophthora genomes reveals oomycete pathogenesis in crops” [1]. These data contain the description of genomic structure of the two plant pathogens, P. fragariae and P. rubi and characterize several gene families associated with pathogenicity of them: P450, ACX gene families, CAZymes and effector. This data presents the relevant results of two newly sequenced P. fragariae and P. rubi, so as to provide data for further studies by researchers.

Genomic signatures of a major adaptive event in the pathogenic fungus Melampsora larici-populina

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host–pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype–phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

Identification and Comparison of Colletotrichum Secreted Effector Candidates Reveal Two Independent Lineages Pathogenic to Soybean

Colletotrichum is one of the most important plant pathogenic genus of fungi due to its scientific and economic impact. A wide range of hosts can be infected by Colletotrichum spp., which causes losses in crops of major importance worldwide, such as soybean. Soybean anthracnose is mainly caused by C. truncatum, but other species have been identified at an increasing rate during the last decade, becoming one of the most important limiting factors to soybean production in several regions. To gain a better understanding of the evolutionary origin of soybean anthracnose, we compared the repertoire of effector candidates of four Colletotrichum species pathogenic to soybean and eight species not pathogenic. Our results show that the four species infecting soybean belong to two lineages and do not share any effector candidates. These results strongly suggest that two Colletotrichum lineages have acquired the capability to infect soybean independently. This study also provides, for each lineage, a set of candidate effectors encoding genes that may have important roles in pathogenicity towards soybean offering a new resource useful for further research on soybean anthracnose management.

The genus Entomophthora: bringing the insect destroyers into the twenty-first century

The fungal genus Entomophthora consists of highly host-specific pathogens that cause deadly epizootics in their various insect hosts. The most well-known among these is the "zombie fly" fungus E. muscae, which, like other Entomophthora species, elicits a series of dramatic behaviors in infected hosts to promote optimal spore dispersal. Despite having been first described more than 160 years ago, there are still many open questions about Entomophthora biology, including the molecular underpinnings of host behavior manipulation and host specificity. This review provides a comprehensive overview of our current understanding of the biology of Entomophthora fungi and enumerates the most pressing outstanding questions that should be addressed in the field. We briefly review the discovery of Entomophthora and provide a summary of the 21 recognized Entomophthora species, including their type hosts, methods of transmission (ejection of spores after or before host death), and for which molecular data are available. Further, we argue that this genus is globally distributed, based on a compilation of Entomophthora records in the literature and in online naturalist databases, and likely to contain additional species. Evidence for strain-level specificity of hosts is summarized and directly compared to phylogenies of Entomophthora and the class Insecta. A detailed description of Entomophthora's life-cycle and observed manipulated behaviors is provided and used to summarize a consensus for ideal growth conditions. We discuss evidence for Entomophthora's adaptation to growth exclusively inside insects, such as producing wall-less hyphal bodies and a unique set of subtilisin-like proteases to penetrate the insect cuticle. However, we are only starting to understand the functions of unusual molecular and genomic characteristics, such as having large > 1 Gb genomes full of repetitive elements and potential functional diploidy. We argue that the high host-specificity and obligate life-style of most Entomophthora species provides ample scope for having been shaped by close coevolution with insects despite the current general lack of such evidence. Finally, we propose six major directions for future Entomophthora research and in doing so hope to provide a foundation for future studies of these fungi and their interaction with insects.

Chromosome-level genomes provide insights into genome evolution, organization and size in Epichloe fungi

Epichloe fungi are endophytes of cool season grasses, both wild species and commercial cultivars, where they may exhibit mutualistic or pathogenic lifestyles. The Epichloe-grass symbiosis is of great interest to agricultural research for the fungal bioprotective properties conferred to host grasses but also serves as an ideal system to study the evolution of fungal plant-pathogens in natural environments. Here, we assembled and annotated gapless chromosome-level genomes of two pathogenic Epichloe sibling species. Both genomes have a bipartite genome organization, with blocks of highly syntenic gene-rich regions separated by blocks of AT-rich DNA. The AT-rich regions show an extensive signature of RIP (repeat-induced point mutation) and the expansion of this compartment accounts for the large difference in genome size between the two species. This study reveals how the rapid evolution of repeat structure can drive divergence between closely related taxa and highlights the evolutionary role of dynamic compartments in fungal genomes.

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.

Genomic Sequence Resource of Kabatiella zeae, the Causative Pathogen of Corn Eyespot Disease

Kabatiella zeae is the causative pathogen of corn eyespot disease, which is an important leaf disease that damages corn (Zea mays L.) worldwide. In this study, we provided an annotated draft of the assembled genome of the K. zeae field strain KZ1 through PacBio and Illumina sequencing. The assembled KZ1 genome size is 23,602,820 bp, and its GC content is 50.71%. The completeness of the assembled genome is 97.6% in this study. The assembly obtained in this study has 94 contigs and the length of N₅₀ is 720,243 bp. This study is the first report of the K. zeae genome, which contributes to further research on the genetic variation and pathogenic mechanism of this important fungal pathogen.[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.

Molecular interaction of charcoal rot pathogenesis in soybean: a complex interaction

Charcoal rot (CR) is a major disease of soybean, which is caused by Macrophomina phaseolina (Mp). Increasing temperatures and low rainfall in recent years have immensely benefitted the pathogen. Hence, the search for genetically acquired resistance to this pathogen is essential. The pathogen is a hemibiotroph, which germinates on the root surface and colonizes epidermal tissue. Several surface receptors initiate pathogenesis, followed by the secretion of various enzymes that provide entry to host tissue. Several enzymes and other converging cascades in the pathogen participate against host defensive responses. β-glucan of the fungal cell wall is recognized as MAMPs (microbe-associated molecular patterns) in plants, which trigger host immune responses. Kinase receptors, resistance, and pathogenesis-related genes correspond to host defense response. They work in conjunction with hormone-mediated defense pathway especially, the systemic acquired resistance, calcium-signaling, and production of phytoalexins. Due to its quantitative nature, limited QTLs have been identified in soybean for CR resistance. The present review attempts to provide a functional link between M. phaseolina pathogenicity and soybean responses. Elucidation of CR resistance responses would facilitate improved designing of breeding programs, and may help in the selection of corresponding genes to introgress CR resistant traits.

Cross-species analysis between the maize smut fungi Ustilago maydis and Sporisorium reilianum highlights the role of transcriptional change of effector orthologs for virulence and disease

The constitution and regulation of effector repertoires shape host-microbe interactions. Ustilago maydis and Sporisorium reilianum are two closely related smut fungi, which both infect maize but cause distinct disease symptoms. Understanding how effector orthologs are regulated in these two pathogens can therefore provide insights into the evolution of different infection strategies. We tracked the infection progress of U. maydis and S. reilianum in maize leaves and used two distinct infection stages for cross-species RNA-sequencing analyses. We identified 207 of 335 one-to-one effector orthologs as differentially regulated during host colonization, which might reflect the distinct disease development strategies. Using CRISPR-Cas9-mediated gene conversion, we identified two differentially expressed effector orthologs with conserved function between two pathogens. Thus, differential expression of functionally conserved genes might contribute to species-specific adaptation and symptom development. Interestingly, another differentially expressed orthogroup (UMAG_05318/Sr10075) showed divergent protein function, providing a possible case for neofunctionalization. Collectively, we demonstrated that the diversification of effector genes in related pathogens can be caused both by alteration on the transcriptional level and through functional diversification of the encoded effector proteins.

Multiple Mechanisms Drive the Evolutionary Adaptation of Phytophthora infestans Effector Avr1 to Host Resistance

Effectors, a group of small proteins secreted by pathogens, play a central role in antagonistic interactions between plant hosts and pathogens. The evolution of effector genes threatens plant disease management and sustainable food production, but population genetic analyses to understand evolutionary mechanisms of effector genes are limited compared to molecular and functional studies. Here we investigated the evolution of the Avr1 effector gene from 111 Phytophthora infestans isolates collected from six areas covering three potato cropping regions in China using a population genetic approach. High genetic variation of the effector gene resulted from diverse mechanisms including base substitution, pre-termination, intragenic recombination and diversifying selection. Nearly 80% of the 111 sequences had a point mutation in the 512th nucleotide (T512G), which generated a pre-termination stop codon truncating 38 amino acids in the C-terminal, suggesting that the C-terminal may not be essential to ecological and biological functions of P. infestans. A significant correlation between the frequency of Avr1 sequences with the pre-termination and annual mean temperature in the collection sites suggests that thermal heterogeneity might be one of contributors to the diversifying selection, although biological and biochemical mechanisms of the likely thermal adaptation are not known currently. Our results highlight the risk of rapid adaptation of P. infestans and possibly other pathogens as well to host resistance, and the application of eco-evolutionary principles is necessary for sustainable disease management in agricultural ecosystems.

A population-level invasion by transposable elements triggers genome expansion in a fungal pathogen

Genome evolution is driven by the activity of transposable elements (TEs). The spread of TEs can have deleterious effects including the destabilization of genome integrity and expansions. However, the precise triggers of genome expansions remain poorly understood because genome size evolution is typically investigated only among deeply divergent lineages. Here, we use a large population genomics dataset of 284 individuals from populations across the globe of Zymoseptoria tritici, a major fungal wheat pathogen. We built a robust map of genome-wide TE insertions and deletions to track a total of 2456 polymorphic loci within the species. We show that purifying selection substantially depressed TE frequencies in most populations, but some rare TEs have recently risen in frequency and likely confer benefits. We found that specific TE families have undergone a substantial genome-wide expansion from the pathogen's center of origin to more recently founded populations. The most dramatic increase in TE insertions occurred between a pair of North American populations collected in the same field at an interval of 25 years. We find that both genome-wide counts of TE insertions and genome size have increased with colonization bottlenecks. Hence, the demographic history likely played a major role in shaping genome evolution within the species. We show that both the activation of specific TEs and relaxed purifying selection underpin this incipient expansion of the genome. Our study establishes a model to recapitulate TE-driven genome evolution over deeper evolutionary timescales.

The stem rust fungus Puccinia graminis f. sp. tritici induces centromeric small RNAs during late infection that are associated with genome-wide DNA methylation

Background

Silencing of transposable elements (TEs) is essential for maintaining genome stability. Plants use small RNAs (sRNAs) to direct DNA methylation to TEs (RNA-directed DNA methylation; RdDM). Similar mechanisms of epigenetic silencing in the fungal kingdom have remained elusive.

Results

We use sRNA sequencing and methylation data to gain insight into epigenetics in the dikaryotic fungus Puccinia graminis f. sp. tritici (Pgt), which causes the devastating stem rust disease on wheat. We use Hi-C data to define the Pgt centromeres and show that they are repeat-rich regions (~250 kb) that are highly diverse in sequence between haplotypes and, like in plants, are enriched for young TEs. DNA cytosine methylation is particularly active at centromeres but also associated with genome-wide control of young TE insertions. Strikingly, over 90% of Pgt sRNAs and several RNAi genes are differentially expressed during infection. Pgt induces waves of functionally diversified sRNAs during infection. The early wave sRNAs are predominantly 21 nts with a 5′ uracil derived from genes. In contrast, the late wave sRNAs are mainly 22-nt sRNAs with a 5′ adenine and are strongly induced from centromeric regions. TEs that overlap with late wave sRNAs are more likely to be methylated, both inside and outside the centromeres, and methylated TEs exhibit a silencing effect on nearby genes.

Conclusions

We conclude that rust fungi use an epigenetic silencing pathway that might have similarity with RdDM in plants. The Pgt RNAi machinery and sRNAs are under tight temporal control throughout infection and might ensure genome stability during sporulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01123-z.

Comparative Genomic Analysis Reveals Genetic Variation and Adaptive Evolution in the Pathogenicity-Related Genes of Phytophthora capsici

Phytophthora capsici is an oomycete pathogen responsible for damping off, root rot, fruit rot, and foliar blight in popular vegetable and legume crops. The existence of distinct aggressiveness levels and physiological races among the P. capsici population is a major constraint to developing resistant varieties of host crops. In the present study, we compared the genomes of three P. capsici isolates with different aggressiveness levels to reveal their genomic differences. We obtained genome sequences using short-read and long-read technologies, which yielded an average genome size of 76 Mbp comprising 514 contigs and 15,076 predicted genes. A comparative genomic analysis uncovered the signatures of accelerated evolution, gene family expansions in the pathogenicity-related genes among the three isolates. Resequencing two additional P. capsici isolates enabled the identification of average 1,023,437 SNPs, revealing the frequent accumulation of non-synonymous substitutions in pathogenicity-related gene families. Furthermore, pathogenicity-related gene families, cytoplasmic effectors and ATP binding cassette (ABC) transporters, showed expansion signals in the more aggressive isolates, with a greater number of non-synonymous SNPs. This genomic information explains the plasticity, difference in aggressiveness levels, and genome structural variation among the P. capsici isolates, providing insight into the genomic features related to the evolution and pathogenicity of this oomycete pathogen.

A devil's bargain with transposable elements in plant pathogens

Transposable elements (TEs) spread in genomes through self-copying mechanisms and are a major cause of genome expansions. Plant pathogens have finely tuned the expression of virulence factors to rely on epigenetic control targeted at nearby TEs. Stress experienced during the plant infection process leads to derepression of TEs and concurrently allows the expression of virulence factors. We argue that the derepression of TEs elements causes an evolutionary conflict by favoring TEs that can be reactivated. Active TEs and recent genome size expansions indicate that plant pathogens could face long-term consequences from the short-term benefit of fine-tuning the infection process. Hence, encoding key virulence factors close to TEs under epigenetic control constitutes a devil's bargain for pathogens.

Cryptic genetic structure and copy-number variation in the ubiquitous forest symbiotic fungus Cenococcum geophilum

Ectomycorrhizal (ECM) fungi associated with plants constitute one of the most successful symbiotic interactions in forest ecosystems. ECM support trophic exchanges with host plants and are important factors for the survival and stress resilience of trees. However, ECM clades often harbour morpho-species and cryptic lineages, with weak morphological differentiation. How this relates to intraspecific genome variability and ecological functioning is poorly known. Here, we analysed 16 European isolates of the ascomycete Cenococcum geophilum, an extremely ubiquitous forest symbiotic fungus with no known sexual or asexual spore-forming structures but with a massively enlarged genome. We carried out whole-genome sequencing to identify single-nucleotide polymorphisms. We found no geographic structure at the European scale but divergent lineages within sampling sites. Evidence for recombination was restricted to specific cryptic lineages. Lineage differentiation was supported by extensive copy-number variation. Finally, we confirmed heterothallism with a single MAT1 idiomorph per genome. Synteny analyses of the MAT1 locus revealed substantial rearrangements and a pseudogene of the opposite MAT1 idiomorph. Our study provides the first evidence for substantial genome-wide structural variation, lineage-specific recombination and low continent-wide genetic differentiation in C. geophilum. Our study provides a foundation for targeted analyses of intra-specific functional variation in this major symbiosis.

Increased abundance of secreted hydrolytic enzymes and secondary metabolite gene clusters define the genomes of latent plant pathogens in the Botryosphaeriaceae

BACKGROUND

The Botryosphaeriaceae are important plant pathogens, but also have the ability to establish asymptomatic infections that persist for extended periods in a latent state. In this study, we used comparative genome analyses to shed light on the genetic basis of the interactions of these fungi with their plant hosts. For this purpose, we characterised secreted hydrolytic enzymes, secondary metabolite biosynthetic gene clusters and general trends in genomic architecture using all available Botryosphaeriaceae genomes, and selected Dothideomycetes genomes.

RESULTS

The Botryosphaeriaceae genomes were rich in carbohydrate-active enzymes (CAZymes), proteases, lipases and secondary metabolic biosynthetic gene clusters (BGCs) compared to other Dothideomycete genomes. The genomes of Botryosphaeria, Macrophomina, Lasiodiplodia and Neofusicoccum, in particular, had gene expansions of the major constituents of the secretome, notably CAZymes involved in plant cell wall degradation. The Botryosphaeriaceae genomes were shown to have moderate to high GC contents and most had low levels of repetitive DNA. The genomes were not compartmentalized based on gene and repeat densities, but genes of secreted enzymes were slightly more abundant in gene-sparse regions.

CONCLUSION

The abundance of secreted hydrolytic enzymes and secondary metabolite BGCs in the genomes of Botryosphaeria, Macrophomina, Lasiodiplodia, and Neofusicoccum were similar to those in necrotrophic plant pathogens and some endophytes of woody plants. The results provide a foundation for comparative genomic analyses and hypotheses to explore the mechanisms underlying Botryosphaeriaceae host-plant interactions.

Population genomics of transposable element activation in the highly repressive genome of an agricultural pathogen

The activity of transposable elements (TEs) can be an important driver of genetic diversity with TE-mediated mutations having a wide range of fitness consequences. To avoid deleterious effects of TE activity, some fungi have evolved highly sophisticated genomic defences to reduce TE proliferation across the genome. Repeat-induced point mutation (RIP) is a fungal-specific TE defence mechanism efficiently targeting duplicated sequences. The rapid accumulation of RIPs is expected to deactivate TEs over the course of a few generations. The evolutionary dynamics of TEs at the population level in a species with highly repressive genome defences is poorly understood. Here, we analyse 366 whole-genome sequences of Parastagonospora nodorum, a fungal pathogen of wheat with efficient RIP. A global population genomics analysis revealed high levels of genetic diversity and signs of frequent sexual recombination. Contrary to expectations for a species with RIP, we identified recent TE activity in multiple populations. The TE composition and copy numbers showed little divergence among global populations regardless of the demographic history. Miniature inverted-repeat transposable elements (MITEs) and terminal repeat retrotransposons in miniature (TRIMs) were largely underlying recent intra-species TE expansions. We inferred RIP footprints in individual TE families and found that recently active, high-copy TEs have possibly evaded genomic defences. We find no evidence that recent positive selection acted on TE-mediated mutations rather that purifying selection maintained new TE insertions at low insertion frequencies in populations. Our findings highlight the complex evolutionary equilibria established by the joint action of TE activity, selection and genomic repression.

Acautalides A-C, Neuroprotective Diels-Alder Adducts from Solid-State Cultivated Acaulium sp. H-JQSF

The solid-state cultivation of Acaulium sp. H-JQSF isolated from Armadillidium vulgare produces acautalides A-C (1-3) as skeletally unprecedented Diels-Alder adducts of a 14-membered macrodiolide to an octadeca-9,11,13-trienoic acid. The acautalide structures, along with the intramolecular transesterifications of 1-acylglycerols, were elucidated by mass spectrometry, nuclear magnetic resonance, chemical transformation, and single-crystal X-ray diffraction. Compounds 1-3 were found to be neuroprotective with antiparkinsonic potential in the 1-methyl-4-phenylpyridinium-challenged nematode model, with the magnitude impacted by the glycerol esterification.

Dynamics of transposable elements in recently diverged fungal pathogens: lineage-specific transposable element content and efficiency of genome defenses

Transposable elements (TEs) impact genome plasticity, architecture, and evolution in fungal plant pathogens. The wide range of TE content observed in fungal genomes reflects diverse efficacy of host-genome defense mechanisms that can counter-balance TE expansion and spread. Closely related species can harbor drastically different TE repertoires. The evolution of fungal effectors, which are crucial determinants of pathogenicity, has been linked to the activity of TEs in pathogen genomes. Here, we describe how TEs have shaped genome evolution of the fungal wheat pathogen Zymoseptoria tritici and four closely related species. We compared de novo TE annotations and repeat-induced point mutation signatures in 26 genomes from the Zymoseptoria species-complex. Then, we assessed the relative insertion ages of TEs using a comparative genomics approach. Finally, we explored the impact of TE insertions on genome architecture and plasticity. The 26 genomes of Zymoseptoria species reflect different TE dynamics with a majority of recent insertions. TEs associate with accessory genome compartments, with chromosomal rearrangements, with gene presence/absence variation, and with effectors in all Zymoseptoria species. We find that the extent of RIP-like signatures varies among Z. tritici genomes compared to genomes of the sister species. The detection of a reduction of RIP-like signatures and TE recent insertions in Z. tritici reflects ongoing but still moderate TE mobility.

Reference genome assembly for Australian Ascochyta lentis isolate Al4

Ascochyta lentis causes ascochyta blight in lentil (Lens culinaris Medik.) and yield loss can be as high as 50%. With careful agronomic management practices, fungicide use, and advances in breeding resistant lentil varieties, disease severity and impact to farmers have been largely controlled. However, evidence from major lentil producing countries, Canada and Australia, suggests that A. lentis isolates can change their virulence profile and level of aggressiveness over time and under different selection pressures. In this paper, we describe the first genome assembly for A. lentis for the Australian isolate Al4, through the integration of data from Illumina and PacBio SMRT sequencing. The Al4 reference genome assembly is almost 42 Mb in size and encodes 11,638 predicted genes. The Al4 genome comprises 21 full-length and gapless chromosomal contigs and two partial chromosome contigs each with one telomere. We predicted 31 secondary metabolite clusters, and 38 putative protein effectors, many of which were classified as having an unknown function. Comparison of A. lentis genome features with the recently published reference assembly for closely related A. rabiei show that genome synteny between these species is highly conserved. However, there are several translocations and inversions of genome sequence. The location of secondary metabolite clusters near transposable element and repeat-rich genomic regions was common for A. lentis as has been reported for other fungal plant pathogens.

Biological Characteristics of Verticillium dahliae MAT1-1 and MAT1-2 Strains

Verticillium dahliae is a soil-borne plant pathogenic fungus that causes Verticillium wilt on hundreds of dicotyledonous plant species. V. dahliae is considered an asexually (clonal) reproducing fungus, although both mating type idiomorphs (MAT1-1 and MAT1-2) are present, and is heterothallic. Most of the available information on V. dahliae strains, including their biology, pathology, and genomics comes from studies on isolates with the MAT1-2 idiomorph, and thus little information is available on the MAT1-1 V. dahliae strains in the literature. We therefore evaluated the growth responses of MAT1-1 and MAT1-2 V. dahliae strains to various stimuli. Growth rates and melanin production in response to increased temperature, alkaline pH, light, and H₂O₂ stress were higher in the MAT1-2 strains than in the MAT1-1 strains. In addition, the MAT1-2 strains showed an enhanced ability to degrade complex polysaccharides, especially starch, pectin, and cellulose. Furthermore, several MAT1-2 strains from both potato and sunflower showed increased virulence on their original hosts, relative to their MAT1-1 counterparts. Thus, compared to MAT1-1 strains, MAT1-2 strains derive their potentially greater fitness from an increased capacity to adapt to their environment and exhibit higher virulence. These competitive advantages might explain the current abundance of MAT1-2 strains relative to MAT1-1 strains in the agricultural and sylvicultural ecosystems, and this study provides the baseline information on the two mating idiomorphs to study sexual reproduction in V. dahliae under natural and laboratory conditions.

Genome Sequence Resource for Nigrospora oryzae, an Important Pathogenic Fungus Threatening Crop Production

Nigrospora oryzae is an important phytopathogenic fungus with a broad host range. Here, we report an annotated draft of the genome of N. oryzae field strain GZL1 collected from maize assembled from PacBio and Illumina sequencing reads. The assembly we obtained has 15 scaffolds with an N₅₀ length of 4,037,616 bp. The resulting GZL1 draft genome is 43,214,190 bp, with GC content of 58.19%. The completeness of GZL1 genome assembly is 99.30%. This study is the first report of the genome sequence of N. oryzae, which can facilitate future study of the genetic variation and pathogenic mechanism of this important fungal pathogen.[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.

Genome reduction and relaxed selection is associated with the transition to symbiosis in the basidiomycete genus Podaxis

Insights into the genomic consequences of symbiosis for basidiomycete fungi associated with social insects remain sparse. Capitalizing on viability of spores from centuries-old herbarium specimens of free-living, facultative, and specialist termite-associated Podaxis fungi, we obtained genomes of 10 specimens, including two type species described by Linnaeus >240 years ago. We document that the transition to termite association was accompanied by significant reductions in genome size and gene content, accelerated evolution in protein-coding genes, and reduced functional capacities for oxidative stress responses and lignin degradation. Functional testing confirmed that termite specialists perform worse under oxidative stress, while all lineages retained some capacity to cleave lignin. Mitochondrial genomes of termite associates were significantly larger; possibly driven by smaller population sizes or reduced competition, supported by apparent loss of certain biosynthetic gene clusters. Our findings point to relaxed selection that mirrors genome traits observed among obligate endosymbiotic bacteria of many insects.