Accessories Make the Outfit: Accessory Chromosomes and Other Dispensable DNA Regions in Plant-Pathogenic Fungi

Multiple Horizontal Mini-chromosome Transfers Drive Genome Evolution of Clonal Blast Fungus Lineages

Crop disease pandemics are often driven by asexually reproducing clonal lineages of plant pathogens that reproduce asexually. How these clonal pathogens continuously adapt to their hosts despite harboring limited genetic variation, and in absence of sexual recombination remains elusive. Here, we reveal multiple instances of horizontal chromosome transfer within pandemic clonal lineages of the blast fungus Magnaporthe (Syn. Pyricularia) oryzae. We identified a horizontally transferred 1.2Mb accessory mini-chromosome which is remarkably conserved between M. oryzae isolates from both the rice blast fungus lineage and the lineage infecting Indian goosegrass (Eleusine indica), a wild grass that often grows in the proximity of cultivated cereal crops. Furthermore, we show that this mini-chromosome was horizontally acquired by clonal rice blast isolates through at least nine distinct transfer events over the past three centuries. These findings establish horizontal mini-chromosome transfer as a mechanism facilitating genetic exchange among different host-associated blast fungus lineages. We propose that blast fungus populations infecting wild grasses act as genetic reservoirs that drive genome evolution of pandemic clonal lineages that afflict cereal crops.

Complete genome of the Medicago anthracnose fungus, Colletotrichum destructivum, reveals a mini-chromosome-like region within a core chromosome

Colletotrichum destructivum (Cd) is a phytopathogenic fungus causing significant economic losses on forage legume crops (Medicago and Trifolium species) worldwide. To gain insights into the genetic basis of fungal virulence and host specificity, we sequenced the genome of an isolate from Medicago sativa using long-read (PacBio) technology. The resulting genome assembly has a total length of 51.7 Mb and comprises ten core chromosomes and two accessory chromosomes, all of which were sequenced from telomere to telomere. A total of 15, 631 gene models were predicted, including genes encoding potentially pathogenicity-related proteins such as candidate-secreted effectors (484), secondary metabolism key enzymes (110) and carbohydrate-active enzymes (619). Synteny analysis revealed extensive structural rearrangements in the genome of Cd relative to the closely related Brassicaceae pathogen, Colletotrichum higginsianum. In addition, a 1.2 Mb species-specific region was detected within the largest core chromosome of Cd that has all the characteristics of fungal accessory chromosomes (transposon-rich, gene-poor, distinct codon usage), providing evidence for exchange between these two genomic compartments. This region was also unique in having undergone extensive intra-chromosomal segmental duplications. Our findings provide insights into the evolution of accessory regions and possible mechanisms for generating genetic diversity in this asexual fungal pathogen.

Cations impact the biodegradation of iodosulfuron-methyl herbicidal ionic liquids by fungi

In the framework of this study, six fungal isolates which demonstrated a high capability for biodegrading iodosulphuron-methyl sodium as well as herbicidal ionic liquids based on this herbicide were isolated from different soil samples. The isolates were identified based on the ITS region, whereas biodegradation residues were determined based on LC-MS/MS. Depending on the isolate, the half-lives values of the biodegraded herbicide or herbicidal ionic liquid ranged significantly from just 1.25 days to more than 40 days. The research findings unveiled that the structure of cations is a central limiting factor affecting fungal growth and herbicide transformation in case of ionic liquids. The length of the alkyl chain has been identified as the primary driver of herbicide toxicity, emphasizing the importance of structural factors in herbicide design. In cases when dodecyl(2-hydroxyethyl)dimethyl cation was used, its biodegradation ranged from 0 to approx. 20% and the biodegradability of the iodosulfuron-methyl was notably limited for the majority of the studied isolates. This knowledge provides guidance for development and selection of herbicides with reduced environmental impact. This study highlights the ecological importance of soil fungi, their potential role in herbicide biodegradation, the influence of cations on fungal growth and herbicide transformation, and the structural factors governing herbicide toxicity. Further research in these areas may lead to more efficient and environmentally friendly approaches to herbicide management.

Template-Based Modelling of the Structure of Fungal Effector Proteins

The discovery of new fungal effector proteins is necessary to enable the screening of cultivars for disease resistance. Sequence-based bioinformatics methods have been used for this purpose, but only a limited number of functional effector proteins have been successfully predicted and subsequently validated experimentally. A significant obstacle is that many fungal effector proteins discovered so far lack sequence similarity or conserved sequence motifs. The availability of experimentally determined three-dimensional (3D) structures of a number of effector proteins has recently highlighted structural similarities amongst groups of sequence-dissimilar fungal effectors, enabling the search for similar structural folds amongst effector sequence candidates. We have applied template-based modelling to predict the 3D structures of candidate effector sequences obtained from bioinformatics predictions and the PHI-BASE database. Structural matches were found not only with ToxA- and MAX-like effector candidates but also with non-fungal effector-like proteins-including plant defensins and animal venoms-suggesting the broad conservation of ancestral structural folds amongst cytotoxic peptides from a diverse range of distant species. Accurate modelling of fungal effectors were achieved using RaptorX. The utility of predicted structures of effector proteins lies in the prediction of their interactions with plant receptors through molecular docking, which will improve the understanding of effector-plant interactions.

Genome comparisons reveal accessory genes crucial for the evolution of apple Glomerella leaf spot pathogenicity in Colletotrichum fungi

Apple Glomerella leaf spot (GLS) is an emerging fungal disease caused by Colletotrichum fructicola and other Colletotrichum species. These species are polyphyletic and it is currently unknown how these pathogens convergently evolved to infect apple. We generated chromosome-level genome assemblies of a GLS-adapted isolate and a non-adapted isolate in C. fructicola using long-read sequencing. Additionally, we resequenced 17 C. fructicola and C. aenigma isolates varying in GLS pathogenicity using short-read sequencing. Genome comparisons revealed a conserved bipartite genome architecture involving minichromosomes (accessory chromosomes) shared by C. fructicola and other closely related species within the C. gloeosporioides species complex. Moreover, two repeat-rich genomic regions (1.61 Mb in total) were specifically conserved among GLS-pathogenic isolates in C. fructicola and C. aenigma. Single-gene deletion of 10 accessory genes within the GLS-specific regions of C. fructicola identified three that were essential for GLS pathogenicity. These genes encoded a putative non-ribosomal peptide synthetase, a flavin-binding monooxygenase and a small protein with unknown function. These results highlight the crucial role accessory genes play in the evolution of Colletotrichum pathogenicity and imply the significance of an unidentified secondary metabolite in GLS pathogenesis.

Frequent horizontal chromosome transfer between asexual fungal insect pathogens

Entire chromosomes are typically only transmitted vertically from one generation to the next. The horizontal transfer of such chromosomes has long been considered improbable, yet gained recent support in several pathogenic fungi where it may affect the fitness or host specificity. To date, it is unknown how these transfers occur, how common they are, and whether they can occur between different species. In this study, we show multiple independent instances of horizontal transfers of the same accessory chromosome between two distinct strains of the asexual entomopathogenic fungus Metarhizium robertsii during experimental co-infection of its insect host, the Argentine ant. Notably, only the one chromosome-but no other-was transferred from the donor to the recipient strain. The recipient strain, now harboring the accessory chromosome, exhibited a competitive advantage under certain host conditions. By phylogenetic analysis, we further demonstrate that the same accessory chromosome was horizontally transferred in a natural environment between M. robertsii and another congeneric insect pathogen, Metarhizium guizhouense. Hence, horizontal chromosome transfer is not limited to the observed frequent events within species during experimental infections but also occurs naturally across species. The accessory chromosome that was transferred contains genes that may be involved in its preferential horizontal transfer or support its establishment. These genes encode putative histones and histone-modifying enzymes, as well as putative virulence factors. Our study reveals that both intra- and interspecies horizontal transfer of entire chromosomes is more frequent than previously assumed, likely representing a not uncommon mechanism for gene exchange.

Rapid mini-chromosome divergence among fungal isolates causing wheat blast outbreaks in Bangladesh and Zambia

The fungal pathogen, Magnaporthe oryzae Triticum pathotype, causing wheat blast disease was first identified in South America and recently spread across continents to South Asia and Africa. Here, we studied the genetic relationship among isolates found on the three continents. Magnaporthe oryzae strains closely related to a South American field isolate B71 were found to have caused the wheat blast outbreaks in South Asia and Africa. Genomic variation among isolates from the three continents was examined using an improved B71 reference genome and whole-genome sequences. We found strong evidence to support that the outbreaks in Bangladesh and Zambia were caused by the introductions of genetically separated isolates, although they were all close to B71 and, therefore, collectively referred to as the B71 branch. In addition, B71 branch strains carried at least one supernumerary mini-chromosome. Genome assembly of a Zambian strain revealed that its mini-chromosome was similar to the B71 mini-chromosome but with a high level of structural variation. Our findings show that while core genomes of the multiple introductions are highly similar, the mini-chromosomes have undergone marked diversification. The maintenance of the mini-chromosome and rapid genomic changes suggest the mini-chromosomes may serve important virulence or niche adaptation roles under diverse environmental conditions.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Meiosis in an asymmetric dikaryotic genome of Tremella fuciformis Tr01 facilitates new chromosome formation

BACKGROUND

The dikaryotic stage dominates most of the life cycle in basidiomycetes, and each cell carries two different haploid nuclei. Accurate phasing of these two nuclear genomes and their interactions have long been of interest.

RESULTS

We combine PacBio HiFi reads, Nanopore ultra-long reads, and Hi-C data to generate a complete, high-quality asymmetric dikaryotic genome of Tremella fuciformis Tr01, including Haplotypes A and B genomes. We assemble a meiotic haploid DBZ04 genome and detect three recombination events in these two haplotypes. We identify several chromosomal rearrangements that lead to differences in chromosome number, length, content, and sequence arrangement between these two haplotypes. Each nucleus contains a two-speed genome, harboring three accessory chromosomes and two accessory compartments that affect horizontal chromatin transfer between nuclei. We find few basidiospores are ejected from fruiting bodies of Tr01. Most monospore isolates sequenced belong to Tr01-Haplotype A genome architecture. More than one-third of monospore isolates carry one or two extra chromosomes including Chr12B and two new chromosomes ChrN1 and ChrN2. We hypothesize that homologous regions of seven sister chromatids pair into a large complex during meiosis, followed by inter-chromosomal recombination at physical contact sites and formation of new chromosomes.

CONCLUSION

We assemble two haplotype genomes of T. fuciformis Tr01 and provide the first overview of basidiomycetous genomes with discrete genomic architecture. Meiotic activities of asymmetric dikaryotic genomes result in formation of new chromosomes, aneuploidy of some daughter cells, and inviability of most other daughter cells. We propose a new approach for breeding of sporeless mushroom.

Population genomic analyses suggest recent dispersal events of the pathogen Cercospora zeina into East and Southern African maize cropping systems

A serious factor hampering global maize production is gray leaf spot disease. Cercospora zeina is one of the causative pathogens, but population genomics analysis of C. zeina is lacking. We conducted whole-genome Illumina sequencing of a representative set of 30 C. zeina isolates from Kenya and Uganda (East Africa) and Zambia, Zimbabwe, and South Africa (Southern Africa). Selection of the diverse set was based on microsatellite data from a larger collection of the pathogen. Pangenome analysis of the C. zeina isolates was done by (1) de novo assembly of the reads with SPAdes, (2) annotation with BRAKER, and (3) protein clustering with OrthoFinder. A published long-read assembly of C. zeina (CMW25467) from Zambia was included and annotated using the same pipeline. This analysis revealed 790 non-shared accessory and 10,677 shared core orthogroups (genes) between the 31 isolates. Accessory gene content was largely shared between isolates from all countries, with a few genes unique to populations from Southern Africa (32) or East Africa (6). There was a significantly higher proportion of effector genes in the accessory secretome (44%) compared to the core secretome (24%). PCA, ADMIXTURE, and phylogenetic analysis using a neighbor-net network indicated a population structure with a geographical subdivision between the East African isolates and the Southern African isolates, although gene flow was also evident. The small pangenome and partial population differentiation indicated recent dispersal of C. zeina into Africa, possibly from 2 regional founder populations, followed by recurrent gene flow owing to widespread maize production across sub-Saharan Africa.

Population genomic analysis reveals geographic structure and climatic diversification for Macrophomina phaseolina isolated from soybean and dry bean across the United States, Puerto Rico, and Colombia

Macrophomina phaseolina causes charcoal rot, which can significantly reduce yield and seed quality of soybean and dry bean resulting from primarily environmental stressors. Although charcoal rot has been recognized as a warm climate-driven disease of increasing concern under global climate change, knowledge regarding population genetics and climatic variables contributing to the genetic diversity of M. phaseolina is limited. This study conducted genome sequencing for 95 M. phaseolina isolates from soybean and dry bean across the continental United States, Puerto Rico, and Colombia. Inference on the population structure using 76,981 single nucleotide polymorphisms (SNPs) revealed that the isolates exhibited a discrete genetic clustering at the continental level and a continuous genetic differentiation regionally. A majority of isolates from the United States (96%) grouped in a clade with a predominantly clonal genetic structure, while 88% of Puerto Rican and Colombian isolates from dry bean were assigned to a separate clade with higher genetic diversity. A redundancy analysis (RDA) was used to estimate the contributions of climate and spatial structure to genomic variation (11,421 unlinked SNPs). Climate significantly contributed to genomic variation at a continental level with temperature seasonality explaining the most variation while precipitation of warmest quarter explaining the most when spatial structure was accounted for. The loci significantly associated with multivariate climate were found closely to the genes related to fungal stress responses, including transmembrane transport, glycoside hydrolase activity and a heat-shock protein, which may mediate climatic adaptation for M. phaseolina. On the contrary, limited genome-wide differentiation among populations by hosts was observed. These findings highlight the importance of population genetics and identify candidate genes of M. phaseolina that can be used to elucidate the molecular mechanisms that underly climatic adaptation to the changing climate.

Attack of the clones: Population genetics reveals clonality of Colletotrichum lupini, the causal agent of lupin anthracnose

Colletotrichum lupini, the causative agent of lupin anthracnose, affects lupin cultivation worldwide. Understanding its population structure and evolutionary potential is crucial to design successful disease management strategies. The objective of this study was to employ population genetics to investigate the diversity, evolutionary dynamics, and molecular basis of the interaction of this notorious lupin pathogen with its host. A collection of globally representative C. lupini isolates was genotyped through triple digest restriction site-associated DNA sequencing, resulting in a data set of unparalleled resolution. Phylogenetic and structural analysis could distinguish four independent lineages (I-IV). The strong population structure and high overall standardized index of association (r̅_d ) indicates that C. lupini reproduces clonally. Different morphologies and virulence patterns on white lupin (Lupinus albus) and Andean lupin (Lupinus mutabilis) were observed between and within clonal lineages. Isolates belonging to lineage II were shown to have a minichromosome that was also partly present in lineage III and IV, but not in lineage I isolates. Variation in the presence of this minichromosome could imply a role in host-pathogen interaction. All four lineages were present in the South American Andes region, which is suggested to be the centre of origin of this species. Only members of lineage II have been found outside South America since the 1990s, indicating it as the current pandemic population. As a seedborne pathogen, C. lupini has mainly spread through infected but symptomless seeds, stressing the importance of phytosanitary measures to prevent future outbreaks of strains that are yet confined to South America.

Combined pangenomics and transcriptomics reveals core and redundant virulence processes in a rapidly evolving fungal plant pathogen

BACKGROUND

Studying genomic variation in rapidly evolving pathogens potentially enables identification of genes supporting their "core biology", being present, functional and expressed by all strains or "flexible biology", varying between strains. Genes supporting flexible biology may be considered to be "accessory", whilst the "core" gene set is likely to be important for common features of a pathogen species biology, including virulence on all host genotypes. The wheat-pathogenic fungus Zymoseptoria tritici represents one of the most rapidly evolving threats to global food security and was the focus of this study.

RESULTS

We constructed a pangenome of 18 European field isolates, with 12 also subjected to RNAseq transcription profiling during infection. Combining this data, we predicted a "core" gene set comprising 9807 sequences which were (1) present in all isolates, (2) lacking inactivating polymorphisms and (3) expressed by all isolates. A large accessory genome, consisting of 45% of the total genes, was also defined. We classified genetic and genomic polymorphism at both chromosomal and individual gene scales. Proteins required for essential functions including virulence had lower-than average sequence variability amongst core genes. Both core and accessory genomes encoded many small, secreted candidate effector proteins that likely interact with plant immunity. Viral vector-mediated transient in planta overexpression of 88 candidates failed to identify any which induced leaf necrosis characteristic of disease. However, functional complementation of a non-pathogenic deletion mutant lacking five core genes demonstrated that full virulence was restored by re-introduction of the single gene exhibiting least sequence polymorphism and highest expression.

CONCLUSIONS

These data support the combined use of pangenomics and transcriptomics for defining genes which represent core, and potentially exploitable, weaknesses in rapidly evolving pathogens.

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

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

Complete Genome Sequences of Septoria linicola: A Resource for Studying a Damaging Flax Pathogen

Fungal genus Septoria causes diseases in a wide range of plants. Here, we report the first genome sequences of two strains of Septoria linicola, the causal agent of the pasmo disease of flax (Linum usitatissimum). The genome of the first strain, SE15195, was fully assembled in 16 chromosomes, while 35 unitigs were obtained for a second strain, SE14017. Structural annotations predicted 13,096 and 13,085 protein-encoding genes and transposable elements content of 19.0 and 18.1% of the genome for SE15195 and SE14017, respectively. The four smaller chromosomes 13 to 16 show genomics features of potential accessory chromosomes. The assembly of these two genomes is a new resource for studying S. linicola and improving management of pasmo. [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.

Species Diversity in the Fusarium solani (Neocosmospora) Complex and Their Pathogenicity for Plants and Humans

The Fusarium solani species complex is a large group of soil saprotrophs with a broad adaptive potential, which allows the fungi to exist under various conditions and to parasitize on different hosts. The review analyzes the modern data concerning the genetic peculiarities of species from this complex by the example of F. solani f. sp. pisi and generalizes the data on the most widespread species pathogenic for both plants and humans. The enhanced resistance of the F. solani species complex to the most of modern antifungal agents and the need for novel therapeutic agents against fusariosis has been considered.

The extrachromosomal circular DNAs of the rice blast pathogen Magnaporthe oryzae contain a wide variety of LTR retrotransposons, genes, and effectors

BACKGROUND

One of the ways genomes respond to stress is by producing extrachromosomal circular DNAs (eccDNAs). EccDNAs can contain genes and dramatically increase their copy number. They can also reinsert into the genome, generating structural variation. They have been shown to provide a source of phenotypic and genotypic plasticity in several species. However, whole circularome studies have so far been limited to a few model organisms. Fungal plant pathogens are a serious threat to global food security in part because of their rapid adaptation to disease prevention strategies. Understanding the mechanisms fungal pathogens use to escape disease control is paramount to curbing their threat.

RESULTS

We present a whole circularome sequencing study of the rice blast pathogen, Magnaporthe oryzae. We find that M. oryzae has a highly diverse circularome that contains many genes and shows evidence of large LTR retrotransposon activity. We find that genes enriched on eccDNAs in M. oryzae occur in genomic regions prone to presence-absence variation and that disease-associated genes are frequently on eccDNAs. Finally, we find that a subset of genes is never present on eccDNAs in our data, which indicates that the presence of these genes on eccDNAs is selected against.

CONCLUSIONS

Our study paves the way to understanding how eccDNAs contribute to adaptation in M. oryzae. Our analysis also reveals how M. oryzae eccDNAs differ from those of other species and highlights the need for further comparative characterization of eccDNAs across species to gain a better understanding of these molecules.

Fungal Effectoromics: A World in Constant Evolution

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

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.

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.

Chromosomal polymorphism of the Ceratocystis fimbriata species complex in Brazil

Ceratocystis fimbriata is an important pathogen that causes wilt in several plant species. Despite the importance of this pathogen, knowledge about its karyotypic polymorphism and genomic architecture is limited. The main objective of this study was to investigate the karyotype of isolates of the C. fimbriata species complex from different host plants and geographical origins in Brazil. First, the identity of the isolates was confirmed conducting multilocus sequence analysis (MLSA) phylogeny using β-tubulin (TUBB), translation elongation factor 1α (TEF-1α) and mating-type (MAT1 and MAT2) gene sequences. To investigate the chromosomal polymorphism, two conditions of pulsed-field gel electrophoresis (PFGE) were used and the karyotypes of the isolates obtained. The retrotransposon-microsatellite amplified polymorphism (REMAP) molecular marker was utilized to assess the genetic variability among isolates. In the MLSA utilizing the concatenated gene sequences, Ceratocystis cacaofunesta and C. fimbriata formed separate clades, but considerable variation among C. fimbriata isolates was observed. Polymorphism in chromosome number and size was found, indicating the existence of genomic differences among isolates and occurrence of chromosomal rearrangements in the species complex. The number of chromosomes varied from seven to nine and the estimated minimum chromosome sizes were estimated to be between 2.7 to 6.0 Mbp. Small polymorphic chromosomes ranging from 1.2 to 1.8 Mbp were observed in all isolates, raising the hypothesis that they could be supernumerary chromosomes. REMAP analysis revealed a high genetic variability and that isolates from the same host tend to group together in a same cluster. Our results bring new insights into the chromosomal diversity and genome organization of the C. fimbriata complex.

Non-Mendelian transmission of accessory chromosomes in fungi

Non-Mendelian transmission has been reported for various genetic elements, ranging from small transposons to entire chromosomes. One prime example of such a transmission pattern are B chromosomes in plants and animals. Accessory chromosomes in fungi are similar to B chromosomes in showing presence/absence polymorphism and being non-essential. How these chromosomes are transmitted during meiosis is however poorly understood—despite their often high impact on the fitness of the host. For several fungal organisms, a non-Mendelian transmission or a mechanistically unique meiotic drive of accessory chromosomes have been reported. In this review, we provide an overview of the possible mechanisms that can cause the non-Mendelian transmission or meiotic drives of fungal accessory chromosomes. We compare processes responsible for the non-Mendelian transmission of accessory chromosomes for different fungal eukaryotes and discuss the structural traits of fungal accessory chromosomes affecting their meiotic transmission. We conclude that research on fungal accessory chromosomes, due to their small size, ease of sequencing, and epigenetic profiling, can complement the study of B chromosomes in deciphering factors that influence and regulate the non-Mendelian transmission of entire chromosomes.

Using genomics to understand the mechanisms of virulence and drug resistance in fungal pathogens

Fungal pathogens pose an increasingly worrying threat to human health, food security and ecosystem diversity. To tackle fungal infections and improve current diagnostic and therapeutic tools it is necessary to understand virulence and antifungal drug resistance mechanisms in diverse species. Recent advances in genomics approaches have provided a suitable framework to understand these phenotypes, which ultimately depend on genetically encoded determinants. In this work, we review how the study of genome sequences has been key to ascertain the bases of virulence and drug resistance traits. We focus on the contribution of comparative genomics, population genomics and directed evolution studies. In addition, we discuss how different types of genomic mutations (small or structural variants) contribute to intraspecific differences in virulence or drug resistance. Finally, we review current challenges in the field and anticipate future directions to solve them. In summary, this work provides a short overview of how genomics can be used to understand virulence and drug resistance in fungal pathogens.

Lifestyle Transitions in Fusarioid Fungi are Frequent and Lack Clear Genomic Signatures

The fungal genus Fusarium (Ascomycota) includes well-known plant pathogens that are implicated in diseases worldwide, and many of which have been genome sequenced. The genus also encompasses other diverse lifestyles, including species found ubiquitously as asymptomatic-plant inhabitants (endophytes). Here, we produced structurally annotated genome assemblies for five endophytic Fusarium strains, including the first whole-genome data for Fusarium chuoi. Phylogenomic reconstruction of Fusarium and closely related genera revealed multiple and frequent lifestyle transitions, the major exception being a monophyletic clade of mutualist insect symbionts. Differential codon usage bias and increased codon optimisation separated Fusarium sensu stricto from allied genera. We performed computational prediction of candidate secreted effector proteins (CSEPs) and carbohydrate-active enzymes (CAZymes)-both likely to be involved in the host-fungal interaction-and sought evidence that their frequencies could predict lifestyle. However, phylogenetic distance described gene variance better than lifestyle did. There was no significant difference in CSEP, CAZyme, or gene repertoires between phytopathogenic and endophytic strains, although we did find some evidence that gene copy number variation may be contributing to pathogenicity. Large numbers of accessory CSEPs (i.e., present in more than one taxon but not all) and a comparatively low number of strain-specific CSEPs suggested there is a limited specialisation among plant associated Fusarium species. We also found half of the core genes to be under positive selection and identified specific CSEPs and CAZymes predicted to be positively selected on certain lineages. Our results depict fusarioid fungi as prolific generalists and highlight the difficulty in predicting pathogenic potential in the group.

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.

Long reads and Hi-C sequencing illuminate the two-compartment genome of the model arbuscular mycorrhizal symbiont Rhizophagus irregularis

Chromosome folding links genome structure with gene function by generating distinct nuclear compartments and topologically associating domains. In mammals, these undergo preferential interactions and regulate gene expression. However, their role in fungal genome biology is unclear. Here, we combine Nanopore (ONT) sequencing with chromatin conformation capture sequencing (Hi-C) to reveal chromosome and epigenetic diversity in a group of obligate plant symbionts: the arbuscular mycorrhizal fungi (AMF). We find that five phylogenetically distinct strains of the model AMF Rhizophagus irregularis carry 33 chromosomes with substantial within-species variability in size, as well as in gene and repeat content. Strain-specific Hi-C contact maps reveal a 'checkerboard' pattern that underline two dominant euchromatin (A) and heterochromatin (B) compartments. Each compartment differs in the level of gene transcription, regulation of candidate effectors and methylation frequencies. The A-compartment is more gene-dense and contains most core genes, while the B-compartment is more repeat-rich and has higher rates of chromosomal rearrangement. While the B-compartment is transcriptionally repressed, it has significantly more secreted proteins and in planta upregulated candidate effectors, suggesting a possible host-induced change in chromosome conformation. Overall, this study provides a fine-scale view into the genome biology and evolution of model plant symbionts, and opens avenues to study the epigenetic mechanisms that modify chromosome folding during host-microbe interactions.

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.

An automated and combinative method for the predictive ranking of candidate effector proteins of fungal plant pathogens

Fungal plant-pathogens promote infection of their hosts through the release of 'effectors'-a broad class of cytotoxic or virulence-promoting molecules. Effectors may be recognised by resistance or sensitivity receptors in the host, which can determine disease outcomes. Accurate prediction of effectors remains a major challenge in plant pathology, but if achieved will facilitate rapid improvements to host disease resistance. This study presents a novel tool and pipeline for the ranking of predicted effector candidates-Predector-which interfaces with multiple software tools and methods, aggregates disparate features that are relevant to fungal effector proteins, and applies a pairwise learning to rank approach. Predector outperformed a typical combination of secretion and effector prediction methods in terms of ranking performance when applied to a curated set of confirmed effectors derived from multiple species. We present Predector ( https://github.com/ccdmb/predector ) as a useful tool for the ranking of predicted effector candidates, which also aggregates and reports additional supporting information relevant to effector and secretome prediction in a simple, efficient, and reproducible manner.

Remote homology clustering identifies lowly conserved families of effector proteins in plant-pathogenic fungi

Plant diseases caused by fungal pathogens are typically initiated by molecular interactions between 'effector' molecules released by a pathogen and receptor molecules on or within the plant host cell. In many cases these effector-receptor interactions directly determine host resistance or susceptibility. The search for fungal effector proteins is a developing area in fungal-plant pathology, with more than 165 distinct confirmed fungal effector proteins in the public domain. For a small number of these, novel effectors can be rapidly discovered across multiple fungal species through the identification of known effector homologues. However, many have no detectable homology by standard sequence-based search methods. This study employs a novel comparison method (RemEff) that is capable of identifying protein families with greater sensitivity than traditional homology-inference methods, leveraging a growing pool of confirmed fungal effector data to enable the prediction of novel fungal effector candidates by protein family association. Resources relating to the RemEff method and data used in this study are available from https://figshare.com/projects/Effector_protein_remote_homology/87965.

Quantitative interactions: the disease outcome of Botrytis cinerea across the plant kingdom

Botrytis cinerea is a fungal pathogen that causes necrotic disease on more than a thousand known hosts widely spread across the plant kingdom. How B. cinerea interacts with such extensive host diversity remains largely unknown. To address this question, we generated an infectivity matrix of 98 strains of B. cinerea on 90 genotypes representing eight host plants. This experimental infectivity matrix revealed that the disease outcome is largely explained by variations in either the host resistance or pathogen virulence. However, the specific interactions between host and pathogen account for 16% of the disease outcome. Furthermore, the disease outcomes cluster among genotypes of a species but are independent of the relatedness between hosts. When analyzing the host specificity and virulence of B. cinerea, generalist strains are predominant. In this fungal necrotroph, specialization may happen by a loss in virulence on most hosts rather than an increase of virulence on a specific host. To uncover the genetic architecture of Botrytis host specificity and virulence, a genome-wide association study (GWAS) was performed and revealed up to 1492 genes of interest. The genetic architecture of these traits is widespread across the B. cinerea genome. The complexity of the disease outcome might be explained by hundreds of functionally diverse genes putatively involved in adjusting the infection to diverse hosts.

Apicidin biosynthesis is linked to accessory chromosomes in Fusarium poae isolates

BACKGROUND

Fusarium head blight is a disease of global concern that reduces crop yields and renders grains unfit for consumption due to mycotoxin contamination. Fusarium poae is frequently associated with cereal crops showing symptoms of Fusarium head blight. While previous studies have shown F. poae isolates produce a range of known mycotoxins, including type A and B trichothecenes, fusarins and beauvericin, genomic analysis suggests that this species may have lineage-specific accessory chromosomes with secondary metabolite biosynthetic gene clusters awaiting description.

METHODS

We examined the biosynthetic potential of 38 F. poae isolates from Eastern Canada using a combination of long-read and short-read genome sequencing and untargeted, high resolution mass spectrometry metabolome analysis of extracts from isolates cultured in multiple media conditions.

RESULTS

A high-quality assembly of isolate DAOMC 252244 (Fp157) contained four core chromosomes as well as seven additional contigs with traits associated with accessory chromosomes. One of the predicted accessory contigs harbours a functional biosynthetic gene cluster containing homologs of all genes associated with the production of apicidins. Metabolomic and genomic analyses confirm apicidins are produced in 4 of the 38 isolates investigated and genomic PCR screening detected the apicidin synthetase gene APS1 in approximately 7% of Eastern Canadian isolates surveyed.

CONCLUSIONS

Apicidin biosynthesis is linked to isolate-specific putative accessory chromosomes in F. poae. The data produced here are an important resource for furthering our understanding of accessory chromosome evolution and the biosynthetic potential of F. poae.

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.

Identification of rye B chromosome-associated peptides by mass spectrometry

B chromosomes (Bs) are supernumerary dispensable components of the standard genome (A chromosomes, As) that have been found in many eukaryotes. So far, it is unkown whether the B-derived transcripts translate to proteins or if the host proteome is changed due to the presence of Bs. Comparative mass spectrometry was performed using the protein samples isolated from shoots of rye plants with and without Bs. We aimed to identify B-associated peptides and analyzed the effects of Bs on the total proteome. Our comparative proteome analysis demonstrates that the presence of rye Bs affects the total proteome including different biological function processes. We found 319 of 16 776 quantified features in at least three out of five +B plants but not in 0B plants; 31 of 319 features were identified as B-associated peptide features. According to our data mining, one B-specific protein fragment showed similarity to a glycine-rich RNA binding protein which differed from its A-paralogue by two amino acid insertions. Our result represents a milestone in B chromosome research, because this is the first report to demonstrate the existence of Bs changing the proteome of the host.

The Plot Thickens: Haploid and Triploid-Like Thalli, Hybridization, and Biased Mating Type Ratios in Letharia

The study of the reproductive biology of lichen fungal symbionts has been traditionally challenging due to their complex lifestyles. Against the common belief of haploidy, a recent genomic study found a triploid-like signal in Letharia. Here, we infer the genome organization and reproduction in Letharia by analyzing genomic data from a pure culture and from thalli, and performing a PCR survey of the MAT locus in natural populations. We found that the read count variation in the four Letharia specimens, including the pure culture derived from a single sexual spore of L. lupina, is consistent with haploidy. By contrast, the L. lupina read counts from a thallus' metagenome are triploid-like. Characterization of the mating-type locus revealed a conserved heterothallic configuration across the genus, along with auxiliary genes that we identified. We found that the mating-type distributions are balanced in North America for L. vulpina and L. lupina, suggesting widespread sexual reproduction, but highly skewed in Europe for L. vulpina, consistent with predominant asexuality. Taken together, we propose that Letharia fungi are heterothallic and typically haploid, and provide evidence that triploid-like individuals are hybrids between L. lupina and an unknown Letharia lineage, reconciling classic systematic and genetic studies with recent genomic observations.

Chromosome-level genome assembly and manually-curated proteome of model necrotroph Parastagonospora nodorum Sn15 reveals a genome-wide trove of candidate effector homologs, and redundancy of virulence-related functions within an accessory chromosome

Background

The fungus Parastagonospora nodorum causes septoria nodorum blotch (SNB) of wheat (Triticum aestivum) and is a model species for necrotrophic plant pathogens. The genome assembly of reference isolate Sn15 was first reported in 2007. P. nodorum infection is promoted by its production of proteinaceous necrotrophic effectors, three of which are characterised – ToxA, Tox1 and Tox3.

Results

A chromosome-scale genome assembly of P. nodorum Australian reference isolate Sn15, which combined long read sequencing, optical mapping and manual curation, produced 23 chromosomes with 21 chromosomes possessing both telomeres. New transcriptome data were combined with fungal-specific gene prediction techniques and manual curation to produce a high-quality predicted gene annotation dataset, which comprises 13,869 high confidence genes, and an additional 2534 lower confidence genes retained to assist pathogenicity effector discovery. Comparison to a panel of 31 internationally-sourced isolates identified multiple hotspots within the Sn15 genome for mutation or presence-absence variation, which was used to enhance subsequent effector prediction. Effector prediction resulted in 257 candidates, of which 98 higher-ranked candidates were selected for in-depth analysis and revealed a wealth of functions related to pathogenicity. Additionally, 11 out of the 98 candidates also exhibited orthology conservation patterns that suggested lateral gene transfer with other cereal-pathogenic fungal species. Analysis of the pan-genome indicated the smallest chromosome of 0.4 Mbp length to be an accessory chromosome (AC23). AC23 was notably absent from an avirulent isolate and is predominated by mutation hotspots with an increase in non-synonymous mutations relative to other chromosomes. Surprisingly, AC23 was deficient in effector candidates, but contained several predicted genes with redundant pathogenicity-related functions.

Conclusions

We present an updated series of genomic resources for P. nodorum Sn15 – an important reference isolate and model necrotroph – with a comprehensive survey of its predicted pathogenicity content.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07699-8.

Pan-Mitogenomics Approach Discovers Diversity and Dynamism in the Prominent Brown Rot Fungal Pathogens

Monilinia fructicola and Monilinia laxa species are the most destructive and economically devastating fungal plant pathogens causing brown rot disease on stone and pome fruits worldwide. Mitochondrial genomes (mitogenomes) play critical roles influencing the mechanisms and directions of the evolution of fungal pathogens. The pan-mitogenomics approach predicts core and accessory regions of the mitochondrial genomes and explains the gain or loss of variation within and between species. The present study is a fungal pan-mitogenome of M. fructicola (N = 8) and M. laxa (N = 8) species. The completely sequenced and annotated mitogenomes showed high variability in size within and between the species. The mitogenomes of M. laxa were larger, ranging from 178,351 to 179,780bp, than the mitogenomes of M. fructicola, ranging from 158,607 to 167,838bp. However, size variation within the species showed that M. fructicola isolates were more variable in the size range than M. laxa isolates. All the mitogenomes included conserved mitochondrial genes, as well as variable regions including different mobile introns encoding homing endonucleases or maturase, non-coding introns, and repetitive elements. The linear model analysis supported the hypothesis that the mitogenome size expansion is due to presence of variable (accessory) regions. Gene synteny was mostly conserved among all samples, with the exception for order of the rps3 in the mitogenome of one isolate. The mitogenomes presented AT richness; however, A/T and G/C skew varied among the mitochondrial genes. The purifying selection was detected in almost all the protein-coding genes (PCGs) between the species. However, cytochrome b was the only gene showing a positive selection signal among the total samples. Combined datasets of amino acid sequences of 14 core mitochondrial PCGs and rps3 obtained from this study together with published mitochondrial genome sequences from some other species from Heliotales were used to infer a maximum likelihood (ML) phylogenetic tree. ML tree indicated that both Monilinia species highly diverged from each other as well as some other fungal species from the same order. Mitogenomes harbor much information about the evolution of fungal plant pathogens, which could be useful to predict pathogenic life strategies.

Accessory Chromosome-Acquired Secondary Metabolism in Plant Pathogenic Fungi: The Evolution of Biotrophs Into Host-Specific Pathogens

Accessory chromosomes are strain- or pathotype-specific chromosomes that exist in addition to the core chromosomes of a species and are generally not considered essential to the survival of the organism. Among pathogenic fungal species, accessory chromosomes harbor pathogenicity or virulence factor genes, several of which are known to encode for secondary metabolites that are involved in plant tissue invasion. Accessory chromosomes are of particular interest due to their capacity for horizontal transfer between strains and their dynamic "crosstalk" with core chromosomes. This review focuses exclusively on secondary metabolism (including mycotoxin biosynthesis) associated with accessory chromosomes in filamentous fungi and the role accessory chromosomes play in the evolution of secondary metabolite gene clusters. Untargeted metabolomics profiling in conjunction with genome sequencing provides an effective means of linking secondary metabolite products with their respective biosynthetic gene clusters that reside on accessory chromosomes. While the majority of literature describing accessory chromosome-associated toxin biosynthesis comes from studies of Alternaria pathotypes, the recent discovery of accessory chromosome-associated biosynthetic genes in Fusarium species offer fresh insights into the evolution of biosynthetic enzymes such as non-ribosomal peptide synthetases (NRPSs), polyketide synthases (PKSs) and regulatory mechanisms governing their expression.

Molecular and Environmental Triggering Factors of Pathogenicity of Fusarium oxysporum and F. solani Isolates Involved in the Coffee Corky-Root Disease

Coffee corky-root disease causes serious damages to coffee crop and is linked to combined infection of Fusarium spp. and root-knot nematodes Meloidogyne spp. In this study, 70 Fusarium isolates were collected from both roots of healthy coffee plants and with corky-root disease symptoms. A phylogenetic analysis, and the detection of pathogenicity SIX genes and toxigenicity Fum genes was performed for 59 F. oxysporum and 11 F. solani isolates. Based on the molecular characterization, seven F. oxysporum and three F. solani isolates were assessed for their pathogenicity on coffee seedlings under optimal watering and water stress miming root-knot nematode effect on plants. Our results revealed that a drastic increment of plant colonization capacity and pathogenicity on coffee plants of some Fusarium isolates was caused by water stress. The pathogenicity on coffee of F. solani linked to coffee corky-root disease and the presence of SIX genes in this species were demonstrated for the first time. Our study provides evidence for understanding the pathogenic basis of F. oxysporum and F. solani isolates on coffee and revealed the presence of SIX and Fum genes as one of their pathogenicity-related mechanisms. We also highlight the relevance of chlorophyll, a fluorescence as an early and high-throughput phenotyping tool in Fusarium pathogenicity studies on coffee.

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

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

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

Dynamics in Secondary Metabolite Gene Clusters in Otherwise Highly Syntenic and Stable Genomes in the Fungal Genus Botrytis

Fungi of the genus Botrytis infect >1,400 plant species and cause losses in many crops. Besides the broad host range pathogen Botrytis cinerea, most other species are restricted to a single host. Long-read technology was used to sequence genomes of eight Botrytis species, mostly pathogenic on Allium species, and the related onion white rot fungus, Sclerotium cepivorum. Most assemblies contained <100 contigs, with the Botrytis aclada genome assembled in 16 gapless chromosomes. The core genome and pan-genome of 16 Botrytis species were defined and the secretome, effector, and secondary metabolite repertoires analyzed. Among those genes, none is shared among all Allium pathogens and absent from non-Allium pathogens. The genome of each of the Allium pathogens contains 8-39 predicted effector genes that are unique for that single species, none stood out as potential determinant for host specificity. Chromosome configurations of common ancestors of the genus Botrytis and family Sclerotiniaceae were reconstructed. The genomes of B. cinerea and B. aclada were highly syntenic with only 19 rearrangements between them. Genomes of Allium pathogens were compared with ten other Botrytis species (nonpathogenic on Allium) and with 25 Leotiomycetes for their repertoire of secondary metabolite gene clusters. The pattern was complex, with several clusters displaying patchy distribution. Two clusters involved in the synthesis of phytotoxic metabolites are at distinct genomic locations in different Botrytis species. We provide evidence that the clusters for botcinic acid production in B. cinerea and Botrytis sinoallii were acquired by horizontal transfer from taxa within the same genus.

Partial pathogenicity chromosomes in Fusarium oxysporum are sufficient to cause disease and can be horizontally transferred

In Fusarium oxysporum f.sp. lycopersici, all effector genes reported so far - also called SIX genes - are located on a single accessory chromosome which is required for pathogenicity and can also be horizontally transferred to another strain. To narrow down the minimal region required for virulence, we selected partial pathogenicity chromosome deletion strains by fluorescence-assisted cell sorting of a strain in which the two arms of the pathogenicity chromosome were labelled with GFP and RFP respectively. By testing the virulence of these deletion mutants, we show that the complete long arm and part of the short arm of the pathogenicity chromosome are not required for virulence. In addition, we demonstrate that smaller versions of the pathogenicity chromosome can also be transferred to a non-pathogenic strain and they are sufficient to turn the non-pathogen into a pathogen. Surprisingly, originally non-pathogenic strains that had received a smaller version of the pathogenicity chromosome were much more aggressive than recipients with a complete pathogenicity chromosome. Whole genome sequencing analysis revealed that partial deletions of the pathogenicity chromosome occurred mainly close to repeats, and that spontaneous duplication of sequences in accessory regions is frequent both in chromosome deletion strains and in horizontal transfer strains.

What are the Top 10 Unanswered Questions in Molecular Plant-Microbe Interactions?

This article is part of the Top 10 Unanswered Questions in MPMI invited review series.The past few decades have seen major discoveries in the field of molecular plant-microbe interactions. As the result of technological and intellectual advances, we are now able to answer questions at a level of mechanistic detail that we could not have imagined possible 20 years ago. The MPMI Editorial Board felt it was time to take stock and reassess. What big questions remain unanswered? We knew that to identify the fundamental, overarching questions that drive our research, we needed to do this as a community. To reach a diverse audience of people with different backgrounds and perspectives, working in different areas of plant-microbe interactions, we queried the more than 1,400 participants at the 2019 International Congress on Molecular Plant-Microbe Interactions meeting in Glasgow. This group effort resulted in a list of ten, broad-reaching, fundamental questions that influence and inform our research. Here, we introduce these Top 10 unanswered questions, giving context and a brief description of the issues. Each of these questions will be the subject of a detailed review in the coming months. We hope that this process of reflecting on what is known and unknown and identifying the themes that underlie our research will provide a framework to use going forward, giving newcomers a sense of the mystery of the big questions and inspiring new avenues and novel insights.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Whole genome analysis of the koa wilt pathogen (Fusarium oxysporum f. sp. koae) and the development of molecular tools for early detection and monitoring

BACKGROUND

Development and application of DNA-based methods to distinguish highly virulent isolates of Fusarium oxysporum f. sp. koae [Fo koae; cause of koa wilt disease on Acacia koa (koa)] will help disease management through early detection, enhanced monitoring, and improved disease resistance-breeding programs.

RESULTS

This study presents whole genome analyses of one highly virulent Fo koae isolate and one non-pathogenic F. oxysporum (Fo) isolate. These analyses allowed for the identification of putative lineage-specific DNA and predicted genes necessary for disease development on koa. Using putative chromosomes and predicted gene comparisons, Fo koae-exclusive, virulence genes were identified. The putative lineage-specific DNA included identified genes encoding products secreted in xylem (e. g., SIX1 and SIX6) that may be necessary for disease development on koa. Unique genes from Fo koae were used to develop pathogen-specific PCR primers. These diagnostic primers allowed target amplification in the characterized highly virulent Fo koae isolates but did not allow product amplification in low-virulence or non-pathogenic isolates of Fo. Thus, primers developed in this study will be useful for early detection and monitoring of highly virulent strains of Fo koae. Isolate verification is also important for disease resistance-breeding programs that require a diverse set of highly virulent Fo koae isolates for their disease-screening assays to develop disease-resistant koa.

CONCLUSIONS

These results provide the framework for understanding the pathogen genes necessary for koa wilt disease and the genetic variation of Fo koae populations across the Hawaiian Islands.

Structural specificity in plant-filamentous pathogen interactions

Plant diseases bear names such as leaf blights, root rots, sheath blights, tuber scabs, and stem cankers, indicating that symptoms occur preferentially on specific parts of host plants. Accordingly, many plant pathogens are specialized to infect and cause disease in specific tissues and organs. Conversely, others are able to infect a range of tissues, albeit often disease symptoms fluctuate in different organs infected by the same pathogen. The structural specificity of a pathogen defines the degree to which it is reliant on a given tissue, organ, or host developmental stage. It is influenced by both the microbe and the host but the processes shaping it are not well established. Here we review the current status on structural specificity of plant-filamentous pathogen interactions and highlight important research questions. Notably, this review addresses how constitutive defence and induced immunity as well as virulence processes vary across plant organs, tissues, and even cells. A better understanding of the mechanisms underlying structural specificity will aid targeted approaches for plant health, for instance by considering the variation in the nature and the amplitude of defence responses across distinct plant organs and tissues when performing selective breeding.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

PacBio genome sequencing reveals new insights into the genomic organisation of the multi-copy ToxB gene of the wheat fungal pathogen Pyrenophora tritici-repentis

Background

Necrotrophic effector proteins secreted by fungal pathogens are important virulence factors that mediate the development of disease in wheat. Pyrenophora tritici-repentis (Ptr), the causal agent of wheat tan spot, has a race structure dependent on the combination of effectors. In Ptr, ToxA and ToxB are known proteinaceous effectors responsible for necrosis and chlorosis respectively. While Ptr ToxA is encoded by the single gene ToxA, ToxB has multiple loci in the Ptr genome, which is postulated to be directly related to the level of ToxB production and leaf chlorosis. Although previous analysis has indicated that the majority of the ToxB loci lie on a single chromosome, the exact number and chromosomal locations for all the ToxB loci have not been fully identified.

Results

In this study, we have sequenced the genome of a race 5 ToxB-producing isolate (DW5), using PacBio long read technology, and found that ToxB duplications are nested in the complex subtelomeric chromosomal regions. A total of ten identical ToxB gene copies were identified and based on flanking sequence identity, nine loci appeared associated with chromosome 11 and a single copy with chromosome 5. Chromosome 11 multiple ToxB gene loci were separated by large sequence regions between 31 and 66 kb within larger segmental duplications in an alternating pattern related to loci strand, and flanked by transposable elements.

Conclusion

This work provides for the first time the full accompaniment of ToxB loci and surrounding regions, and identifies the organization and distribution of ten ToxB loci to subtelomeric regions. To our knowledge, this is the first report of an interwoven strand-related duplication pattern event. This study further highlights the importance of resolving the highly complex distal chromosomal regions, that remain difficult to assemble, and can harbour important effectors and virulence factors.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12864-020-07029-4.

Accessory Chromosomes in Fusarium oxysporum

Most genomes within the species complex of Fusarium oxysporum are organized into two compartments: the core chromosomes (CCs) and accessory chromosomes (ACs). As opposed to CCs, which are conserved and vertically transmitted to carry out essential housekeeping functions, lineage- or strain-specific ACs are believed to be initially horizontally acquired through unclear mechanisms. These two genomic compartments are different in terms of gene density, the distribution of transposable elements, and epigenetic markers. Although common in eukaryotes, the functional importance of ACs is uniquely emphasized among fungal species, specifically in relationship to fungal pathogenicity and their adaptation to diverse hosts. With a focus on the cross-kingdom fungal pathogen F. oxysporum, this review provides a summary of the differences between CCs and ACs based on current knowledge of gene functions, genome structures, and epigenetic signatures, and explores the transcriptional crosstalk between the core and accessory genomes.

The rise and fall of genes: origins and functions of plant pathogen pangenomes

Plant pathogens can rapidly overcome resistance of their hosts by mutating key pathogenicity genes encoding for effectors. Pathogen adaptation is fuelled by extensive genetic variability in populations and different strains may not share the same set of genes. Recently, such an intra-specific variation in gene content became formalized as pangenomes distinguishing core genes (i.e. shared) and accessory genes (i.e. lineage or strain-specific). Across pathogens species, key effectors tend to be part of the rapidly evolving accessory genome. Here, we show how the construction and analysis of pathogen pangenomes provide deep insights into the dynamic host adaptation process. We also discuss how pangenomes should ideally be built and how geography, niche and lifestyle likely determine pangenome sizes.

Advances in understanding the evolution of fungal genome architecture

Diversity within the fungal kingdom is evident from the wide range of morphologies fungi display as well as the various ecological roles and industrial purposes they serve. Technological advances, particularly in long-read sequencing, coupled with the increasing efficiency and decreasing costs across sequencing platforms have enabled robust characterization of fungal genomes. These sequencing efforts continue to reveal the rampant diversity in fungi at the genome level. Here, we discuss studies that have furthered our understanding of fungal genetic diversity and genomic evolution. These studies revealed the presence of both small-scale and large-scale genomic changes. In fungi, research has recently focused on many small-scale changes, such as how hypermutation and allelic transmission impact genome evolution as well as how and why a few specific genomic regions are more susceptible to rapid evolution than others. High-throughput sequencing of a diverse set of fungal genomes has also illuminated the frequency, mechanisms, and impacts of large-scale changes, which include chromosome structural variation and changes in chromosome number, such as aneuploidy, polyploidy, and the presence of supernumerary chromosomes. The studies discussed herein have provided great insight into how the architecture of the fungal genome varies within species and across the kingdom and how modern fungi may have evolved from the last common fungal ancestor and might also pave the way for understanding how genomic diversity has evolved in all domains of life.