Skip to main content
  • Research article
  • Open access
  • Published:

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



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.


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 (pG), however, four SNPs were associated with major quantitative trait loci (QTLs) with pG ranging from 12 to 48%. The QTL with the highest pG was involved in DON production and associated with a SNP most probably located within the Tri4 gene.


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.


Fusarium belongs to one of the most diverse and widely distributed fungal genera producing a broad variety of secondary metabolites [1]. It causes many economically important diseases in crops, the most devastating is Fusarium head blight (FHB). In Europe, FHB is caused by about 15 Fusarium species [2], where F. graminearum and F. culmorum are frequently isolated members [3,4,5]. They both cause stalk and ear rot in maize [6, 7] and seedling blight, foot rot, and head blight in small-grain cereals including wheat, barley, rye and triticale [8]. A common symptom of FHB is the pre-mature bleaching of infected spikelets. Kernels can be aborted completely when the ears are colonized during flowering. Even with low levels of infection, FHB diminishes grain quality and can reduce yield by nearly 50% [9]. In addition to yield loss, F. culmorum infection can be even more problematic since the fungus produces an array of mycotoxins, including type B trichothecenes, which are inhibitors of protein synthesis in plants and animals. The main trichothecenes reported in F. culmorum-infected cereals include deoxynivalenol (DON), nivalenol (NIV), and their acetylated derivatives [10], but in addition to these, each isolate can produce numerous other mycotoxins. For example, a single F. culmorum isolate was described that produced nine additional mycotoxins on wheat, one of which was zearalenone [11]. Crop mycotoxin contamination by Fusarium spp. is a worldwide problem and therefore strict regulatory limits are established in many countries regarding the amount of mycotoxins legally permitted in grains for food/feed usage and in processed food in the EU, USA, and Canada [10]. Consequently, Fusarium infection of cereals results in decreased market value or even the rejection of the grain.

Pathogenicity is the ability of the pathogen to infect a susceptible host [12]. If infection is caused by a gene-for-gene relationship, then the observed pathogenicity is a qualitative trait called ‘virulence’, determined by the presence or absence of resistance genes that give rise to plant immunity. Conversely, if the phenotype of pathogenicity represents a continuous distribution from low to high infection rate, it is a quantitative trait denoted as ‘aggressiveness’ [13]. In order to assess quantitative pathogenicity, symptom development and toxin production are often used [14, 15]. As a quantitative trait, aggressiveness is highly influenced by the environment [15], controlled by multiple genes and expected to be slowly evolving [14, 16]. Variation in aggressiveness within pathogen populations could lead to a gradual adaptation of the pathogen to resistant hosts in the long term [16, 17]. Since F. culmorum has a wide host range, it offers an ideal opportunity to unravel the hidden genetic basis of quantitative pathogenicity. Previous studies consistently reported a large variation in aggressiveness of F. culmorum [18,19,20,21]. Older studies reported several genes to be associated with aggressiveness on the basis of replacement mutants, among them are a secreted lipase FgFGL1 [22], other hydrolyzing enzymes [23] and ABC transporters in F. culmorum and F. graminearum [24, 25]. However, due to the large progress in transcriptomics and genomics in F. graminearum, many other potential genes and genomic regions have been identified [23, 26,27,28,29]. About 70 h after infection, mycotoxins trigger plant tissue necrosis [30]. Accordingly, the trichothecene pathway is a major factor that has been definitively and repeatedly attributed with F. graminearum aggressiveness, in particular the Tri5 gene encoding the first pathway enzyme was reported in multiple studies to be linked with aggressiveness [31,32,33,34]. DON-nonproducing mutants of F. culmorum and F. graminearum cause reduced symptoms compared with wild-type isolates [31, 33, 35]. In either species, while they are still able to infect and cause symptoms on the inoculated spikelets, they do not spread within the head. Therefore, DON is a factor of aggressiveness rather than a pathogenicity factor. More recently, the trichothecene toxin gene cluster was found among 14 regions with distinct signatures of selection in North American F. graminearum populations [29] providing another evidence for the importance of the trichothecenes.

Twelve of the trichothecene genes in F. graminearum are clustered within a 25-kb region on chromosome II that is inherited as one core TRI cluster [36, 37]. Two other loci are involved: the single gene Tri101 and the two-gene Tri1-Tri16 locus, both located in unlinked chromosomal regions [10]. The presence of functional Tri7 and Tri13 genes is required for NIV-producing chemotypes in both F. graminearum [38] and F. culmorum [39]. The trichothecene cluster of the F. culmorum genome is structured similarly to that of F. graminearum, as shown by Schmidt et al. [1] in their analysis of a DON-producing F. culmorum isolate. In the meantime, however, NIV-producing isolates are increasing in frequency in several geographic regions, as has been reported for different Fusarium species [7, 40, 41]. For example, in the UK they are already dominating the F. culmorum population with a frequency of 75.1% compared to 15.0% in F. graminearum at least in maize [7]. Although DON seems to play an important role in the development of disease and aggressiveness, there appears to be a multitude of genes responsible for the broad genetic variation in aggressiveness of Fusarium spp. Indeed, several other genes have been described as aggressiveness factors in candidate-gene association studies in F. culmorum [21, 42] and in a whole-genome based association mapping study of F. graminearum [43].

The understanding of the physiological function of genes, chromosomal organization and related gene products or metabolites has been boosted by the arrival of next generation sequencing technologies, assisting to elucidate the genetic and molecular features that facilitate success of some fungal species as plant pathogens [29, 44, 45]. F. graminearum was the first member of cereal-infecting Fusaria whose complete genome assembly was provided [23], which has since been refined several times [46,47,48]. For F. culmorum, there are two genome sequences available, one from the Australian strain CS7071 [(PRJEB1738) [49], GenBank accession CBMH010000000], whose genome sequence is a fragmented assembly that belongs to the “Wheat Pathogen Genomes” project of the Bioplatforms Australia ( and the other from a United Kingdom field strain, UK99 [50]. The latter is a genome sequence of nearly 42 Mb, after removal of unknown bases just over 39 Mb, which is close to the genome size of F. graminearum strain PH-1 with 36.6 Mb [47]. The draft genome assembly of F. culmorum revealed four chromosomes using alignments to F. graminearum PH-1, but also yielded two small pseudomolecules concatenated by unplaced contigs > 2 kb (“chromosome 5”) and < 2 kb (“chromosome 6”) [50].

Although the full genome sequence of F. culmorum is available, to date only two genes have been associated with pathogenicity or aggressiveness according to the Pathogen-Host-Interaction database ( an ABC transporter and a ROGDI like leucine zipper domain regulating Trichothecene type B biosynthesis [24, 51]. Different methods are now available to further analyze the underlying mechanisms of the large quantitative variation in aggressiveness among strains. Comparative genomics enables the detection of genes responsible for basic pathogenicity across and within Fusarium species [29, 44, 45, 52]. Wang et al. [52] postulated a conserved two-speed genome from the resequencing of three Chinese strains of F. graminearum with the fast sub genome relating to host-pathogen interaction. High-throughput gene deletion studies have also proven useful in identifying genes essential for pathogenicity, e.g., [53, 54]. Neither of these approaches, however, are designed to search for genes specifically responsible for quantitative, within-species differences in aggressiveness or mycotoxin production. Therefore, we applied an association mapping that is widely used for studies of quantitative inheritance of key traits in crops [55]. Candidate-gene based association mapping revealed a few genes that could be associated with differences of aggressiveness in F. culmorum [21, 42], but only a small number of genes with annotated functions and identified by gene-deletion studies in the greenhouse are accessible.

To obtain a broader picture of the genetics of aggressiveness and DON production of this important fungal pathogen in its natural environment, we re-sequenced the genome of 92 isolates of F. culmorum and used 134,789 high quality single-nucleotide polymorphisms (SNPs) to conduct a genome-wide association study (GWAS). We thereby took advantage of the natural diversity produced by historical mutation events among field populations of F. culmorum sampled from different countries and phenotyped them in multi-environment trials. This study is unique since, until now, no resequencing study, e.g., [1, 29, 45] comprised phenotypic data on aggressiveness and DON content from field environments.


Whole-genome data analysis identified 134,789 SNPs

After mapping the sequences from the 92 isolates against the reference genome, we identified 134,789 SNPs. The density of SNPs on F. culmorum chromosomes was 3.66 per kb on average (S2 Table) and varied from 2.96 to 5.19 SNPs per kb. SNP density for each chromosome is shown in S1 Fig. From these, 103,167 SNPs are transitions and 31,622 are transversions, representing a Ts/Tv ratio of 3.26 (S2 Fig). Based on predictions by SnpEff, most of the SNPs were found in upstream or downstream gene regions (38.7 and 38.5%, respectively), while 12.3% are located in transcripts, 4.8% in exons, 4.2% in intergenic regions and the remaining 1.4% in other regions such as splice sites and UTR regions. Taking into account the annotation of the genes, it was possible to predict that 56.8% of the identified SNPs in coding regions are silent, 42.5% of them are missense and 0.7% are nonsense changes. After removing SNPs with missing values > 5% and a minor allele frequency < 5%, 131,605 SNPs were used for further analyses.

Phenotypic data analysis revealed large genotypic variance in aggressiveness

For DON content in wheat (DON-WH), the lowest concentration detected was 0.44 mg kg− 1 for FC50, which is in fact a NIV chemotype and the DON value might be a contamination, and the highest value of 20.7 mg kg− 1 was detected for isolate 9D22 (Table 1, S3 Table). The mean value for DON content was 12.3 mg kg− 1. Regarding aggressiveness on wheat (AGG-WH), the mean value was 15.5% with a minimum of 6.9% and a maximum of 19.5% for isolates S129 and S109, respectively. Values of aggressiveness in rye (AGG-RYE) ranged from 3.4% for isolate S129 to 27.9% for isolate FC95 with a mean of 15.7%. The entry-mean heritabilities were high for all traits, ranging from 0.85 to 0.89, thus providing a good starting point for association analysis.

Table 1 Phenotypic data and variance components for genotype (\( {\sigma}_G^2 \)), genotype-by-environment interaction (\( {\sigma}_{G\times E}^2 \)), error (\( {\sigma}_e^2 \)), and entry-mean heritabilities for aggressiveness in wheat (AGG-WH) and rye (AGG-RYE) and deoxynivalenol content in wheat (DON-WH). Estimates were calculated with arcsine square root transformed data

When color-coding the subpopulations in a scatterplot for aggressiveness on wheat and rye (Fig. 1), no clustering based on origin occurs. For example, the aggressiveness of the Syrian isolates spans the entire range of aggressiveness for both cereal species. In addition, the correlation of aggressiveness between rye and wheat is significant, but not complete across four environments.

Fig. 1
figure 1

Relationship between wheat and rye aggressiveness calculated as BLUES with mean Fusarium head blight (FHB) rating of 92 Fusarium culmorum isolates across four environments; r = coefficient of correlation, p = probability of error

The genotypic variance (\( {\sigma}_G^2\Big) \) as well as the genotype-by-environment interaction variance (\( {\sigma}_{G\times E}^2\Big) \) were significantly (P < 0.001) different from zero for all traits. The correlation between aggressiveness in wheat and DON production was significant (r = 0.77, p < 0.001) (S3 Fig). Furthermore, the correlation between aggressiveness in rye and aggressiveness or DON production in wheat were moderate and significant (r2 = 0.62 and r2 = 0.54, respectively, p < 0.001).

Detection of population structure

Classical multidimensional scaling (MDS) of the distance matrix was performed and the resulting principal components revealed some degree of population structure (S4 Fig). Resulting clusters represent the geographical origin of the populations and therefore the population membership of the isolates. Particularly the Syrian population is grouped apart from the German and the Russian clusters. The first principal component explained 25.93% of the total genetic variance, the second, third and fourth an additional 8.67, 5.72 and 4.65%, respectively. Consequently, principal components were added to the mixed linear model in order to correct for population structure (S5 Fig).

Identification of SNPs associated with DON production and aggressiveness

The association analysis resulted in the identification of 27 SNPs that were associated with DON production (Table 2).

Table 2 Detected QTL for deoxynivalenol production in wheat (DON-WH), aggressiveness in wheat (AGG-WH) and rye (AGG-RYE)

Two SNPs are significantly associated with the aggressiveness of F. culmorum in rye and 20 SNPs with aggressiveness in wheat. Interestingly, 10 of the SNPs that are significantly associated with DON-WH are also significantly associated with AGG-WH (Fig. 2) and the two SNPs that are associated with AGG-RYE are also associated with DON-WH. Surprisingly, there are no shared SNPs between aggressiveness in both hosts.

Fig. 2
figure 2

Venn diagram for the number of significantly associated SNPs for deoxynivalenol (DON) content in wheat and aggressiveness (AGG) in wheat and in rye

The largest QTL identified for aggressiveness in rye, located on chromosome I, explains 17.2% of the genotypic variance. For aggressiveness in wheat, the largest QTL on chromosome II explains 33.4% of the genotypic variance and another, but linked QTL is the largest for DON concentration, explaining 48.4% of the genotypic variance (Table 2). When the allele frequency of the SNPs with the highest proportion of explained variance in the whole population was assessed, we found that most of the isolates carry the high aggressiveness allele (Fig. 3).

Fig. 3
figure 3

Boxplots of the SNPs explaining the largest proportion of genotypic variance for deoxynivalenol production in wheat and aggressiveness in wheat and rye with the two SNP alleles and the number of isolates bearing the respective allele (count); horizontal line within boxes = median, and  = outliers

Two SNPs, FC2_6904196 and FC4_3857293, associated with aggressiveness in wheat are located within two genes, FCUL_06122 and FCUL_10689, and cause amino acid changes (Table 3). Furthermore, the SNP FC1_7462132 associated with aggressiveness and DON production in wheat is a splice region variant (intron variant). All mentioned protein sequences have an unknown function.

Table 3 SNPs that are significantly associated with aggressiveness (AGG-WH) and deoxynivalenol content (DON-WH) in wheat, explaining more than 2% of genotypic variance and located in gene sequences

We also identified a SNP (FC2_5593616) that is significantly associated with DON production and explains 48.4% of the genotypic variance. Although this SNP explains just 1.33% of the genotypic variance of aggressiveness in wheat, it is the only SNP with a significance value above the Bonferroni threshold (Fig. 4). This SNP is located in the gene FCUL_05634 that encodes for a Cytochrome P450 [56, 57] with a length of 520 amino acids. The translated sequence of this gene shares 95.7% of identity with the sequence of Tri4 from F. graminearum and 96.9% with the respective protein from F. pseudograminearum, as confirmed by reciprocal blastp. Additionally, we identified the gene cluster for biosynthesis of trichothecenes in F. culmorum using the antiSMASH pipeline, and each of the core cluster candidate genes was confirmed by blastp with the F. graminearum sequences as query. The FCUL_05634 gene is the third gene of the core cluster and based on these data, we conclude that the gene FCUL_05634 is orthologous to the Tri4 gene previously reported in other species [58]. Given that the gene cluster where the DON biosynthesis genes are located is in a linkage disequilibrium (LD) block (Fig. 5), the associated SNP could either represent an effect of the Tri4 gene or be linked to any other gene within the cluster. It thus represents a QTL that contributes to aggressiveness and DON content of F. culmorum.

Fig. 4
figure 4

Manhattan plots of the genome-wide association study of (A) deoxynivalenol content in wheat, (B) aggressiveness in wheat, (C) aggressiveness in rye. The dotted line represents the Bonferroni significance threshold of p < 3,8 × 10−7 and the horizontal blue line the explorative significance threshold of P-value < 0.0001

Fig. 5
figure 5

Manhattan plots for deoxynivalenol (DON) production and aggressiveness in wheat and of the genomic segment from 5.570.000 bp to 5.620.000 bp on chromosome 2 (49.145 bp, 467 SNPs) and results from the genome browser with the proposed genes in this region and LD block


Fusarium culmorum is, besides F. graminearum, one of the most frequently occurring Fusarium species associated with Fusarium head blight disease in Europe [4, 5]. However, not much is known about the basis of its pathogenicity or aggressiveness. The availability of an annotated whole-genome assembly [50] opens new routes towards an understanding of aggressiveness. Association mapping is based on a diversity panel of individuals with complex relatedness that have captured historical recombination events during evolution [55]. To identify significant associations between genotypes and phenotypes, high-density genetic markers are necessary, because the method is based on linkage disequilibrium between markers and causative nucleotides. For the first time in F. culmorum, we have analyzed the whole-genome sequences of 92 isolates aiming for an association study with two important phenotypic features, aggressiveness and DON production.

Large genetic variation of aggressiveness and DON production

The results of this study support the outcome of previous studies revealing large genetic variation of aggressiveness and DON production among F. culmorum isolates [18,19,20, 59]. The genetic variation was estimated statistically by partitioning the phenotypic variance measured in our field experiment across four location × year combinations (=environments) into the components of genotype, genotype-by-environment interaction, and error. The genetic component was the largest for all three traits (Table 1). The genotype-by-environment interaction, however, had a similar impact for aggressiveness in rye, but was smaller for the traits in wheat. The genetic causes for this large genetic variation in F. culmorum are not clear, because there is no sexual stage known in this fungal species. However, neutral molecular markers also revealed a large variation among 186 isolates [60]. This points to a recombining structure of F. culmorum [61], which may be explained by one or any of the following hypotheses: (1) cryptic variation by a sexual stage; (2) the teleomorph has become extinct; or (3) asexual variation, e.g. by heterokaryosis. The first two hypotheses are supported by the fact that F. culmorum has two transcribed mating type gene homologues [62, 63] occurring at similar frequencies in natural populations [61].

In our study, the aggressiveness in rye reached higher values than that in wheat. This does not reflect the reality when testing larger arrays of host genotypes [64]. From our previous knowledge that rye is generally less susceptible to FHB than wheat, we have chosen a highly susceptible rye genotype to get an optimal differentiation of the isolates.

Several QTLs identified for aggressiveness and DON production

Because aggressiveness is inherited quantitatively, each gene contributing to this complex trait is expected to explain only a small proportion of the genotypic variance. This has been recently shown by candidate-gene based association mapping [21, 42] where the associated candidate genes had contributions to the genotypic variation varying from 6.05 to 16.05%. Moreover, Castiblanco et al. [42] attributed different candidate genes to aggressiveness of wheat and rye. Especially the gene cutinase (CUT) was significantly associated with aggressiveness towards rye but not towards wheat. This was attributed by the authors to the different composition of the cuticles of wheat and rye. Moreover, rye and wheat were overlapping only in 1 year (2015), in 2014 only wheat and in 2016 only rye was analyzed. This might have caused differences due to the highly significant genotype × environment interaction observed also in this study that is typical for the pathosystem. The complex genetic control underlying aggressiveness is corroborated by the finding of this study, that no QTL for aggressiveness were common between wheat and rye.

When the allele frequency of the SNPs was assessed in the whole population, we found that most of the isolates carry the high aggressiveness allele (Fig. 3). This might be related to the fact that the isolates used in this study were sampled from symptomatic heads of commercial wheat fields where low aggressiveness isolates are hardly found.

Two SNPs significantly associated with aggressiveness in wheat were found to be the proteins FCUL_06122 and FCUL_10689 with no predicted function in the F. culmorum genome. A blast of the genomic and the peptide sequences of both genes against the NCBI database and the UniProt database did not reveal for FCUL_06122 any more specific functions. The best hits, all on Fusarium species, also resulted in “uncharacterized protein”. For FCUL_10689 the best hits on UniProts were also unknown proteins, but in Fusarium longipes, Fusarium sporotrichioides, and Fusarium vanettenii (syn. Fusarium solani f. sp. pisi) peptide sequences were annotated as Transaldolase. Transaldolase (EC was associated with the sugar and energy metabolism and one of the most abundant protein spots identified in B. cinerea mycelium [65]. Interestingly, a transaldolase-like protein (TAL 1) of Cochliobolus heterostrophus was identified as a Hog1-dependent gene playing an important role in protection against oxidative stress [66]. And Hog-1 was significantly associated with aggressiveness and DON production in our candidate-gene association mapping of 100 F. culmorum isolates [21] of which 92 were included in this study.

A large part of aggressiveness is explained by QTL responsible for DON production

The correlation between aggressiveness and DON production in wheat was significant, corroborating previous studies [21, 42, 67, 68]. Interestingly, all QTL associated with aggressiveness in rye and 10 out of 20 QTL associated with aggressiveness in wheat were also associated with DON content in wheat (Fig. 2). This points to a strong effect of DON on the infection process in the head. To date, however, cause and effect could not been separated. Gang et al. [18], for example, proposed that higher DON levels might just be the consequence of higher fungal biomass. When additionally estimating fungal biomass by measuring exoantigen content, the correlation between aggressiveness and DON content per unit fungal protein was non-significant [69]. On the other hand, strains with higher DON production were found to be more aggressive in some wheat genotypes [70, 71]. Whether this is related to the level of DON production or other genetic factors remains unknown. Importantly, while DON-nonproducing isolates are not impaired in very early infection events and form infection structures and induce the development of necrotic lesions [72], they are slower at developing the disease and are not able to grow outside the inoculation site [31, 68]. Thus, DON content might be important for host tissue proliferation [73, 74] until a threshold is reached in the plant tissue, where more DON does not necessarily lead to more symptoms.

Results of this association study strongly promote an important role of DON content for aggressiveness. The SNP FC2_5593616, explaining 48.4% of the total genotypic variation in DON content and as the only SNP exceeding the conservative Bonferroni threshold, is located in a gene with a very high identity (95.7%) with the published sequence of the Tri4 gene of F. graminearum. This gene encodes a multifunctional cytochrome P450 monooxygenase for different oxygenation and hydroxylation steps in the trichothecene pathway leading to DON [58] and is, therefore, an indispensable component for trichothecene biosynthesis [75]. Alternatively, our most prominent SNP FC2_5593616 could also be in linkage disequilibrium with another gene from this cluster and not directly represent Tri4, but it is nonetheless an important factor in DON production having a high impact on aggressiveness. However, our study also shows that the large differences in aggressiveness among F. culmorum isolates most probably involve a multitude of quantitative loci of which Tri4 is only one, although an important one. Accordingly, Talas et al. [43] found 26 and 17 genes significantly associated with aggressiveness and DON production, respectively, among 119 isolates of F. graminearum, each explaining 8 to 24% of the genotypic variance. As in our study, most genes detected did not have a known function.


The high phenotypic variation of aggressiveness and DON production of F. culmorum is most likely caused by a multitude of genes of whom those of the DON pathway have an important function. Among them, we could show here that the Tri4 gene or a gene in close linkage to Tri4 explained nearly half of the total genotypic variation. More functional studies will be needed in future to further decipher the molecular mechanisms underlying their associations with fungal aggressiveness.


Fusarium culmorum isolates

A total of 92 isolates of F. culmorum collected from three different geographic origins, were used (S1 table). A group of 33 isolates were collected in three commercial fields in Germany, 18 isolates originated from a Russian field population and a transect population of 25 isolates was collected in Syria. Additionally, 19 isolates from the international collection of F. culmorum sampled from 12 different countries were included. These isolates were available in the rye group of the State Breeding Institute at the University of Hohenheim and came from single-spore isolation. Isolates were stored as colonized agar plugs in sterile water in 2 ml Eppendorf tubes at + 6 °C.

Phenotypic data analyses of F. culmorum aggressiveness

Aggressiveness of F. culmorum in wheat and rye, and DON content in wheat were assessed. Raw phenotypic data from previous field studies were used as described in detail by Castiblanco et al. [21, 42]. Briefly, two host species with one cultivar each were used for inoculation: a winter wheat (T. aestivum L.) cultivar that is moderately susceptible (“Inspiration”, KWS LOCHOW SE, Bergen, Germany) and a highly susceptible, homogeneous rye single cross (S. cereale L., “L2177-P × L2184-N”, HYBRO Saatzucht GMBH & Co., KG, Schenkenberg, Germany). Inoculations and phenotyping were done at two locations: Hohenheim (HOH, 48° 42′ 54″ N, 9° 12′ 58″ E, 400 m) and Oberer Lindenhof (OLI, 48° 28′ 25″ N, 9° 18′ 12″E, altitude 700 m) in Germany in 2 years (2014 and 2015 for wheat as host and 2015 and 2016 for rye), which is equivalent to four environments (=location-year combinations).

Seeds were grown in two-row plots (rye) or three-row plots (wheat), respectively, with 1 m length and a row distance of 0.21 m. To decrease the drifting or secondary spore dispersal and avoid possible interference among plots, a chessboard-like design was used to arrange the entry plots that were bordered by long-strawed rye in the rye experiment and triticale (xTriticosecale) in the wheat experiment. Plots were sown with 220 kernels m− 2 resulting in a homogeneous dense plant stand. Experiments were designed following an incomplete block design (alpha design) with two replications per environment and an incomplete block size of 10.

Inoculum of all isolates was prepared in shaking liquid cultures according to the procedure of Reid et al. [76] and aliquots frozen at -80 °C. Before application, the samples were thawed in water at about 20 °C. Inoculation was performed by conidiospores in a dose of approximately 100 ml suspension per square meter and a concentration of 2 × 105 conidiospores ml− 1. Inoculum for each isolate was sprayed onto the heads using a hand atomizer with constant air pressure of 3 bar produced by a tractor to ensure full coverage of all heads of the plot with the same dosage. All plots flowered simultaneously, because only one homogeneous rye or wheat genotype was used. This allowed inoculation and ratings for all plots at the same dates per location.

Aggressiveness values came from the arithmetic mean from at least three symptom ratings taken over time. The ratings started with the beginning of symptom development and were continued at 3- to 5-day intervals until the beginning of yellow ripening stages. Ratings represented the percentage of infected spikelets per ear (severity) and the percentage of infected ears per plot (incidence) in a single value from 0 to 100% [21, 42]. For determining DON content, wheat plots were harvested by hand at full ripening, carefully threshed in a single-head thresher (Walter-Wintersteiger, Austria) and cleaned with reduced wind speed. Cleaned wheat grain was ground in a commercial laboratory mill with a sieve size of 1 mm. The coarse meal was analyzed by a commercially available immunotest (Ridascreen FAST DON®, R-Biopharm AG, Darmstadt, Germany).

Statistical analyses of phenotypic data were done using the software R [77], outliers were detected in each environmental dataset with the BH-MADR method (Bonferroni–Holm with re-scaled MAD standardized residuals) suggested by Bernal-Vasquez et al. [78] and deleted accordingly. All data was arcsin transformed to meet the prerequisite of normality. Entry-mean heritabilities were computed as suggested by Piepho and Möhring [79]. The following mixed model was used to obtain the BLUEs (best linear unbiased estimators) for each isolate:

$$ \mathrm{G}+\mathrm{E}+\mathrm{G}\times \mathrm{E}+\mathrm{R}+\mathrm{B} $$

where G, E, R and B denote genotype, environment, replication and block, respectively, and G, E and R were modeled as fixed effects and G × E and B as random effects. R was modelled as nested within E, and B nested within R. Host (wheat or rye) was not considered in the model, given that pathogen fitness could depend on the host and since plant pathogens may use different strategies to infect different hosts [31]. Instead, we compared the outcome of our association mapping from the two hosts. The variance components were estimated using restricted maximum likelihood (REML) method and their significance was assessed by comparing the models using Likelihood Ratio test [80]. Pearson’s correlation (r2) between traits was calculated based on the BLUEs. The software package ASReml-R 3.0 [81] was used to fit the linear mixed models. All plots, except the local GWAS plot (see below), were created using basic R functionalities [77].

DNA extraction, whole-genome resequencing, annotation and SNP calling

Plugs from each isolate were cultivated on PDA plates during 5 days at 20 °C. From each isolate, 200 mg of mycelia was scraped off from several plates with a sterile spatula and dried out with silica gel for 1 week. DNA extraction was performed following the cetyltrimethylammonium bromide (CTAB) method. DNA samples were run on a 1% agarose gel to check DNA stability and integrity with a O’GeneRuler 1 kb DNA ladder (Thermo Scientific) and the 23 kb DNA Molecular Weight Marker II (Sigma–Aldrich). DNA quality was measured with nanodrop, and quantity was measured with QuantiFluor® dsDNA System (Promega). Finally, the DNA was concentrated to 50 ng/μl. DNA samples were sent to Novogene Corporation (Beijing, China) where Illumina libraries were prepared for whole-genome sequencing using a Illumina HiSeq X with 150 bp paired-end sequencing. The resulting sequencing libraries had an average insert size of 300 bp.

Pair-end reads from all isolates were mapped to the reference genome with the Burrows-Wheeler Aligner software package (BWA, release 0.7.17) [82] using the MEM algorithm. As a reference genome, we used the draft assembly of F. culmorum strain UK99 that has almost 42 Mbp and consists of the 4 F. culmorum chromosomes plus two scaffolds of concatenated sequences that could not be assigned to the chromosomes [50]. The NGSEP software package [83] was used to identity SNPs, small indels, and structural variants in the assembled sequences, to make the functional annotation of each variant and to genotype each sample. Only bi-allelic SNPs with a quality call score larger than 40 (Phred scale), a minor allele frequency (MAF) of ≥0.05 and a minimum coverage of 95% of genotyped individuals were kept after filtering. Descriptive statistics such as SNPs density and position in F. culmorum chromosomes were calculated using the SnpEff tool [84] and vcftools [85] and mapping metrics were calculated with BBMap aligner [86]. Upstream and downstream gene region size was set to 5 kb. Around 9 millions of pair-end reads on average per each isolate were obtained after sequence quality check. From these reads, 91.6% were successfully aligned to the reference genome of F. culmorum strain UK99 covering an average of 87.7% of it, with a minimum of 81% and a maximum of 91.6%. The average sequencing depth was 33X and the mean insert length was 287.8 bp and ranged from 2 up to 636.3 kb. During the filtering process, SNPs that were mapped to repetitive sequences and copy number variations were excluded. Genotypic quality control of the final SNPs was performed using the R package GenABEL [87].

Genome-wide association analysis

Given that isolates coming from different geographic regions were used for this study, stratification was assessed by principal component analysis with the R package GenABEL [87]. For this, a genomic relationship matrix between all pairs of individuals was computed and from this a distance matrix was constructed. Subsequently, a classical multidimensional scaling of the distance matrix was performed to obtain the principal components analysis of genetic variation.

Association analysis between SNP markers and traits was done with GenABEL [87], using a linear mixed model that included four principal coordinates for wheat and two for rye and a kinship matrix to adjust for possible population stratification. To control for multiple testing, Bonferroni-corrected significance threshold and, in addition, an explorative threshold of p < 0.0001 were used. All significantly associated markers of a trait were fitted into a linear model in the order of the strength of their association, i.e. with the most strongly associated marker fitted first [88]. The proportion of genotypic variance explained by all significantly associated SNPs was estimated as the ratio of the adjusted R2 of this linear model to the heritability of the trait \( {p}_G=\frac{R_{adj}^2}{h^2} \) [89]. The proportion of genotypic variance explained by each single SNP was likewise derived from the sums of squares from this model (pG-combined). In addition, we also calculated the proportion of explained genotypic variance of each marker when fitted singly in a linear model (pG-single). The combination of the two values allows to draw conclusions whether significantly associated markers from one genomic region identify the same QTL, as the joint fit of all significantly associated markers corrects for collinearity among them. The additive effect (α-effect) of each individual significant SNP was derived from a linear model in which just the respective SNP was considered. The local GWAS plot was created using the R package LDheatmap [90].

To analyze the sequences that harbor the significantly associated SNPs, sequences of 100 bp around the SNPs were extracted with BEDTools [91] and genomic position and allele variation were confirmed in the F. culmorum genome. Translated sequences were compared with the protein sequences of closely related species. For the identification of secondary metabolite biosynthesis gene clusters in the whole genome of F. culmorum the web server antiSMASH 3.0 was used [92].

Availability of data and materials

All details on the analyzed isolates can be found in S1 Table, all phenotypic data are available as best linear unbiased estimates (BLUES) in the S3 Table (Supplementary Information).

Sequencing data generated and analyzed during the current study are available in the NCBI (Bethesda, Maryland USA) database as BioProject ID PRJNA749912 [, release at Aug 08, 2021].





Best linear unbiased estimator


Deoxyribonucleic acid




Fusarium culmorum


Fusarium head blight


Genome-wide association study




Potato dextrose agar


Quantitative trait locus


Single nucleotide polymorphisms

Tri :

Trichothecene gene




  1. Schmidt R, Durling MB, de Jager V, Menezes RC, Nordkvist E, Svatoš A, et al. Deciphering the genome and secondary metabolome of the plant pathogen Fusarium culmorum. FEMS Microbiol Ecol. 2018;94:fiy078.

    Article  Google Scholar 

  2. Parry DW, Jenkinson P, McLeod L. Fusarium ear blight (scab) in small grain cereals-a review. Plant Pathol. 1995;44(2):207–38.

    Article  Google Scholar 

  3. Xu X-M, Parry DW, Nicholson P, Thomsett MA, Simpson D, Edwards SG, et al. Predominance and association of pathogenic fungi causing Fusarium ear blight in wheat in four European countries. Eur J Plant Pathol. 2005;112(2):143–54.

    Article  Google Scholar 

  4. Pasquali M, Beyer M, Logrieco A, Audenaert K, Balmas V, Basler R, et al. European database of Fusarium graminearum and Fusarium culmorum trichothecene genotypes. Front Microbiol. 2016;7:406.

    Article  Google Scholar 

  5. Birr T, Hasler M, Verreet JA, Klink H. Composition and predominance of Fusarium species causing Fusarium head blight in winter wheat grain depending on cultivar susceptibility and meteorological factors. Microorganisms. 2020;8(4):617.

    Article  CAS  PubMed Central  Google Scholar 

  6. Bottalico A. Fusarium diseases of cereals: species complex and related mycotoxin profiles, in Europe. J Plant Pathol. 1998;80:85–103.

    CAS  Google Scholar 

  7. Basler R. Diversity of Fusarium species isolated from UK forage maize and the population structure of F. graminearum from maize and wheat. PeerJ. 2016;4:e2143.

    Article  Google Scholar 

  8. Scherm B, Balmas V, Spanu F, Pani G, Delogu G, Pasquali M, et al. Fusarium culmorum: causal agent of foot and root rot and head blight on wheat. Mol Plant Pathol. 2013;14(4):323–41.

    Article  CAS  PubMed  Google Scholar 

  9. Becher R, Miedaner T, Wirsel SGR. Biology, diversity, and management of FHB-causing Fusarium species in small-grain cereals. In: Kempken F, editor. Agricultural applications. The Mycota XI. 2nd ed. Berlin: Springer; 2013. p. 199–241.

    Chapter  Google Scholar 

  10. Torres AM, Palacios SA, Yerkovich N, Palazzini JM, Battilani P, Leslie JF, et al. Fusarium head blight and mycotoxins in wheat: prevention and control strategies across the food chain. World Mycotoxin J. 2019;12(4):333–55.

    Article  CAS  Google Scholar 

  11. Suchowilska E, Kandler W, Sulyok M, Wiwart M, Krska R. Mycotoxin profiles in the grain of Triticum monococcum, Triticum dicoccum and Triticum spelta after head infection with Fusarium culmorum. J Sci Food Agric. 2010;90(4):556–65.

    Article  CAS  PubMed  Google Scholar 

  12. Agrios BN. Plant pathology. 3rd ed. San Diego: Academic; 1988.

    Google Scholar 

  13. Van der Plank J. Plant diseases: epidemics and control. New York: Academic; 1963.

    Google Scholar 

  14. Pariaud B, Ravigné V, Halkett F, Goyeau H, Carlier J, Lannou C. Aggressiveness and its role in the adaptation of plant pathogens. Plant Pathol. 2009;58(3):409–24.

    Article  Google Scholar 

  15. Lannou C. Variation and selection of quantitative traits in plant pathogens. Annu Rev Phytopathol. 2012;50(1):319–38.

    Article  CAS  PubMed  Google Scholar 

  16. Pariaud B, Robert C, Goyeau H, Lannou C. Aggressiveness components and adaptation to a host cultivar in wheat leaf rust. Phytopathology. 2009;99(7):869–78.

    Article  PubMed  Google Scholar 

  17. McDonald BA, Linde C. Pathogen population genetics, evolutionary potential, and durable resistance. Annu Rev Phytopathol. 2002;40:349–79.

    Article  CAS  PubMed  Google Scholar 

  18. Gang G, Miedaner T, Schuhmacher U, Schollenberger M, Geiger HH. Deoxynivalenol and nivalenol production by Fusarium culmorum isolates differing in aggressiveness toward winter rye. Phytopathology. 1998;88(9):879–84.

    Article  CAS  PubMed  Google Scholar 

  19. Muthomi JW, Schütze A, Dehne HW, Mutitu EW, Oerke EC. Characterization of Fusarium culmorum isolates by mycotoxin production and aggressiveness to winter wheat. J Plant Dis Prot. 2000;107:113–23.

    CAS  Google Scholar 

  20. Miedaner T, Schilling AG, Geiger HH. Molecular genetic diversity and variation for aggressiveness in populations of Fusarium graminearum and Fusarium culmorum sampled from wheat fields in different countries. J Phytopathol. 2001;149(11-12):641–8.

    Article  CAS  Google Scholar 

  21. Castiblanco V, Marulanda JJ, Würschum T, Miedaner T. Candidate gene based association mapping in Fusarium culmorum for field quantitative pathogenicity and mycotoxin production in wheat. BMC Genet. 2017;18(1):49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Voigt CA, Schafer W, Salomon S. A secreted lipase of Fusarium graminearum is a virulence factor required for infection of cereals. Plant J. 2005;42(3):364–75.

    Article  CAS  PubMed  Google Scholar 

  23. Cuomo CA, Güldener U, Xu JR, Trail F, Turgeon BG, Di Pietro A, et al. The Fusarium graminearum genome reveals a link between localized polymorphism and pathogen specialization. Science. 2007;317(5843):1400–2.

    Article  CAS  PubMed  Google Scholar 

  24. Skov J, Lemmens M, Giese H. Role of a Fusarium culmorum ABC transporter (FcABC1) during infection of wheat and barley. Physiol Mol Plant Pathol. 2004;64(5):245–54.

    Article  CAS  Google Scholar 

  25. Abou Ammar G, Tryono R, Döll K, Karlovsky P, Deising HB, Wirsel SG. Identification of ABC transporter genes of Fusarium graminearum with roles in azole tolerance and/or virulence. PloS One. 2013;8:e79042.

    Article  CAS  Google Scholar 

  26. Walter S, Nicholson P, Doohan FM. Action and reaction of host and pathogen during Fusarium head blight disease. New Phytol. 2010;185(1):54–66.

    Article  CAS  PubMed  Google Scholar 

  27. Urban M, King R, Hassani-Pak K, Hammond-Kosack KE. Whole-genome analysis of Fusarium graminearum insertional mutants identifies virulence associated genes and unmasks untagged chromosomal deletions. BMC Genomics. 2015;16(1):261.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Boedi S, Berger H, Sieber C, Münsterkötter M, Maloku I, Warth B, et al. Comparison of Fusarium graminearum transcriptomes on living or dead wheat differentiates substrate-responsive and defense-responsive genes. Front Microbiol. 2016;7:1113.

    Article  Google Scholar 

  29. Kelly AC, Ward TJ. Population genomics of Fusarium graminearum reveals signatures of divergent evolution within a major cereal pathogen. PLoS One. 2018;13(3):e0194616.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Kang Z, Buchenauer H. Immunocytochemical localization of Fusarium toxins in infected wheat spikes by Fusarium culmorum. Physiol Mol Plant Pathol. 1999;55(5):275–88.

    Article  CAS  Google Scholar 

  31. Proctor RH, Hohn TM, McCormick SP. Reduced virulence of Gibberella zeae caused by disruption of a trichothecene toxin biosynthetic gene. Mol Plant-Microbe Interact. 1995;8(4):593–601.

    Article  CAS  PubMed  Google Scholar 

  32. Proctor RH, Hohn TM, McCormick SP. Restoration of wild-type virulence to Tri5 disruption mutants of Gibberella zeae via gene reversion and mutant complementation. Microbiology. 1997;143(8):2583–91.

    Article  CAS  PubMed  Google Scholar 

  33. Desjardins AE, Proctor RH, Bai G, McCormick SP, Shaner G, Buechley G, et al. Reduced virulence of trichothecene-nonproducing mutants of Gibberella zeae in wheat field tests. Mol Plant-Microbe Interact. 1996;9(9):775–81.

    Article  CAS  Google Scholar 

  34. Maier FJ, Miedaner T, Hadeler B, Felk A, Salomon S, Lemmens M, et al. Involvement of trichothecenes in fusarioses of wheat, barley and maize evaluated by gene disruption of the trichodiene synthase (Tri5) gene in three field isolates of different chemotype and virulence. Mol Plant Pathol. 2006;7(6):449–61.

    Article  CAS  PubMed  Google Scholar 

  35. Scherm B, Orru M, Balmas V, Spanu F, Azara E, Delogu G, et al. Altered trichothecene biosynthesis in Tri6-silenced transformants of Fusarium culmorum influences the severity of crown and foot rot on durum wheat seedlings. Mol Plant Pathol. 2011;12(8):759–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Kimura M, Tokai T, O’Donnell K, Ward TJ, Fujimura M, Hamamoto H, et al. The trichothecene biosynthesis gene cluster of Fusarium graminearum F15 contains a limited number of essential pathway genes and expressed non-essential genes. FEBS Lett. 2003;539(1-3):105–10.

    Article  CAS  PubMed  Google Scholar 

  37. Brown DW, Dyer RB, McCormick SP, Kendra DF, Plattner RD. Functional demarcation of the Fusarium core trichothecene gene cluster. Fungal Genet Biol. 2004;41(4):454–62.

    Article  CAS  PubMed  Google Scholar 

  38. Lee T, Han YK, Kim KH, Yun SH, Lee YW. Tri13 and Tri7 determine deoxynivalenol-and nivalenol-producing chemotypes of Gibberella zeae. Appl Environ Microbiol. 2002;68(5):2148–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Jennings P, Coates ME, Turner JA, Chandler EA, Nicholson P. Determination of deoxynivalenol and nivalenol chemotypes of Fusarium culmorum isolates from England and Wales by PCR assay. Plant Pathol. 2004;53(2):182–90.

    Article  CAS  Google Scholar 

  40. Gale LR, Harrison SA, Ward TJ, O'Donnell K, Milus EA, Gale SW, et al. Nivalenol-type populations of Fusarium graminearum and F. asiaticum are prevalent on wheat in southern Louisiana. Phytopathology. 2011;101(1):124–34.

    Article  CAS  PubMed  Google Scholar 

  41. Malihipour A, Gilbert J, Piercey-Normore M, Cloutier S. Molecular phylogenetic analysis, trichothecene chemotype patterns, and variation in aggressiveness of Fusarium isolates causing head blight in wheat. Plant Dis. 2012;96(7):1016–25.

    Article  CAS  PubMed  Google Scholar 

  42. Castiblanco V, Castillo HE, Miedaner T. Candidate genes for aggressiveness in a natural Fusarium culmorum population greatly differ between wheat and rye head blight. J Fungi. 2018;4:14.

    Article  Google Scholar 

  43. Talas F, Kalih R, Miedaner T, McDonald BA. Genome-wide association study identifies novel candidate genes for aggressiveness, deoxynivalenol production and azole sensitivity in natural field populations of Fusarium graminearum. Mol Plant-Microbe Interact. 2016;29(5):417–30.

    Article  CAS  PubMed  Google Scholar 

  44. Rep M, Kistler HC. The genomic organization of plant pathogenicity in Fusarium species. Curr Opin Plant Biol. 2010;13(4):420–6.

    Article  CAS  PubMed  Google Scholar 

  45. Walkowiak S, Rowland O, Rodrigue N, Subramaniam R. Whole genome sequencing and comparative genomics of closely related Fusarium head blight fungi: Fusarium graminearum, F meridionale and F asiaticum. BMC Genomics. 2016;17(1):1014.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Wong P, Walter M, Lee W, Mannhaupt G, Münsterkötter M, Mewes HW, et al. FGDB: revisiting the genome annotation of the plant pathogen Fusarium graminearum. Nucleic Acids Res. 2010;39:D637–9.

    Article  Google Scholar 

  47. King R, Urban M, Hammond-Kosack MCU, Hassani-Pak K, Hammond-Kosack KE. The completed genome sequence of the pathogenic ascomycete fungus Fusarium graminearum. BMC Genomics. 2015;16(1):544.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. King R, Urban M, Hammond-Kosack KE. Annotation of Fusarium graminearum (PH-1) version 5.0. Genome Announc. 2017;5:e01479–16.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Sperschneider J, Gardiner DM, Thatcher LF, Lyons R, Singh KB, Manners JM, et al. Genome-wide analysis in three Fusarium pathogens identifies rapidly evolving chromosomes and genes associated with pathogenicity. Genome Biol Evol. 2015;7(6):1613–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Urban M, King R, Andongabo A, Maheswari U, Pedro H, Kersey P, et al. First draft genome sequence of a UK strain (UK99) of Fusarium culmorum. Genome Announc. 2016;4:e00771–16.

    Article  Google Scholar 

  51. Spanu F, Scherm B, Camboni I, Balmas,V, Pani, G, Oufensou S, et al. FcRav2, a gene with a ROGDI domain involved in Fusarium head blight and crown rot on durum wheat caused by Fusarium culmorum. Mol Plant Pathol. 2018;19:677–88, 3,

  52. Wang Q, Jiang C, Wang C, Chen C, Xu JR, Liu H. Characterization of the two-speed subgenomes of Fusarium graminearum reveals the fast-speed subgenome specialized for adaption and infection. Front Plant Sci. 2017;8:140.

    PubMed  PubMed Central  Google Scholar 

  53. Son H, Seo YS, Min K, Park AR, Lee J, Jin JM, et al. A phenome-based functional analysis of transcription factors in the cereal head blight fungus, Fusarium graminearum. PLoS Pathog. 2011;7(10):e1002310.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Wang C, Zhang S, Hou R, Zhao Z, Zheng Q, Xu Q, et al. Functional analysis of the kinome of the wheat scab fungus Fusarium graminearum. PLoS Pathog. 2011;7(12):e1002460.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Lipka AE, Kandianis CB, Hudson ME, Yu J, Drnevich J, Bradbury PJ, et al. From association to prediction: statistical methods for the dissection and selection of complex traits in plants. Curr Opin Plant Biol. 2015;24:110–8.

    Article  PubMed  Google Scholar 

  56. Alexander NJ, Hohn TM, McCormick SP. The Tri11 gene of Fusarium sporotrichioides encodes a cytochrome P450 monooxygenase required for C-15 hydroxylation in trichothecene biosynthesis. Appl Environ Microbiol. 1998;64(1):221–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Hohn TM, Desjardins AE, McCormick SP. The Tri4 gene of Fusarium sporotrichioides encodes a cytochrome P450 monooxygenase involved in trichothecene biosynthesis. Mol Gen Genet. 1995;248(1):95–102.

    Article  CAS  PubMed  Google Scholar 

  58. Tokai T, Koshino H, Takahashi-Ando N, Sato M, Fujimura M, Kimura M. Fusarium Tri4 encodes a key multifunctional cytochrome P450 monooxygenase for four consecutive oxygenation steps in trichothecene biosynthesis. Biochem Biophys Res Commun. 2007;353(2):412–7.

    Article  CAS  PubMed  Google Scholar 

  59. Miedaner T, Gang G, Geiger HH. Quantitative-genetic basis of aggressiveness of 42 isolates of Fusarium culmorum for winter rye head blight. Plant Dis. 1996;80(5):500–4.

    Article  Google Scholar 

  60. Miedaner T, Caixeta F, Talas F. Head-blighting populations of Fusarium culmorum from Germany, Russia, and Syria analyzed by microsatellite markers show a recombining structure. Eur J Plant Pathol. 2013;137(4):743–52.

    Article  Google Scholar 

  61. Obanor F, Erginbas-Orakci G, Tunali B, Nicol JM, Chakraborty S. Fusarium culmorum is a single phylogenetic species based on multilocus sequence analysis. Fungal Biol. 2010;114(9):753–65.

    Article  PubMed  Google Scholar 

  62. Mishra PK, Fox RT, Culham A. Inter-simple sequence repeat and aggressiveness analyses revealed high genetic diversity, recombination and long-range dispersal in Fusarium culmorum. Ann Appl Biol. 2003;143(3):291–301.

    Article  CAS  Google Scholar 

  63. Tóth B, Mesterházy Á, Nicholson P, Téren J, Varga J. Mycotoxin production and molecular variability of European and American isolates of Fusarium culmorum. In: Mulè G, Bailey JA, Cooke BM, Logrieco A, editors. Molecular diversity and PCR-detection of toxigenic Fusarium species and ochratoxigenic fungi. Dordrecht: Springer; 2004. p. 587–99.

    Chapter  Google Scholar 

  64. Gaikpa DS, Lieberherr B, Maurer HP, Longin CFH, Miedaner T. Comparison of rye, triticale, durum, and bread wheat genotypes for Fusarium head blight resistance and deoxynivalenol contamination. Plant Breed. 2020;139(2):251–62.

    Article  CAS  Google Scholar 

  65. González-Fernández R, Aloria K, Valero-Galván J, Redondo I, Arizmendi JM, Jorrín-Novo JV. Proteomic analysis of mycelium and secretome of different Botrytis cinerea wild-type strains. J Proteome. 2014;97:195–221.

    Article  CAS  Google Scholar 

  66. Igbaria A, Lev S, Rose MS, Lee BN, Hadar R, Degani O, et al. Distinct and combined roles of the MAP kinases of Cochliobolus heterostrophus in virulence and stress responses. Mol Plant Microbe Int. 2008;21:769–80.

    Article  CAS  Google Scholar 

  67. Mesterházy Á. Role of deoxynivalenol in aggressiveness of Fusarium graminearum and F. culmorum and in resistance to Fusarium head blight. In: Logrieco A, Bailey JA, Corazza L, Cooke BM, editors. Mycotoxins in plant disease. Dordrecht: Springer; 2002. p. 675–84.

    Chapter  Google Scholar 

  68. Hestbjerg H, Felding G, Elmholt S. Fusarium culmorum infection of barley seedlings: correlation between aggressiveness and deoxynivalenol content. J Phytopathol. 2002;150(6):308–12.

    Article  CAS  Google Scholar 

  69. Cumagun CJR, Rabenstein F, Miedaner T. Genetic variation and covariation for aggressiveness, deoxynivalenol production and fungal colonization among progeny of Gibberella zeae in wheat. Plant Pathol. 2004;53(4):446–53.

    Article  CAS  Google Scholar 

  70. Foroud NA, McCormick SP, MacMillan T, Badea A, Kendra DF, Ellis BE, et al. Greenhouse studies reveal increased aggressiveness of emergent Canadian Fusarium graminearum chemotypes in wheat. Plant Dis. 2012;96(9):1271–9.

    Article  CAS  PubMed  Google Scholar 

  71. Puri KD, Zhong S. The 3ADON population of Fusarium graminearum found in North Dakota is more aggressive and produces a higher level of DON than the prevalent 15ADON population in spring wheat. Phytopathology. 2010;100(10):1007–14.

    Article  CAS  PubMed  Google Scholar 

  72. Boenisch MJ, Schäfer W. Fusarium graminearum forms mycotoxin producing infection structures on wheat. BMC Plant Biol. 2011;11(1):110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Eudes F, Comeau A, Rioux S, Collin J. Impact of trichothecenes on Fusarium head blight [Fusarium graminearum] development in spring wheat (Triticum aestivum). Can J Plant Pathol. 2001;23(3):318–22.

    Article  Google Scholar 

  74. Bai GH, Plattner R, Desjardins A, Kolb F, McIntosh RA. Resistance to Fusarium head blight and deoxynivalenol accumulation in wheat. Plant Breed. 2001;120(1):1–6.

    Article  CAS  Google Scholar 

  75. McCormick SP, Harris LJ, Alexander NJ, Ouellet T, Saparno A, Allard S, et al. Tri1 in Fusarium graminearum encodes a P450 oxygenase. Appl Environ Microbiol. 2004;70(4):2044–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Reid LM, Mather DE, Hamilton RI, Bolton AT. Genotypic differences in the resistance of maize silk to Fusarium graminearum. Can J Plant Pathol. 1992;14(3):211–4.

    Article  Google Scholar 

  77. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013. Available from:

    Google Scholar 

  78. Bernal-Vasquez A-M, Utz H-F, Piepho H-P. Outlier detection methods for generalized lattices: a case study on the transition from ANOVA to REML. Theor Appl Genet. 2016;129(4):787–804.

    Article  PubMed  Google Scholar 

  79. Piepho H-P, Möhring J. Computing heritability and selection response from unbalanced plant breeding trials. Genetics. 2007;177(3):1881–8.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Stram DO, Lee JW. Variance components testing in the longitudinal mixed effects model. Biometrics. 1994;50(4):1171–7.

    Article  CAS  PubMed  Google Scholar 

  81. Gilmour AR, Gogel BJ, Cullis BR, Thompson R. ASReml user guide release 3.0. VSN International Ltd: Hemel Hempstead; 2009.

    Google Scholar 

  82. Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Duitama J, Quintero JC, Cruz DF, Quintero C, Hubmann G, Foulquié-Moreno MR, et al. An integrated framework for discovery and genotyping of genomic variants from high-throughput sequencing experiments. Nucleic Acids Res. 2014;42:e44.

    Article  Google Scholar 

  84. Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff. Fly. 2012;6(2):80–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Bushnell B. BBMap short read aligner. Berkeley: University of California; 2016. Available from:

    Google Scholar 

  87. Aulchenko YS, Ripke S, Isaacs A, van Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007;23(10):1294–6.

    Article  CAS  PubMed  Google Scholar 

  88. Würschum T, Langer SM, Longin CFL, Tucker MR, Leiser WL. A modern green revolution gene for reduced height in wheat. Plant J. 2017;92(5):892–903.

    Article  CAS  PubMed  Google Scholar 

  89. Utz HF, Melchinger AE, Schön CC. Bias and sampling error of the estimated proportion of genotypic variance explained by quantitative trait loci determined from experimental data in maize using cross validation and validation with independent samples. Genetics. 2000;154(4):1839–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Shin JH, Blay S, McNeney B, Graham J. LDheatmap: an R function for graphical display of pairwise linkage disequilibria between single nucleotide polymorphisms. J Stat Software. 2006;16(3):1–10

    Google Scholar 

  91. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  92. Weber T, Blin K, Duddela S, Krug D, Kim HU, Bruccoleri R, et al. AntiSMASH 3.0-a comprehensive resource for the genome mining of biosynthetic gene clusters. Nucleic Acids Res. 2015;43(W1):W237–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We highly acknowledge the work of our technicians on both experimental stations.


This research was funded by the German Federal Ministry of Education and Research (BMBF), Bonn, in the framework of Bioeconomy International (FusResist) to T.M. (No. 031B0011A). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The responsibility of the content of this publication rests with the authors. Open Access funding enabled and organized by Projekt DEAL.

Author information

Authors and Affiliations



TM designed the experiments. AV analyzed the sequence data and calculated all molecular data with support from TW and WL; VC and HEC performed the field phenotyping of all isolates in four environments and made the biometrical analyses of the phenotypic data. TM and AV wrote the paper that was edited and refined by TW, WL and NF. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Thomas Miedaner.

Ethics declarations

Ethics approval and consent to participate

Our experiments with cultivated plants are in consent with the relevant institutional and national legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Fusarium culmorum isolates used in this study.

Additional file 2: Table S2.

Number of SNPs identified and SNP density in each chromosome.

Additional file 3: Table S3.

Best linear unbiased estimates (BLUES) for aggressiveness in wheat (AGG.WHEAT), aggressiveness in rye (AGG.RYE), and deoxynivalenol content in wheat (DON-WHEAT) for all tested isolates across locations and years.

Additional file 4: Figure S1.

SNP density along the four chromosomes of F. culmorum. Histogram of variants per 100 Kb.

Additional file 5: Figure S2.

SNP count across all the samples of F. culmorum classified by the mutation type.

Additional file 6: Figure S3.

Frequency distritbution, phenotypic correlation coefficients and scatter plots for and between the evaluated traits for 92 isolates: Aggressiveness in wheat (AGG-WH) and in rye (AGG-RYE) and deoxynivalenol production in wheat (DON-WH). *** Significantly different from zero at 0.001 level of probability.

Additional file 7: Figure S4.

Scatterplot of the first two principal components (PC) of 92 isolates of F. culmorum originating from Germany, Syria, Russia and isolates from the international collection. The proportion of explained variance is shown in brackets at the corresponding axes.

Additional file 8: Figure S5.

Quantile-quantile plots (QQ-plots) for the aggressiveness in wheat (AGG-WH) and in rye (AGG-RYE) and deoxynivalenol production in wheat (DON-WH) showing the relation between expected and observed P values (−log10) for all SNPs after GWAS analysis adjusted for population stratification; Genomic inflation factor λ is given. The straight black diagonal lines represent the values under the null hypothesis of no association.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Miedaner, T., Vasquez, A., Castiblanco, V. et al. Genome-wide association study for deoxynivalenol production and aggressiveness in wheat and rye head blight by resequencing 92 isolates of Fusarium culmorum. BMC Genomics 22, 630 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: