Gapless genome assembly of Colletotrichum higginsianum reveals chromosome structure and association of transposable elements with secondary metabolite gene clusters

The Molecular Genetics and Genomics of Plant-Pathogen Interactions

Plants have evolved an intricate immune system to protect themselves from potential pathogens [...].

Guanosine-specific single-stranded ribonuclease effectors of a phytopathogenic fungus potentiate host immune responses

Plants activate immunity upon recognition of pathogen-associated molecular patterns. Although phytopathogens have evolved a set of effector proteins to counteract plant immunity, some effectors are perceived by hosts and induce immune responses. Here, we show that two secreted ribonuclease effectors, SRN1 and SRN2, encoded in a phytopathogenic fungus, Colletotrichum orbiculare, induce cell death in a signal peptide- and catalytic residue-dependent manner, when transiently expressed in Nicotiana benthamiana. The pervasive presence of SRN genes across Colletotrichum species suggested the conserved roles. Using a transient gene expression system in cucumber (Cucumis sativus), an original host of C. orbiculare, we show that SRN1 and SRN2 potentiate host pattern-triggered immunity responses. Consistent with this, C. orbiculare SRN1 and SRN2 deletion mutants exhibited increased virulence on the host. In vitro analysis revealed that SRN1 specifically cleaves single-stranded RNAs at guanosine, leaving a 3'-end phosphate. Importantly, the potentiation of C. sativus responses by SRN1 and SRN2, present in the apoplast, depends on ribonuclease catalytic residues. We propose that the pathogen-derived apoplastic guanosine-specific single-stranded endoribonucleases lead to immunity potentiation in plants.

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.

Assaying Effector Cell-to-Cell Mobility in Plant Tissues Identifies Hypermobility and Indirect Manipulation of Plasmodesmata

In plants, plasmodesmata establish cytoplasmic continuity between cells to allow for communication and resource exchange across the cell wall. While plant pathogens use plasmodesmata as a pathway for both molecular and physical invasion, the benefits of molecular invasion (cell-to-cell movement of pathogen effectors) are poorly understood. To establish a methodology for identification and characterization of the cell-to-cell mobility of effectors, we performed a quantitative live imaging-based screen of candidate effectors of the fungal pathogen Colletotrichum higginsianum. We predicted C. higginsianum effectors by their expression profiles, the presence of a secretion signal, and their predicted and in planta localization when fused to green fluorescent protein. We assayed for cell-to-cell mobility of nucleocytosolic effectors and identified 14 that are cell-to-cell mobile. We identified that three of these effectors are "hypermobile," showing cell-to-cell mobility greater than expected for a protein of that size. To explore the mechanism of hypermobility, we chose two hypermobile effectors and measured their impact on plasmodesmata function and found that even though they show no direct association with plasmodesmata, each increases the transport capacity of plasmodesmata. Thus, our methods for quantitative analysis of cell-to-cell mobility of candidate microbe-derived effectors, or any suite of host proteins, can identify cell-to-cell hypermobility and offer greater understanding of how proteins affect plasmodesmal function and intercellular connectivity. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

A gln-tRNA-based CRISPR/Cas9 knockout system enables the functional characterization of genes in the genetically recalcitrant brassica anthracnose fungus Colletotrichum higginsianum

Colletotrichum higginsianum causes anthracnose disease in brassicas. The availability of the C. higginsianum genome has paved the way for the genome-wide exploration of genes associated with virulence/pathogenicity. However, delimiting the biological functions of these genes remains an arduous task due to the recalcitrance of C. higginsianum to genetic manipulations. Here, we report a CRISPR/Cas9-based system that can knock out the genes in C. higginsianum with a staggering 100% homologous recombination frequency (HRF). The system comprises two vectors: pCas9-Ch_tRp-sgRNA, in which a C. higginsianum glutaminyl-tRNA drives the expression of sgRNA, and pCE-Zero-HPT carrying a donor DNA cassette containing the marker gene HPT flanked by homology arms. Upon co-transformation of the C. higginsianum protoplasts, pCas9-Ch_tRp-sgRNA causes a DNA double-strand break in the targeted gene, followed by homology-directed replacement of the gene with HPT by pCE-Zero-HPT, thereby generating loss-of-function mutants. Using the system, we generated the knockout mutants of two effector candidates (ChBas3 and OBR06881) with a 100% HRF. Interestingly, the ΔChBas3 and ΔOBR06881 mutants did not seem to affect the C. higginsianum infection of Arabidopsis thaliana. Altogether, the CRISPR/Cas9 system developed in the study enables the targeted deletion of genes, including effectors, in C. higginsianum, thus determining their biological functions.

Comparative Genomic Analysis of Colletotrichum lini Strains with Different Virulence on Flax

Colletotrichum lini is a flax fungal pathogen. The genus comprises differently virulent strains, leading to significant yield losses. However, there were no attempts to investigate the molecular mechanisms of C. lini pathogenicity from high-quality genome assemblies until this study. In this work, we sequenced the genomes of three C. lini strains of high (#390-1), medium (#757), and low (#771) virulence. We obtained more than 100× genome coverage with Oxford Nanopore Technologies reads (N50 = 12.1, 6.1, 5.0 kb) and more than 50× genome coverage with Illumina data (150 + 150 bp). Several assembly strategies were tested. The final assemblies were obtained using the Canu-Racon ×2-Medaka-Polca scheme. The assembled genomes had a size of 54.0-55.3 Mb, 26-32 contigs, N50 values > 5 Mb, and BUSCO completeness > 96%. A comparative genomic analysis showed high similarity among mitochondrial and nuclear genomes. However, a rearrangement event and the loss of a 0.7 Mb contig were revealed. After genome annotation with Funannotate, secreting proteins were selected using SignalP, and candidate effectors were predicted among them using EffectorP. The analysis of the InterPro annotations of predicted effectors revealed unique protein categories in each strain. The assembled genomes and the conducted comparative analysis extend the knowledge of the genetic diversity of C. lini and form the basis for establishing the molecular mechanisms of its pathogenicity.

The Hidden Truths of Fungal Virulence and Adaptation on Hosts: Unraveling the Conditional Dispensability of Minichromosomes in the Hemibiotrophic Colletotrichum Pathogens

Colletotrichum spp. are ascomycete fungi and cause anthracnose disease in numerous crops of economic significance. The genomes of these fungi are distributed among ten core chromosomes and two to three minichromosomes. While the core chromosomes regulate fungal growth, development and virulence, the extent to which the minichromosomes are involved in these processes is still uncertain. Here, we discuss the minichromosomes of three hemibiotrophic Colletotrichum pathogens, i.e., C. graminicola, C. higginsianum and C. lentis. These minichromosomes are typically less than one megabase in length, characterized by containing higher repetitive DNA elements, lower GC content, higher frequency of repeat-induced point mutations (RIPMs) and sparse gene distribution. Molecular genetics and functional analyses have revealed that these pathogens harbor one conditionally dispensable minichromosome, which is dispensable for fungal growth and development but indispensable for fungal virulence on hosts. They appear to be strain-specific innovations and are highly compartmentalized into AT-rich and GC-rich blocks, resulting from RIPMs, which may help protect the conditionally dispensable minichromosomes from erosion of already scarce genes, thereby helping the Colletotrichum pathogens maintain adaptability on hosts. Overall, understanding the mechanisms underlying the conditional dispensability of these minichromosomes could lead to new strategies for controlling anthracnose disease in crops.

Efficient multiple gene knockout in Colletotrichum higginsianum via CRISPR/Cas9 ribonucleoprotein and URA3-based marker recycling

Colletotrichum higginsianum is a hemibiotrophic pathogen that causes anthracnose disease on crucifer hosts, including Arabidopsis thaliana. Despite the availability of genomic and transcriptomic information and the ability to transform both organisms, identifying C. higginsianum genes involved in virulence has been challenging due to recalcitrance to gene targeting and redundancy of virulence factors. To overcome these obstacles, we developed an efficient method for multiple gene disruption in C. higginsianum by combining CRISPR/Cas9 and a URA3-based marker recycling system. Our method significantly increased the efficiency of gene knockout via homologous recombination by introducing genomic DNA double-strand breaks. We demonstrated the applicability of the URA3-based marker recycling system for multiple gene targeting in the same strain. Using our technology, we successfully targeted two melanin biosynthesis genes, SCD1 and PKS1, which resulted in deficiency in melanization and loss of pathogenicity in the mutants. Our findings demonstrate the effectiveness of our methods in analysing virulence factors in C. higginsianum, thus accelerating research on plant-fungus interactions.

Yeast-based heterologous production of the Colletochlorin family of fungal secondary metabolites

Transcriptomic studies have revealed that fungal pathogens of plants activate the expression of numerous biosynthetic gene clusters (BGC) exclusively when in presence of a living host plant. The identification and structural elucidation of the corresponding secondary metabolites remain challenging. The aim was to develop a polycistronic system for heterologous expression of fungal BGCs in Saccharomyces cerevisiae. Here we adapted a polycistronic vector for efficient, seamless and cost-effective cloning of biosynthetic genes using in vivo assembly (also called transformation-assisted recombination) directly in Escherichia coli followed by heterologous expression in S. cerevisiae. Two vectors were generated with different auto-inducible yeast promoters and selection markers. The effectiveness of these vectors was validated with fluorescent proteins. As a proof-of-principle, we applied our approach to the Colletochlorin family of molecules. These polyketide secondary metabolites were known from the phytopathogenic fungus Colletotrichum higginsianum but had never been linked to their biosynthetic genes. Considering the requirement for a halogenase, and by applying comparative genomics, we identified a BGC putatively involved in the biosynthesis of Colletochlorins in C. higginsianum. Following the expression of those genes in S. cerevisiae, we could identify the presence of the precursor Orsellinic acid, Colletochlorins and their non-chlorinated counterparts, the Colletorins. In conclusion, the polycistronic vectors described herein were adapted for the host S. cerevisiae and allowed to link the Colletochlorin compound family to their corresponding biosynthetic genes. This system will now enable the production and purification of infection-specific secondary metabolites of fungal phytopathogens. More widely, this system could be applied to any fungal BGC of interest.

Dynamics of transposable element accumulation in the non-recombining regions of mating-type chromosomes in anther-smut fungi

In the absence of recombination, the number of transposable elements (TEs) increases due to less efficient selection, but the dynamics of such TE accumulations are not well characterized. Leveraging a dataset of 21 independent events of recombination cessation of different ages in mating-type chromosomes of Microbotryum fungi, we show that TEs rapidly accumulated in regions lacking recombination, but that TE content reached a plateau at ca. 50% of occupied base pairs by 1.5 million years following recombination suppression. The same TE superfamilies have expanded in independently evolved non-recombining regions, in particular rolling-circle replication elements (Helitrons). Long-terminal repeat (LTR) retrotransposons of the Copia and Ty3 superfamilies also expanded, through transposition bursts (distinguished from gene conversion based on LTR divergence), with both non-recombining regions and autosomes affected, suggesting that non-recombining regions constitute TE reservoirs. This study improves our knowledge of genome evolution by showing that TEs can accumulate through bursts, following non-linear decelerating dynamics.

The evolution of mini-chromosomes in the fungal genus Colletotrichum

Anthracnose diseases caused by Colletotrichum species are among the most common fungal diseases. These symptoms typically manifest as dark, sunken lesions on leaves, stems, and fruit. In China, mango anthracnose seriously affects fruit yield and quality. Genome sequencing of several species shows the presence of mini-chromosomes. These are thought to contribute to virulence, but their formation and activity remain to be fully elucidated. Here, we assembled 17 Colletotrichum genomes (16 isolated from mango plus one from persimmon) through PacBio long-read sequencing. Half of the assembled scaffolds had telomeric repeats at both ends indicating full-length chromosomes. Based on comparative genomics analysis at interspecies and intraspecies levels, we identified extensive chromosomal rearrangements events. We analyzed mini-chromosomes of Colletotrichum spp. and found large variation among close relatives. In C. fructicola, homology between core chromosomes and mini-chromosomes suggested that some mini-chromosomes were generated by recombination of core chromosomes. In C. musae GZ23-3, we found 26 horizontally transferred genes arranged in clusters on mini-chromosomes. In C. asianum FJ11-1, several potential pathogenesis-related genes on mini-chromosomes were upregulated, especially in strains with highly pathogenic phenotypes. Mutants of these upregulated genes showed obvious defects in virulence. Our findings provide insights into the evolution and potential relationships to virulence associated with mini-chromosomes. IMPORTANCE Colletotrichum is a cosmopolitan fungal genus that seriously affects fruit yield and quality of many plant species. Mini-chromosomes have been found to be related to virulence in Colletotrichum. Further examination of mini-chromosomes can help us elucidate some pathogenic mechanisms of Colletotrichum. In this study, we generated novel assemblies of several Colletotrichum strains. Comparative genomic analyses within and between Colletotrichum species were conducted. We then identified mini-chromosomes in our sequenced strains systematically. The characteristics and generation of mini-chromosomes were investigated. Transcriptome analysis and gene knockout revealed pathogenesis-related genes located on mini-chromosomes of C. asianum FJ11-1. This study represents the most comprehensive investigation of chromosome evolution and potential pathogenicity of mini-chromosomes in the Colletotrichum genus.

CRISPR/Cas9-mediated deletion of large chromosomal segments identifies a minichromosome modulating the Colletotrichum graminicola virulence on maize

Colletotrichum graminicola causes anthracnose on maize, an economically significant disease worldwide. To decipher how the pathogen controls its virulence/pathogenicity on maize at the minichromosomal level, we sequenced the genome and transcriptome of the C. graminicola strain T1-3-3. The 61.91 Mb genome contains three transcriptionally repressed, full-length strain-specific minichromosomes (<1 Mb; Chr11 through Chr13). A CRISPR/Cas9-based system was developed to knock out large chromosomal segments; it involved the generation of multiple simultaneous DNA double-strand breaks across a targeted genomic region, followed by homology-directed replacement thereof with a donor DNA template carrying the selectable marker hygromycin phosphotransferase gene flanked by homologous sequence arms of the targeted region. Using this system, we obtained distinct mutants functionally nullisomic for individual minichromosomes. Only the ΔChr12 mutant lacking the 498.44 Kb genomic region carrying all of the 31 genes of Chr12 exhibited attenuated virulence on maize and was indistinguishable from T1-3-3 in fungal growth and conidiation, indicating that Chr12 is a conditionally dispensable minichromosome and imparts full virulence to C. graminicola on maize. The CRISPR/Cas9-mediated genome editing system developed in this study will enable the determination of the biological functions of minichromosomes or large chromosomal segments in fungal plant pathogens.

Quantitative trait locus for calving traits on Bos taurus autosome 18 in Holstein cattle is embedded in a complex genomic region

Although the quantitative trait locus (QTL) on chromosome 18 (BTA18) associated with paternal calving ease and stillbirth in Holstein Friesian cattle and its cross has been known for over 20 years, to our knowledge, the exact causal genetic sequence has yet escaped identification. The aim of this study was to re-examine the region of the published QTL on BTA18 and to investigate the possible reasons behind this elusiveness. For this purpose, we carried out a combined linkage disequilibrium and linkage analysis using genotyping data of 2,697 German Holstein Friesian (HF) animals and subsequent whole-genome sequencing (WGS) data analyses and genome assembly of HF samples. We confirmed the known QTL in the 95% confidence interval of 1.089 Mbp between 58.34 and 59.43 Mbp on BTA18. Additionally, these 4 SNPs in the near-perfect linkage disequilibrium with the QTL haplotype were identified: rs381577268 (on 57,816,137 bp, C/T), rs381878735 (on 59,574,329 bp, A/T), rs464221818 (on 59,329,176 bp, C/T), and rs472502785 (on 59,345,689 bp, T/C). Search for the causal mutation using short and long-read sequences, and methylation data of the BTA18 QTL region did not reveal any candidates though. The assembly showed problems in the region, as well as an abundance of segmental duplications within and around the region. Taking the QTL of BTA18 in Holstein cattle as an example, the data presented in this study comprehensively characterize the genomic features that could also be relevant for other such elusive QTL in various other cattle breeds and livestock species as well.

CusProSe: a customizable protein annotation software with an application to the prediction of fungal secondary metabolism genes

We report here a new application, CustomProteinSearch (CusProSe), whose purpose is to help users to search for proteins of interest based on their domain composition. The application is customizable. It consists of two independent tools, IterHMMBuild and ProSeCDA. IterHMMBuild allows the iterative construction of Hidden Markov Model (HMM) profiles for conserved domains of selected protein sequences, while ProSeCDA scans a proteome of interest against an HMM profile database, and annotates identified proteins using user-defined rules. CusProSe was successfully used to identify, in fungal genomes, genes encoding key enzyme families involved in secondary metabolism, such as polyketide synthases (PKS), non-ribosomal peptide synthetases (NRPS), hybrid PKS-NRPS and dimethylallyl tryptophan synthases (DMATS), as well as to characterize distinct terpene synthases (TS) sub-families. The highly configurable characteristics of this application makes it a generic tool, which allows the user to refine the function of predicted proteins, to extend detection to new enzymes families, and may also be applied to biological systems other than fungi and to other proteins than those involved in secondary metabolism.

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.

Fungal Secondary Metabolites and Small RNAs Enhance Pathogenicity during Plant-Fungal Pathogen Interactions

Fungal plant pathogens use proteinaceous effectors as well as newly identified secondary metabolites (SMs) and small non-coding RNA (sRNA) effectors to manipulate the host plant's defense system via diverse plant cell compartments, distinct organelles, and many host genes. However, most molecular studies of plant-fungal interactions have focused on secreted effector proteins without exploring the possibly equivalent functions performed by fungal (SMs) and sRNAs, which are collectively known as "non-proteinaceous effectors". Fungal SMs have been shown to be generated throughout the plant colonization process, particularly in the early biotrophic stages of infection. The fungal repertoire of non-proteinaceous effectors has been broadened by the discovery of fungal sRNAs that specifically target plant genes involved in resistance and defense responses. Many RNAs, particularly sRNAs involved in gene silencing, have been shown to transmit bidirectionally between fungal pathogens and their hosts. However, there are no clear functional approaches to study the role of these SM and sRNA effectors. Undoubtedly, fungal SM and sRNA effectors are now a treasured land to seek. Therefore, understanding the role of fungal SM and sRNA effectors may provide insights into the infection process and identification of the interacting host genes that are targeted by these effectors. This review discusses the role of fungal SMs and sRNAs during plant-fungal interactions. It will also focus on the translocation of sRNA effectors across kingdoms, the application of cross-kingdom RNA interference in managing plant diseases and the tools that can be used to predict and study these non-proteinaceous effectors.

Proteinaceous Effector Discovery and Characterization in Plant Pathogenic Colletotrichum Fungi

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

Fungal Pathogenesis-Related Cell Wall Biogenesis, with Emphasis on the Maize Anthracnose Fungus Colletotrichum graminicola

The genus Colletotrichum harbors many plant pathogenic species, several of which cause significant yield losses in the field and post harvest. Typically, in order to infect their host plants, spores germinate, differentiate a pressurized infection cell, and display a hemibiotrophic lifestyle after plant invasion. Several factors required for virulence or pathogenicity have been identified in different Colletotrichum species, and adaptation of cell wall biogenesis to distinct stages of pathogenesis has been identified as a major pre-requisite for the establishment of a compatible parasitic fungus-plant interaction. Here, we highlight aspects of fungal cell wall biogenesis during plant infection, with emphasis on the maize leaf anthracnose and stalk rot fungus, Colletotrichum graminicola.

Fast characterization of segmental duplication structure in multiple genome assemblies

Motivation

The increasing availability of high-quality genome assemblies raised interest in the characterization of genomic architecture. Major architectural elements, such as common repeats and segmental duplications (SDs), increase genome plasticity that stimulates further evolution by changing the genomic structure and inventing new genes. Optimal computation of SDs within a genome requires quadratic-time local alignment algorithms that are impractical due to the size of most genomes. Additionally, to perform evolutionary analysis, one needs to characterize SDs in multiple genomes and find relations between those SDs and unique (non-duplicated) segments in other genomes. A naïve approach consisting of multiple sequence alignment would make the optimal solution to this problem even more impractical. Thus there is a need for fast and accurate algorithms to characterize SD structure in multiple genome assemblies to better understand the evolutionary forces that shaped the genomes of today.

Results

Here we introduce a new approach, BISER, to quickly detect SDs in multiple genomes and identify elementary SDs and core duplicons that drive the formation of such SDs. BISER improves earlier tools by (i) scaling the detection of SDs with low homology to multiple genomes while introducing further 7–33\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}× speed-ups over the existing tools, and by (ii) characterizing elementary SDs and detecting core duplicons to help trace the evolutionary history of duplications to as far as 300 million years.

Availability and implementation

BISER is implemented in Seq programming language and is publicly available at https://github.com/0xTCG/biser.

A Comprehensive Assessment of the Secretome Responsible for Host Adaptation of the Legume Root Pathogen Aphanomyces euteiches

The soil-borne oomycete pathogen Aphanomyces euteiches causes devastating root rot diseases in legumes such as pea and alfalfa. The different pathotypes of A. euteiches have been shown to exhibit differential quantitative virulence, but the molecular basis of host adaptation has not yet been clarified. Here, we re-sequenced a pea field reference strain of A. euteiches ATCC201684 with PacBio long-reads and took advantage of the technology to generate the mitochondrial genome. We identified that the secretome of A. euteiches is characterized by a large portfolio of secreted proteases and carbohydrate-active enzymes (CAZymes). We performed Illumina sequencing of four strains of A. euteiches with contrasted specificity to pea or alfalfa and found in different geographical areas. Comparative analysis showed that the core secretome is largely represented by CAZymes and proteases. The specific secretome is mainly composed of a large set of small, secreted proteins (SSP) without any predicted functional domain, suggesting that the legume preference of the pathogen is probably associated with unknown functions. This study forms the basis for further investigations into the mechanisms of interaction of A. euteiches with legumes.

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

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

The Conserved Colletotrichum spp. Effector Candidate CEC3 Induces Nuclear Expansion and Cell Death in Plants

Plant pathogens secrete proteins, known as effectors, that promote infection by manipulating host cells. Members of the phytopathogenic fungal genus Colletotrichum collectively have a broad host range and generally adopt a hemibiotrophic lifestyle that includes an initial biotrophic phase and a later necrotrophic phase. We hypothesized that Colletotrichum fungi use a set of conserved effectors during infection to support the two phases of their hemibiotrophic lifestyle. This study aimed to examine this hypothesis by identifying and characterizing conserved effectors among Colletotrichum fungi. Comparative genomic analyses using genomes of ascomycete fungi with different lifestyles identified seven effector candidates that are conserved across the genus Colletotrichum. Transient expression assays showed that one of these putative conserved effectors, CEC3, induces nuclear expansion and cell death in Nicotiana benthamiana, suggesting that CEC3 is involved in promoting host cell death during infection. Nuclear expansion and cell death induction were commonly observed in CEC3 homologs from four different Colletotrichum species that vary in host specificity. Thus, CEC3 proteins could represent a novel class of core effectors with functional conservation in the genus Colletotrichum.

Ornithine decarboxylase of the fungal pathogen Colletotrichum higginsianum plays an important role in regulating global metabolic pathways and virulence

Colletotrichum higginsianum is an important fungal pathogen causing anthracnose disease of cruciferous plants. In this study, we characterized a putative orthologue of yeast SPE1 in C. higginsianum, named ChODC. Deletion mutants of ChODC were defective in hyphal and conidial development. Importantly, deletion of ChODC significantly affected appressorium-mediated penetration in C. higginsianum. However, polyamines partially restore appressorium function and virulence indicating that loss of ChODC caused significantly decreased virulence by the crosstalk between polyamines and other metabolic pathways. Subsequently, transcriptomic and metabolomic analyses demonstrated that ChODC played an important role in metabolism of various carbon and nitrogen compounds including amino acids, carbohydrates and lipids. Along with these clues, we found deletion of ChODC affected glycogen and lipid metabolism, which were important for conidial storage utilization and functional appressorium formation. Loss of ChODC affected the mTOR signalling pathway via modulation of autophagy. Interestingly, cAMP treatment restored functional appressoria to the ΔChODC mutant, and rapamycin treatment also stimulated formation of functional appressoria in the ΔChODC mutant. Overall, ChODC was associated with the polyamine biosynthesis pathway, as a mediator of cAMP and mTOR signalling pathways to regulate appressorium function. Our study provides evidence of a link between ChODC and the cAMP signalling pathway and defines a novel mechanism by which ChODC regulates infection-associated autophagy and plant infection by fungi.

Use of Telomere Fingerprinting to Identify Clonal Lineages of Colletotrichum fioriniae in Kentucky Mixed-Fruit Orchards

Multiple species in the fungal genus Colletotrichum cause anthracnose fruit rot diseases that are responsible for major yield losses of as much as 100%. Individual species of Colletotrichum typically have broad host ranges and can infect multiple fruit species. Colletotrichum fioriniae causes anthracnose fruit rots of apples, blueberries, and strawberries in Kentucky orchards where these fruits grow in close proximity. This raises the possibility of cross-infection, which may have significant management implications. The potential occurrence of cross-infection was investigated by using telomere fingerprinting to identify C. fioriniae clones in several mixed-fruit orchards. Telomere fingerprints were highly polymorphic among a test group of C. fioriniae strains and effectively defined clonal lineages. Fingerprints were compared among apple, blueberry, and strawberry isolates of C. fioriniae from three different orchards and similarity matrices were calculated to build phylograms for each orchard group. Multiple clonal lineages of C. fioriniae were identified within each orchard on the same fruit host. Related lineages were found among isolates from different hosts, but the results did not provide direct evidence for cross-infection of different fruit species by the same clones. Recovery of the same clonal lineages within orchards across multiple years suggested that local dispersal was important in pathogen population structure and that C. fioriniae strains persisted within orchards over time. Isolates from blueberry were less diverse than isolates from apple, perhaps related to more intensive anthracnose management protocols on apple versus blueberry. Telomere fingerprinting is a valuable tool for understanding population dynamics of Colletotrichum fruit rot fungi.

Genome-Wide Analysis of Nutrient Signaling Pathways Conserved in Arbuscular Mycorrhizal Fungi

Arbuscular mycorrhizal (AM) fungi form a mutualistic symbiosis with a majority of terrestrial vascular plants. To achieve an efficient nutrient trade with their hosts, AM fungi sense external and internal nutrients, and integrate different hierarchic regulations to optimize nutrient acquisition and homeostasis during mycorrhization. However, the underlying molecular networks in AM fungi orchestrating the nutrient sensing and signaling remain elusive. Based on homology search, we here found that at least 72 gene components involved in four nutrient sensing and signaling pathways, including cAMP-dependent protein kinase A (cAMP-PKA), sucrose non-fermenting 1 (SNF1) protein kinase, target of rapamycin kinase (TOR) and phosphate (PHO) signaling cascades, are well conserved in AM fungi. Based on the knowledge known in model yeast and filamentous fungi, we outlined the possible gene networks functioning in AM fungi. These pathways may regulate the expression of downstream genes involved in nutrient transport, lipid metabolism, trehalase activity, stress resistance and autophagy. The RNA-seq analysis and qRT-PCR results of some core genes further indicate that these pathways may play important roles in spore germination, appressorium formation, arbuscule longevity and sporulation of AM fungi. We hope to inspire further studies on the roles of these candidate genes involved in these nutrient sensing and signaling pathways in AM fungi and AM symbiosis.

Genome-wide mapping of histone modifications during axenic growth in two species of Leptosphaeria maculans showing contrasting genomic organization

Leptosphaeria maculans 'brassicae' (Lmb) and Leptosphaeria maculans 'lepidii' (Lml) are closely related phytopathogenic species that exhibit a large macrosynteny but contrasting genome structure. Lmb has more than 30% of repeats clustered in large repeat-rich regions, while the Lml genome has only a small amount of evenly distributed repeats. Repeat-rich regions of Lmb are enriched in effector genes, expressed during plant infection. The distinct genome structures of Lmb and Lml provide an excellent model for comparing the organization of pathogenicity genes in relation to the chromatin landscape in two closely related phytopathogenic fungi. Here, we performed chromatin immunoprecipitation (ChIP) during axenic culture, targeting histone modifications typical for heterochromatin or euchromatin, combined with transcriptomic analysis to analyze the influence of chromatin organization on gene expression. In both species, we found that facultative heterochromatin is enriched with genes lacking functional annotation, including numerous effector and species-specific genes. Notably, orthologous genes located in H3K27me3 domains are enriched with effector genes. Compared to other fungal species, including Lml, Lmb is distinct in having large H3K9me3 domains associated with repeat-rich regions that contain numerous species-specific effector genes. Discovery of these two distinctive heterochromatin landscapes now raises questions about their involvement in the regulation of pathogenicity, the dynamics of these domains during plant infection and the selective advantage to the fungus to host effector genes in H3K9me3 or H3K27me3 domains.

Construction of a chromosome-level genome and variation map for the Pacific oyster Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely distributed marine bivalve of great ecological and economic importance. In this study, we provide a high-quality chromosome-level genome assembled using Pacific Bioscience long reads and Hi-C-based and linkage-map-based scaffolding technologies and a high-resolution variation map constructed using large-scale resequencing analysis. The 586.8 Mb genome consists of 10 pseudochromosome sequences ranging from 38.6 to 78.9 Mb, containing 301 contigs with an N50 size of 3.1 Mb. A total of 30,078 protein-coding genes were predicted, of which 22,757 (75.7%) were high-reliability annotations supported by a homologous match to a curated protein in the SWISS-PROT database or transcript expression. Although a medium level of repeat components (57.2%) was detected, the genomic content of the segmental duplications reached 26.2%, which is the highest among the reported genomes. By whole genome resequencing analysis of 495 Pacific oysters, a comprehensive variation map was built, comprised of 4.78 million single nucleotide polymorphisms, 0.60 million short insertions and deletions, and 49,333 copy number variation regions. The structural variations can lead to an average interindividual genomic divergence of 0.21, indicating their crucial role in shaping the Pacific oyster genome diversity. The large amount of mosaic distributed repeat elements, small variations, and copy number variations indicate that the Pacific oyster is a diploid organism with an extremely high genomic complexity at the intra- and interindividual level. The genome and variation maps can improve our understanding of oyster genome diversity and enrich the resources for oyster molecular evolution, comparative genomics, and genetic research.

Telomeres and a repeat-rich chromosome encode effector gene clusters in plant pathogenic Colletotrichum fungi

Members of the Colletotrichum gloeosporioides species complex are causal agents of anthracnose in many commercially important plants. Closely related strains have different levels of pathogenicity on hosts despite their close phylogenetic relationship. To gain insight into the genetics underlying these differences, we generated and annotated whole-genome assemblies of multiple isolates of C. fructicola (Cf) and C. siamense (Cs), as well as three previously unsequenced species, C. aenigma (Ca), C. tropicale and C. viniferum with different pathogenicity on strawberry. Based on comparative genomics, we identified accessory regions with a high degree of conservation in strawberry-pathogenic Cf, Cs and Ca strains. These regions encode homologs of pathogenicity-related genes known as effectors, organized in syntenic gene clusters, with copy number variations in different strains of Cf, Cs and Ca. Analysis of highly contiguous assemblies of Cf, Cs and Ca revealed the association of related accessory effector gene clusters with telomeres and repeat-rich chromosomes and provided evidence of exchange between these two genomic compartments. In addition, expression analysis indicated that orthologues in syntenic gene clusters showed a tendency for correlated gene expression during infection. These data provide insight into mechanisms by which Colletotrichum genomes evolve, acquire and organize effectors.

A highly contiguous reference genome assembly for Colletotrichum falcatum pathotype Cf08 causing red rot disease in sugarcane

Among the biotic factors, which affect the productivity and quality of sugarcane, red rot disease caused by the fungal pathogen, Colletotrichum falcatum is the most devastating that cause enormous loss to millers as well as cane growers. We present a highly contiguous genome assembly of C. falcatum pathotype Cf08 which is virulent to popular sugarcane varieties grown in more than 3 million hectares in sub-tropical India. By performing long read sequencing on PacBio RSII system, 56.06 Mb assemblies with 238 contigs having N50 of 0.51 Mb and L50 of 34 was produced. A BUSCO completeness score of 97.24% (including 4.1% fragmented) of the entire C. falcatum Cf08 nuclear genome, greatly improved contiguity compared to an existing highly fragmented draft of C. falcatum Cf671 genome (48.13 Mb) was obtained. This Cf08 assembly had 54.14% GC content and possessed < 1% repetitive elements. A total of 18,635 protein-coding genes were predicted compared with 12,270 for Cf671. Among 617 CAZymes predicted, glycoside hydrolases were the predominant (298), and among 7264 genes associated with pathogenicity/virulence, 77 genes having effector functions were identified. The assembled genome showed its similarity with the genome of C. graminicola and C. higginsianum, the causal organisms of anthracnose in maize and in members of Brassicaceae, respectively. A total of 94 large sequences (> 100 kb) of Cf08 were mapped over C. higginsianum 10 of 12 chromosomes with 106 synteny blocks. Results discussed here would provide an important tool for future studies of evolutionary and functional genomics in C. falcatum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02695-x.

Plants under the Attack of Allies: Moving towards the Plant Pathobiome Paradigm

Plants are functional macrobes living in a close association with diverse communities of microbes and viruses as complex systems that continuously interact with the surrounding environment. The microbiota within the plant holobiont serves various essential and beneficial roles, such as in plant growth at different stages, starting from seed germination. Meanwhile, pathogenic microbes—differentiated from the rest of the plant microbiome based on their ability to damage the plant tissues through transient blooming under specific conditions—are also a part of the plant microbiome. Recent advances in multi-omics have furthered our understanding of the structure and functions of plant-associated microbes, and a pathobiome paradigm has emerged as a set of organisms (i.e., complex eukaryotic, microbial, and viral communities) within the plant’s biotic environment which interact with the host to deteriorate its health status. Recent studies have demonstrated that the one pathogen–one disease hypothesis is insufficient to describe the disease process in many cases, particularly when complex organismic communities are involved. The present review discusses the plant holobiont and covers the steady transition of plant pathology from the one pathogen–one disease hypothesis to the pathobiome paradigm. Moreover, previous reports on model plant diseases, in which more than one pathogen or co-operative interaction amongst pathogenic microbes is implicated, are reviewed and discussed.

Chromatin Dynamics Contribute to the Spatiotemporal Expression Pattern of Virulence Genes in a Fungal Plant Pathogen

Fungal plant pathogens possess a large repertoire of genes encoding putative effectors, which are crucial for infection. Many of these genes are expressed at low levels in the absence of the host but are strongly induced at specific stages of the infection. The mechanisms underlying this transcriptional reprogramming remain largely unknown. We investigated the role of the genomic environment and associated chromatin modifications of effector genes in controlling their expression pattern in the fungal wheat pathogen Zymoseptoria tritici. Depending on their genomic location, effector genes are epigenetically repressed in the absence of the host and during the initial stages of infection. Derepression of effector genes occurs mainly during and after penetration of plant leaves and is associated with changes in histone modifications. Our work demonstrates the role of chromatin in shaping the expression of virulence components and, thereby, the interaction between fungal pathogens and their plant hosts.

Dynamic changes in transcription profiles are key for the success of pathogens in colonizing their hosts. In many pathogens, genes associated with virulence, such as effector genes, are located in regions of the genome that are rich in transposable elements and heterochromatin. The contribution of chromatin modifications to gene expression in pathogens remains largely unknown. Using a combination of a reporter gene-based approach and chromatin immunoprecipitation, we show that the heterochromatic environment of effector genes in the fungal plant pathogen Zymoseptoria tritici is a key regulator of their specific spatiotemporal expression patterns. Enrichment in trimethylated lysine 27 of histone H3 dictates the repression of effector genes in the absence of the host. Chromatin decondensation during host colonization, featuring a reduction in this repressive modification, indicates a major role for epigenetics in effector gene induction. Our results illustrate that chromatin modifications triggered during host colonization determine the specific expression profile of effector genes at the cellular level and, hence, provide new insights into the regulation of virulence in fungal plant pathogens.

Effector Biology of Biotrophic Plant Fungal Pathogens: Current Advances and Future Prospects

The interaction of fungal pathogens with their host requires a novel invading mechanism and the presence of various virulence-associated components responsible for promoting the infection. The small secretory proteins, explicitly known as effector proteins, are one of the prime mechanisms of host manipulation utilized by the pathogen to disarm the host. Several effector proteins are known to translocate from fungus to the plant cell for host manipulation. Many fungal effectors have been identified using genomic, transcriptomic, and bioinformatics approaches. Most of the effector proteins are devoid of any conserved signatures, and their prediction based on sequence homology is very challenging, therefore by combining the sequence consensus based upon machine learning features, multiple tools have also been developed for predicting apoplastic and cytoplasmic effectors. Various post-genomics approaches like transcriptomics of virulent isolates have also been utilized for identifying active consortia of effectors. Significant progress has been made in understanding biotrophic effectors; however, most of it is underway due to their complex interaction with host and complicated recognition and signaling networks. This review discusses advances, and challenges in effector identification and highlighted various features of the potential effector proteins and approaches for understanding their genetics and strategies for regulation.

Inhibition of jasmonate-mediated plant defences by the fungal metabolite higginsianin B

Infection of Arabidopsis thaliana by the ascomycete fungus Colletotrichum higginsianum is characterized by an early symptomless biotrophic phase followed by a destructive necrotrophic phase. The fungal genome contains 77 secondary metabolism-related biosynthetic gene clusters, whose expression during the infection process is tightly regulated. Deleting CclA, a chromatin regulator involved in the repression of some biosynthetic gene clusters through H3K4 trimethylation, allowed overproduction of three families of terpenoids and isolation of 12 different molecules. These natural products were tested in combination with methyl jasmonate, an elicitor of jasmonate responses, for their capacity to alter defence gene induction in Arabidopsis. Higginsianin B inhibited methyl jasmonate-triggered expression of the defence reporter VSP1p:GUS, suggesting it may block bioactive jasmonoyl isoleucine (JA-Ile) synthesis or signalling in planta. Using the JA-Ile sensor Jas9-VENUS, we found that higginsianin B, but not three other structurally related molecules, suppressed JA-Ile signalling by preventing the degradation of JAZ proteins, the repressors of jasmonate responses. Higginsianin B likely blocks the 26S proteasome-dependent degradation of JAZ proteins because it inhibited chymotrypsin- and caspase-like protease activities. The inhibition of target degradation by higginsianin B also extended to auxin signalling, as higginsianin B treatment reduced auxin-dependent expression of DR5p:GUS. Overall, our data indicate that specific fungal secondary metabolites can act similarly to protein effectors to subvert plant immune and developmental responses.

The Many Questions about Mini Chromosomes in Colletotrichum spp

Many fungal pathogens carry accessory regions in their genome, which are not required for vegetative fitness. Often, although not always, these regions occur as relatively small chromosomes in different species. Such mini chromosomes appear to be a typical feature of many filamentous plant pathogens. Since these regions often carry genes coding for effectors or toxin-producing enzymes, they may be directly related to virulence of the respective pathogen. In this review, we outline the situation of small accessory chromosomes in the genus Colletotrichum, which accounts for ecologically important plant diseases. We summarize which species carry accessory chromosomes, their gene content, and chromosomal makeup. We discuss the large variation in size and number even between different isolates of the same species, their potential roles in host range, and possible mechanisms for intra- and interspecies exchange of these interesting genetic elements.

Colletotrichum: species complexes, lifestyle, and peculiarities of some sources of genetic variability

The genus Colletotrichum comprises species with different lifestyles but is mainly known for phytopathogenic species that infect crops of agronomic relevance causing considerable losses. The fungi of the genus Colletotrichum are distributed in species complexes and within each complex some species have particularities regarding their lifestyle. The most commonly found and described lifestyles in Colletotrichum are endophytic and hemibiotrophic phytopathogenic. Several of these phytopathogenic species show wide genetic variability, which makes long-term maintenance of resistance in plants difficult. Different mechanisms may play an important role in the emergence of genetic variants but are not yet fully understood in this genus. These mechanisms include heterokaryosis, a parasexual cycle, sexual cycle, transposable element activity, and repeat-induced point mutations. This review provides an overview of the genus Colletotrichum, the species complexes described so far and the most common lifestyles in the genus, with a special emphasis on the mechanisms that may be responsible, at least in part, for the emergence of new genotypes under field conditions.

A detailed in silico analysis of secondary metabolite biosynthesis clusters in the genome of the broad host range plant pathogenic fungus Sclerotinia sclerotiorum

BACKGROUND

The broad host range pathogen Sclerotinia sclerotiorum infects over 400 plant species and causes substantial yield losses in crops worldwide. Secondary metabolites are known to play important roles in the virulence of plant pathogens, but little is known about the secondary metabolite repertoire of S. sclerotiorum. In this study, we predicted secondary metabolite biosynthetic gene clusters in the genome of S. sclerotiorum and analysed their expression during infection of Brassica napus using an existing transcriptome data set. We also investigated their sequence diversity among a panel of 25 previously published S. sclerotiorum isolate genomes.

RESULTS

We identified 80 putative secondary metabolite clusters. Over half of the clusters contained at least three transcriptionally coregulated genes. Comparative genomics revealed clusters homologous to clusters in the closely related plant pathogen Botrytis cinerea for production of carotenoids, hydroxamate siderophores, DHN melanin and botcinic acid. We also identified putative phytotoxin clusters that can potentially produce the polyketide sclerin and an epipolythiodioxopiperazine. Secondary metabolite clusters were enriched in subtelomeric genomic regions, and those containing paralogues showed a particularly strong association with repeats. The positional bias we identified was borne out by intraspecific comparisons that revealed putative secondary metabolite genes suffered more presence / absence polymorphisms and exhibited a significantly higher sequence diversity than other genes.

CONCLUSIONS

These data suggest that S. sclerotiorum produces numerous secondary metabolites during plant infection and that their gene clusters undergo enhanced rates of mutation, duplication and recombination in subtelomeric regions. The microevolutionary regimes leading to S. sclerotiorum secondary metabolite diversity have yet to be elucidated. Several potential phytotoxins documented in this study provide the basis for future functional analyses.

A draft genome of a field-collected Steinernema feltiae strain NW

Advances in sequencing technologies have accelerated our understanding of the complex genetic network of organisms and genomic divergences that are linked to evolutionary processes. While many model organisms and laboratory strains have been sequenced, wild populations are underrepresented in the growing list of sequenced genomes. Here, we present a de novo assembly of Steinernema feltiae, strain NW, collected from a working agricultural field in south central Washington, USA. Leveraging Pacific Biosciences (PacBio) long reads, we sequenced strain NW to a high depth (99×). The resulting de novo assembly is significantly larger than the previous assembly generated from the laboratory strain SN, with a noticeable improvement in continuity and completeness. Comparative analysis of two assemblies revealed numerous single nucleotide polymorphisms (SNPs), breakpoints, and indels present between the two genomes. This alternative genome resource and annotation could benefit the research community to examine the genetic foundation of evolutionary processes as well as genomic variation among conspecific populations.

Assembly, annotation, and comparison of Macrophomina phaseolina isolates from strawberry and other hosts

Background

Macrophomina phaseolina is a fungal plant pathogen with a broad host range, but one genotype was shown to exhibit host preference/specificity on strawberry. This pathogen lacked a high-quality genome assembly and annotation, and little was known about genomic differences among isolates from different hosts.

Results

We used PacBio sequencing and Hi-C scaffolding to provide nearly complete genome assemblies for M. phaseolina isolates representing the strawberry-specific genotype and another genotype recovered from alfalfa. The strawberry isolate had 59 contigs/scaffolds with an N50 of 4.3 Mb. The isolate from alfalfa had an N50 of 5.0 Mb and 14 nuclear contigs with half including telomeres. Both genomes were annotated with MAKER using transcript evidence generated in this study with over 13,000 protein-coding genes predicted. Unique groups of genes for each isolate were identified when compared to closely related fungal species. Structural comparisons between the isolates reveal large-scale rearrangements including chromosomal inversions and translocations. To include isolates representing a range of pathogen genotypes, an additional 30 isolates were sequenced with Illumina, assembled, and compared to the strawberry genotype assembly. Within the limits of comparing Illumina and PacBio assemblies, no conserved structural rearrangements were identified among the isolates from the strawberry genotype compared to those from other hosts, but some candidate genes were identified that were largely present in isolates of the strawberry genotype and absent in other genotypes.

Conclusions

High-quality reference genomes of M. phaseolina have allowed for the identification of structural changes associated with a genotype that has a host preference toward strawberry and will enable future comparative genomics studies. Having more complete assemblies allows for structural rearrangements to be more fully assessed and ensures a greater representation of all the genes. Work with Illumina data from additional isolates suggests that some genes are predominately present in isolates of the strawberry genotype, but additional work is needed to confirm the role of these genes in pathogenesis. Additional work is also needed to complete the scaffolding of smaller contigs identified in the strawberry genotype assembly and to determine if unique genes in the strawberry genotype play a role in pathogenicity.

Colletotrichum shisoi sp. nov., an anthracnose pathogen of Perilla frutescens in Japan: molecular phylogenetic, morphological and genomic evidence

Species of the fungal genus Colletotrichum are among the most devastating pathogens of agricultural crops in the world. Based on DNA sequence data (ITS, GAPDH, CHS-1, ACT, TUB2) and morphology, we revealed Colletotrichum isolates infecting the oil crop Perilla frutescens, commonly known as shiso, to represent a previously unknown species of the C. destructivum species complex and described it as C. shisoi. We found that C. shisoi appears to be able to adopt a hemibiotrophic lifestyle, characterised by the formation of biotrophic hyphae followed by severe necrotic lesions on P. frutescens, but is less virulent on Arabidopsis, compared to its close relative C. higginsianum which also belongs to the C. destructivum species complex. The genome of C. shisoi was sequenced, annotated and its predicted proteome compared with four other Colletotrichum species. The predicted proteomes of C. shisoi and C. higginsianum, share many candidate effectors, which are small, secreted proteins that may contribute to infection. Interestingly, C. destructivum species complex-specific secreted proteins showed evidence of increased diversifying selection which may be related to their host specificities.

Small RNAs from the plant pathogenic fungus Sclerotinia sclerotiorum highlight host candidate genes associated with quantitative disease resistance

Fungal plant pathogens secrete effector proteins and metabolites to cause disease. Additionally, some species transfer small RNAs (sRNAs) into plant cells to silence host mRNAs through complementary base pairing and suppress plant immunity. The fungus Sclerotinia sclerotiorum infects over 600 plant species, but little is known about the molecular processes that govern interactions with its many hosts. In particular, evidence for the production of sRNAs by S. sclerotiorum during infection is lacking. We sequenced sRNAs produced by S. sclerotiorum in vitro and during infection of two host species, Arabidopsis thaliana and Phaseolus vulgaris. We found that S. sclerotiorum produces at least 374 distinct highly abundant sRNAs during infection, mostly originating from repeat-rich plastic genomic regions. We predicted the targets of these sRNAs in A. thaliana and found that these genes were significantly more down-regulated during infection than the rest of the genome. Predicted targets of S. sclerotiorum sRNAs in A. thaliana were enriched for functional domains associated with plant immunity and were more strongly associated with quantitative disease resistance in a genome-wide association study (GWAS) than the rest of the genome. Mutants in A. thaliana predicted sRNA target genes SERK2 and SNAK2 were more susceptible to S. sclerotiorum than wild-type, suggesting that S. sclerotiorum sRNAs may contribute to the silencing of immune components in plants. The prediction of fungal sRNA targets in plant genomes can be combined with other global approaches, such as GWAS, to assist in the identification of plant genes involved in quantitative disease resistance.

Genomic Plasticity Mediated by Transposable Elements in the Plant Pathogenic Fungus Colletotrichum higginsianum

Phytopathogen genomes are under constant pressure to change, as pathogens are locked in an evolutionary arms race with their hosts, where pathogens evolve effector genes to manipulate their hosts, whereas the hosts evolve immune components to recognize the products of these genes. Colletotrichum higginsianum (Ch), a fungal pathogen with no known sexual morph, infects Brassicaceae plants including Arabidopsis thaliana. Previous studies revealed that Ch differs in its virulence toward various Arabidopsis thaliana ecotypes, indicating the existence of coevolutionary selective pressures. However, between-strain genomic variations in Ch have not been studied. Here, we sequenced and assembled the genome of a Ch strain, resulting in a highly contiguous genome assembly, which was compared with the chromosome-level genome assembly of another strain to identify genomic variations between strains. We found that the two closely related strains vary in terms of large-scale rearrangements, the existence of strain-specific regions, and effector candidate gene sets and that these variations are frequently associated with transposable elements (TEs). Ch has a compartmentalized genome consisting of gene-sparse, TE-dense regions with more effector candidate genes and gene-dense, TE-sparse regions harboring conserved genes. Additionally, analysis of the conservation patterns and syntenic regions of effector candidate genes indicated that the two strains vary in their effector candidate gene sets because of de novo evolution, horizontal gene transfer, or gene loss after divergence. Our results reveal mechanisms for generating genomic diversity in this asexual pathogen, which are important for understanding its adaption to hosts.

Near-Complete Genome Assembly of Alternaria brassicae-A Necrotrophic Pathogen of Brassica Crops

Alternaria brassicae, a necrotrophic fungal pathogen, causes Alternaria blight, an important disease of brassica crops. Although many Alternaria spp. have been sequenced, no genome information is available for A. brassicae, a monotypic lineage within the Alternaria genus. A highly contiguous genome assembly of A. brassicae has been generated using Nanopore MinION sequencing with an N50 of 2.98 Mb, yielding nine full chromosome-level sequences. This study adds to the current genomic resources available for the genus Alternaria and will provide opportunities for further analyses to unravel the mechanisms underlying pathogenicity of this important necrotrophic pathogen.

Nonproteinaceous effectors: the terra incognita of plant-fungal interactions

Molecular plant-fungal interaction studies have mainly focused on small secreted protein effectors. However, accumulating evidence shows that numerous fungal secondary metabolites are produced at all stages of plant colonization, especially during early asymptomatic/biotrophic phases. The discovery of fungal small RNAs targeting plant transcripts has expanded the fungal repertoire of nonproteinaceous effectors even further. The challenge now is to develop specific functional methods to fully understand the biological roles of these effectors. Studies on fungal extracellular vesicles are also needed because they could be the universal carriers of all kinds of fungal effectors. With this review, we aim to stimulate the nonproteinaceous effector research field to move from descriptive to functional studies, which should bring a paradigm shift in plant-fungal interactions.

H3K4 trimethylation by CclA regulates pathogenicity and the production of three families of terpenoid secondary metabolites in Colletotrichum higginsianum

The role of histone 3 lysine 4 (H3K4) methylation is poorly understood in plant pathogenic fungi. Here, we analysed the function of CclA, a subunit of the COMPASS complex mediating H3K4 methylation, in the brassica anthracnose pathogen Colletotrichum higginsianum. We show that CclA is required for full genome-wide H3K4 trimethylation. The deletion of cclA strongly reduced mycelial growth, asexual sporulation and spore germination but did not impair the morphogenesis of specialized infection structures (appressoria and biotrophic hyphae). Virulence of the ΔcclA mutant on plants was strongly attenuated, associated with a marked reduction in appressorial penetration ability on both plants and inert cellophane membranes. The secondary metabolite profile of the ΔcclA mutant was greatly enriched compared to that of the wild type, with three different families of terpenoid compounds being overproduced by the mutant, namely the colletochlorins, higginsianins and sclerosporide. These included five novel molecules that were produced exclusively by the ΔcclA mutant: colletorin D, colletorin D acid, higginsianin C, 13-epi-higginsianin C and sclerosporide. Taken together, our findings indicate that H3K4 trimethylation plays a critical role in regulating fungal growth, development, pathogenicity and secondary metabolism in C. higginsianum.

Comparative genome analysis indicates high evolutionary potential of pathogenicity genes in Colletotrichum tanaceti

Colletotrichum tanaceti is an emerging foliar fungal pathogen of commercially grown pyrethrum (Tanacetum cinerariifolium). Despite being reported consistently from field surveys in Australia, the molecular basis of pathogenicity of C. tanaceti on pyrethrum is unknown. Herein, the genome of C. tanaceti (isolate BRIP57314) was assembled de novo and annotated using transcriptomic evidence. The inferred putative pathogenicity gene suite of C. tanaceti comprised a large array of genes encoding secreted effectors, proteases, CAZymes and secondary metabolites. Comparative analysis of its putative pathogenicity gene profiles with those of closely related species suggested that C. tanaceti likely has additional hosts to pyrethrum. The genome of C. tanaceti had a high repeat content and repetitive elements were located significantly closer to genes inferred to influence pathogenicity than other genes. These repeats are likely to have accelerated mutational and transposition rates in the genome, resulting in a rapid evolution of certain CAZyme families in this species. The C. tanaceti genome showed strong signals of Repeat Induced Point (RIP) mutation which likely caused its bipartite nature consisting of distinct gene-sparse, repeat and A-T rich regions. Pathogenicity genes within these RIP affected regions were likely to have a higher evolutionary rate than the rest of the genome. This "two-speed" genome phenomenon in certain Colletotrichum spp. was hypothesized to have caused the clustering of species based on the pathogenicity genes, to deviate from taxonomic relationships. The large repertoire of pathogenicity factors that potentially evolve rapidly due to the plasticity of the genome, indicated that C. tanaceti has a high evolutionary potential. Therefore, C. tanaceti poses a high-risk to the pyrethrum industry. Knowledge of the evolution and diversity of the putative pathogenicity genes will facilitate future research in disease management of C. tanaceti and other Colletotrichum spp.

Deleting a Chromatin Remodeling Gene Increases the Diversity of Secondary Metabolites Produced by Colletotrichum higginsianum

Colletotrichum higginsianum is the causal agent of crucifer anthracnose disease, responsible for important economic losses in Brassica crops. A mutant lacking the CclA subunit of the COMPASS complex was expected to undergo chromatin decondensation and the activation of cryptic secondary metabolite biosynthetic gene clusters. Liquid-state fermentation of the Δ cclA mutant coupled with in situ solid-phase extraction led to the production of three families of compounds, namely, colletorin and colletochlorin derivatives with two new representatives, colletorin D (1) and colletorin D acid (2), the diterpenoid α-pyrone higginsianin family with two new analogues, higginsianin C (3) and 13- epi-higginsianin C (4), and sclerosporide (5) coupling a sclerosporin moiety with dimethoxy inositol.

Giant Transposons in Eukaryotes: Is Bigger Better?

Transposable elements (TEs) are ubiquitous in both prokaryotes and eukaryotes, and the dynamic character of their interaction with host genomes brings about numerous evolutionary innovations and shapes genome structure and function in a multitude of ways. In traditional classification systems, TEs are often being depicted in simplistic ways, based primarily on the key enzymes required for transposition, such as transposases/recombinases and reverse transcriptases. Recent progress in whole-genome sequencing and long-read assembly, combined with expansion of the familiar range of model organisms, resulted in identification of unprecedentedly long transposable units spanning dozens or even hundreds of kilobases, initially in prokaryotic and more recently in eukaryotic systems. Here, we focus on such oversized eukaryotic TEs, including retrotransposons and DNA transposons, outline their complex and often combinatorial nature and closely intertwined relationship with viruses, and discuss their potential for participating in transfer of long stretches of DNA in eukaryotes.

Establishment of a selection marker recycling system for sequential transformation of the plant-pathogenic fungus Colletotrichum orbiculare

Genome sequencing of pathogenic fungi has revealed the presence of various effectors that aid pathogen invasion by the manipulation of plant immunity. Effectors are often individually dispensable because of duplication and functional redundancy as a result of the arms race between host plants and pathogens. To study effectors that have functional redundancy, multiple gene disruption is often required. However, the number of selection markers that can be used for gene targeting is limited. Here, we established a marker recycling system that allows the use of the same selection marker in successive transformations in the model fungal pathogen Colletotrichum orbiculare, a causal agent of anthracnose disease in plants belonging to the Cucurbitaceae. We identified two C. orbiculare homologues of yeast URA3/pyrG, designated as URA3A and URA3B, which can be used as selection markers on medium with no uridine. The gene can then be removed from the genome via homologous recombination when the fungus is grown in the presence of 5-fluoroorotic acid (5-FOA), a chemical that is converted into a toxin by URA3 activity. The ura3a/b double mutants showed auxotrophy for uridine and insensitivity to 5-FOA. Using the ura3a/b mutants, transformation with the URA3B marker and its removal were successfully applied to disrupt the virulence-related gene, PKS1. The pks1 mutants showed a reduction in virulence, demonstrating that the method can be used to study virulence-related genes in C. orbiculare. The establishment of a URA3-based marker recycling system in plant-pathogenic fungi enables the genetic analysis of multiple genes that have redundant functions, including effector genes.

The Modern View of B Chromosomes Under the Impact of High Scale Omics Analyses

Supernumerary B chromosomes (Bs) are extra karyotype units in addition to A chromosomes, and are found in some fungi and thousands of animals and plant species. Bs are uniquely characterized due to their non-Mendelian inheritance, and represent one of the best examples of genomic conflict. Over the last decades, their genetic composition, function and evolution have remained an unresolved query, although a few successful attempts have been made to address these phenomena. A classical concept based on cytogenetics and genetics is that Bs are selfish and abundant with DNA repeats and transposons, and in most cases, they do not carry any function. However, recently, the modern quantum development of high scale multi-omics techniques has shifted B research towards a new-born field that we call "B-omics". We review the recent literature and add novel perspectives to the B research, discussing the role of new technologies to understand the mechanistic perspectives of the molecular evolution and function of Bs. The modern view states that B chromosomes are enriched with genes for many significant biological functions, including but not limited to the interesting set of genes related to cell cycle and chromosome structure. Furthermore, the presence of B chromosomes could favor genomic rearrangements and influence the nuclear environment affecting the function of other chromatin regions. We hypothesize that B chromosomes might play a key function in driving their transmission and maintenance inside the cell, as well as offer an extra genomic compartment for evolution.

Chromosome rearrangements shape the diversification of secondary metabolism in the cyclosporin producing fungus Tolypocladium inflatum

BACKGROUND

Genes involved in production of secondary metabolites (SMs) in fungi are exceptionally diverse. Even strains of the same species may exhibit differences in metabolite production, a finding that has important implications for drug discovery. Unlike in other eukaryotes, genes producing SMs are often clustered and co-expressed in fungal genomes, but the genetic mechanisms involved in the creation and maintenance of these secondary metabolite biosynthetic gene clusters (SMBGCs) remains poorly understood.

RESULTS

In order to address the role of genome architecture and chromosome scale structural variation in generating diversity of SMBGCs, we generated chromosome scale assemblies of six geographically diverse isolates of the insect pathogenic fungus Tolypocladium inflatum, producer of the multi-billion dollar lifesaving immunosuppressant drug cyclosporin, and utilized a Hi-C chromosome conformation capture approach to address the role of genome architecture and structural variation in generating intraspecific diversity in SMBGCs. Our results demonstrate that the exchange of DNA between heterologous chromosomes plays an important role in generating novelty in SMBGCs in fungi. In particular, we demonstrate movement of a polyketide synthase (PKS) and several adjacent genes by translocation to a new chromosome and genomic context, potentially generating a novel PKS cluster. We also provide evidence for inter-chromosomal recombination between nonribosomal peptide synthetases located within subtelomeres and uncover a polymorphic cluster present in only two strains that is closely related to the cluster responsible for biosynthesis of the mycotoxin aflatoxin (AF), a highly carcinogenic compound that is a major public health concern worldwide. In contrast, the cyclosporin cluster, located internally on chromosomes, was conserved across strains, suggesting selective maintenance of this important virulence factor for infection of insects.

CONCLUSIONS

This research places the evolution of SMBGCs within the context of whole genome evolution and suggests a role for recombination between chromosomes in generating novel SMBGCs in the medicinal fungus Tolypocladium inflatum.