The two-speed genomes of filamentous pathogens: waltz with plants

Regressive evolution of an effector following a host jump in the Irish potato famine pathogen lineage

In order to infect a new host species, the pathogen must evolve to enhance infection and transmission in the novel environment. Although we often think of evolution as a process of accumulation, it is also a process of loss. Here, we document an example of regressive evolution of an effector activity in the Irish potato famine pathogen (Phytophthora infestans) lineage, providing evidence that a key sequence motif in the effector PexRD54 has degenerated following a host jump. We began by looking at PexRD54 and PexRD54-like sequences from across Phytophthora species. We found that PexRD54 emerged in the common ancestor of Phytophthora clade 1b and 1c species, and further sequence analysis showed that a key functional motif, the C-terminal ATG8-interacting motif (AIM), was also acquired at this point in the lineage. A closer analysis showed that the P. mirabilis PexRD54 (PmPexRD54) AIM is atypical, the otherwise-conserved central residue mutated from a glutamate to a lysine. We aimed to determine whether this PmPexRD54 AIM polymorphism represented an adaptation to the Mirabilis jalapa host environment. We began by characterizing the M. jalapa ATG8 family, finding that they have a unique evolutionary history compared to previously characterized ATG8s. Then, using co-immunoprecipitation and isothermal titration calorimetry assays, we showed that both full-length PmPexRD54 and the PmPexRD54 AIM peptide bind weakly to the M. jalapa ATG8s. Through a combination of binding assays and structural modelling, we showed that the identity of the residue at the position of the PmPexRD54 AIM polymorphism can underpin high-affinity binding to plant ATG8s. Finally, we conclude that the functionality of the PexRD54 AIM was lost in the P. mirabilis lineage, perhaps owing to as-yet-unknown selection pressure on this effector in the new host environment.

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.

The pangenome of the wheat pathogen Pyrenophora tritici-repentis reveals novel transposons associated with necrotrophic effectors ToxA and ToxB

BACKGROUND

In fungal plant pathogens, genome rearrangements followed by selection pressure for adaptive traits have facilitated the co-evolutionary arms race between hosts and their pathogens. Pyrenophora tritici-repentis (Ptr) has emerged recently as a foliar pathogen of wheat worldwide and its populations consist of isolates that vary in their ability to produce combinations of different necrotrophic effectors. These effectors play vital roles in disease development. Here, we sequenced the genomes of a global collection (40 isolates) of Ptr to gain insights into its gene content and genome rearrangements.

RESULTS

A comparative genome analysis revealed an open pangenome, with an abundance of accessory genes (~ 57%) reflecting Ptr's adaptability. A clear distinction between pathogenic and non-pathogenic genomes was observed in size, gene content, and phylogenetic relatedness. Chromosomal rearrangements and structural organization, specifically around effector coding genes, were detailed using long-read assemblies (PacBio RS II) generated in this work in addition to previously assembled genomes. We also discovered the involvement of large mobile elements associated with Ptr's effectors: ToxA, the gene encoding for the necrosis effector, was found as a single copy within a 143-kb 'Starship' transposon (dubbed 'Horizon') with a clearly defined target site and target site duplications. 'Horizon' was located on different chromosomes in different isolates, indicating mobility, and the previously described ToxhAT transposon (responsible for horizontal transfer of ToxA) was nested within this newly identified Starship. Additionally, ToxB, the gene encoding the chlorosis effector, was clustered as three copies on a 294-kb element, which is likely a different putative 'Starship' (dubbed 'Icarus') in a ToxB-producing isolate. ToxB and its putative transposon were missing from the ToxB non-coding reference isolate, but the homolog toxb and 'Icarus' were both present in a different non-coding isolate. This suggests that ToxB may have been mobile at some point during the evolution of the Ptr genome which is contradictory to the current assumption of ToxB vertical inheritance. Finally, the genome architecture of Ptr was defined as 'one-compartment' based on calculated gene distances and evolutionary rates.

CONCLUSIONS

These findings together reflect on the highly plastic nature of the Ptr genome which has likely helped to drive its worldwide adaptation and has illuminated the involvement of giant transposons in facilitating the evolution of virulence in Ptr.

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.

Comparative genomic and transcriptomic analyses of trans-kingdom pathogen Fusarium solani species complex reveal degrees of compartmentalization

BACKGROUND

The Fusarium solani species complex (FSSC) comprises fungal pathogens responsible for mortality in a diverse range of animals and plants, but their genome diversity and transcriptome responses in animal pathogenicity remain to be elucidated. We sequenced, assembled and annotated six chromosome-level FSSC clade 3 genomes of aquatic animal and plant host origins. We established a pathosystem and investigated the expression data of F. falciforme and F. keratoplasticum in Chinese softshell turtle (Pelodiscus sinensis) host.

RESULTS

Comparative analyses between the FSSC genomes revealed a spectrum of conservation patterns in chromosomes categorised into three compartments: core, fast-core (FC), and lineage-specific (LS). LS chromosomes contribute to variations in genomes size, with up to 42.2% of variations between F. vanettenii strains. Each chromosome compartment varied in structural architectures, with FC and LS chromosomes contain higher proportions of repetitive elements with genes enriched in functions related to pathogenicity and niche expansion. We identified differences in both selection in the coding sequences and DNA methylation levels between genome features and chromosome compartments which suggest a multi-speed evolution that can be traced back to the last common ancestor of Fusarium. We further demonstrated that F. falciforme and F. keratoplasticum are opportunistic pathogens by inoculating P. sinensis eggs and identified differentially expressed genes also associated with plant pathogenicity. These included the most upregulated genes encoding the CFEM (Common in Fungal Extracellular Membrane) domain.

CONCLUSIONS

The high-quality genome assemblies provided new insights into the evolution of FSSC chromosomes, which also serve as a resource for studies of fungal genome evolution and pathogenesis. This study also establishes an animal model for fungal pathogens of trans-kingdom hosts.

Karyotype variation, spontaneous genome rearrangements affecting chemical insensitivity, and expression level polymorphisms in the plant pathogen Phytophthora infestans revealed using its first chromosome-scale assembly

Natural isolates of the potato and tomato pathogen Phytophthora infestans exhibit substantial variation in virulence, chemical sensitivity, ploidy, and other traits. A chromosome-scale assembly was developed to expand genomic resources for this oomyceteous microbe, and used to explore the basis of variation. Using PacBio and Illumina data, a long-range linking library, and an optical map, an assembly was created and coalesced into 15 pseudochromosomes spanning 219 Mb using SNP-based genetic linkage data. De novo gene prediction combined with transcript evidence identified 19,981 protein-coding genes, plus about eight thousand tRNA genes. The chromosomes were comprised of a mosaic of gene-rich and gene-sparse regions plus very long centromeres. Genes exhibited a biased distribution across chromosomes, especially members of families encoding RXLR and CRN effectors which clustered on certain chromosomes. Strikingly, half of F1 progeny of diploid parents were polyploid or aneuploid. Substantial expression level polymorphisms between strains were identified, much of which could be attributed to differences in chromosome dosage, transposable element insertions, and adjacency to repetitive DNA. QTL analysis identified a locus on the right arm of chromosome 3 governing sensitivity to the crop protection chemical metalaxyl. Strains heterozygous for resistance often experienced megabase-sized deletions of that part of the chromosome when cultured on metalaxyl, increasing resistance due to loss of the sensitive allele. This study sheds light on diverse phenomena affecting variation in P. infestans and relatives, helps explain the prevalence of polyploidy in natural populations, and provides a new foundation for biologic and genetic investigations.

A global pangenome for the wheat fungal pathogen Pyrenophora tritici-repentis and prediction of effector protein structural homology

The adaptive potential of plant fungal pathogens is largely governed by the gene content of a species, consisting of core and accessory genes across the pathogen isolate repertoire. To approximate the complete gene repertoire of a globally significant crop fungal pathogen, a pan genomic analysis was undertaken for Pyrenophora tritici-repentis (Ptr), the causal agent of tan (or yellow) spot disease in wheat. In this study, 15 new Ptr genomes were sequenced, assembled and annotated, including isolates from three races not previously sequenced. Together with 11 previously published Ptr genomes, a pangenome for 26 Ptr isolates from Australia, Europe, North Africa and America, representing nearly all known races, revealed a conserved core-gene content of 57 % and presents a new Ptr resource for searching natural homologues (orthologues not acquired by horizontal transfer from another species) using remote protein structural homology. Here, we identify for the first time a non-synonymous mutation in the Ptr necrotrophic effector gene ToxB, multiple copies of the inactive toxb within an isolate, a distant natural Pyrenophora homologue of a known Parastagonopora nodorum necrotrophic effector (SnTox3), and clear genomic break points for the ToxA effector horizontal transfer region. This comprehensive genomic analysis of Ptr races includes nine isolates sequenced via long read technologies. Accordingly, these resources provide a more complete representation of the species, and serve as a resource to monitor variations potentially involved in pathogenicity.

Chromosome-level assemblies from diverse clades reveal limited structural and gene content variation in the genome of Candida glabrata

Background

Candida glabrata is an opportunistic yeast pathogen thought to have a large genetic and phenotypic diversity and a highly plastic genome. However, the lack of chromosome-level genome assemblies representing this diversity limits our ability to accurately establish how chromosomal structure and gene content vary across strains.

Results

Here, we expanded publicly available assemblies by using long-read sequencing technologies in twelve diverse strains, obtaining a final set of twenty-one chromosome-level genomes spanning the known C. glabrata diversity. Using comparative approaches, we inferred variation in chromosome structure and determined the pan-genome, including an analysis of the adhesin gene repertoire. Our analysis uncovered four new adhesin orthogroups and inferred a rich ancestral adhesion repertoire, which was subsequently shaped through a still ongoing process of gene loss, gene duplication, and gene conversion.

Conclusions

C. glabrata has a largely stable pan-genome except for a highly variable subset of genes encoding cell wall-associated functions. Adhesin repertoire was established for each strain and showed variability among clades.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01412-1.

Comparative profiling of canonical and non-canonical small RNAs in the rice blast fungus, Magnaporthe oryzae

RNA interference (RNAi) is divided into canonical, Dicer-dependent and non-canonical, Dicer-independent pathways according to Dicer protein dependency. However, sRNAs processed in a Dicer-independent manner have not been reported in plant pathogenic fungi, including Magnaporthe oryzae. We comparatively profiled the Dicer-dependent and -independent sRNAs of M. oryzae. Dicer-dependent sRNAs were 19–24-nt in length, had low strand-specificity, and showed a preference for uracil at the 5′-end. By contrast, Dicer-independent sRNAs presented irregular patterns in length distribution, high strand-specificity, and a preference for cytosine at the penultimate position. Dicer-dependent sRNA loci were mainly associated with LTR-transposons, while Dicer-independent sRNAs were associated with protein-coding genes and transposons. We identified MoERI-1, a non-canonical RNAi component, and profiled the sRNA and mRNA transcriptomes of ΔMoeri-1 at the mycelia and conidiation stages, as the mutant showed increased conidiation. We found that genes involved in conidiation and cell cycle were upregulated by MoERI-1 deletion. Furthermore, a comparison between sRNA and mRNA transcriptome revealed that MoERI-1-dependent sRNAs mediate the regulation of gene expression. Overall, these results showed that M. oryzae has non-canonical RNAi pathways distinct to the Dicer-dependent manner and exploits MoERI-1-dependent sRNAs to regulate the conidiation process.

Homokaryotic High-Quality Genome Assembly of Medicinal Fungi Wolfiporia hoelen Reveals Auto-Regulation and High-Temperature Adaption of Probable Two-Speed Genome

Sclerotia of Wolfiporia hoelen are one of the most important traditional Chinese medicines and are commonly used in China, Japan, Korea, and other Asian countries. In the present study, we presented the first high-quality homokaryotic genome of W. hoelen with 14 chromosomes which was evaluated with assembly index, telomere position detection, and whole-genome collinearity. A 64.44 Mb genome was assembled with a Contig N50 length of 3.76 Mb. The imbalanced distribution of transposons and chromosome characters revealed the probable two-speed genome of W. hoelen. High consistency between methylation and transposon conserved the genome stability. The expansion of the gene family about signal transduction and nutritional transport has intimate relationships with sclerotial formation. Up-regulation of expression for distinctive decomposition enzymes, ROS clearance genes, biosynthesis of unsaturated fatty acids, and change of the cell wall components maintained high-speed growth of mycelia that may be the high-temperature adaption strategy of W. hoelen. Further, the analysis of mating-control genes demonstrated that HD3 probably had no function on mating recognition, with the HD protein in a distant genetic with known species. Overall, the high-quality genome of W. hoelen provided crucial information for genome structure and stability, high-temperature adaption, and sexual and asexual process.

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.

A potato late blight resistance gene protects against multiple Phytophthora species by recognizing a broadly conserved RXLR-WY effector

Species of the genus Phytophthora, the plant killer, cause disease and reduce yields in many crop plants. Although many Resistance to Phytophthora infestans (Rpi) genes effective against potato late blight have been cloned, few have been cloned against other Phytophthora species. Most Rpi genes encode nucleotide-binding domain, leucine-rich repeat-containing (NLR) immune receptor proteins that recognize RXLR (Arg-X-Leu-Arg) effectors. However, whether NLR proteins can recognize RXLR effectors from multiple Phytophthora species has rarely been investigated. Here, we identified a new RXLR-WY effector AVRamr3 from P. infestans that is recognized by Rpi-amr3 from a wild Solanaceae species Solanum americanum. Rpi-amr3 associates with AVRamr3 in planta. AVRamr3 is broadly conserved in many different Phytophthora species, and the recognition of AVRamr3 homologs by Rpi-amr3 activates resistance against multiple Phytophthora pathogens, including the tobacco black shank disease and cacao black pod disease pathogens P. parasitica and P. palmivora. Rpi-amr3 is thus the first characterized resistance gene that acts against P. parasitica or P. palmivora. These findings suggest a novel path to redeploy known R genes against different important plant pathogens.

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.

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

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

Unveiling the Core Effector Proteins of Oil Palm Pathogen Ganoderma boninense via Pan-Secretome Analysis

Ganoderma boninense is the major causal agent of basal stem rot (BSR) disease in oil palm, causing the progressive rot of the basal part of the stem. Despite its prominence, the key pathogenicity determinants for the aggressive nature of hemibiotrophic infection remain unknown. In this study, genome sequencing and the annotation of G. boninense T10 were carried out using the Illumina sequencing platform, and comparative genome analysis was performed with previously reported G. boninense strains (NJ3 and G3). The pan-secretome of G. boninense was constructed and comprised 937 core orthogroups, 243 accessory orthogroups, and 84 strain-specific orthogroups. In total, 320 core orthogroups were enriched with candidate effector proteins (CEPs) that could be classified as carbohydrate-active enzymes, hydrolases, and non-catalytic proteins. Differential expression analysis revealed an upregulation of five CEP genes that was linked to the suppression of PTI signaling cascade, while the downregulation of four CEP genes was linked to the inhibition of PTI by preventing host defense elicitation. Genome architecture analysis revealed the one-speed architecture of the G. boninense genome and the lack of preferential association of CEP genes to transposable elements. The findings obtained from this study aid in the characterization of pathogenicity determinants and molecular biomarkers of BSR disease.

Nutrient Availability Does Not Affect Community Assembly in Root-Associated Fungi but Determines Fungal Effects on Plant Growth

Nonmycorrhizal root-colonizing fungi are key determinants of plant growth, driving processes ranging from pathogenesis to stress alleviation. Evidence suggests that they might also facilitate host access to soil nutrients in a mycorrhiza-like manner, but the extent of their direct contribution to plant nutrition is unknown. To study how widespread such capacity is across root-colonizing fungi, we surveyed soils in nutrient-limiting habitats using plant baits to look for fungal community changes in response to nutrient conditions. We established a fungal culture collection and used Arabidopsis thaliana inoculation bioassays to assess the ability of fungi to facilitate host's growth in the presence of organic nutrients unavailable to plants. Plant baits captured a representation of fungal communities extant in natural habitats and showed that nutrient limitation has little influence on community assembly. Arabidopsis thaliana inoculated with 31 phylogenetically diverse fungi exhibited a consistent fungus-driven growth promotion when supplied with organic nutrients compared to untreated plants. However, direct phosphorus measurement and RNA-seq data did not support enhanced nutrient uptake but rather that growth effects may result from changes in the plant's immune response to colonization. The widespread and consistent host responses to fungal colonization suggest that distinct, locally adapted nonmycorrhizal fungi affect plant performance across habitats. IMPORTANCE Recent studies have shown that root-associated fungi that do not engage in classical mycorrhizal associations can facilitate the hosts' access to nutrients in a mycorrhiza-like manner. However, the generality of this capacity remains to be tested. Root-associated fungi are frequently deemed major determinants of plant diversity and performance, but in the vast majority of cases their ecological roles in nature remain unknown. Assessing how these plant symbionts affect plant productivity, diversity, and fitness is important to understanding how plant communities function. Recent years have seen important advances in the understanding of the main drivers of the diversity and structure of plant microbiomes, but a major challenge is still linking community properties with function. This study contributes to the understanding of the cryptic function of root-associated fungi by testing their ability to participate in a specific process: nutrient acquisition by plants.

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.

Ancestral Chromosomes for Family Peronosporaceae Inferred from a Telomere-to-Telomere Genome Assembly of Peronospora effusa

Downy mildew disease of spinach, caused by the oomycete Peronospora effusa, causes major losses to spinach production. In this study, the 17 chromosomes of P. effusa were assembled telomere-to-telomere, using Pacific Biosciences high-fidelity reads. Of these, 16 chromosomes are complete and gapless; chromosome 15 contains one gap bridging the nucleolus organizer region. This is the first telomere-to-telomere genome assembly for an oomycete. Putative centromeric regions were identified on all chromosomes. This new assembly enables a reevaluation of the genomic composition of Peronospora spp.; the assembly was almost double the size and contained more repeat sequences than previously reported for any Peronospora species. Genome fragments consistently underrepresented in six previously reported assemblies of P. effusa typically encoded repeats. Some genes annotated as encoding effectors were organized into multigene clusters on several chromosomes. Putative effectors were annotated on 16 of the 17 chromosomes. The intergenic distances between annotated genes were consistent with compartmentalization of the genome into gene-dense and gene-sparse regions. Genes encoding putative effectors were enriched in gene-sparse regions. The near-gapless assembly revealed apparent horizontal gene transfer from Ascomycete fungi. Gene order was highly conserved between P. effusa and the genetically oriented assembly of the oomycete Bremia lactucae; high levels of synteny were also detected with Phytophthora sojae. Extensive synteny between phylogenetically distant species suggests that many other oomycete species may have similar chromosome organization. Therefore, this assembly provides the foundation for genomic analyses of diverse oomycetes.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Evolutionary Significance of Fungal Hypermutators: Lessons Learned from Clinical Strains and Implications for Fungal Plant Pathogens

Rapid evolution of fungal pathogens poses a serious threat to medicine and agriculture. The mutation rate determines the pace of evolution of a fungal pathogen. Hypermutator fungal strains have an elevated mutation rate owing to certain defects such as those in the DNA mismatch repair system. Studies in Saccharomyces cerevisiae show that hypermutators expedite evolution by generating beneficial alleles at a faster pace than the wild-type strains. However, an accumulation of deleterious alleles in a hypermutator may reduce its fitness. The balance between fitness cost and mutation benefit determines the prevalence of hypermutators in a population. This balance is affected by a complex interaction of ploidy, mode of reproduction, population size, and recent population history. Studies in human fungal pathogens like Aspergillus fumigatus, Candida albicans, Candida glabrata, Cryptococcus deuterogattii, and Cryptococcus neoformans have highlighted the importance of hypermutators in host adaptation and development of antifungal resistance. However, a critical examination of hypermutator biology, experimental evolution studies, and epidemiological studies suggests that hypermutators may impact evolutionary investigations. This review aims to integrate the knowledge about biology, experimental evolution, and dynamics of fungal hypermutators to critically examine the evolutionary role of hypermutators in fungal pathogen populations and project implications of hypermutators in the evolution of fungal plant pathogen populations. Understanding the factors determining the emergence and evolution of fungal hypermutators can open a novel avenue of managing rapidly evolving fungal pathogens in medicine and agriculture.

De Novo Assembly of Plasmodium knowlesi Genomes From Clinical Samples Explains the Counterintuitive Intrachromosomal Organization of Variant SICAvar and kir Multiple Gene Family Members

Plasmodium knowlesi, a malaria parasite of Old World macaque monkeys, is used extensively to model Plasmodium biology. Recently, P. knowlesi was found in the human population of Southeast Asia, particularly Malaysia. P. knowlesi causes uncomplicated to severe and fatal malaria in the human host with features in common with the more prevalent and virulent malaria caused by Plasmodium falciparum. As such, P. knowlesi presents a unique opportunity to develop experimental translational model systems for malaria pathophysiology informed by clinical data from same-species human infections. Experimental lines of P. knowlesi represent well-characterized genetically stable parasites, and to maximize their utility as a backdrop for understanding malaria pathophysiology, genetically diverse contemporary clinical isolates, essentially wild-type, require comparable characterization. The Oxford Nanopore PCR-free long-read sequencing platform was used to sequence and de novo assemble P. knowlesi genomes from frozen clinical samples. The sequencing platform and assembly pipelines were designed to facilitate capturing data and describing, for the first time, P. knowlesi schizont-infected cell agglutination (SICA) var and Knowlesi-Interspersed Repeats (kir) multiple gene families in parasites acquired from nature. The SICAvar gene family members code for antigenically variant proteins analogous to the virulence-associated P. falciparum erythrocyte membrane protein (PfEMP1) multiple var gene family. Evidence presented here suggests that the SICAvar family members have arisen through a process of gene duplication, selection pressure, and variation. Highly evolving genes including PfEMP1family members tend to be restricted to relatively unstable sub-telomeric regions that drive change with core genes protected in genetically stable intrachromosomal locations. The comparable SICAvar and kir gene family members are counter-intuitively located across chromosomes. Here, we demonstrate that, in contrast to conserved core genes, SICAvar and kir genes occupy otherwise gene-sparse chromosomal locations that accommodate rapid evolution and change. The novel methods presented here offer the malaria research community not only new tools to generate comprehensive genome sequence data from small clinical samples but also new insight into the complexity of clinically important real-world parasites.

Enemy or ally: a genomic approach to elucidate the lifestyle of Phyllosticta citrichinaensis

Members of the fungal genus Phyllosticta can colonize a variety of plant hosts, including several Citrus species such as Citrus sinensis (orange), Citrus limon (lemon), and Citrus maxima (pomelo). Some Phyllosticta species have the capacity to cause disease, such as Citrus Black Spot, while others have only been observed as endophytes. Thus far, genomic differences underlying lifestyle adaptations of Phyllosticta species have not yet been studied. Furthermore, the lifestyle of Phyllosticta citrichinaensis is ambiguous, as it has been described as a weak pathogen but Koch's postulates may not have been established and the presence of this species was never reported to cause any crop or economic losses. Here, we examined the genomic differences between pathogenic and endophytic Phyllosticta spp. colonizing Citrus and specifically aimed to elucidate the lifestyle of Phyllosticta citrichinaensis. We found several genomic differences between species of different lifestyles, including groups of genes that were only present in pathogens or endophytes. We also observed that species, based on their carbohydrate active enzymes, group independent of their phylogenetic association, and this clustering correlated with trophy prediction. Phyllosticta citrichinaensis shows an intermediate lifestyle, sharing genomic and phenotypic attributes of both pathogens and endophytes. We thus present the first genomic comparison of multiple citrus-colonizing pathogens and endophytes of the genus Phyllosticta, and therefore provide the basis for further comparative studies into the lifestyle adaptations within this genus.

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.

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.

Tracking of Diversity and Evolution in the Brown Rot Fungi Monilinia fructicola, Monilinia fructigena, and Monilinia laxa

Monilinia species are among the most devastating fungi worldwide as they cause brown rot and blossom blight on fruit trees. To understand the molecular bases of their pathogenic lifestyles, we compared the newly assembled genomes of single strains of Monilinia fructicola, M. fructigena and M. laxa, with those of Botrytis cinerea and Sclerotinia sclerotiorum, as the closest species within Sclerotiniaceae. Phylogenomic analysis of orthologous proteins and syntenic investigation suggest that M. laxa is closer to M. fructigena than M. fructicola, and is closest to the other investigated Sclerotiniaceae species. This indicates that M. laxa was the earliest result of the speciation process. Distinct evolutionary profiles were observed for transposable elements (TEs). M. fructicola and M. laxa showed older bursts of TE insertions, which were affected (mainly in M. fructicola) by repeat-induced point (RIP) mutation gene silencing mechanisms. These suggested frequent occurrence of the sexual process in M. fructicola. More recent TE expansion linked with low RIP action was observed in M. fructigena, with very little in S. sclerotiorum and B. cinerea. The detection of active non-syntenic TEs is indicative of horizontal gene transfer and has resulted in alterations in specific gene functions. Analysis of candidate effectors, biosynthetic gene clusters for secondary metabolites and carbohydrate-active enzymes, indicated that Monilinia genus has multiple virulence mechanisms to infect host plants, including toxins, cell-death elicitor, putative virulence factors and cell-wall-degrading enzymes. Some species-specific pathogenic factors might explain differences in terms of host plant and organ preferences between M. fructigena and the other two Monilinia species.

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.

Chromosome-Level Genome Assembly of a Human Fungal Pathogen Reveals Synteny among Geographically Distinct Species

Histoplasma capsulatum, a dimorphic fungal pathogen, is the most common cause of fungal respiratory infections in immunocompetent hosts. Histoplasma is endemic in the Ohio and Mississippi River Valleys in the United States and is also distributed worldwide. Previous studies have revealed at least eight clades, each specific to a geographic location: North American classes 1 and 2 (NAm 1 and NAm 2), Latin American groups A and B (LAm A and LAm B), Eurasian, Netherlands, Australian and African, and an additional distinct lineage (H81) comprised of Panamanian isolates. Previously assembled Histoplasma genomes are highly fragmented, with the highly repetitive G217B (NAm 2) strain, which has been used for most whole-genome-scale transcriptome studies, assembled into over 250 contigs. In this study, we set out to fully assemble the repeat regions and characterize the large-scale genome architecture of Histoplasma species. We resequenced five Histoplasma strains (WU24 [NAm 1], G217B [NAm 2], H88 [African], G186AR [Panama], and G184AR [Panama]) using Oxford Nanopore Technologies long-read sequencing technology. Here, we report chromosomal-level assemblies for all five strains, which exhibit extensive synteny among the geographically distant Histoplasma isolates. The new assemblies revealed that RYP2, a major regulator of morphology and virulence, is duplicated in G186AR. In addition, we mapped previously generated transcriptome data sets onto the newly assembled chromosomes. Our analyses revealed that the expression of transposons and transposon-embedded genes are upregulated in yeast phase compared to mycelial phase in the G217B and H88 strains. This study provides an important resource for fungal researchers and further highlights the importance of chromosomal-level assemblies in analyzing high-throughput data sets. IMPORTANCE Histoplasma species are dimorphic fungi causing significant morbidity and mortality worldwide. These fungi grow as mold in the soil and as budding yeast within the human host. Histoplasma can be isolated from soil in diverse regions, including North America, South America, Africa, and Europe. Phylogenetically distinct species of Histoplasma have been isolated and sequenced. However, for the commonly used strains, genome assemblies have been fragmented, leading to underutilization of genome-scale data. This study provides chromosome-level assemblies of the commonly used Histoplasma strains using long-read sequencing technology. Comparative analysis of these genomes shows largely conserved gene order within the chromosomes. Mapping existing transcriptome data on these new assemblies reveals clustering of transcriptionally coregulated genes. The results of this study highlight the importance of obtaining chromosome-level assemblies in understanding the biology of human fungal pathogens.

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.

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.

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 whole genome duplication drives the genome evolution of Phytophthora betacei, a closely related species to Phytophthora infestans

BACKGROUND

Pathogens of the genus Phytophthora are the etiological agents of many devastating diseases in several high-value crops and forestry species such as potato, tomato, cocoa, and oak, among many others. Phytophthora betacei is a recently described species that causes late blight almost exclusively in tree tomatoes, and it is closely related to Phytophthora infestans that causes the disease in potato crops and other Solanaceae. This study reports the assembly and annotation of the genomes of P. betacei P8084, the first of its species, and P. infestans RC1-10, a Colombian strain from the EC-1 lineage, using long-read SMRT sequencing technology.

RESULTS

Our results show that P. betacei has the largest sequenced genome size of the Phytophthora genus so far with 270 Mb. A moderate transposable element invasion and a whole genome duplication likely explain its genome size expansion when compared to P. infestans, whereas P. infestans RC1-10 has expanded its genome under the activity of transposable elements. The high diversity and abundance (in terms of copy number) of classified and unclassified transposable elements in P. infestans RC1-10 relative to P. betacei bears testimony of the power of long-read technologies to discover novel repetitive elements in the genomes of organisms. Our data also provides support for the phylogenetic placement of P. betacei as a standalone species and as a sister group of P. infestans. Finally, we found no evidence to support the idea that the genome of P. betacei P8084 follows the same gene-dense/gense-sparse architecture proposed for P. infestans and other filamentous plant pathogens.

CONCLUSIONS

This study provides the first genome-wide picture of P. betacei and expands the genomic resources available for P. infestans. This is a contribution towards the understanding of the genome biology and evolutionary history of Phytophthora species belonging to the subclade 1c.

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.

Sugarcane Smut, Caused by Sporisorium scitamineum, a Major Disease of Sugarcane: A Contemporary Review

Sugarcane smut caused by the fungus Sporisorium scitamineum is one of the major diseases of sugarcane worldwide, causing significant losses in productivity and profitability of this perennial crop. Teliospores of this fungus are airborne, can travel long distances, and remain viable in hot and dry conditions for >6 months. The disease is easily recognized by its long whiplike sorus produced on the apex or side shoots of sugarcane stalks. Each sorus can release ≤100 million teliospores in a day; the spores are small (≤7.5 µ) and light and can survive in harsh environmental conditions. The airborne teliospores are the primary mode of smut spread around the world and across cane-growing regions. The most effective method of managing this disease is via resistant varieties. Because of the complex genomic makeup of sugarcane, selection for resistant traits is difficult in sugarcane breeding programs. In recent times, the application of molecular markers as a rapid tool of discarding susceptible genotypes early in the selection program has been investigated. Large effect resistance loci have been identified and have the potential to be used for marker-assisted selection to increase the frequency of resistant breeding lines in breeding programs. Recent developments in omics technologies (genomics, transcriptomics, proteomics, and metabolomics) have contributed to our understanding and provided insights into the mechanism of resistance and susceptibility. This knowledge will further our understanding of smut and its interactions with sugarcane genotypes and aid in the development of durable resistant varieties.

Epigenetic modifications affect the rate of spontaneous mutations in a pathogenic fungus

Mutations are the source of genetic variation and the substrate for evolution. Genome-wide mutation rates appear to be affected by selection and are probably adaptive. Mutation rates are also known to vary along genomes, possibly in response to epigenetic modifications, but causality is only assumed. In this study we determine the direct impact of epigenetic modifications and temperature stress on mitotic mutation rates in a fungal pathogen using a mutation accumulation approach. Deletion mutants lacking epigenetic modifications confirm that histone mark H3K27me3 increases whereas H3K9me3 decreases the mutation rate. Furthermore, cytosine methylation in transposable elements (TE) increases the mutation rate 15-fold resulting in significantly less TE mobilization. Also accessory chromosomes have significantly higher mutation rates. Finally, we find that temperature stress substantially elevates the mutation rate. Taken together, we find that epigenetic modifications and environmental conditions modify the rate and the location of spontaneous mutations in the genome and alter its evolutionary trajectory.

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.

Metatranscriptomic Comparison of Endophytic and Pathogenic Fusarium-Arabidopsis Interactions Reveals Plant Transcriptional Plasticity

Plants are continuously exposed to beneficial and pathogenic microbes, but how plants recognize and respond to friends versus foes remains poorly understood. Here, we compared the molecular response of Arabidopsis thaliana independently challenged with a Fusarium oxysporum endophyte Fo47 versus a pathogen Fo5176. These two F. oxysporum strains share a core genome of about 46 Mb, in addition to 1,229 and 5,415 unique accessory genes. Metatranscriptomic data reveal a shared pattern of expression for most plant genes (about 80%) in responding to both fungal inoculums at all timepoints from 12 to 96 h postinoculation (HPI). However, the distinct responding genes depict transcriptional plasticity, as the pathogenic interaction activates plant stress responses and suppresses functions related to plant growth and development, while the endophytic interaction attenuates host immunity but activates plant nitrogen assimilation. The differences in reprogramming of the plant transcriptome are most obvious in 12 HPI, the earliest timepoint sampled, and are linked to accessory genes in both fungal genomes. Collectively, our results indicate that the A. thaliana and F. oxysporum interaction displays both transcriptome conservation and plasticity in the early stages of infection, providing insights into the fine-tuning of gene regulation underlying plant differential responses to fungal endophytes and pathogens.[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.

Remarkable recent changes in the genetic diversity of the avirulence gene AvrStb6 in global populations of the wheat pathogen Zymoseptoria tritici

Septoria tritici blotch (STB), caused by the fungus Zymoseptoria tritici, is one of the most economically important diseases of wheat. Recently, both factors of a gene-for-gene interaction between Z. tritici and wheat, the wheat receptor-like kinase Stb6 and the Z. tritici secreted effector protein AvrStb6, have been identified. Previous analyses revealed a high diversity of AvrStb6 haplotypes present in earlier Z. tritici isolate collections, with up to c.18% of analysed isolates possessing the avirulence isoform of AvrStb6 identical to that originally identified in the reference isolate IPO323. With Stb6 present in many commercial wheat cultivars globally, we aimed to assess potential changes in AvrStb6 genetic diversity and the incidence of haplotypes allowing evasion of Stb6-mediated resistance in more recent Z. tritici populations. Here we show, using targeted resequencing of AvrStb6, that this gene is universally present in field isolates sampled from major wheat-growing regions of the world in 2013-2017. However, in contrast to the data from previous AvrStb6 population studies, we report a complete absence of the originally described avirulence isoform of AvrStb6 amongst modern Z. tritici isolates. Moreover, a remarkably small number of haplotypes, each encoding AvrStb6 protein isoforms conditioning virulence on Stb6-containing wheat, were found to predominate among modern Z. tritici isolates. A single virulence isoform of AvrStb6 was found to be particularly abundant throughout the global population. These findings indicate that, despite the ability of Z. tritici to sexually reproduce on resistant hosts, AvrStb6 avirulence haplotypes tend to be eliminated in subsequent populations.

Small RNAs in Plant Immunity and Virulence of Filamentous Pathogens

Gene silencing guided by small RNAs governs a broad range of cellular processes in eukaryotes. Small RNAs are important components of plant immunity because they contribute to pathogen-triggered transcription reprogramming and directly target pathogen RNAs. Recent research suggests that silencing of pathogen genes by plant small RNAs occurs not only during viral infection but also in nonviral pathogens through a process termed host-induced gene silencing, which involves trans-species small RNA trafficking. Similarly, small RNAs are also produced by eukaryotic pathogens and regulate virulence. This review summarizes the small RNA pathways in both plants and filamentous pathogens, including fungi and oomycetes, and discusses their role in host-pathogen interactions. We highlight secondarysmall interfering RNAs of plants as regulators of immune receptor gene expression and executors of host-induced gene silencing in invading pathogens. The current status and prospects of trans-species gene silencing at the host-pathogen interface are discussed.

Plant and Fungal Genome Editing to Enhance Plant Disease Resistance Using the CRISPR/Cas9 System

Crop production has been substantially reduced by devastating fungal and oomycete pathogens, and these pathogens continue to threaten global food security. Although chemical and cultural controls have been used for crop protection, these involve continuous costs and time and fungicide resistance among plant pathogens has been increasingly reported. The most efficient way to protect crops from plant pathogens is cultivation of disease-resistant cultivars. However, traditional breeding approaches are laborious and time intensive. Recently, the CRISPR/Cas9 system has been utilized to enhance disease resistance among different crops such as rice, cacao, wheat, tomato, and grape. This system allows for precise genome editing of various organisms via RNA-guided DNA endonuclease activity. Beyond genome editing in crops, editing the genomes of fungal and oomycete pathogens can also provide new strategies for plant disease management. This review focuses on the recent studies of plant disease resistance against fungal and oomycete pathogens using the CRISPR/Cas9 system. For long-term plant disease management, the targeting of multiple plant disease resistance mechanisms with CRISPR/Cas9 and insights gained by probing fungal and oomycete genomes with this system will be powerful approaches.

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.

The "Bipartite" Structure of the First Genome of Ampelomyces quisqualis, a Common Hyperparasite and Biocontrol Agent of Powdery Mildews, May Point to Its Evolutionary Origin from Plant Pathogenic Fungi

Powdery mildews are among the most important plant pathogens worldwide, which are often attacked in the field by mycoparasitic fungi belonging to the genus Ampelomyces. The taxonomy of the genus Ampelomyces is unresolved, but well-supported molecular operational taxonomic units were repeatedly defined suggesting that the genus may include at least four to seven species. Some Ampelomyces strains were commercialized as biocontrol agents of crop pathogenic powdery mildews. However, the genomic mechanisms underlying their mycoparasitism are still poorly understood. To date, the draft genome of a single Ampelomyces strain, designated as HMLAC 05119, has been released. We report a high-quality, annotated hybrid draft genome assembly of A. quisqualis strain BRIP 72107, which, based on phylogenetic analyses, is not conspecific with HMLAC 05119. The constructed genome is 40.38 Mb in size, consisting of 24 scaffolds with an N50 of 2.99 Mb and 96.2% completeness. Our analyses revealed "bipartite" structure of Ampelomyces genomes, where GC-balanced genomic regions are interspersed by longer or shorter stretches of AT-rich regions. This is also a hallmark of many plant pathogenic fungi and provides further evidence for evolutionary affinity of Ampelomyces species to plant pathogenic fungi. The high-quality genome and annotation produced here provide an important resource for future genomic studies of mycoparasitisim to decipher molecular mechanisms underlying biocontrol processes and natural tritrophic interactions.

A 20-kb lineage-specific genomic region tames virulence in pathogenic amphidiploid Verticillium longisporum

Amphidiploid fungal Verticillium longisporum strains Vl43 and Vl32 colonize the plant host Brassica napus but differ in their ability to cause disease symptoms. These strains represent two V. longisporum lineages derived from different hybridization events of haploid parental Verticillium strains. Vl32 and Vl43 carry same-sex mating-type genes derived from both parental lineages. Vl32 and Vl43 similarly colonize and penetrate plant roots, but asymptomatic Vl32 proliferation in planta is lower than virulent Vl43. The highly conserved Vl43 and Vl32 genomes include less than 1% unique genes, and the karyotypes of 15 or 16 chromosomes display changed genetic synteny due to substantial genomic reshuffling. A 20 kb Vl43 lineage-specific (LS) region apparently originating from the Verticillium dahliae-related ancestor is specific for symptomatic Vl43 and encodes seven genes, including two putative transcription factors. Either partial or complete deletion of this LS region in Vl43 did not reduce virulence but led to induction of even more severe disease symptoms in rapeseed. This suggests that the LS insertion in the genome of symptomatic V. longisporum Vl43 mediates virulence-reducing functions, limits damage on the host plant, and therefore tames Vl43 from being even more virulent.

Selection for Plastic, Pathogen-Inducible Recombination in a Red Queen Model with Diploid Antagonists

Antagonistic interactions and co-evolution between a host and its parasite are known to cause oscillations in the population genetic structure of both species (Red Queen dynamics). Potentially, such oscillations may select for increased sex and recombination in the host, although theoretical models suggest that this happens under rather restricted values of selection intensity, epistasis, and other parameters. Here, we explore a model in which the diploid parasite succeeds to infect the diploid host only if their phenotypes at the interaction-mediating loci match. Whenever regular oscillations emerge in this system, we test whether plastic, pathogen-inducible recombination in the host can be favored over the optimal constant recombination. Two forms of the host recombination dependence on the parasite pressure were considered: either proportionally to the risk of infection (prevention strategy) or upon the fact of infection (remediation strategy). We show that both forms of plastic recombination can be favored, although relatively infrequently (up to 11% of all regimes with regular oscillations, and up to 20% of regimes with obligate parasitism). This happens under either strong overall selection and high recombination rate in the host, or weak overall selection and low recombination rate in the host. In the latter case, the system’s dynamics are considerably more complex. The prevention strategy is favored more often than the remediation one. It is noteworthy that plastic recombination can be favored even when any constant recombination is rejected, making plasticity an evolutionary mechanism for the rescue of host recombination.

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.

Impact of Repetitive DNA Elements on Snake Genome Biology and Evolution

The distinctive biology and unique evolutionary features of snakes make them fascinating model systems to elucidate how genomes evolve and how variation at the genomic level is interlinked with phenotypic-level evolution. Similar to other eukaryotic genomes, large proportions of snake genomes contain repetitive DNA, including transposable elements (TEs) and satellite repeats. The importance of repetitive DNA and its structural and functional role in the snake genome, remain unclear. This review highlights the major types of repeats and their proportions in snake genomes, reflecting the high diversity and composition of snake repeats. We present snakes as an emerging and important model system for the study of repetitive DNA under the impact of sex and microchromosome evolution. We assemble evidence to show that certain repetitive elements in snakes are transcriptionally active and demonstrate highly dynamic lineage-specific patterns as repeat sequences. We hypothesize that particular TEs can trigger different genomic mechanisms that might contribute to driving adaptive evolution in snakes. Finally, we review emerging approaches that may be used to study the expression of repetitive elements in complex genomes, such as snakes. The specific aspects presented here will stimulate further discussion on the role of genomic repeats in shaping snake evolution.

Transposable Elements Contribute to Genome Dynamics and Gene Expression Variation in the Fungal Plant Pathogen Verticillium dahliae

Transposable elements (TEs) are a major source of genetic and regulatory variation in their host genome and are consequently thought to play important roles in evolution. Many fungal and oomycete plant pathogens have evolved dynamic and TE-rich genomic regions containing genes that are implicated in host colonization and adaptation. TEs embedded in these regions have typically been thought to accelerate the evolution of these genomic compartments, but little is known about their dynamics in strains that harbor them. Here, we used whole-genome sequencing data of 42 strains of the fungal plant pathogen Verticillium dahliae to systematically identify polymorphic TEs that may be implicated in genomic as well as in gene expression variation. We identified 2,523 TE polymorphisms and characterize a subset of 8% of the TEs as polymorphic elements that are evolutionary younger, less methylated, and more highly expressed when compared with the remaining 92% of the total TE complement. As expected, the polyrmorphic TEs are enriched in the adaptive genomic regions. Besides, we observed an association of polymorphic TEs with pathogenicity-related genes that localize nearby and that display high expression levels. Collectively, our analyses demonstrate that TE dynamics in V. dahliae contributes to genomic variation, correlates with expression of pathogenicity-related genes, and potentially impacts the evolution of adaptive genomic regions.

Comparative analyses of the Hymenoscyphus fraxineus and Hymenoscyphus albidus genomes reveals potentially adaptive differences in secondary metabolite and transposable element repertoires

Background

The dieback epidemic decimating common ash (Fraxinus excelsior) in Europe is caused by the invasive fungus Hymenoscyphus fraxineus. In this study we analyzed the genomes of H. fraxineus and H. albidus, its native but, now essentially displaced, non-pathogenic sister species, and compared them with several other members of Helotiales. The focus of the analyses was to identify signals in the genome that may explain the rapid establishment of H. fraxineus and displacement of H. albidus.

Results

The genomes of H. fraxineus and H. albidus showed a high level of synteny and identity. The assembly of H. fraxineus is 13 Mb longer than that of H. albidus’, most of this difference can be attributed to higher dispersed repeat content (i.e. transposable elements [TEs]) in H. fraxineus. In general, TE families in H. fraxineus showed more signals of repeat-induced point mutations (RIP) than in H. albidus, especially in Long-terminal repeat (LTR)/Copia and LTR/Gypsy elements.

Comparing gene family expansions and 1:1 orthologs, relatively few genes show signs of positive selection between species. However, several of those did appeared to be associated with secondary metabolite genes families, including gene families containing two of the genes in the H. fraxineus-specific, hymenosetin biosynthetic gene cluster (BGC).

Conclusion

The genomes of H. fraxineus and H. albidus show a high degree of synteny, and are rich in both TEs and BGCs, but the genomic signatures also indicated that H. albidus may be less well equipped to adapt and maintain its ecological niche in a rapidly changing environment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07837-2.

Diverse Trajectories Drive the Expression of a Giant Virus in the Oomycete Plant Pathogen Phytophthora parasitica

Giant viruses of amoebas, recently classified in the class Megaviricetes, are a group of viruses that can infect major eukaryotic lineages. We previously identified a set of giant virus sequences in the genome of Phytophthora parasitica, an oomycete and a devastating major plant pathogen. How viral insertions shape the structure and evolution of the invaded genomes is unclear, but it is known that the unprecedented functional potential of giant viruses is the result of an intense genetic interplay with their hosts. We previously identified a set of giant virus sequences in the genome of P. parasitica, an oomycete and a devastating major plant pathogen. Here, we show that viral pieces are found in a 550-kb locus and are organized in three main clusters. Viral sequences, namely RNA polymerases I and II and a major capsid protein, were identified, along with orphan sequences, as a hallmark of giant viruses insertions. Mining of public databases and phylogenetic reconstructions suggest an ancient association of oomycetes and giant viruses of amoeba, including faustoviruses, African swine fever virus (ASFV) and pandoraviruses, and that a single viral insertion occurred early in the evolutionary history of oomycetes prior to the Phytophthora–Pythium radiation, estimated at ∼80 million years ago. Functional annotation reveals that the viral insertions are located in a gene sparse region of the Phytophthora genome, characterized by a plethora of transposable elements (TEs), effectors and other genes potentially involved in virulence. Transcription of viral genes was investigated through analysis of RNA-Seq data and qPCR experiments. We show that most viral genes are not expressed, and that a variety of mechanisms, including deletions, TEs insertions and RNA interference may contribute to transcriptional repression. However, a gene coding a truncated copy of RNA polymerase II along a set of neighboring sequences have been shown to be expressed in a wide range of physiological conditions, including responses to stress. These results, which describe for the first time the endogenization of a giant virus in an oomycete, contribute to challenge our view of Phytophthora evolution.