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

Reference-free identification and pangenome analysis of accessory chromosomes in a major fungal plant pathogen

Accessory chromosomes, found in some but not all individuals of a species, play an important role in pathogenicity and host specificity in fungal plant pathogens. However, their variability complicates reference-based analysis, especially when these chromosomes are missing in the reference genome. Pangenome variation graphs offer a reference-free alternative for studying these chromosomes. Here, we constructed a pangenome variation graph for 73 diverse Fusarium oxysporum genomes, a major fungal plant pathogen with a compartmentalized genome that includes conserved core as well as variable accessory chromosomes. To obtain insights into accessory chromosome dynamics, we first constructed a chromosome similarity network using all-vs-all similarity mapping. We identified eleven core chromosomes conserved across all strains and a substantial number of highly variable accessory chromosomes. Some of these accessory chromosomes are host-specific and likely play a role in determining host range. Using a k-mer based approach, we further identified the presence of these accessory chromosomes in all available (581) F. oxysporum assemblies and corroborated the occurrence of host-specific accessory chromosomes. To further analyze the evolution of chromosomes in F. oxysporum, we constructed a pangenome variation graph per group of homologous chromosomes. This reveals that accessory chromosomes are composed of different stretches of accessory regions, and possibly rearrangements between accessory regions gave rise to these mosaic accessory chromosomes. Furthermore, we show that accessory chromosomes are likely horizontally transferred in natural populations. Our findings demonstrate that a pangenome variation graph is a powerful approach to elucidate the evolutionary dynamics of accessory chromosomes in F. oxysporum, which is not only a useful resource for Fusarium but also provides a framework for similar analyses in other species containing accessory chromosomes.

Comparative genomic analysis of Fusarium oxysporum f. sp. lycopersici reveals telomeric duplications of a lineage-specific region carrying SIX8 and PSL1 and genome-wide expansion of Foxy transposable elements

Fusarium oxysporum f. sp. lycopersici (Fol), the causal agent of tomato wilt disease, is a soil-borne, vascular-colonizing fungal pathogen that severely impacts tomato production in most growing regions worldwide. Despite the availability of over thirty Fol genome sequences in public databases, only one chromosome-scale assembly exists, comprising the low sequence-coverage Fol4287 reference genome generated using Sanger sequencing. Thus, genome structural variation and comparative genomics analyses of Fol remain largely unexplored. Here, using third generation Nanopore long-read sequencing in combination with high-throughput chromosome conformation capture (Hi-C) data, we have independently constructed a high-quality assembly and annotation of the Fol007 race 2 genome that consists of 15 pseudochromosomes with 25 telomeres, including 10 telomere-to-telomere chromosomes. The Fol007 genome is predicted to contain 29,148 protein-encoding genes, including all known SIX genes except for SIX4, which is absent in Fol race 2. Compared to the genome of Fusarium verticillioides, the Fol007 genome contains four complete lineage-specific (LS) chromosomes, including chromosome 5 (Chr5), chromosome 13 (Chr13), and two small chromosomes, Chr14 and Chr15. Comparative genomic analysis between the newly assembled Fol007 genome and Fol4287_2010 deposited in the GenBank database reveals that the completely assembled Chr13 of Fol007 carrying all known SIX genes (except SIX4 and SIX8) and three novel candidate effector genes is relatively stable and corresponds to pathogenicity chromosome Chr14 in Fol4287_2010. It also reveals a telomeric duplication of an LS region carrying SIX8 and PSL1 on core chromosomes Chr6, Chr7 and Chr11, and LS chromosome Chr15, as well as the absence of segmental duplications on LS chromosomes of the Fol007 genome. Furthermore, compared to the genomes of Fol race 1 and non-pathogenic Fo strains, an active and specific Foxy transposable element, responsible for the inactivation of a second copy of SIX13, was identified and found to have expanded in the genomes of Fol race 2 and 3 strains. This element may contribute to Fusarium oxysporum genome evolution and has potential as a genetic marker for studying phylogenetic relationships among formae speciales of Fusarium oxysporum. These findings provide a basis for further genetic and genomic understanding of Fol evolution and virulence mechanisms employed by Fol and contribute another reference genome for Fusarium study more broadly.

Frequent genetic exchanges revealed by a pan-mitogenome graph of a fungal plant pathogen

Mitochondria are present in almost all eukaryotic lineages. The mitochondrial genomes (mitogenomes) evolve separately from nuclear genomes, and they can therefore provide relevant insights into the evolution of their host species. Fusarium oxysporum is a major fungal plant pathogen that is assumed to reproduce clonally. However, horizontal chromosome transfer between strains can occur through heterokaryon formation, and recently, signs of sexual recombination have been observed. Similarly, signs of recombination in F. oxysporum mitogenomes challenged the prevailing assumption of clonal reproduction in this species. Here, we construct, to our knowledge, the first fungal pan-mitogenome graph of nearly 500 F. oxysporum mitogenome assemblies to uncover the variation and evolution. In general, the gene order of fungal mitogenomes is not well conserved, yet the mitogenome of F. oxysporum and related species are highly colinear. We observed two strikingly contrasting regions in the F. oxysporum pan-mitogenome, comprising a highly conserved core mitogenome and a long variable region (6-16 kb in size), of which we identified three distinct types. The pan-mitogenome graph reveals that only five intron insertions occurred in the core mitogenome and that the long variable regions drive the difference between mitogenomes. Moreover, we observed that their evolution is neither concurrent with the core mitogenome nor with the nuclear genome. Our large-scale analysis of long variable regions uncovers frequent recombination between mitogenomes, even between strains that belong to different taxonomic clades. This challenges the common assumption of incompatibility between genetically diverse F. oxysporum strains and provides new insights into the evolution of this fungal species.IMPORTANCEInsights into plant pathogen evolution is essential for the understanding and management of disease. Fusarium oxysporum is a major fungal pathogen that can infect many economically important crops. Pathogenicity can be transferred between strains by the horizontal transfer of pathogenicity chromosomes. The fungus has been thought to evolve clonally, yet recent evidence suggests active sexual recombination between related isolates, which could at least partially explain the horizontal transfer of pathogenicity chromosomes. By constructing a pan-genome graph of nearly 500 mitochondrial genomes, we describe the genetic variation of mitochondria in unprecedented detail and demonstrate frequent mitochondrial recombination. Importantly, recombination can occur between genetically diverse isolates from distinct taxonomic clades and thus can shed light on genetic exchange between fungal strains.

Ectopic expression of a truncated NLR gene from wild Arachis enhances resistance to Fusarium oxysporum

Fusarium oxysporum causes devastating vascular wilt diseases in numerous crop species, resulting in substantial yield losses. The Arabidopsis thaliana-F. oxysporum f.sp. conglutinans (FOC) model system enables the identification of meaningful genotype-phenotype correlations and was applied in this study to evaluate the effects of overexpressing an NLR gene (AsTIR19) from Arachis stenosperma against pathogen infection. AsTIR19 overexpression (OE) lines exhibited enhanced resistance to FOC without any discernible phenotype penalties. To elucidate the underlying resistance mechanisms mediated by AsTIR19 overexpression, we conducted whole transcriptome sequencing of an AsTIR19-OE line and non-transgenic wild-type (WT) plants inoculated and non-inoculated with FOC using Illumina HiSeq4000. Comparative analysis revealed 778 differentially expressed genes (DEGs) attributed to transgene overexpression, while fungal inoculation induced 434 DEGs in the OE line, with many falling into defense-related Gene Ontology (GO) categories. GO and KEGG enrichment analysis showed that DEGs were enriched in the phenylpropanoid and flavonoid pathways in the OE plants. This comprehensive transcriptomic analysis underscores how AsTIR19 overexpression reprograms transcriptional networks, modulating the expression of stress-responsive genes across diverse metabolic pathways. These findings provide valuable insights into the molecular mechanisms underlying the role of this NLR gene under stress conditions, highlighting its potential to enhance resistance to Fusarium oxysporum.

Comparative genomics and transcriptomics reveal differences in effector complement and expression between races of Fusarium oxysporum f.sp. lactucae

This study presents the first genome and transcriptome analyses for Fusarium oxysporum f. sp. lactucae (Fola) which causes Fusarium wilt disease of lettuce. Long-read genome sequencing of three race 1 (Fola1) and three race 4 (Fola4) isolates revealed key differences in putative effector complement between races and with other F. oxysporum ff. spp. following mimp-based bioinformatic analyses. Notably, homologues of Secreted in Xylem (SIX) genes, also present in many other F. oxysporum ff. spp, were identified in Fola, with both SIX9 and SIX14 (multiple copies with sequence variants) present in both Fola1 and Fola4. All Fola4 isolates also contained an additional single copy of SIX8. RNAseq of lettuce following infection with Fola1 and Fola4 isolates identified highly expressed effectors, some of which were homologues of those reported in other F. oxysporum ff. spp. including several in F. oxysporum f. sp. apii. Although SIX8, SIX9 and SIX14 were all highly expressed in Fola4, of the two SIX genes present in Fola1, only SIX9 was expressed as further analysis revealed that SIX14 gene copies were disrupted by insertion of a transposable element. Two variants of Fola4 were also identified based on different genome and effector-based analyses. This included two different SIX8 sequence variants which were divergently transcribed from a shared promoter with either PSE1 or PSL1 respectively. In addition, there was evidence of two independent instances of HCT in the different Fola4 variants. The involvement of helitrons in Fola genome rearrangement and gene expression is discussed.

Evaluation of the Genetic Diversity, Haplotype, and Virulence of Fusarium oxysporum f. sp. vasinfectum Field Isolates from Alabama

The United States is the third largest producer of cotton and the largest exporter of cotton globally. Fusarium wilt, caused by the soilborne fungal pathogen Fusarium oxysporum f. sp. vasinfectum (Fov), was estimated to cause a $21 million cotton yield loss in 2022. Historically, Alabama was an important producer of cotton in the Southeastern United States and was the first state in which Fusarium wilt on cotton was described. To assess the genetic diversity of Fov field isolates in Alabama, 118 field isolates were collected from six counties across the state from 2014 to 2016. Phylogenetic analysis using TEF1 and RPB2 placed the Fov field isolates into 18 haplotypes. Upon profiling the Tfo1 transposon insertion in the NAT gene, it was determined that no race 4 isolates were recovered in Alabama. Representatives of all field isolate haplotypes caused disease on Upland cotton variety Rowden in a hydroponic test tube assay. Two haplotype A isolates were the most aggressive isolates recovered, and haplotype A isolate TF1 was more aggressive than the race 4 isolate 89-1A on Upland cotton and had similar symptom severity on Pima cotton. Karyotype profiling indicted an abundance of small chromosomes characteristic of karyotypes that include accessory chromosomes, with considerable variability between isolates. Collectively, our study indicates that Fov isolates from Alabama are genetically diverse, which may have been promoted by its persistence in cotton fields.

Pangenomics of flax fungal parasite Fusarium oxysporum f. sp. lini

To assess the genomic diversity of Fusarium oxysporum f. sp. lini strains and compile a comprehensive gene repertoire, we constructed a pangenome using 13 isolates from four different clonal lineages, each exhibiting distinct levels of virulence. Syntenic analyses of two selected genomes revealed significant chromosomal rearrangements unique to each genome. A comprehensive examination of both core and accessory pangenome content and diversity points at an open genome state. Additionally, Gene Ontology (GO) enrichment analysis indicated that non-core pangenome genes are associated with pathogen recognition and immune signaling. Furthermore, the Folini pansecterome, encompassing secreted proteins critical for fungal pathogenicity, primarily consists of three functional classes: effector proteins, CAZYmes, and proteases. These three classes account for approximately 3.5% of the pangenome. Each functional class within the pansecterome was meticulously annotated and characterized with respect to pangenome category distribution, PFAM domain frequency, and strain virulence assessment. This analysis revealed that highly virulent isolates have specific types of PFAM domains that are exclusive to them. Upon examining the repertoire of SIX genes known for virulence in other formae speciales, it was found that all isolates had a similar gene content except for two, which lacked SIX genes entirely.

In silico identification of a promising inhibitor of Fusarium oxysporum f. sp. Lycopersici, Secreted in Xylem 1 protein

Fusarium oxysporum f. sp. Lycopersici (FOL) is a soilborne pathogen that infects tomato plants and inflicts severe damage, resulting in heavy yield losses worldwide, causing Fusarium wilt disease. FOL encodes several pathogenicity factors necessary for colonizing and invading the host plants. Secreted in Xylem (SIX), a pathogenicity factor, is a small cysteine-rich fungal protein found in the xylem sap of FOL-infected tomato plants, which plays a major role in determining host specificity and in contributing to pathogenicity/virulence. However, the structure of SIX1 has not been modeled yet. Therefore, this study aimed to elucidate the structure of SIX1 by comparative modeling using Robetta server. The best possible structures obtained were then refined, validated, and utilized for subsequent analysis. An antifungal library comprising 16,824 compounds was screened to determine small molecules that can interact with SIX1. Five antifungal compounds were identified from the library. Further analyses revealed that, of the five ligands, 4-[(2-(3-methoxyphenoxy)acetyl)amino] benzamide exhibited the capacity to stably interact with SIX1. This shows that 4-[[2-(3-methoxyphenoxy)acetyl]amino] benzamide can be used as a potential candidate in the prevention of FOL infection. In summary, small-molecule inhibitors such as 4-[[2-(3-methoxyphenoxy)acetyl]amino] benzamide could be highly effective in combating FOL infection, along with biocontrol methods and strategies that use transgenic plants overexpressing resistance genes.

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

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

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.

The structural repertoire of Fusarium oxysporum f. sp. lycopersici effectors revealed by experimental and computational studies

Plant pathogens secrete proteins, known as effectors, that function in the apoplast or inside plant cells to promote virulence. Effector recognition by cell-surface or cytosolic receptors results in the activation of defence pathways and plant immunity. Despite their importance, our general understanding of fungal effector function and recognition by immunity receptors remains poor. One complication often associated with effectors is their high sequence diversity and lack of identifiable sequence motifs precluding prediction of structure or function. In recent years, several studies have demonstrated that fungal effectors can be grouped into structural classes, despite significant sequence variation and existence across taxonomic groups. Using protein X-ray crystallography, we identify a new structural class of effectors hidden within the secreted in xylem (SIX) effectors from Fusarium oxysporum f. sp. lycopersici (Fol). The recognised effectors Avr1 (SIX4) and Avr3 (SIX1) represent the founding members of the Fol dual-domain (FOLD) effector class, with members containing two distinct domains. Using AlphaFold2, we predicted the full SIX effector repertoire of Fol and show that SIX6 and SIX13 are also FOLD effectors, which we validated experimentally for SIX6. Based on structural prediction and comparisons, we show that FOLD effectors are present within three divisions of fungi and are expanded in pathogens and symbionts. Further structural comparisons demonstrate that Fol secretes effectors that adopt a limited number of structural folds during infection of tomato. This analysis also revealed a structural relationship between transcriptionally co-regulated effector pairs. We make use of the Avr1 structure to understand its recognition by the I receptor, which leads to disease resistance in tomato. This study represents an important advance in our understanding of Fol-tomato, and by extension plant-fungal interactions, which will assist in the development of novel control and engineering strategies to combat plant pathogens.

Pathogenicity chromosome of Fusarium oxysporum f. sp. cepae

Fusarium oxysporum f. sp. cepae (Foc) is the causative agent of Fusarium basal rot disease in onions, which leads to catastrophic global crop production losses. Therefore, the interaction of Foc with its host has been actively investigated, and the pathogen-specific (PS) regions of the British strain Foc_FUS2 have been identified. However, it has not been experimentally determined whether the identified PS region plays a role in pathogenicity. To identify the pathogenicity chromosome in the Japanese strain Foc_TA, we initially screened effector candidates, defined as small proteins with a signal peptide that contain two or more cysteines, from genome sequence data. Twenty-one candidate effectors were identified, five of which were expressed during infection. Of the expressed effector candidates, four were located on the 4-Mb-sized chromosome in Foc_TA. To clarify the relationship between pathogenicity and the 4-Mb-sized chromosome in Foc_TA, nine putative 4-Mb-sized chromosome loss strains were generated by treatment with benomyl (a mitotic inhibitor drug). A pathogenicity test with putative 4-Mb-sized chromosome loss strains showed that these strains were impaired in their pathogenicity toward onions. Genome analysis of three putative 4-Mb-sized chromosome loss strains revealed that two strains lost a 4-Mb-sized chromosome in common, and another strain maintained a 0.9-Mb region of the 4-Mb-sized chromosome. Our findings show that the 4-Mb-sized chromosome is the pathogenicity chromosome in Foc_TA, and the 3.1-Mb region within the 4-Mb-sized chromosome is required for full pathogenicity toward onion.

Real-Time PCR Assays for Races of the Spinach Fusarium Wilt Pathogen, Fusarium oxysporum f. sp. spinaciae

Fusarium wilt of spinach, caused by Fusarium oxysporum f. sp. spinaciae, is a significant limitation for producers of vegetative spinach and spinach seed crops during warm temperatures and/or on acid soils. Identification of isolates of F. oxysporum f. sp. spinaciae, and distinction of isolates of the two known races, entails time-intensive pathogenicity tests. In this study, two real-time PCR assays were developed: one for a candidate effector gene common to both races of F. oxysporum f. sp. spinaciae, and another for a candidate effector gene unique to isolates of race 2. The assays were specific to isolates of F. oxysporum f. sp. spinaciae (n = 44) and isolates of race 2 (n = 23), respectively. Neither assay amplified DNA from 10 avirulent isolates of F. oxysporum associated with spinach, 57 isolates of other formae speciales and Fusarium spp., or 7 isolates of other spinach pathogens. When the assays were used to detect DNA extracted from spinach plants infected with an isolate of race 1, race 2, or a 1:1 mixture of both races, the amount of target DNA detected increased with increasing severity of wilt. Plants infected with one or both isolates could be distinguished based on the ratio in copy number for each target locus. The real-time PCR assays enable rapid diagnosis of Fusarium wilt of spinach and will facilitate research on the epidemiology and management of this disease, as well as surveys on the prevalence of this understudied pathogen in regions of spinach and/or spinach seed production.

Identifying and tracking mobile elements in evolving compost communities yields insights into the nanobiome

Microbial evolution is driven by rapid changes in gene content mediated by horizontal gene transfer (HGT). While mobile genetic elements (MGEs) are important drivers of gene flux, the nanobiome-the zoo of Darwinian replicators that depend on microbial hosts-remains poorly characterised. New approaches are necessary to increase our understanding beyond MGEs shaping individual populations, towards their impacts on complex microbial communities. A bioinformatic pipeline (xenoseq) was developed to cross-compare metagenomic samples from microbial consortia evolving in parallel, aimed at identifying MGE dissemination, which was applied to compost communities which underwent periodic mixing of MGEs. We show that xenoseq can distinguish movement of MGEs from demographic changes in community composition that otherwise confounds identification, and furthermore demonstrate the discovery of various unexpected entities. Of particular interest was a nanobacterium of the candidate phylum radiation (CPR) which is closely related to a species identified in groundwater ecosystems (Candidatus Saccharibacterium), and appears to have a parasitic lifestyle. We also highlight another prolific mobile element, a 313 kb plasmid hosted by a Cellvibrio lineage. The host was predicted to be capable of nitrogen fixation, and acquisition of the plasmid coincides with increased ammonia production. Taken together, our data show that new experimental strategies combined with bioinformatic analyses of metagenomic data stand to provide insight into the nanobiome as a driver of microbial community evolution.

Major proliferation of transposable elements shaped the genome of the soybean rust pathogen Phakopsora pachyrhizi

With >7000 species the order of rust fungi has a disproportionately large impact on agriculture, horticulture, forestry and foreign ecosystems. The infectious spores are typically dikaryotic, a feature unique to fungi in which two haploid nuclei reside in the same cell. A key example is Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease, one of the world's most economically damaging agricultural diseases. Despite P. pachyrhizi's impact, the exceptional size and complexity of its genome prevented generation of an accurate genome assembly. Here, we sequence three independent P. pachyrhizi genomes and uncover a genome up to 1.25 Gb comprising two haplotypes with a transposable element (TE) content of ~93%. We study the incursion and dominant impact of these TEs on the genome and show how they have a key impact on various processes such as host range adaptation, stress responses and genetic plasticity.

The Botrytis cinerea Gene Expression Browser

For comprehensive gene expression analyses of the phytopathogenic fungus Botrytis cinerea, which infects a number of plant taxa and is a cause of substantial agricultural losses worldwide, we developed BEB, a web-based B. cinerea gene Expression Browser. This computationally inexpensive web-based application and its associated database contain manually curated RNA-Seq data for B. cinerea. BEB enables expression analyses of genes of interest under different culture conditions by providing publication-ready heatmaps depicting transcript levels, without requiring advanced computational skills. BEB also provides details of each experiment and user-defined gene expression clustering and visualization options. If needed, tables of gene expression values can be downloaded for further exploration, including, for instance, the determination of differentially expressed genes. The BEB implementation is based on open-source computational technologies that can be deployed for other organisms. In this case, the new implementation will be limited only by the number of transcriptomic experiments that are incorporated into the platform. To demonstrate the usability and value of BEB, we analyzed gene expression patterns across different conditions, with a focus on secondary metabolite gene clusters, chromosome-wide gene expression, previously described virulence factors, and reference genes, providing the first comprehensive expression overview of these groups of genes in this relevant fungal phytopathogen. We expect this tool to be broadly useful in B. cinerea research, providing a basis for comparative transcriptomics and candidate gene identification for functional assays.

The contribution of DNA repair pathways to genome editing and evolution in filamentous pathogens

DNA double-strand breaks require repair or risk corrupting the language of life. To ensure genome integrity and viability, multiple DNA double-strand break repair pathways function in eukaryotes. Two such repair pathways, canonical non-homologous end joining and homologous recombination, have been extensively studied, while other pathways such as microhomology-mediated end joint and single-strand annealing, once thought to serve as back-ups, now appear to play a fundamental role in DNA repair. Here, we review the molecular details and hierarchy of these four DNA repair pathways, and where possible, a comparison for what is known between animal and fungal models. We address the factors contributing to break repair pathway choice, and aim to explore our understanding and knowledge gaps regarding mechanisms and regulation in filamentous pathogens. We additionally discuss how DNA double-strand break repair pathways influence genome engineering results, including unexpected mutation outcomes. Finally, we review the concept of biased genome evolution in filamentous pathogens, and provide a model, termed Biased Variation, that links DNA double-strand break repair pathways with properties of genome evolution. Despite our extensive knowledge for this universal process, there remain many unanswered questions, for which the answers may improve genome engineering and our understanding of genome evolution.

Fusarium oxysporum effector clustering version 2: An updated pipeline to infer host range

The fungus Fusarium oxysporum is infamous for its devastating effects on economically important crops worldwide. F. oxysporum isolates are grouped into formae speciales based on their ability to cause disease on different hosts. Assigning F. oxysporum strains to formae speciales using non-experimental procedures has proven to be challenging due to their genetic heterogeneity and polyphyletic nature. However, genetically diverse isolates of the same forma specialis encode similar repertoires of effectors, proteins that are secreted by the fungus and contribute to the establishment of compatibility with the host. Based on this observation, we previously designed the F. oxysporum Effector Clustering (FoEC) pipeline which is able to classify F. oxysporum strains by forma specialis based on hierarchical clustering of the presence of predicted putative effector sequences, solely using genome assemblies as input. Here we present the updated FoEC2 pipeline which is more user friendly, customizable and, due to multithreading, has improved scalability. It is designed as a Snakemake pipeline and incorporates a new interactive visualization app. We showcase FoEC2 by clustering 537 publicly available F. oxysporum genomes and further analysis of putative effector families as multiple sequence alignments. We confirm classification of isolates into formae speciales and are able to further identify their subtypes. The pipeline is available on github: https://github.com/pvdam3/FoEC2.

Pathotypes of Fusarium oxysporum f. sp. fragariae express discrete repertoires of accessory genes and induce distinct host transcriptional responses during root infection

Convergent evolution of phytopathogenicity is poorly described, especially among multiple strains of a single microbial species. We investigated this phenomenon with genetically diverse isolates of Fusarium oxysporum f. sp. fragariae (Fof) that cause one of two syndromes: chlorosis and wilting (the 'yellows-fragariae' pathotype), or only wilting (the 'wilt-fragariae' pathotype). We challenged strawberry (Fragaria × ananassa) plants to root infection by five fungal isolates: three yellows-fragariae, one wilt-fragariae and one that is not pathogenic to strawberry. All Fof isolates had chromosome-level assemblies; three were newly generated. The two pathotypes triggered distinct host responses, especially among phytohormone-associated genes; yellows-fragariae isolates strongly induced jasmonic acid-associated genes, whereas the wilt-fragariae isolate primarily induced ethylene biosynthesis and signalling. The differentially expressed genes on fungal accessory chromosomes were almost entirely distinct between pathotypes. We identified an ~150 kbp 'pathogenicity island' that was horizontally transferred between wilt-fragariae strains. This predicted pathogenicity island was enriched with differentially expressed genes whose predicted functions were related to plant infection, and only one of these genes was also upregulated in planta by yellows-fragariae isolates. These results support the conclusion that wilt- and yellows-fragariae cause physiologically distinct syndromes by the expression of discrete repertoires of genes on accessory chromosomes.

Genome editing using a versatile vector-based CRISPR/Cas9 system in Fusarium species

Fusarium species include important filamentous fungal pathogens that can infect plants, animals, and humans. Meanwhile, some nonpathogenic Fusarium species are promising biocontrol agents against plant pathogens. Here, we developed a genome editing technology using a vector-based CRISPR/Cas9 system for Fusarium oxysporum f. sp. lycopersici (Fol). This optimized CRISPR/Cas9 system, harboring an endogenous U6 small nuclear RNA promoter for the expression of single-guide RNA and an endogenous H2B nuclear localization signal for the localization of Cas9, enabled efficient targeted gene knock-out, including in the accessory chromosomal regions in Fol. We further demonstrated single crossover-mediated targeted base editing and endogenous gene tagging. This system was also applicable for genome editing in F. oxysporum f. sp. spinaciae and F. commune without any modifications, suggesting that this CRISPR/Cas9 vector has a potential application for a broad range of researches on other Fusarium species.

Transcriptome Analysis of Fusarium–Tomato Interaction Based on an Updated Genome Annotation of Fusarium oxysporum f. sp. lycopersici Identifies Novel Effector Candidates That Suppress or Induce Cell Death in Nicotiana benthamiana

Fusarium oxysporum f. sp. lycopersici (Fol) causes vascular wilt disease in tomato. Upon colonization of the host, Fol secretes many small effector proteins into the xylem sap to facilitate infection. Besides known SIX (secreted in xylem) proteins, the identity of additional effectors that contribute to Fol pathogenicity remains largely unexplored. We performed a deep RNA-sequencing analysis of Fol race 2-infected tomato, used the sequence data to annotate a published genome assembly generated via PacBio SMRT sequencing of the Fol race 2 reference strain Fol4287, and analysed the resulting transcriptome to identify Fol effector candidates among the newly annotated genes. We examined the Fol-infection expression profiles of all 13 SIX genes present in Fol race 2 and identified 27 new candidate effector genes that were likewise significantly upregulated upon Fol infection. Using Agrobacterium-mediated transformation, we tested the ability of 22 of the new candidate effector genes to suppress or induce cell death in leaves of Nicotiana benthamiana. One effector candidate designated Fol-EC19, encoding a secreted guanyl-specific ribonuclease, was found to trigger cell death and two effector candidates designated Fol-EC14 and Fol-EC20, encoding a glucanase and a secreted trypsin, respectively, were identified that can suppress Bax-mediated cell death. Remarkably, Fol-EC14 and Fol-EC20 were also found to suppress I-2/Avr2- and I/Avr1-mediated cell death. Using the yeast secretion trap screening system, we showed that these three biologically-active effector candidates each contain a functional signal peptide for protein secretion. Our findings provide a basis for further understanding the virulence functions of Fol effectors.

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.

Breeding for Resistance to Fusarium Wilt of Tomato: A Review

For over a century, breeders have worked to develop tomato (Solanum lycopersicum) cultivars with resistance to Fusarium wilt (Fol) caused by the soilborne fungus Fusarium oxysporum f. sp. lycopersici. Host resistance is the most effective strategy for the management of this disease. For each of the three Fol races, resistance has been introgressed from wild tomato species, predominately in the form of R genes. The I, I-2, I-3, and I-7 R genes have each been identified, as well as the corresponding Avr effectors in the fungus with the exception of Avr7. The mechanisms by which the R gene protein products recognize these effectors, however, has not been elucidated. Extensive genetic mapping, gene cloning, and genome sequencing efforts support the development of tightly-linked molecular markers, which greatly expedite tomato breeding and the development of elite, Fol resistant cultivars. These resources also provide important tools for pyramiding resistance genes and should support the durability of host resistance.

Homeobox Transcription Factors Are Required for Fungal Development and the Suppression of Host Defense Mechanisms in the Colletotrichum scovillei-Pepper Pathosystem

Colletotrichum scovillei, an ascomycete phytopathogenic fungus, is the main causal agent of serious yield losses of economic crops worldwide. The fungus causes anthracnose disease on several fruits, including peppers. However, little is known regarding the underlying molecular mechanisms involved in the development of anthracnose caused by this fungus. In an initial step toward understanding the development of anthracnose on pepper fruits, we retrieved 624 transcription factors (TFs) from the whole genome of C. scovillei and comparatively analyzed the entire repertoire of TFs among phytopathogenic fungi. Evolution and proliferation of members of the homeobox-like superfamily, including homeobox (HOX) TFs that regulate the development of eukaryotic organisms, were demonstrated in the genus Colletotrichum. C. scovillei was found to contain 10 HOX TF genes (CsHOX1 to CsHOX10), which were functionally characterized using deletion mutants of each CsHOX gene. Notably, CsHOX1 was identified as a pathogenicity factor required for the suppression of host defense mechanisms, which represents a new role for HOX TFs in pathogenic fungi. CsHOX2 and CsHOX7 were found to play essential roles in conidiation and appressorium development, respectively, in a stage-specific manner in C. scovillei. Our study provides a molecular basis for understanding the mechanisms associated with the development of anthracnose on fruits caused by C. scovillei, which will aid in the development of novel approaches for disease management. IMPORTANCE The ascomycete phytopathogenic fungus, Colletotrichum scovillei, causes serious yield loss on peppers. However, little is known about molecular mechanisms involved in the development of anthracnose caused by this fungus. We analyzed whole-genome sequences of C. scovillei and isolated 624 putative TFs, revealing the existence of 10 homeobox (HOX) transcription factor (TF) genes. We found that CsHOX1 is a pathogenicity factor required for the suppression of host defense mechanism, which represents a new role for HOX TFs in pathogenic fungi. We also found that CsHOX2 and CsHOX7 play essential roles in conidiation and appressorium development, respectively, in a stage-specific manner in C. scovillei. Our study contributes to understanding the mechanisms associated with the development of anthracnose on fruits caused by C. scovillei, which will aid for initiating novel approaches for disease management.

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.

A Genomic Blueprint of Flax Fungal Parasite Fusarium oxysporum f. sp. lini

Fusarium wilt of flax is an aggressive disease caused by the soil-borne fungal pathogen Fusarium oxysporum f. sp. lini. It is a challenging pathogen presenting a constant threat to flax production industry worldwide. Previously, we reported chromosome-level assemblies of 5 highly pathogenic F. oxysporum f. sp. lini strains. We sought to characterize the genomic architecture of the fungus and outline evolutionary mechanisms shaping the pathogen genome. Here, we reveal the complex multi-compartmentalized genome organization and uncover its diverse evolutionary dynamics, which boosts genetic diversity and facilitates host adaptation. In addition, our results suggest that host of functions implicated in the life cycle of mobile genetic elements are main contributors to dissimilarity between proteomes of different Fusaria. Finally, our experiments demonstrate that mobile genetics elements are expressed in planta upon infection, alluding to their role in pathogenicity. On the whole, these results pave the way for further in-depth studies of evolutionary forces shaping the host–pathogen interaction.

Putative Novel Effector Genes Revealed by the Genomic Analysis of the Phytopathogenic Fungus Fusarium oxysporum f. sp. physali (Foph) That Infects Cape Gooseberry Plants

The vascular wilt disease caused by the fungus Fusarium oxysporum f. sp. physali (Foph) is one of the most limiting factors for the production and export of cape gooseberry (Physalis peruviana) in Colombia. A transcriptomic analysis of a highly virulent strain of F. oxysporum in cape gooseberry plants, revealed the presence of secreted in the xylem (SIX) effector genes, known to be involved in the pathogenicity of other formae speciales (ff. spp.) of F. oxysporum. This pathogenic strain was classified as a new f. sp. named Foph, due to its specificity for cape gooseberry hosts. Here, we sequenced and assembled the genome of five strains of F. oxysporum from a fungal collection associated to the cape gooseberry crop (including Foph), focusing on the validation of the presence of SIX homologous and on the identification of putative effectors unique to Foph. By comparative and phylogenomic analyses based on single-copy orthologous, we found that Foph is closely related to F. oxysporum ff. spp., associated with solanaceous hosts. We confirmed the presence of highly identical homologous genomic regions between Foph and Fol that contain effector genes and identified six new putative effector genes, specific to Foph pathogenic strains. We also conducted a molecular characterization using this set of putative novel effectors in a panel of 36 additional stains of F. oxysporum including two of the four sequenced strains, from the fungal collection mentioned above. These results suggest the polyphyletic origin of Foph and the putative independent acquisition of new candidate effectors in different clades of related strains. The novel effector candidates identified in this genomic analysis, represent new sources involved in the interaction between Foph and cape gooseberry, that could be implemented to develop appropriate management strategies of the wilt disease caused by Foph in the cape gooseberry crop.