Comparative genomics of geographically distant Fusarium fujikuroi isolates revealed two distinct pathotypes correlating with secondary metabolite profiles

Fusarium fujikuroi causes bakanae (“foolish seedling”) disease of rice which is characterized by hyper-elongation of seedlings resulting from production of gibberellic acids (GAs) by the fungus. This plant pathogen is also known for production of harmful mycotoxins, such as fusarins, fusaric acid, apicidin F and beauvericin. Recently, we generated the first de novo genome sequence of F. fujikuroi strain IMI 58289 combined with extensive transcriptional, epigenetic, proteomic and chemical product analyses. GA production was shown to provide a selective advantage during infection of the preferred host plant rice. Here, we provide genome sequences of eight additional F. fujikuroi isolates from distant geographic regions. The isolates differ in the size of chromosomes, most likely due to variability of subtelomeric regions, the type of asexual spores (microconidia and/or macroconidia), and the number and expression of secondary metabolite gene clusters. Whilst most of the isolates caused the typical bakanae symptoms, one isolate, B14, caused stunting and early withering of infected seedlings. In contrast to the other isolates, B14 produced no GAs but high amounts of fumonisins during infection on rice. Furthermore, it differed from the other isolates by the presence of three additional polyketide synthase (PKS) genes (PKS40, PKS43, PKS51) and the absence of the F. fujikuroi-specific apicidin F (NRPS31) gene cluster. Analysis of additional field isolates confirmed the strong correlation between the pathotype (bakanae or stunting/withering), and the ability to produce either GAs or fumonisins. Deletion of the fumonisin and fusaric acid-specific PKS genes in B14 reduced the stunting/withering symptoms, whereas deletion of the PKS51 gene resulted in elevated symptom development. Phylogenetic analyses revealed two subclades of F. fujikuroi strains according to their pathotype and secondary metabolite profiles.


Introduction
The heterothallic ascomycete Fusarium fujikuroi Nirenberg is a member of the Fusarium fujikuroi species complex (FFC), a monophyletic lineage which includes at least eleven mating populations (MPs A-K) that are sexually infertile with one another, and numerous distinct anamorphic species [1].
F. fujikuroi (MP-C) is the causal agent of the rice disease bakanae ("foolish seedlings"), one of the most notorious seed-borne diseases with increasing economic importance in the major rice-growing countries in the world, including all rice-growing Asian and African countries, California, and more recently, Italy [1][2][3]. The fungus was one of the first fungal pathogens to be described, and bakanae is one of the oldest known diseases of rice being reported more than 100 years ago by Japanese scientists [4,5]. The most prominent symptoms of the disease are chlorotic, elongated and thin seedlings that are often several inches taller than healthy plants, and empty panicles leading to yield losses ranging from ca. 30-95% [6][7][8][9]. However, not all infected seedlings show the bakanae symptoms: sometimes they may be stunted or appear symptomless [10]. The incidence and severity of the bakanae or stunting disease symptoms varies with regions and isolate. The pathogen is dispersed predominantly with infected seeds, infected crop residues from the previous season in the soil, or by conidia on diseased stems which can be transmitted by rain and wind [6]. Disease control has become increasingly difficult due to rapidly developing fungicide resistance in the fungal population [4].
Although all newly sequenced strains except for V64-1 clearly belong to the species F. fujikuroi, visualization of chromosome content of the ten strains by pulse field gel electrophoresis (PFGE) combined with clamped homogeneous electric fields (CHEF) indicated that all strains contained multiple chromosomes of varying sizes ( Fig 1C). However, the precise number of chromosomes in each strain could not be determined because several chromosomes had a similar size and could not be distinguished. Table 2 summarizes physical genome features of the newly and previously (IMI 58289) sequenced strains. The genome size for the F. fujikuroi strains is in the range of 43.9 Mb (IMI 58289) to 46.1 Mb (E282 and FSU48). The outgroup genome F. oxysporum V64-1 is smaller (49.1 Mb) than the reference F. oxysporum strain 4287 (61.4 Mb). Most of these differences are likely due to a different read coverage and different completeness of assemblies for repetitive regions. Despite the varying number of protein-encoding genes (14,817 for IMI 58289 to 16,088 for E282) which partially may be due to a manual gene structure validation, most key genome features of the newly sequenced strains are similar to those of F. fujikuroi IMI 58289 [15].
The completeness of the new draft genomes was explored by comparing each predicted proteome to two different, highly conserved eukaryote protein sets by BLAST [20,21]. Orthologs for all conserved proteins were found in all proteomes with the exception of a missing ortholog to 'T-complex protein 1 subunit theta' and an ortholog to 'translation factor eIF6' in F. fujikuroi B20. In F. oxysporum V64-1, an ortholog to 'MTO1 -mitochondrial translation optimization' is missing suggesting that also these genomes are more than 99% complete. In addition, the protein sets were subjected to BUSCO analysis and subsequently compared to published Fusarium proteome sets [22]. Of the 3,725 single-copy orthologs searched (library Sordariomyceta_odb9), 97.9%-99.2% were detected as complete and single-copy in the F. fujikuroi strains and F. oxysporum V64-1 which even exceeds the completeness of F. graminearum PH1 (94.7%) and F. oxysporum (94.2%) (S1 Table).
Species of the FFC are heterothallic [23]. To determine the mating type of the nine F. fujikuroi isolates, we searched the genomes for the presence of genes either of the MAT1-1 or MAT1-2 mating types. Five of the strains contain the genes MAT1-1-1, MAT1-1-2 and MAT1-1-3, all belonging to the MAT1-1 mating type, while the other four contain the MAT1-2 mating type locus with the HMG-box type transcription factor (TF)-encoding gene MAT1-2-1 and an additional gene, MAT1-2-3, that has been recently identified as specific to the Hypocreales [23] (S2 Table). Differences between the isolates were found regarding the size of gene families prevalent in the genome of B14, which is the only one of the ten examined strains causing stunting of rice seedlings (see below). The number of transporter-encoding genes (954), Zn(II) 2 Cys 6 fungal type TFs (567), cytochrome P450-(173) and dehydrogenase-encoding genes (404) were higher when compared to the other F. fujikuroi strains. Also the number of polyketide synthase (PKS) genes (18 plus two truncated pseudogenes) was larger than in the genomes of the other strains ( Table 3).
The strain NCIM 1100 encodes less PKS (eight plus one truncated pseudogene) and nine terpene cyclase (TC) genes which results in the inability to produce gibepyrones, fujikurin and certain terpenes besides other unknown products (S3 Table and below).
To gain a deeper insight into the variation at the chromosome ends between the ten analyzed strains, we used a PCR approach with primer pairs from the two ends of each chromosome based on the genome sequence of strain IMI 58289 [15]. This analysis revealed differences between the isolates either on one side or on both sides of the chromosomes (S1 Fig) variations were found for subtelomeric regions of chromosomes 3, 5, 7, 10 and 11. Two of the key enzyme-encoding genes of SM gene clusters, NRPS31 (apicidin F) and PKS16 (unknown product), are located at the far end of chromosome 1 and 11, respectively, and the PCR analysis clearly showed differences in the presence of these SM genes (S1 Fig).
Variation in growth, colony morphology and asexual spore formation All isolates showed differences in colony morphology and pigmentation on solidified and in liquid media, respectively, indicating metabolic variation between them (Fig 2A-2E). There were differences in the ability to produce the red pigments, bikaverin and fusarubins, despite the presence of the respective gene clusters in all F. fujikuroi isolates (Fig 2C-2E). Previously, we have shown that both PKS-derived pigments are only produced under low nitrogen conditions. However, while bikaverin biosynthesis is induced at acidic pH, the perithecial pigments fusarubins are produced only under alkaline pH conditions [24,25]. Strain C1995 showed no coloration under bikaverin-production conditions (Fig 2C), and C1995, B14, FSU48 and NCIM 1100 showed no pigmentation under fusarubin production conditions ( Fig 2E). To verify that the biosynthetic genes are not expressed in these isolates, we performed Northern blot analyses using the bikaverin and fusarubin biosynthetic genes as probes (S2 Fig). FSR2 (encoding an O-methyltransferase in the fusarubin gene cluster) was expressed in IMI 58289, but not in the remaining strains suggesting that the fusarubin genes are only slightly expressed in m567, MRC2276, B20, E282 and F. oxysporum V64-1, or alternatively that the red pigmentation in these strains under fusarubin production conditions ( Fig 2F) might be due to the expression of the bikaverin biosynthetic genes. The species F. fujikuroi is described to produce slender macroconidia with three to five septa and oval or club-shaped microconidia, mostly without or with one septum [1]. However, the strains used in this work differ in their ability to produce asexual spores, i.e. micro-and macroconidia ( Fig 2F). While most of the strains produce both types of conidia, others produce predominantly (B14 and NCIM 1100) or exclusively (IMI 58289) microconidia or predominantly macroconidia (E282), respectively, on V8 agar under light conditions. Strain m567 hardly forms any spores. In A. nidulans, asexual reproduction has been extensively studied for several decades [26][27][28]. Sequential activation of three major regulators, BrlA, AbaA and WetA, is necessary for the fungus to undergo asexual development. In addition, there are a number of fluffy regulatory genes (FLBA-FLBE) which regulate BRLA expression. Recently, it has been shown that an AbaA-WetA pathway is conserved in the distantly related species Fusarium graminearum [29][30][31]. There are highly conserved homologs for FlbB, FlbC, FlbD, FlbE and the central regulators WetA and AbaA also in the genomes of the F. fujikuroi isolates. However, no close relative to BrlA, which is responsible for the vesicle formation during conidiogenesis in A. nidulans, has been found in any Fusarium genome.
Bioinformatic analysis of the nine F. fujikuroi strains revealed several differences in the presence of PKS, NRPS, dimethylallyltryptophan synthase (DMATS) and TC gene clusters between them and also compared to closely related FFC members, such as F. mangiferae and F. proliferatum [43] (S3 Table). Altogether, there are 17 NRPS-, 23 PKS-, three DMATS-, and twelve TC-encoding genes present in the nine F. fujikuroi strains, i.e. 55 unique core-enzymeencoding genes that could give rise to 54 distinct SMs (the fusaric acid cluster encodes two core enzymes, PKS6 and NRPS34). Only twelve of the 23 PKS-, 13 of the 18 NRPS-, eight of the twelve TC-and two of the three DMATS-encoding genes are present in the genomes of all F. fujikuroi isolates. The Korean strain B14 has three additional PKS genes (PKS40, 43 and 51) not present in any other analyzed F. fujikuroi strain (S3 Table, Fig 3A and 3B). PKS51 is also not present in any other member of the FFC, while PKS40 is present in F. proliferatum ET1 and F. verticillioides, and PKS43 in F. mangiferae, respectively (S3 Table). In addition, B14 is the only F. fujikuroi isolate with a complete PKS5 cluster. This yet uncharacterized cluster is either absent or contains pseudogenes in the other isolates of this species, but seems to be functional in F. proliferatum and F. mangiferae ( Fig 4A). The fujikurin gene cluster (PKS19) has been previously described as F. fujikuroi-specific [15], but was recently discovered in two newly sequenced strains of F. proliferatum [43]. However, the fujikurin cluster is absent in three of the nine analyzed F. fujikuroi isolates (S3 Table, Fig 4B). The apicidin F gene cluster (NRPS31) which is located at the far end of chromosome 1 in strain IMI 58289 [15] was shown to be unique for F. fujikuroi [32]. This cluster is present in all isolates but one, B14, most likely due to chromosome rearrangements in the subtelomere regions (S1 and S3A Figs). While the entire fumonisin gene cluster (PKS11) is present in most members of the FFC except for F. mangiferae (S3 Table), nine genes of the cluster, the homologs of FFUJ_09248 -FFUJ_09256, are missing in F. fujikuroi B20 (S3B Fig) and C1995.
The Indian isolate NCIM 1100 has only 14 PKS-encoding genes (S3 Table). The absence of five PKS genes/clusters is due to large parts of chromosome 11 missing at both ends in this strain compared to IMI 58289 (Figs 3C and 4C, S1  Table). Also the discontinuous distribution of some other SM gene clusters among the F. fujikuroi isolates is likely the result of gene loss and gene gain events, respectively, at the chromosome ends. However, gene clusters located in central parts of chromosomes, e.g. PKS19 (fujikurin) and PKS11 (fumonisin) clusters, are also only present in some and absent in other strains, while the genes at the right and left borders of the clusters show colinearity between all isolates. It is not clear whether these gene clusters are the result of horizontal gene transfer [44] or cluster duplication and loss (birth and death) as shown for the fumonisin gene cluster [45].
Thanks to the genome sequencing and subsequent identification of putative gene clusters [15], SM products have now been assigned to 22 of the 54 predicted PKS-, NRPS-, DMATS-, and TC-derived SM gene clusters in isolates of F. fujikuroi, and an additional two (ferrirhodin and depudecin) in members of the FFC. One reason for the inability to identify more products is probably due to the fact that many of the SM genes are not or only minimally expressed under laboratory conditions [46]. Only 15 of these predicted clusters with known products are present in the genomes of all nine F. fujikuroi isolates indicating the genetic diversity among isolates of one species. The genes which are commonly present in all F. fujikuroi genomes are four NRPS genes required for synthesis of the siderophores ferricrocin and fusarinine and the mycotoxins beauvericin and fusaric acid (together with PKS6); six PKS genes required for Two pathotypes in Fusarium fujikuroi synthesis of fusarubins, bikaverin, fusaric acid (together with NRPS8), fusarins, trichosetin, and fumonisins; six TCs required for synthesis of GAs, phytoene, eremophilene, (-)-α-acorenol, (-)-guaia-6,10(14)-diene, and (+)-koraiol, and one DMATS responsible for biosynthesis of r-N-DMAT (S3 Table).
To study the variability in the production of the most prominent SMs, we cultivated the strains under three standardized culture conditions (6 mM and 60 mM glutamine, 6 mM NaNO 3 ). The optimal conditions for the production of the different metabolites were previously shown to vary considerably regarding nitrogen availability and pH in strain IMI 58289 [15]. Therefore, we analysed these SMs in the three media by high-performance liquid chromatography coupled to a Fourier transform mass spectrometer (HPLC-FTMS) or by HPLC with a diode array detector (HPLC-DAD) in the case of GAs. We compared the production levels ( Table 4) with transcriptome profiles of the SM biosynthetic genes generated by RNAseq for two of these conditions (6 mM and 60 mM glutamine) (Table 5A and 5B). For most of the SMs, we observed similar regulatory mechanisms regarding nitrogen availability as previously described for strain IMI 58289. However, there were strain-specific differences.
The most prominent and species-specific SMs are the GAs causing the F. fujikuroi-specific bakanae disease of rice. Although some recently sequenced species of the FFC such as F. mangiferae and F. proliferatum contain one or even two GA gene clusters, these species produce either no or only very small amounts of GAs [43]. In this work, we examined whether all analyzed F. fujikuroi strains produce GAs, and whether the GA levels can be correlated to the virulence on rice.
Previously, the regulation of GA biosynthesis has been extensively studied for strain IMI 58289. It has been shown that GA gene expression is strictly regulated by nitrogen availability in an AreA-and AreB-dependent manner [47][48][49]. To examine whether low nitrogen conditions are optimal for GA production also in the other isolates, we performed HPLC-DAD analysis of all strains and analysed the expression by Northern blot analysis in addition to RNAseq data. Accordingly, we observed the highest GA yields and GA gene expression levels under low nitrogen conditions for eight of the nine F. fujikuroi strains. The only exception among the F. fujikuroi isolates is strain B14 which showed no production of GAs and no visible expression of the GA genes despite the presence of the complete GA gene cluster in the genome (Table 5A and 5B; Fig 5A and 5B). Besides the very low expression of GA genes, strain B14 differs from the other strains by high expression of fumonisin biosynthetic genes under low nitrogen conditions (Table 5A, Fig 5A and 5C). Recently, we have shown that the fumonisin gene cluster in F. fujikuroi IMI 58289 is almost silent. Consequently, only very low amounts of fumonisins are produced in comparison to F. verticillioides. In F. fujikuroi IMI58289, the genes are only expressed and fumonisins are only produced, when the cluster-specific TF gene FUM21 is constitutively and strongly expressed [50]. Here, similar results were observed for eight of the nine F. fujikuroi isolates (except for B14) which all produce either no or only very small amounts of fumonisins and showed very low expression levels.
In addition, B14 is the only isolate producing fujikurins under alkaline and low nitrogen conditions (Table 4). Previously, the fujikurin gene cluster (PKS19) was shown to be silent under all conditions tested in strain IMI 58289, and the products have been identified only gene cluster is present in some and totally missing in other isolates. (C) PKS16 at the end of chromosome 11 in F. fujikuroi IMI 58289 is truncated in stains C1995 and E282, and missing in strain NCIM 1100. Arrows in blue represent genes belonging to a specific gene cluster. Green arrows represent genes that have closely related homologs in two or more isolates while light-gray arrows represent genes that do not have closely related homologs in other isolates. Ψ indicates a pseudogene. https://doi.org/10.1371/journal.ppat.1006670.g004 Two pathotypes in Fusarium fujikuroi after simultaneous over-expression of both the cluster-specific TF gene and the PKS19 gene itself [15,34].
In most of the strains, fusarins, fusaric acid, and apicidin F (except for B14) were only produced under high nitrogen conditions (60 mM glutamine) in accordance with the higher expression levels of the corresponding biosynthetic genes. The genes of the beauvericin cluster were recently shown to be silent under all conditions tested and activated only after deletion of the histone deacetylase gene HDA1 and knock-down of the histone methyltransferase gene KMT6 in strain IMI 58289 [33,38]. Most of the other strains analyzed here showed no or very low expression levels and no beauvericin production. The only exception was strain B20 that      Two pathotypes in Fusarium fujikuroi highly expressed the biosynthetic genes under both nitrogen conditions (low and high amounts of glutamine) and produced high levels of beauvericin at high nitrogen (Table 4;  Table 5A and 5B) suggesting that the chromatin status around the beauvericin cluster differs in B20.

Virulence assay on rice
To examine whether the isolates differ in their virulence on rice, and whether the different levels of GAs or the strain-specific production of other SMs (e.g. fumonisins in B14) correleate with the extent of symptom development, we performed assays with both germinating seeds and rice seedlings.
To determine the ability of all ten strains to impair the seed germination, rice seeds were soaked for 18 h in spore suspensions of the respective strains. The percentage of germinated seeds was counted after 14 days post inoculation (dpi). Strains B14 and B20 induced a high percentage of seedling death indicating a high potential to kill the host seedlings (S4 Table). As B14 does not produce much GAs, the aggressiveness of this strain cannot be caused by these phytohormones.
Next, we performed a rice seedling assay to compare the virulence of the different isolates to assess their ability to induce bakanae symptoms. The surface-sterilized seeds were first germinated, and the young healthy seedlings were then inoculated with spores. All strains except for B14 and F. oxysporum V64-1 (outgroup) induced the formation of slender, elongated and yellowish stems (Fig 6). Wheras V64-1-infected seedlings behaved like the water control, B14-infected seedlings were stunted and showed withering instead of the typical bakanae symptoms. The total seedling length and the length of the internodes were smaller than, or comparable to uninfected seedlings (water control). Previously, it has been reported that heavily infected seedlings can also be stunted and can show severe crown and root rot [51]. The type of symptoms and severity of disease depends on the fungal isolate and is thought to be affected by the proportions of GA and fusaric acid produced by the fungus, which potentially cause elongation of the plants or stunting, respectively [52,53]. However, this assumption has never been proved experimentally.

In planta expression of SM gene clusters and product analysis
To examine which of the SM gene clusters are similarly expressed in rice in all isolates, and which of them are specifically expressed or not expressed in strain B14, and therefore could be relevant for disease symptoms (especially stunting), we performed RNA-seq for all isolates also from rice seedlings at 7 dpi (Table 5C). The expression patterns were compared with those under in vitro conditions at high (60 mM glutamine) and low (6 mM glutamine) nitrogen (Table 5A and 5B).
In general, the GA biosynthetic genes (shown for the key enzyme gene CPS/KS) are the most highly expressed genes in rice roots except for B14 in which they are hardly expressed (Table 5C).
To check whether the expression profiles of SM genes correspond to the in planta production levels, we also performed SM analysis of infected rice plants by HPLS-FTMS at 7 dpi. In accordance with gene expression, the in planta GA production levels also differed between the isolates: B14-infected rice plants contained no detectable GA amounts in rice compared to the bakanae strains ( Table 6). The most obvious difference of strain B14 compared to the other isolates was the high expression of the fumonisin genes resulting in significant levels of fumonisins (FB 1 and FB 2 ) in planta, similar to those observed in vitro (Table 5C; Table 6). Furthermore, while all isolates showed similar expression patterns for fusaric acid biosynthetic genes (PKS6; NRPS34) and similar fusaric acid production levels (Table 5C, Table 6), B14 was the only strain with low expression of the fusarin C genes. As expected, no fusarins were detectable in rice roots inoculated with this strain (Table 6).

Two pathotypes in Fusarium fujikuroi
In addition, the yet uncharacterized PKS51 gene cluster, which is only present in B14, was exclusively expressed in planta, suggesting that the unknown product of this gene cluster might play a role during infection (Table 5C). Besides B20, B14 also gave high expression for beauvericin genes, but only low expression of gibepyrone (PKS13) and acorenol (STC6) biosynthetic genes.
In conclusion, the very low expression of GA genes and the lack of detectable GA levels in B14-infected rice seedlings after 7 dpi are most likely responsible for the absence of bakanae symptoms in this isolate. Instead, the high levels of fumonisins which are exclusively produced only in this isolate, may overrule the growth-stimulating effect of the GAs and cause stunting/ withering. Do the fumonisins, fusaric acid and/or the PKS51 product cause the stunting symptom development of F. fujikuroi B14?
To find out whether one of the gene clusters that are specifically expressed in B14 during infection on rice (Table 5C) might indeed cause the stunting effect of this isolate, the key enzymeencoding genes for fumonisins (FUM1 = PKS11), and the yet uncharacterized B14-specific PKS51 gene were deleted in this strain. In addition, we also deleted the key gene for fusaric acid biosynthesis (FUB1 = PKS6) due to the speculation that fusaric acid production might cause the stunting pathotype [52].
The 5-day-old healthy rice seedlings were soaked in the conidial suspension of B14 or the mutant strains. At 5 dpi, the B14-infected seedlings were already stunted compared to the water control and showed withering symptoms which were even more obvious at 7 dpi and 9 dpi (Fig 7A). In contrast, the Δfum1 and Δfub1 strains of B14 appeared healthy at 7 dpi and caused delayed withering symptoms at 9 dpi only. The roots of all infected seedlings were heavily colonized with fungal mycelia at 9 dpi (Fig 7D). The double deletion mutants of B14 lacking both FUM1 and FUB1 behaved like the mock control and did not induce stunting at 5 dpi (Fig 7A) or occasionally upto 9 dpi. However, in some cases, they caused a similar delayed disease development as the single deletion strains at 9 dpi (Fig 7B and 7C). The add-back strains, which were generated by introducing the native copy of FUM1 and FUB1 into Δfum1  Two pathotypes in Fusarium fujikuroi and Δfub1 mutants, respectively, caused typical stunting/withering symptoms similar to B14 (S4A Fig). The delay of stunting/withering symptom development by Δfum1 and Δfub1 strains indicates that the production of fumonisins and/or fusaric acid, in combination with the nondetectable levels of GAs, play an important role for the development of this specific pathotype in B14. For comparison, the genes FUM1 and FUB1 were also deleted in strain B20 (bakanae strain). Unlike the B14-derived mutants, all of the gene deletion strains of B20 we examined caused typical bakanae symptoms (elongated, slender and chlorotic shoots), as did B20. No withering of the seedlings or mycelial colonization on roots were visible suggesting that the absence of fumonisin production by B20 and B20-1 derived mutants protect the rice seedlings from development of B14-like symptoms (S5 Fig).
Surprisingly, rice seedlings inoculated with the ΔPKS51 mutant strains showed an earlier and more severe symptom development compared to those with B14. The stunting/withering symptoms were clearly shown already at 5 dpi (Fig 7A). The more severe symptoms caused by the ΔPKS51 strains may suggest a possible role of the PKS51 product as an avirulence determinant. A similar role was described for the product of the Magnaporthe grisea PKS-NRPS gene ACE1 (Avirulence Conferring Enzyme1) which is also specifically expressed only on rice. Its yet unknown product is probably recognized by rice cultivars carrying a specific resistance gene [54,55]. However, further investigations will be needed to show whether the product of PKS51 acts in a similar way in F. fujikuroi strain B14.
To further study the impact of these SMs on disease symptom development, fumonisin FB 1 , fusaric acid, or GA 3 were exogenously supplied to rice seedlings infected with wild-type or mutant strains (Fig 7A). Addition of FB 1 and fusaric acid to rice seedlings infected with Δfum1 and Δfub1, respectively, restored the WT phenotype and resulted in stronger withering symptoms compared to the deletion strains without the toxins. However, the even more reduced virulence of the double deletion strain (Δfum1/Δfub1) was not clearly complemented by exogenous supplies of both fumonisin and fusaric acid (Fig 7A). Interestingly, exogenous supply of GA 3 did not cause bakanae symptoms on rice seedlings inoculated with the wild-type B14 strain, while addition of GA 3 to the Δfum1/Δfub1 double deletion strain caused stem elongation of rice seedlings at 5 dpi, similar to the bakanae symptom caused by B20 (Fig 7A). Therefore, conversion of B14 into a bakanae pathotype by addition of GAs was only possible after deleting the key genes for the production of fumonisins and fusaric acid.
Furthermore, exogenous supply of culture filtrate from strain B14 to seedlings infected with B14, B20 or B20 Δcps/ks caused stunting and withering symptoms while the culture filtrate of the B14 Δfum1/Δfub1 mutant caused milder symptoms (S4B and S4C Fig). Addition of culture fluid to B14-infected seedlings resulted in even more severe stunting than B14 alone (S4B Fig). Exogenous supply of B14 culture filtrate to seedlings inoculated with B20 or the GA-deficient B20 Δcps/ks mutant revealed severe stunting symptoms in both cases irrespective of the ability to produce GAs (S4C and S4D Fig). These data support our suggestion that fumonisins and fusaric acid play an important role for symptom development.
To further investigate the role of fumonisins for causing stunting and withering, we performed the rice seedling pathotest with two F. verticillioides strains from corn which were shown to produce high amounts (more than 3,000 μg/g) of fumonisins and no GAs in rice seedlings. Both F. verticilloides isolates (Os35 and Os40) [56] caused withering at 7 dpi although the plants were not stunted compared to those of the mock control (S6A Fig). In addition, mycelia of both F. verticillioides strains colonized the roots of infected rice seedling as much as those of B14 (S6B Fig). This strong root colonization was not observed for roots infected with B20. Previously, it has been already reported that B14 triggered severe inhibition of root growth, and that its own growth rate in rice roots was more than 4 times higher compared to that of B20.
Taking together the results of FUM1 and FUB1 deletion, exogenous addition of pure toxins or culture filtrate of B14 to rice seedlings inoculated with wild-type or mutant strains, and the pathotests with GA-deficient, highly fumonisin-producing F. verticillioides isolates provide strong indications that the biosynthesis of both fumonisins and fusaric acid and the lack of GA biosynthesis in B14 play cruicial roles for causing the stunting/withering phenotype on rice seedlings. However, additional factors are probably involved in inducing the stunting symptoms.

Comparisons of TF-encoding genes between B14 (stunting) and B20 (bakanae) strains of F. fujikuroi
Besides SMs, stunting/withering might be caused also by the different sets of TFs present in the genome of B14 and the bakanae strains, or by different expression levels of TF-encoding genes in rice. Therefore, we compared the expression levels of TF-encoding genes between B14 and the other eight F. fujikuroi isolates (S5A Table). There are 37 genes which are present in most of the nine genomes and which were specifically up-regulated during infection of rice (S5B Table). B14 had slightly higher expression levels only for three of them: FFB14_03090, FFB14_05980 and FFB14_01631. In addition, B14 has 28 strain-specific TFs which are not present in the genomes of the other strains (S5C Table). The most highly expressed gene in planta was FFB14_06367 encoding the pathway-specific TF of the putative PKS51 gene cluster. The high expression of PKS51 and the adjacent genes, including the TF-encoding gene, supports our assumption that this unique gene cluster is involved in determining the severity of disease symptom development.

Population analyses revealed the segregation of field isolates into the two different pathotypes
Because B14 was the only isolate causing the stunting pathotype among the analyzed ten strains, we attempted to determine how often this pathotype can be found in rice fields. Therefore, we inoculated rice seedlings with 15 field isolates, which were collected from rice grains and air above rice paddy fields in Korea between 2014 and 2016. Among the 15 isolates, we identified additional nine field isolates causing stunting and early withering symptoms similar to B14 while six isolates caused typical bakanae symptoms (Fig 8A).
To determine whether these isolates can be phylogenetically distinguished from each other, we generated a phylogenetic tree using the nucleotide sequences of the RPB2 and EF1A genes [57] from the new field isolates, the ten isolates used in this study, and some closely related Fusarium species of the FFC. The F. fujikuroi clade, which was clearly separated from those of other closely related species such as F. proliferatum and F. verticillioides, contained two strongly supported subclades (with 76% bootstrap support: BS). Interestingly, the subclade with B14 consists of all of the field isolates causing stunting symptoms, while the other subclade contained B20 and the other F. fujikuroi strains used in this study and all bakanae-type field Two pathotypes in Fusarium fujikuroi isolates (Fig 8B). This result indicates that the two pathotypes of F. fujikuroi exist as phylogenetically distinct groups within the population.
To examine if there is a correlation between the pathotype and the presence of one of the SM clusters specific to B14 or B20 (PKS51 for B14, unknown product, and NRPS31 for B20, apicidin F) among the field isolates of each pathotype, we performed a PCR amplification using primer sets derived from PKS51 (unknown product) and NRPS31 (apicidin F), respectively. The PKS51-specific primer set amplified a fragment from all of the stunting-type isolates examined including B14, but not from the bakanae-type isolates. Similarly, the NRPS31-specific primer set amplified a fragment only from bakanae-type isolates examined (S7A Fig).
Based on the correlation between the phylogenetically distinct pathotypes on the one hand, and the presence of either the PKS51 or the apicidin F (NRPS31) gene clusters, we examined additional field isolates from Korea using these primer sets. Among a total of 151 isolates examined, the B14-specific PCR amplification pattern was found in 69% of rice isolates, 100% of corn isolates, and 85% of airspora isolates. A phylogenetic tree generated from all of these isolates revealed a clustering of the stunting-type and bakanae-type isolates similar to that for the field isolates examined by pathogenicity tests (S7B Fig). However, it is currently unclear why the stunting-type population is predominant among the field isolates in Korea.
To gain more information on the SM profile of these new stunting-type isolates, we examined the expression levels of the key genes of GA (CPS/KS) and fumonisin (FUM1) biosynthesis in the field isolates grown in liquid culture with 6 mM glutamine by use of quantitative real-time PCR (qPCR). The CPS/KS transcript levels from all the stunting-type isolates examined were similar to or even lower than that in B14, while those from all the bakanae-type isolates, including B20, were 8-to 24-fold higher than those in the stunting-type strains (S8A Fig). Furthermore, all stunting-type isolates showed clear expression of FUM1 in contrast to the bakanae isolates (shown for B20) (S8B Fig). These data indicate that both subgroups of isolates differ in a whole set of characteristic features (presence of either the PKS51 or NRPS31 genes, expression of either GA or fumonisin genes) which correlate with either the stunting or the bakanae pathotype. We also examined the production levels for both SMs in some of the new B14-like strains in vitro and in planta. As expected, the stunting-type isolates produced no or 10 to 15 times less GAs than the bakanae strain IMI 58289, while only the stunting-type strains produced fumonisins in submers cultures (S9A and S9B Fig). The analysis of GAs and fumonisins in rice roots and shoots revealed no GAs at all in seedlings infected with stuntingtype isolates, while only the latter produced significant amounts of fumonisins in planta (S9C and S9D Fig). A similar correlation between low GA and high fumonisin levels on the one hand, and a pathotype called "dwarfism" on the other hand, was also described for the Italian F. fujikuroi isolate CSV1 [58].
In conclusion, we provide genome sequences of eight new F. fujikuroi isolates and one F. oxysporum isolate, all collected in different rice growing regions worldwide. We show that all these strains differ in genome and chromosome size, number of genes in major gene families such as TFs, transporters, SM biosynthetic genes and others. In addition, the isolates differ in colony morphology, pigmentation and the number and type of asexual spores (micro-and/or phylogenetic tree constructed by the NJ method using the nucleotide sequences of combined TEF1 and RPB2 from the two pathotypes of field isolates, determined by pathogenicity test (A), and from the worldwide isolates used in this study. The phylogenetic tree supports the separation of the two pathotypes into two different phylogenetic subclades. https://doi.org/10.1371/journal.ppat.1006670.g008 Two pathotypes in Fusarium fujikuroi macroconidia). Major differences were identified in subtelomeric regions of the chromosomes where several SM gene clusters are located. Besides the differences in the presence of gene clusters, we observed variations in the ability to express SM genes and to produce the respective metabolites under in vitro and in planta (rice) conditions. Among the nine F. fujikuroi isolates analyzed, eight cause typical bakanae symptoms on rice seedlings due to their ability to produce GAs. Only one isolate, B14, does not cause elongation of infected plants, but instead causes stunting combined with early withering of rice seedlings. This isolate is the only one which does not produce GAs under in vitro and in planta conditions. Instead, B14 produces high amounts of fumonisins both under in vitro and in planta conditions. Furthermore, it is the only strain containing a putative gene cluster with PKS51 as key gene which is highly expressed in rice. To examine the determinant(s) of the stunting/withering pathotype, several key enzyme-encoding candidate genes were deleted in B14. The data demonstrate that the formation of fumonisins, and probably also fusaric acid, on the one hand, and the lack of GA production on the other hand, contribute to the stunting/withering pathotype, because the deletion of the respective key genes resulted in reduced virulence. Furthermore, the unknown PKS51 gene cluster seems to produce a SM which acts as an attenuator of disease, because the deletion of the PKS51 gene caused early withering of infected seedlings. Examination of more field isolates from Korean rice fields revealed a correlation between the pathotype and the ability to produce either fumonisins or GAs which is supported by the clear separation into two distinct phylogenetic clades.

Fungal and bacterial strains
The fungal strains used in this work and their origin are shown in Table 1. Strain IMI 58289 was derived from Commonwealth Mycological Institute, Kew, United Kingdom. Strains m567 and FSU48 were provided by the Fungal Stock Center at the University Jena, Germany, C1995 by J.F. Leslie, Kansas State University, E282 by S. Tonti, University Bologna, Italy, MRC2276 by W.C.A. Gelderblom, South Africa, and NCIM 1100 was provided by the National Collection of Industrial Microorganisms, India. Strains B14 and B20 were provided by S.-H. Yun, Korea. Strain V64-1 was kindly provided by T. Kyndt from the University Ghent, Belgium.
Escherichia coli strain Top10 F' (Invitrogen, Groningen, The Netherlands) was used for plasmid propagation. The uracil-auxotrophic Saccharomyces cerevisiae FGSC 9721 (FY 834) was provided by the Fungal Genetics Stock Center (Kansas State University) and used for yeast recombination cloning.
Genome sequencing and assembly F. fujikuroi B20: Illumina TrueSeq genome sequencing by TheragenEtex, Suwon, Korea. The assembly was performed using Celera Assembler version 7.0 [59], 'overlap minimum length' set to 150 bases. The assembly resulted in 318 scaffolds with a 21-fold coverage of the TrueSeq large reads.
For all other strains sequencing was carried out by shot gun sequencing of an 8 kb library with paired end 100 bp read length using Illumina HiSeq 2000 by Eurofins MWG Operon, Germany. The assemblies were performed by ALLPATHS-LG [60] and the scaffolds were error corrected by mapping all Illumina shotgun paired-end data and further scaffolded using SSPACE [61] ( Table 2). The data on the new genomes, including annotation, was submitted to the European Nucleotide Archive, study accession PRJEB14872 available at: http://www.ebi.ac. uk/ena/data/view/PRJEB14872. Sample accession numbers are listed in Table 1. RNA-seq data are available at: https://www.ncbi.nlm.nih.gov/gds/?term=GSE89480.

Structural annotation
Draft gene models for all genomes were generated by three de novo prediction programs: 1) Fgenesh [62] with different matrices (trained on Aspergillus nidulans, Neurospora crassa and a mixed matrix based on different species); 2) GeneMark-ES [63] and 3) Augustus [64] with Fusarium ESTs and RNA-seq based transcripts as training sets. Annotation was aided by exonerate [65] hits of protein sequences from F. fujikuroi IMI 58289 and F. oxysporum 4287 to uncover gene annotation gaps and to validate de novo predictions. Transcripts were assembled on the RNA-seq data sets using Trinity [66]. The different gene structures and evidences (exonerate mapping, RNA-seq reads and transcripts) were visualized in GBrowse [67] allowing manual validation of coding sequences with a focus on SM cluster genes and other genes of interest. The best fitting model per locus was selected manually and gene structures were adjusted by splitting or fusion of gene models and redefining exon-intron boundaries if necessary. tRNAs were predicted using tRNAscan-SE [68]. The predicted protein sets were searched for highly conserved single (low) copy genes to assess the completeness of the genomic sequences and gene predictions. Orthologous genes to all 246 single copy genes were searched for all proteomes by BLASTp comparisons (eVal: 10 −3 ) against the single-copy families from all 21 species available from the FunyBASE [21]. Additionally, the proteomes were searched for the 248 core-genes commonly present in higher eukaryotes (CEGs) by BLASTp comparisons (eVal: 10 −3 ) [20]. We also used BUSCO Version 3.0.1 in the Ubuntu virtual machine with the lineage specific profile library Sordariomyceta_odb9 (3.725 BUSCO groups), downloaded from http://busco.ezlab.org. The analysis was performed in gene set (protein) assessment mode running the python script run_ BUSCO.py [22]. All genomes were analyzed using the PEDANT system [69]. To avoid misleading ortholog information based on similarity and bi-directional best hits, QuartetS [70] was applied to retrieve a reliable ortholog matrix which was used for all comparative representations.

Phylogenetic tree analyses
The phylogenetic tree of Fusarium species was calculated based on the protein sequences of 5,181 single copy genes that are shared among all analyzed species. Orthologs of the sequences were aligned separately using MAFFT [71]. After that, we concatenated the alignments and removed columns with gaps using Gblocks [72]. The evolutionary history was inferred using the Maximum Likelihood method PhyML [73] with default parameters and the amino acid substitution model LG. Branch support was tested using approximate likelihood ratio test (aLRT) based on the Shimodaira-Hasegawa-like (SH-like) procedure [74]. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site.
We calculated single copy genes in clustering proteins of all genomes and selecting clusters with exactly one representative from each genome. Protein clusters were calculated using usearch [75] (e-value cutoff: 0.01) and mcl [76] (inflation value: 2).

Contour-clamped homogeneous electric field (CHEF)-gel analysis
Protoplasts from Fusarium strains were prepared as described previously. The protoplast suspension was mixed with 1.2% InCert agarose (Lonza Group AG, Basel, Switzerland) and then loaded on a CHEF gel as described in [80]. Chromosomes of Schizosaccharomyces pombe and S. cerevisiae were used as a molecular size marker (Bio-Rad, Munich, Germany).

In silico identification of secondary metabolite clusters
To identify SM clusters in each genome, the InterPro scan results of the PEDANT analysis were used as described [81]. Essentially, predicted proteins with homology to a domain of a signature SM core enzyme (e.g. PKS, NRPS, TC or DMATS) were considered a marker for a gene cluster. A cluster was verified if any neighboring genes with homology to typical SM cluster enzymes, like P450 monooxygenases, oxidases, methyltransferases, MFS or ABC transporters or TFs, were identified. The extent of each putative gene cluster was then adjusted by comparison to previously published data and to homologous clusters in other Fusarium species.

Cultivation methods
For SM production experiments and RNA preparation under in vitro conditions, strains were first cultivated for 3 days in 300 mL Erlenmeyer flasks with 100 mL Darken medium [82] on a rotary shaker at 180 rpm at 28˚C. 500 μL of this culture were then used to inoculate 100 mL of ICI (Imperial Chemical Industries, UK) media [83] containing either 6 mM glutamine, 60 mM glutamine, or 6 mM NaNO 3 for 3 days (RNA extraction) or 7 days (SM analysis), respectively. For RNA extraction, mycelia were flash-frozen with liquid nitrogen prior to lyophilization. For in planta expression studies by RNA-seq, lyophilized roots of infected rice plants were used for RNA preparation.

Plant material for RNA sequencing
For the generation of infected rice samples for transcriptome analysis by RNA-seq, rice seeds of the cultivar GSOR 100, Nipponbare, and Dongjin were used. The seeds of the former two were provided by Genetic Stocks-Oryza (GSOR) Collection, USDA ARS Dale Bumpers National Rice Research Center, Hwy, Arkansas, USA.

RNA isolation, RNA sequencing and quantification
Total RNA was extracted from mycelia grown for 3 days in liquid ICI media (containing either 6 or 60 mM glutamine) and from infected rice seedlings after 7 dpi using TRIzol Reagent (Life Technology, Karlsruhe, Germany) and purified using an RNeasy Plant MinElute Cleanup Kit (Qiagen, Hilden, Germany). The quality of DNase-treated RNA (28S:18S > 1.0; RIN! 6.5; OD 260/280 !1.8; OD 260/230 ! 1.8) was determined using an Agilent Bioanalyzer. The high quality RNA was sent to BGI Tech Solutions Co., Limited (Hong Kong) for library construction and sequencing by Illumina HiSeq2000 technology.
RNA-seq reads were mapped on the reference genome using tophat2 (v2.0.8). The interval for allowed intron lengths was set to min 20 nt and max 1 kb [84][85][86]. We used cufflinks to determine the abundance of transcripts in FPKM (fragments per kilobase of exon per million fragments mapped) and calculated differentially expressed genes using cuffdiff [85,86]. The gene models were included as raw junctions. Genes with a minimum of fourfold increase or decrease in expression (|log2 of the FPKM values +1| ! 2) between the two experimental conditions were considered as regulated. The RNA-seq data has been deposited in NCBI's Gene Expression Omnibus [87] and are accessible through GEO Series accession number GSE89480.

Generation of vectors for targeted gene deletions and complementations
The DNA constructs for deletion of FUM1 (FFB14_08440, FFB20_01984), FUB1 (FFB14_016 51FFB20_13404), and PKS51 (FFB14_06372) from the genomes of F. fujikuroi strains B14 or B20 were created using a split-marker recombination procedure as previously described [88]. To delete FUM1, the 5' and 3' flanking regions of the FUM1 ORF were amplified with the primer pairs JFUM1f5/JFUM1rt5 and JFUM1ft3/JFUM1r3, respectively (in the first round of PCR), fused to the hygromycin B resistance gene (hygB) cassette, which was amplified from pBCATPH [88] with the primers HygB-for and HygB-rev (in the second round of PCR), and used as a template to generate split markers with the new nested primer sets, JFUM1fn/ pUH-BC/H3 and JFUM1rn/pUH-BC/H2, respectively (in the third round of PCR) (S6 Table). Similarly, DNA constructs for targeted deletions of the other genes were created using the strategy described above. For the deletion of FUB1, the primer pairs JFUMB1f5/JFUM1Brt5 and JFUB1ft3/JFUB1r3 were used for amplification of 5' and 3' flanking regions of the FUB1 ORF, respectively, and JFUB1fn/pUH-BC/H3 and JFUB1rn/pUH-BC/H2 were used as the nested primer sets, respectively. For the deletion of PKS51, the primer pairs B14_6372For5/ B14_6372rev5t and B14_6372for3t/B14_6372rev3 were used in the first round of PCR, and B14_6372forN/ pUH-BC/H3 and B14_6372revN/ pUH-BC/H2 for the third round of PCR. Additionally, for double deletion of FUM1 and FUB1, we generated a knock-out construct through which the 5 0 -and 3 0 -flanking regions of the FUM1 ORF were fused to a geneticin resistance gene cassette (gen) amplified from pII99 using the primer pair Gen-for and Genrev, as described above. The resulting construct was transformed into the deletion strain of FUB1.
For the complementation of each deletion mutant of FUM1 or FUB1 derived from B14, intact copies of each gene were amplified from the genome of B14 using the primers JFUMf5/ JFUMr3 and JFUB1f5/JFUBr3, respectively, and directly added into protoplasts of each deletion strain along with pSK660 including the geneticin resistance (gen) gene for co-transformation as previously described [89]. All primer sequences are listed in S6 Table. Pathotests on rice The effect of the single isolates on rice seedling germination was studied by infecting the seeds with conidia of each isolate. For spore formation, the fungal strains were cultured on PDA plates for 7 days at 28˚C under light/dark (12 h/12 h) conditions. The plates were flooded with sterile water to obtain a conidial suspension (1 x 10 −6 ). Seeds were soaked in the suspension for 18 h. Inoculated and non-inoculated (control) seeds were sown into 100 mL plastic pots. Eight days after sowing the number of germinated seeds was assessed. The number of dead, chlorotic and elongated seedlings was measured 15 days after sowing.
For pathogenicity assays, healthy rice seeds were surface-sterilized by submersion in 70% of ethanol followed by 1% of sodium hypochlorite. Sterilized seeds were germinated in Murashige and Skoog (MS) agar [7] at 26˚C for 5 days. The Fusarium isolates were first grown on oatmeal agar for one week. The pathogenicity assays were performed as previously described [7]. Agar plugs from the oatmeal plates were placed on top of 3 cm of sterilized vermiculite in a glass tube (18 cm x 1.6 cm). The agar plugs were then covered with 3 cm of vermiculite. Fiveday-old seedlings were transferred to the surface of the vermiculite layer to avoid the direct contact between seedlings and fungal inocula. Before covering the tubes with a cap, 4 mL of Yoshida solution was gently added to each test tube to help retain high humidity [90] and placed at 26˚C for 3, 5, 7 or 9 days. Their heights and internode lengths were measured and photographs of the seedlings and infected roots were taken. The symptom development caused by wild-type and mutant strains was examined in five independent pathogenicity tests. For pathotests with exogenous supply of culture filtrates, the wild-type B14 strain or its mutant lacking FUM1 and FUB1 (Δfum1/Δfub1) was inoculated into 50 mL of PDB (potato dextrose broth) and incubated for 5 days. The fungal liquid cultures were filtered through 2 × cheesecloth followed by filtration through 0.25 μm membranes. The culture fluid was dried to 5 ml by lyophilization (10-fold concentration). For inoculation assay, 500 μl of the concentrated culture filtrate was exogenously supplied to a single rice seedling.

Transformations
F. fujikuroi B14 and B20 strains were transformed as previously described for F. graminearum [91]. Vector integration events were confirmed by diagnostic PCR (S10 Fig) using specific primers as indicated (S6 Table).

PCR
PCR mixtures contained 25 ng of template DNA, 50 ng of each primer (S6 Table), 0.2 mM deoxynucleoside triphosphates, and 1 U of Biotherm Taq polymerase (Genecraft, Lüdinghausen, Germany). The cDNA synthesis was performed using Superscript II (Invitrogen, Groningen, The Netherlands) and 1.5 μg of total RNA as the template, according to the manufacturer's instructions. The qPCR was performed using iTaq Universal SYBR Green Supermix (BioRad) and Superscript II cDNA as template, in a Biorad thermocycler iTaq. In all cases, the qPCR efficiency was between 90-110% and the annealing temperature was 60˚C. Gene expression was measured in three biological replicates from each time point, and the relative expression levels were calculated using the ΔΔCt method [92]. The expression of a translation elongation factor α gene (EF1A), amplified by a primer pair (EF1-PS1 and EF1-PS2) (S6 Table), was used as an endogenous reference for data normalization.

Chemical analysis of secondary metabolites in vitro and in planta
For analyses of the SMs, the strains were grown in submerged cultures as described above. After 7 days, mycelia were removed from the culture by filtration through Mirachloth (Calbiochem, Merck KGaA, Darmstadt, Germany). The culture filtrates were filter-sterilized using 0.45 μm syringe filters (BGB, Schloßböckelheim, Germany).
Liquid culture samples were prepared as following: 10 μL of the culture filtrate and 10 μL of MePa (100 μg/mL) as internal standard were added to 80 μL of water. For in planta analysis, ten rice samples were combined and freeze-dried. The dried samples were treated with liquid nitrogen and pestled simultaneously, larger pieces were cut first with a scalpel. The samples were weighed and extracted with 1.5 mL of the following mixture: ethyl acetate:methanol (MeOH):dichloromethane, 3:2:1. Precellys ceramic beads 1.4/2.8 (Peqlab, Erlangen, Germany) were added to the samples, and the mixture was vortexed for 1 min. Afterwards, the samples were shaken for 1 h on a rotary shaker with 150 rpm. After a short centrifugation step (3 min, 2900 g), 500 μL of the supernatant were transferred to a new vial and the solvent was evaporated to dryness under constant nitrogen flow. The residue was dissolved in 100 μL MeOH and put to an ultrasonic bath for 10 min. After vortexing the samples, they were centrifuged again with 5000 g, and 50 μL of the supernatant were collected. Afterwards, the samples were dried again under nitrogen flow and 1.5 mL of MeOH/water, 3/1 (v/v), + 0.1% formic acid (FA) were added. The extraction procedure described above was repeated. 10 μL of each extract were mixed with 10 μL MePa and 80 μL water for analysis. Some extracts were diluted again 1:10; the corresponding values are labelled.
The liquid culture samples were analyzed as following: A Shimadzu LC-20AD HPLC (Shimadzu, Kyoto, Japan) with a SIL-20ACXR autosampler coupled to a Sciex QTRAP 5500 (Sciex, Darmstadt, Germany) mass spectrometer was used. Separation was carried out on a Synergi Hydro-RP column from Phenomenex (Aschaffenburg, Germany) with 50 × 2.0 mm and 2.5 μm particle size, protected by a KrudKatcher classic filter (Phenomenex), and on a Nucleodur C18 Gravity-SB column from Macherey-Nagel (Düren, Germany) with 50 × 2.0 mm and 3 μm particle size, protected by a KrudKatcher classic filter (Phenomenex), at 50˚C. MeOH + 1% FA + 5 mM NH 4 Ac was used as eluent A, water + 1% FA + 5 mM NH 4 Ac was used as eluent B. A flow volume of 450 μL/min was applied, and the gradient started at 5% A. This condition was held for 1.5 min. The gradient rose up to 98% A in 10.5 min, and finally the column was rinsed for 3 min with 98%. After that, the column was equilibrated for 2.5 min with 5% A. The integrated valco valve was used, discarding the first 2 min of the run, and the samples were cooled to 7˚C. 5 μL of each sample was injected.
Advanced scheduled multiple reaction monitoring (MRM) mode was used for acquisition. Both positive and negative ionization mode were applied. The curtain gas (CUR) was set to 35 psi, the collision gas was set to medium. The temperature of the heater gas (TEM) in the ion source was set to 450˚C. Nebulizer gas (GS1) and auxiliary gas (GS2) were adjusted to 35 and 45 psi, respectively.
In the positive MRM mode, the ion spray voltage (IS) was set to 4500 V, and the target scan time of this experiment was adjusted to 0.3 s, resulting in a cycle time of 0.5 s. The positive ionization mode was applied for the relative quantification of O-methyl-fusarubin, fusaric acid, gibepyrone A, apicidin F, beauvericin, fumonisins (FB 1 +FB 2 ), fusarins, trichosetin and fujikurin A. The cell exit potential (CXP) was set to 11 V, the entrance potential (EP) was set to 10 V. For the negative MRM mode, the IS was set to -4500 V and the target scan time was adjusted to 0.2 s, resulting in a cycle time of 0.5 s. The negative ionization mode was applied for the relative quantification of MePa, and the gibberellins GA 3 , GA 4 and GA 7 . The CXP was set to 11 V, the EP was set to 10 V.
MRM transition for the quantification were as follows: For bikaverin determination in liquid culture samples, HPLC-UV measurements were carried out on a Shimadzu LC-20AT pump system with a Shimadzu SIL autosampler and a photodiode array (PDA). A Gemini 5 u C6-Phenyl 110A, 250 × 4.60 mm, 5 μm column (Phenomenex) was used, with water + 1% FA as eluent A and acetonitril + 1% FA as eluent B. The column oven was set to 40˚C. The gradient started with 10% B with 1.35 μL/min. After 3 min, the gradient rose up to 100% B during 17 min. The column was washed with 100% B for 6 min, and afterwards the column was equilibrated with 10% B again for 4 min. The wavelength for PDA analysis ranged from 220-600 nm. 100 μL of the sterile culture filtrate were injected. Peak areas were determined at 508 nm, and bikaverin obtained from Sigma Aldrich was used as standard substance.
In planta sample analysis was performed with a different HPLC-system but the same mass spectrum. An Agilent 1260 HPLC system (Santa Clara, USA) was used, and the expected retention time for scheduled MRM analysis needed to be adjusted (S7-S10 Tables). Furthermore, bikaverin was analyzed in positive ionization mode.  Table) used in the PCR amplification and expected size of PCR products (designated a or b) are indicated in the deletion schemes and the gels, respectively. In most cases, two independent strains with a gene deletion or a gene add-back (designated 1 and 2) were used in PCR along with their wild-type (WT) progenitor and those carrying a single gene deletion (SD). For gene complementation, we used a co-transformation strategy as previously described [88]. (TIF) S1 Table. BUSCO single-copy analysis, performed in gene set (protein) assessment mode on the library Sordariomyceta_odb9.  Table. MRM transitions for in planta analysis in positive ionization mode, depicting also transition specific variables. (DOCX) S10 Table. MRM transitions for in planta analysis in negative ionization mode, depicting also transition specific variables. (DOCX)