Genetic diversity of the banana Fusarium wilt pathogen in Cuba and across Latin America and the Caribbean.

Fusarium wilt of bananas (FWB) is a severe plant disease that leads to substantial losses in banana production worldwide. It remains a major concern for Cuban banana cultivation. The disease is caused by members of the soil-borne Fusarium oxysporum species complex. However, the genetic diversity among Fusarium species infecting bananas in Cuba has remained largely unexplored. In our comprehensive survey, we examined symptomatic banana plants across all production zones in the country, collecting 170 Fusarium isolates. Leveraging genotyping-by-sequencing and whole-genome comparisons, we investigated the genetic diversity within these isolates and compared it with a global Fusarium panel. Notably, typical FWB symptoms were observed in Bluggoe cooking bananas and Pisang Awak subgroups across 14 provinces. Our phylogenetic analysis revealed that F. purpurascens, F. phialophorum, and F. tardichlamydosporum are responsible for FWB in Cuba, with F. tardichlamydosporum dominating the population. Furthermore, we identified between five and seven distinct genetic clusters, with F. tardichlamydosporum isolates forming at least two subgroups. This finding underscores the high genetic diversity of Fusarium spp. contributing to FWB in the Americas. Our study sheds light on the population genetic structure and diversity of the FWB pathogen in Cuba and the broader Latin American and Caribbean regions.

Receptor-like cytoplasmic kinases of different subfamilies differentially regulate SOBIR1/BAK1-mediated immune responses in Nicotiana benthamiana.

Cell-surface receptors form the front line of plant immunity. The leucine-rich repeat (LRR)-receptor-like kinases SOBIR1 and BAK1 are required for the functionality of the tomato LRR-receptor-like protein Cf-4, which detects the secreted effector Avr4 of the pathogenic fungus Fulvia fulva. Here, we show that the kinase domains of SOBIR1 and BAK1 directly phosphorylate each other and that residues Thr522 and Tyr469 of the kinase domain of Nicotiana benthamiana SOBIR1 are required for its kinase activity and for interacting with signalling partners, respectively. By knocking out multiple genes belonging to different receptor-like cytoplasmic kinase (RLCK)-VII subfamilies in N. benthamiana:Cf-4, we show that members of RLCK-VII-6, -7, and -8 differentially regulate the Avr4/Cf-4-triggered biphasic burst of reactive oxygen species. In addition, members of RLCK-VII-7 play an essential role in resistance against the oomycete pathogen Phytophthora palmivora. Our study provides molecular evidence for the specific roles of RLCKs downstream of SOBIR1/BAK1-containing immune complexes.

Oomycete metabolism is highly dynamic and reflects lifestyle adaptations.

The selective pressure of pathogen-host symbiosis drives adaptations. How these interactions shape the metabolism of pathogens is largely unknown. Here, we use comparative genomics to systematically analyse the metabolic networks of oomycetes, a diverse group of eukaryotes that includes saprotrophs as well as pathogens of animal- and plant pathogens, the latter causing devastating diseases with significant economic and/or ecological impact. In our analyses of 44 oomycete species, we uncover considerable variation in metabolism that can be linked to lifestyle differences. Comparisons of metabolic gene content reveal that plant pathogenic oomycetes have a bipartite metabolism consisting of a conserved core and an accessory set. The accessory set can be associated with the degradation of defence compounds produced by plants when challenged by pathogens. Obligate biotrophic oomycetes have smaller metabolic networks, and taxonomically distantly related biotrophic lineages display convergent evolution by repeated gene losses in both the conserved as well as the accessory set of metabolism. When investigating to what extent the metabolic networks in obligate biotrophs differ from those in hemibiotrophic plant pathogens, we observe that the losses of metabolic enzymes in obligate biotrophs are not random and that gene losses predominantly influence the terminal branches of the metabolic networks. Our analyses represent the first metabolism-focused comparison of oomycetes at this scale and will contribute to a better understanding of the evolution of oomycete metabolism in relation to lifestyle adaptation.

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.

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

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

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

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

A chromosome-level genome assembly of Zasmidium syzygii isolated from banana leaves.

Accurate taxonomic classification of samples from infected host material is essential for disease diagnostics and genome analyses. Despite the importance, diagnosis of fungal pathogens causing banana leaf diseases remains challenging. Foliar diseases of bananas are mainly caused by 3 Pseudocercospora species, of which the most predominant causal agent is Pseudocercospora fijiensis. Here, we sequenced and assembled four fungal isolates obtained from necrotic banana leaves in Bohol (Philippines) and obtained a high-quality genome assembly for one of these isolates. The samples were initially identified as P. fijiensis using PCR diagnostics; however, the assembly size was consistently 30 Mb smaller than expected. Based on the internal transcribed spacer (ITS) sequences, we identified the samples as Zasmidium syzygii (98.7% identity). The high-quality Zasmidium syzygii assembly is 42.5 Mb in size, comprising 16 contigs, of which 11 are most likely complete chromosomes. The genome contains 98.6% of the expected single-copy BUSCO genes and contains 14,789 genes and 10.3% repeats. The 3 short-read assemblies are less continuous but have similar genome sizes (40.4-42.4 Mb) and contain between 96.5 and 98.4% BUSCO genes. All 4 isolates have identical ITS sequences and are distinct from Zasmidium isolates that were previously sampled from banana leaves. We thus report the first continuous genome assembly of a member of the Zasmidium genus, forming an essential resource for further analysis to enhance our understanding of the diversity of pathogenic fungal isolates as well as fungal diversity.

Uncoupled evolution of the Polycomb system and deep origin of non-canonical PRC1.

Polycomb group proteins, as part of the Polycomb repressive complexes, are essential in gene repression through chromatin compaction by canonical PRC1, mono-ubiquitylation of histone H2A by non-canonical PRC1 and tri-methylation of histone H3K27 by PRC2. Despite prevalent models emphasizing tight functional coupling between PRC1 and PRC2, it remains unclear whether this paradigm indeed reflects the evolution and functioning of these complexes. Here, we conduct a comprehensive analysis of the presence or absence of cPRC1, nPRC1 and PRC2 across the entire eukaryotic tree of life, and find that both complexes were present in the Last Eukaryotic Common Ancestor (LECA). Strikingly, ~42% of organisms contain only PRC1 or PRC2, showing that their evolution since LECA is largely uncoupled. The identification of ncPRC1-defining subunits in unicellular relatives of animals and fungi suggests ncPRC1 originated before cPRC1, and we propose a scenario for the evolution of cPRC1 from ncPRC1. Together, our results suggest that crosstalk between these complexes is a secondary development in evolution.

Epigenetic regulation of nuclear processes in fungal plant pathogens.

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

Convergent evolution of plant pattern recognition receptors sensing cysteine-rich patterns from three microbial kingdoms.

The Arabidopsis thaliana Receptor-Like Protein RLP30 contributes to immunity against the fungal pathogen Sclerotinia sclerotiorum. Here we identify the RLP30-ligand as a small cysteine-rich protein (SCP) that occurs in many fungi and oomycetes and is also recognized by the Nicotiana benthamiana RLP RE02. However, RLP30 and RE02 share little sequence similarity and respond to different parts of the native/folded protein. Moreover, some Brassicaceae other than Arabidopsis also respond to a linear SCP peptide instead of the folded protein, suggesting that SCP is an eminent immune target that led to the convergent evolution of distinct immune receptors in plants. Surprisingly, RLP30 shows a second ligand specificity for a SCP-nonhomologous protein secreted by bacterial Pseudomonads. RLP30 expression in N. tabacum results in quantitatively lower susceptibility to bacterial, fungal and oomycete pathogens, thus demonstrating that detection of immunogenic patterns by Arabidopsis RLP30 is involved in defense against pathogens from three microbial kingdoms.

Detection and quantification of the histone code in the fungal genus Aspergillus.

In eukaryotes, the combination of different histone post-translational modifications (PTMs) - the histone code - impacts the chromatin organization as compact and transcriptionally silent heterochromatin or accessible and transcriptionally active euchromatin. Although specific histone PTMs have been studied in fungi, an overview of histone PTMs and their relative abundance is still lacking. Here, we used mass spectrometry to detect and quantify histone PTMs in three fungal species belonging to three distinct taxonomic sections of the genus Aspergillus (Aspergillus niger, Aspergillus nidulans (two strains), and Aspergillus fumigatus). We overall detected 23 different histone PTMs, including a majority of lysine methylations and acetylations, and 23 co-occurrence patterns of multiple histone PTMs. Among those, we report for the first time the detection of H3K79me1, H3K79me2, and H4K31ac in Aspergilli. Although all three species harbour the same PTMs, we found significant differences in the relative abundance of H3K9me1/2/3, H3K14ac, H3K36me1 and H3K79me1, as well as the co-occurrence of acetylation on both K18 and K23 of histone H3 in a strain-specific manner. Our results provide novel insights about the underexplored complexity of the histone code in filamentous fungi, and its functional implications on genome architecture and gene regulation.

Nuclear genome organization in fungi: from gene folding to Rabl chromosomes.

Comparative genomics has recently provided unprecedented insights into the biology and evolution of the fungal lineage. In the postgenomics era, a major research interest focuses now on detailing the functions of fungal genomes, i.e. how genomic information manifests into complex phenotypes. Emerging evidence across diverse eukaryotes has revealed that the organization of DNA within the nucleus is critically important. Here, we discuss the current knowledge on the fungal genome organization, from the association of chromosomes within the nucleus to topological structures at individual genes and the genetic factors required for this hierarchical organization. Chromosome conformation capture followed by high-throughput sequencing (Hi-C) has elucidated how fungal genomes are globally organized in Rabl configuration, in which centromere or telomere bundles are associated with opposite faces of the nuclear envelope. Further, fungal genomes are regionally organized into topologically associated domain-like (TAD-like) chromatin structures. We discuss how chromatin organization impacts the proper function of DNA-templated processes across the fungal genome. Nevertheless, this view is limited to a few fungal taxa given the paucity of fungal Hi-C experiments. We advocate for exploring genome organization across diverse fungal lineages to ensure the future understanding of the impact of nuclear organization on fungal genome function.

A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen.

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

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

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

Dissemination of Fusarium Wilt of Banana in Mozambique Caused by Fusarium odoratissimum Tropical Race 4.

Fusarium wilt of banana (FWB) is a serious soil-borne fungal disease. In the previous century, FWB already destroyed Gros Michel-based banana cultures in Central America, and currently, the disease threatens all major banana-producing regions of the world. The causal agents of these epidemics, however, are diverse. Gros Michel was infected by a wide range of Fusarium species, the so-called Race 1 strains, whereas the contemporary Cavendish-based cultures are affected by Fusarium odoratissimum, colloquially called Tropical Race 4 (TR4). TR4 was reported in Mozambique on two commercial banana farms in 2013, but no incursions were found outside the farm boundaries in 2015, suggesting that the disease was under control. Here we report the presence of TR4 outside of these farm boundaries. We obtained fungal samples from 13 banana plants in smallholder and roadside plantings at various locations throughout northern Mozambique. These samples tested positive for TR4 by molecular diagnostics and in greenhouse pathogenicity assays. The results were confirmed with reisolations, thereby completing Koch's postulates. To study the diversity of TR4 isolates in Mozambique, we selected five samples for whole-genome sequencing. Comparison with a global collection of TR4 samples revealed very little genetic variation, indicating that the fungus is clonally spreading in Mozambique. Furthermore, isolates from Mozambique are clearly genetically separated from other geographic incursions, and thus we cannot trace the origin of TR4 in Mozambique. Nevertheless, our data demonstrates the dissemination of TR4 in Mozambique, underscoring the failure of disease management strategies. This threatens African banana production.

The histone code of the fungal genus Aspergillus uncovered by evolutionary and proteomic analyses.

Chemical modifications of DNA and histone proteins impact the organization of chromatin within the nucleus. Changes in these modifications, catalysed by different chromatin-modifying enzymes, influence chromatin organization, which in turn is thought to impact the spatial and temporal regulation of gene expression. While combinations of different histone modifications, the histone code, have been studied in several model species, we know very little about histone modifications in the fungal genus Aspergillus, whose members are generally well studied due to their importance as models in cell and molecular biology as well as their medical and biotechnological relevance. Here, we used phylogenetic analyses in 94 Aspergilli as well as other fungi to uncover the occurrence and evolutionary trajectories of enzymes and protein complexes with roles in chromatin modifications or regulation. We found that these enzymes and complexes are highly conserved in Aspergilli, pointing towards a complex repertoire of chromatin modifications. Nevertheless, we also observed few recent gene duplications or losses, highlighting Aspergillus species to further study the roles of specific chromatin modifications. SET7 (KMT6) and other components of PRC2 (Polycomb Repressive Complex 2), which is responsible for methylation on histone H3 at lysine 27 in many eukaryotes including fungi, are absent in Aspergilli as well as in closely related Penicillium species, suggesting that these lost the capacity for this histone modification. We corroborated our computational predictions by performing untargeted MS analysis of histone post-translational modifications in Aspergillus nidulans. This systematic analysis will pave the way for future research into the complexity of the histone code and its functional implications on genome architecture and gene regulation in fungi.

Induced resistance to Fusarium wilt of banana caused by Tropical Race 4 in Cavendish cv Grand Naine bananas after challenging with avirulent Fusarium spp.

In the last century, Fusarium wilt of banana (FWB) destroyed the banana cultivar Gros Michel. The Cavendish cultivars saved the global banana industry, and currently they dominate global production (~50%) and the export trade (~95%). However, a new strain called Tropical Race 4 (TR4) surfaced in the late 1960's, spread globally and greatly damages Cavendish plantations as well as manifold local varieties that are primarily grown by small holders. Presently, there is no commercially available replacement for Cavendish and hence control strategies must be developed and implemented to manage FWB. Here, we studied whether it is possible to induce resistance to TR4 by pre-inoculations with different Fusarium spp. Only pre-treatments with an avirulent Race 1 strain significantly reduced disease development of TR4 in a Cavendish genotype and this effect was stable at various nutritional and pH conditions. We then used transcriptome analysis to study the molecular basis of this response. Several genes involved in plant defence responses were up-regulated during the initial stages of individual infections with TR4 and Race 1, as well as in combined treatments. In addition, a number of genes in the ethylene and jasmonate response pathways as well as several gibberellin synthesis associated genes were induced. We observed upregulation of RGA2 like genes in all treatments. Hence, RGA2 could be a key factor involved in both R1 and TR4 resistance. The data support the hypothesis that activating resistance to Race 1 in Cavendish bananas affects TR4 development and provide a first insight of gene expression during the interaction between various Fusarium spp. and banana.

Caenorhabditis elegans MES-3 is a highly divergent ortholog of the canonical PRC2 component SUZ12.

Polycomb Repressive Complex 2 (PRC2) catalyzes the mono-, di-, and trimethylation of histone protein H3 on lysine 27 (H3K27), which is strongly associated with transcriptionally silent chromatin. The functional core of PRC2 is highly conserved in animals and consists of four subunits. One of these, SUZ12, has not been identified in the genetic model Caenorhabditis elegans, whereas C. elegans PRC2 contains the clade-specific MES-3 protein. Through unbiased sensitive sequence similarity searches complemented by high-quality structure predictions of monomers and multimers, we here demonstrate that MES-3 is a highly divergent ortholog of SUZ12. MES-3 shares protein folds and conserved residues of key domains with SUZ12 and is predicted to interact with core PRC2 members similar to SUZ12 in human PRC2. Thus, in agreement with previous genetic and biochemical studies, we provide evidence that C. elegans contains a diverged yet evolutionary conserved core PRC2, like other animals.