First Draft Genome Sequence of a UK Strain (UK99) of Fusarium culmorum

Measuring air metagenomic diversity in an agricultural ecosystem

All species shed DNA during life or in death, providing an opportunity to monitor biodiversity via environmental DNA (eDNA). In recent years, combining eDNA, high-throughput sequencing technologies, bioinformatics, and increasingly complete sequence databases has promised a non-invasive and non-destructive environmental monitoring tool. Modern agricultural systems are often large monocultures and so are highly vulnerable to disease outbreaks. Pest and pathogen monitoring in agricultural ecosystems is key for efficient and early disease prevention, lower pesticide use, and better food security. Although the air is rich in biodiversity, it has the lowest DNA concentration of all environmental media and yet is the route for windborne spread of many damaging crop pathogens. Our work suggests that ecosystems can be monitored efficiently using airborne nucleic acid information. Here, we show that the airborne DNA of microbes can be recovered, shotgun sequenced, and taxonomically classified, including down to the species level. We show that by monitoring a field growing key crops we can identify the presence of agriculturally significant pathogens and quantify their changing abundance over a period of 1.5 months, often correlating with weather variables. We add to the evidence that aerial eDNA can be used as a source for biomonitoring in terrestrial ecosystems, specifically highlighting agriculturally relevant species and how pathogen levels correlate with weather conditions. Our ability to detect dynamically changing levels of species and strains highlights the value of airborne eDNA in agriculture, monitoring biodiversity changes, and tracking taxa of interest.

Barley shows reduced Fusarium head blight under drought and modular expression of differentially expressed genes under combined stress

Plants often face simultaneous abiotic and biotic stress conditions; however, physiological and transcriptional responses under such combined stress conditions are still not fully understood. Spring barley (Hordeum vulgare) is susceptible to Fusarium head blight (FHB), which is strongly affected by weather conditions. We therefore studied the potential influence of drought on FHB severity and plant responses in three varieties of different susceptibility. We found strongly reduced FHB severity in susceptible varieties under drought. The number of differentially expressed genes (DEGs) and strength of transcriptomic regulation reflected the concentrations of physiological stress markers such as abscisic acid or fungal DNA contents. Infection-related gene expression was associated with susceptibility rather than resistance. Weighted gene co-expression network analysis revealed 18 modules of co-expressed genes that reflected the pathogen- or drought-response in the three varieties. A generally infection-related module contained co-expressed genes for defence, programmed cell death, and mycotoxin detoxification, indicating that the diverse genotypes used a similar defence strategy towards FHB, albeit with different degrees of success. Further, DEGs showed co-expression in drought- or genotype-associated modules that correlated with measured phytohormones or the osmolyte proline. The combination of drought stress with infection led to the highest numbers of DEGs and resulted in a modular composition of the single-stress responses rather than a specific transcriptional output.

Determining the biocontrol capacities of Trichoderma spp. originating from Turkey on Fusarium culmorum by transcriptional and antagonistic analyses

In this study aiming to investigate potential fungal biocontrol agents for Fusarium culmorum, several isolates of Trichoderma spp. were evaluated for their antagonistic effects by means of transcriptional analyses. At first, 21 monosporic Trichoderma spp. isolates were obtained from natural wood debris and wood area soils in Manisa, Turkey. Trichoderma spp. Isolates were identified as belonging to four different species (T. atroviride, T. harzianum, T. koningii, and T. brevicompactum) by tef1-α sequencing. Then, the linear growth rate (LGR) of each species was calculated and determined to be in a range between 13.22 ± 0.71 mm/day (T. atroviride TR2) and 25.06 ± 1.45 mm/day (T. harzianum K30). Inter-simple sequence repeat (ISSR) genotyping validated the tef1-α sequencing results by presenting two sub-clusters in the dendrogram. We determined the genetically most similar (TR1 & TR2; 97.77%) and dissimilar (K9 & K17; 40.40%) individuals belonging to the same and different species, respectively. Dual sandwich culture tests (which are useful for antagonism studies) revealed that T. harzianum K21 (the least suppressive) and T. brevicompactum K26 (the most suppressive) isolates suppressed F. culmorum with growth rates of 3% and 46%, respectively. Expressions of genes previously associated with mycoparasitism-plant protection-secondary metabolism (nag1, tgf-1, and tmk-1) were tested by quantitative real-time polymerase chain reaction (qRT-PCR) in both those isolates. While there were no significant differences (p>0.05) in expression that were present in the K21 isolate, those three genes were upregulated with fold change values of 2.69 ± 0.26 (p<0.001), 2.23 ± 0.16 (p<0.001), and 5.38 ± 2.01 (p<0.05) in K26, meaning that the presence of significant alteration in the physiological processes of the fungus. Also, its mycoparasitism potential was tested on Triticum aestivum L. cv Basribey in planta, which was infected with the F. culmorum FcUK99 strain. Results of the trials, including specific plant growth parameters (weight or length of plantlets), confirmed the mycoparasitic potential of the isolate. It can be concluded that (i) nag1, tgf-1, and tmk-1 genes could be approved as reliable markers for evaluation of BCA capacities of Trichoderma spp. and (ii) the T. brevicompactum K26 strain can be suggested as a promising candidate for combating in F. culmorum diseases following the necessary procedures to ensure it is non-hazardous and safe.

Genomic and Metabolomic Analysis of the Endophytic Fungus Fusarium sp. VM-40 Isolated from the Medicinal Plant Vinca minor

The genus Fusarium is well-known to comprise many pathogenic fungi that affect cereal crops worldwide, causing severe damage to agriculture and the economy. In this study, an endophytic fungus designated Fusarium sp. VM-40 was isolated from a healthy specimen of the traditional European medicinal plant Vinca minor. Our morphological characterization and phylogenetic analysis reveal that Fusarium sp. VM-40 is closely related to Fusarium paeoniae, belonging to the F. tricinctum species complex (FTSC), the genomic architecture and secondary metabolite profile of which have not been investigated. Thus, we sequenced the whole genome of Fusarium sp. VM-40 with the new Oxford Nanopore R10.4 flowcells. The assembled genome is 40 Mb in size with a GC content of 47.72%, 15 contigs (≥50,000 bp; N 50~4.3 Mb), and 13,546 protein-coding genes, 691 of which are carbohydrate-active enzyme (CAZyme)-encoding genes. We furthermore predicted a total of 56 biosynthetic gene clusters (BGCs) with antiSMASH, 25 of which showed similarity with known BGCs. In addition, we explored the potential of this fungus to produce secondary metabolites through untargeted metabolomics. Our analyses reveal that this fungus produces structurally diverse secondary metabolites of potential pharmacological relevance (alkaloids, peptides, amides, terpenoids, and quinones). We also employed an epigenetic manipulation method to activate cryptic BGCs, which led to an increased abundance of several known compounds and the identification of several putative new compounds. Taken together, this study provides systematic research on the whole genome sequence, biosynthetic potential, and metabolome of the endophytic fungus Fusarium sp. VM-40.

Independent Evolution Has Led to Distinct Genomic Signatures in Dutch Elm Disease-Causing Fungi and Other Vascular Wilts-Causing Fungal Pathogens

Vascular wilts are important diseases caused by plant pathogenic fungi that result in the rapid death of their plant hosts. This is due to a systemic defense mechanism whereby the plant induces the compartmentalization of the infected vascular system in order to reduce the propagation of the fungus. The ascomycete class Sordariomycetes contains several species that cause vascular wilts in diverse plant hosts, and they can be classified into four taxonomic orders. The genetic mechanisms of pathogenesis have already been investigated in Fusarium and Verticillium species, but they have not yet been compared with other well-known wilt-causing species, especially fungi causing oak wilt or Dutch elm disease (DED). Here we analyzed 20 whole genome assemblies of wilt-causing fungi together with 56 other species using phylogenetic approaches to trace expansions and contractions of orthologous gene families and gene classes related to pathogenicity. We found that the wilt-causing pathogens evolved seven times, experiencing the largest fold changes in different classes of genes almost every time. However, some similarities exist across groups of wilt pathogens, particularly in Microascales and Ophiostomatales, and these include the common gains and losses of genes that make up secondary metabolite clusters (SMC). DED pathogens do not experience large-scale gene expansions, with most of the gene classes, except for some SMC families, reducing in number. We also found that gene family expansions in the most recent common ancestors of wilt pathogen groups are enriched for carbohydrate metabolic processes. Our study shows that wilt-causing species evolve primarily through distinct changes in their repertoires of pathogenicity-related genes and that there is the potential importance of carbohydrate metabolism genes for regulating osmosis in those pathogens that penetrate the plant vascular system.

Insights on KP4 Killer Toxin-like Proteins of Fusarium Species in Interspecific Interactions

KP4 killer toxins are secreted proteins that inhibit cell growth and induce cell death in target organisms. In Fusarium graminearum, KP4-like (KP4L) proteins contribute to fungal virulence in wheat seedling rot and are expressed during Fusarium head blight development. However, fungal KP4L proteins are also hypothesized to support fungal antagonism by permeabilizing cell walls of competing fungi to enable penetration of toxic compounds. Here, we report the differential expression patterns of F. graminearum KP4L genes (Fgkp4l-1, -2, -3 and -4) in a competitive interaction, using Trichoderma gamsii as the antagonist. The results from dual cultures indicate that Fgkp4l-3 and Fgkp4l-4 could participate in the recognition at the distance of the antagonist, while all Fgkp4l genes were highly activated in the pathogen during the physical interaction of both fungi. Only Fgkp4l-4 was up-regulated during the interaction with T. gamsii in wheat spikes. This suggests the KP4L proteins could participate in supporting F. graminearum interspecific interactions, even in living plant tissues. The distribution of KP4L orthologous within the genus Fusarium revealed they are more represented in species with broad host-plant range than in host-specific species. Phylogeny inferred provides evidence that KP4L genes evolved through gene duplications, gene loss and sequence diversification in the genus Fusarium.

Genome-wide association study for deoxynivalenol production and aggressiveness in wheat and rye head blight by resequencing 92 isolates of Fusarium culmorum

Background

Fusarium culmorum is an important pathogen causing head blight of cereals in Europe. This disease is of worldwide importance leading to reduced yield, grain quality, and contamination by mycotoxins. These mycotoxins are harmful for livestock and humans; therefore, many countries have strict regulatory limits for raw materials and processed food. Extensive genetic diversity is described among field populations of F. culmorum isolates for aggressiveness and production of the trichothecene mycotoxin deoxynivalenol (DON). However, the causes for this quantitative variation are not clear, yet. We analyzed 92 isolates sampled from different field populations in Germany, Russia, and Syria together with an international collection for aggressiveness and DON production in replicated field experiments at two locations in two years with two hosts, wheat and rye. The 30x coverage whole-genome resequencing of all isolates resulted in the identification of 130,389 high quality single nucleotide polymorphisms (SNPs) that were used for the first genome-wide association study in this phytopathogenic fungus.

Results

In wheat, 20 and 27 SNPs were detected for aggressiveness and DON content, respectively, of which 10 overlapped. Additionally, two different SNPs were significantly associated with aggressiveness in rye that were among those SNPs being associated with DON production in wheat. Most of the SNPs explained only a small proportion of genotypic variance (p_G), however, four SNPs were associated with major quantitative trait loci (QTLs) with p_G ranging from 12 to 48%. The QTL with the highest p_G was involved in DON production and associated with a SNP most probably located within the Tri4 gene.

Conclusions

The diversity of 92 isolates of F. culmorum were captured using a heuristic approach. Key phenotypic traits, SNPs, and candidate genes underlying aggressiveness and DON production were identified. Clearly, many QTLs are responsible for aggressiveness and DON content in wheat, both traits following a quantitative inheritance. Several SNPs involved in DON metabolism, among them the Tri4 gene of the trichothecene pathway, were inferred as important source of variation in fungal aggressiveness. Using this information underlying the phenotypic variation will be of paramount importance in evaluating strategies for successful resistance breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07931-5.

Mycotoxin Production in Fusarium According to Contemporary Species Concepts

Fusarium is one of the most important genera of plant-pathogenic fungi in the world and arguably the world's most important mycotoxin-producing genus. Fusarium species produce a staggering array of toxic metabolites that contribute to plant disease and mycotoxicoses in humans and other animals. A thorough understanding of the mycotoxin potential of individual species is crucial for assessing the toxicological risks associated with Fusarium diseases. There are thousands of reports of mycotoxin production by various species, and there have been numerous attempts to summarize them. These efforts have been complicated by competing classification systems based on morphology, sexual compatibility, and phylogenetic relationships. The current depth of knowledge of Fusarium genomes and mycotoxin biosynthetic pathways provides insights into how mycotoxin production is distributedamong species and multispecies lineages (species complexes) in the genus as well as opportunities to clarify and predict mycotoxin risks connected with known and newly described species. Here, we summarize mycotoxin production in the genus Fusarium and how mycotoxin risk aligns with current phylogenetic species concepts.

Phylogeny and Mycotoxin Profile of Pathogenic Fusarium Species Isolated from Sudden Decline Syndrome and Leaf Wilt Symptoms on Date Palms (Phoenix dactylifera) in Tunisia

In 2017–2018, extensive symptoms of sudden decline and fruit rot were observed on date palms in southern Tunisia. Samples of diseased plants were randomly collected in six localities. Based on morphological identification, Fusarium was the most frequent fungal genus detected. A sequencing of translation elongation factor, calmodulin, and second largest subunit of RNA polymerase II genes was used to identify 63 representative Fusarium strains at species level and investigate their phylogenetic relationships. The main species detected was Fusarium proliferatum, and at a much lesser extent, Fusarium brachygibbosum, Fusarium caatingaense, Fusarium clavum, Fusarium incarnatum, and Fusarium solani. Pathogenicity on the Deglet Nour variety plantlets and the capability to produce mycotoxins were also assessed. All Fusarium species were pathogenic complying Koch’s postulates. Fusarium proliferatum strains produced mainly fumonisins (FBs), beauvericin (BEA), and, to a lesser extent, enniatins (ENNs) and moniliformin (MON). All F. brachygibbosum strains produced low levels of BEA, diacetoxyscirpenol, and neosolaniol; two strains produced also T-2 toxin, and a single strain produced HT-2 toxin. Fusarium caatingaense, F. clavum, F. incarnatum produced only BEA. Fusarium solani strains produced MON, BEA, and ENNs. This work reports for the first time a comprehensive multidisciplinary study of Fusarium species on date palms, concerning both phytopathological and food safety issues.

The adaptation of Fusarium culmorum to DMI Fungicides Is Mediated by Major Transcriptome Modifications in Response to Azole Fungicide, Including the Overexpression of a PDR Transporter (FcABC1)

Fusarium culmorum is a fungal pathogen causing economically important diseases on a variety of crops. Fungicides can be applied to control this species with triazoles being the most efficient molecules. F. culmorum strains resistant to these molecules have been reported, but the underlying resistance mechanisms remain unknown. In this study, a tebuconazole-adapted F. culmorum strain was developed with a level of fitness similar to its parental strain. The adapted strain showed cross-resistance to all demethylation inhibitors (DMIs), but not to other classes of fungicides tested. RNA-Seq analysis revealed high transcriptomic differences between the resistant strain and its parental strain after tebuconazole treatment. Among these changes, FcABC1 (FCUL_06717), a pleiotropic drug resistance transporter, had a 30-fold higher expression level upon tebuconazole treatment in the adapted strains as compared to the wild-type strain. The implication of this transporter in triazole resistance was subsequently confirmed in field strains harboring distinct levels of sensitivity to triazoles. FcABC1 is present in other species/genera, including F. graminearum in which it is known to be necessary for azole resistance. No difference in FcABC1 sequences, including the surrounding regions, were found when comparing the resistant strain to the wild-type strain. Fusarium culmorum is therefore capable to adapt to triazole pressure by overexpressing a drug resistance transporter when submitted to triazoles and the same mechanism is anticipated to occur in other species.

Development of a Highly Sensitive FcMito qPCR Assay for the Quantification of the Toxigenic Fungal Plant Pathogen Fusarium culmorum

Fusarium culmorum is a ubiquitous, soil-borne fungus (ascomycete) causing foot and root rot and Fusarium head blight on cereals. It is responsible for yield and quality losses as well as grain contamination with mycotoxins, which are a potential health hazard. An extremely sensitive mitochondrial-based qPCR assay (FcMito qPCR) for quantification of F. culmorum was developed in this study. To provide specificity, the FcMito assay was successfully validated against 85 F. culmorum strains and 53 isolates of 30 other fungal species. The assay efficiency and sensitivity were evaluated against different F. culmorum strains with various amounts of pure fungal DNA and in the presence of background wheat DNA. The results demonstrated the high efficiency of the assay (97.2⁻106.0%, R²-values > 0.99). It was also shown that, in the presence of background DNA, 0.01 pg of fungal template could be reliably quantified. The FcMito assay was used to quantify F. culmorum DNA using 108 grain samples with different trichothecene levels. A significant positive correlation was found between fungal DNA quantity and the total trichothecene content. The obtained results showed that the sensitivity of the FcMito assay was much higher than the nuclear-based qPCR assay for F. culmorum.

FcRav2, a gene with a ROGDI domain involved in Fusarium head blight and crown rot on durum wheat caused by Fusarium culmorum

Fusarium culmorum is a soil-borne fungal pathogen which causes foot and root rot and Fusarium head blight on small-grain cereals, in particular wheat and barley. It causes significant yield and quality losses and results in the contamination of kernels with type B trichothecene mycotoxins. Our knowledge of the pathogenicity factors of this fungus is still limited. A transposon tagging approach based on the mimp1/impala double-component system has allowed us to select a mutant altered in multiple metabolic and morphological processes, trichothecene production and virulence. The flanking regions of mimp1 were used to seek homologies in the F. culmorum genome, and revealed that mimp1 had reinserted within the last exon of a gene encoding a hypothetical protein of 318 amino acids which contains a ROGDI-like leucine zipper domain, supposedly playing a protein-protein interaction or regulatory role. By functional complementation and bioinformatic analysis, we characterized the gene as the yeast Rav2 homologue, confirming the high level of divergence in multicellular fungi. Deletion of FcRav2 or its orthologous gene in F. graminearum highlighted its ability to influence a number of functions, including virulence, trichothecene type B biosynthesis, resistance to azoles and resistance to osmotic and oxidative stress. Our results indicate that the FcRav2 protein (and possibly the RAVE complex as a whole) may become a suitable target for new antifungal drug development or the plant-mediated resistance response in filamentous fungi of agricultural interest.

Candidate Genes for Aggressiveness in a Natural Fusarium culmorum Population Greatly Differ between Wheat and Rye Head Blight

Fusarium culmorum is one of the species causing Fusarium head blight (FHB) in cereals in Europe. We aimed to investigate the association between the nucleotide diversity of ten F. culmorum candidate genes and field ratings of aggressiveness in winter rye. A total of 100 F. culmorum isolates collected from natural infections were phenotyped for FHB at two locations and two years. Variance components for aggressiveness showed significant isolate and isolate-by-environment variance, as expected for quantitative host-pathogen interactions. Further analysis of the isolate-by-environment interaction revealed the dominant role of the isolate-by-year over isolate-by-location interaction. One single-nucleotide polymorphism (SNP) in the cutinase (CUT) gene was found to be significantly (p < 0.001) associated with aggressiveness and explained 16.05% of the genotypic variance of this trait in rye. The SNP was located 60 base pairs before the start codon, which suggests a role in transcriptional regulation. Compared to a previous study in winter wheat with the same nucleotide sequences, a larger variation of pathogen aggressiveness on rye was found and a different candidate gene was associated with pathogen aggressiveness. This is the first report on the association of field aggressiveness and a host-specific candidate gene codifying for a protein that belongs to the secretome in F. culmorum.

Candidate gene based association mapping in Fusarium culmorum for field quantitative pathogenicity and mycotoxin production in wheat

Background

Quantitative traits are common in nature, but quantitative pathogenicity has received only little attention in phytopathology. In this study, we used 100 Fusarium culmorum isolates collected from natural field environments to assess their variation for two quantitative traits, aggressiveness and deoxynivalenol (DON) production on wheat plants grown in four different field environments (location-year combinations). Seventeen Fusarium graminearum pathogenicity candidate genes were assessed for their effect on the aggressiveness and DON production of F. culmorum under field conditions.

Results

For both traits, genotypic variance among isolates was high and significant while the isolate-by-environment interaction was also significant, amounting to approximately half of the genotypic variance. Among the studied candidate genes, the mitogen-activated protein kinase (MAPK) HOG1 was found to be significantly associated with aggressiveness and deoxynivalenol (DON) production, explaining 10.29 and 6.05% of the genotypic variance, respectively.

Conclusions

To the best of our knowledge, this is the first report of a protein kinase regulator explaining differences in field aggressiveness and mycotoxin production among individuals from natural populations of a plant pathogen.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0511-9) contains supplementary material, which is available to authorized users.

Annotation of Fusarium graminearum (PH-1) Version 5.0

Fusarium graminearum floral infections are a major risk to the global supply of safe cereal grains. We report updates to the PH-1 reference genome and significant improvements to the annotation. Changes include introduction of legacy annotation identifiers, new gene models, secretome and effectorP predictions, and inclusion of extensive untranslated region (UTR) annotations.