Accurate, high-throughput phenotyping for quantitative traits is a limiting factor for progress in plant breeding. We developed an automated image analysis to measure quantitative resistance to septoria tritici blotch (STB), a globally important wheat disease, enabling identification of small chromosome intervals containing plausible candidate genes for STB resistance. 335 winter wheat cultivars were included in a replicated field experiment that experienced natural epidemic development by a highly diverse but fungicide-resistant pathogen population. More than 5.4 million automatically generated phenotypes were associated with 13,648 SNP markers to perform the GWAS. We identified 26 chromosome intervals explaining 1.9-10.6% of the variance associated with four independent resistance traits. Sixteen of the intervals overlapped with known STB resistance intervals, suggesting that our phenotyping approach can identify simultaneously (i.e., in a single experiment) many previously defined STB resistance intervals. Seventeen of the intervals were less than 5 Mbp in size and encoded only 173 genes, including many genes associated with disease resistance. Five intervals contained four or fewer genes, providing high priority targets for functional validation. Ten chromosome intervals were not previously associated with STB resistance, perhaps representing resistance to pathogen strains that had not been tested in earlier experiments. The SNP markers associated with these chromosome intervals can be used to recombine different forms of quantitative STB resistance that are likely to be more durable than pyramids of major resistance genes. Our experiment illustrates how high-throughput automated phenotyping can accelerate breeding for quantitative disease resistance.
Genome-wide association studies (GWAS) provide a powerful approach to identify genetic markers associated with important quantitative traits in crops (e.g., [1, 2]). The single nucleotide polymorphism (SNP) markers significantly associated with a trait in the GWAS can be directly used in breeding programs for marker-assisted selection or genomic selection and also as tools to enable map-based cloning of the corresponding genes underlying quantitative traits.
An abundant supply of SNP genetic markers is now available for most crops as a result of rapid advances in sequencing technologies. Because phenotyping technologies have not developed as quickly as genotyping technologies, the ability to generate accurate and reproducible phenotypes for quantitative traits is now the primary limitation to progress in breeding for favorable traits [3, 4], including resistance to pests and pathogens . Many research groups are working to develop automated/semiautomated and high-throughput phenotyping of important traits under field conditions, with some reports of success [5–7], but we remain far from the goal of using automated phenotyping to speed progress in plant breeding for useful traits.
Septoria tritici blotch (STB), caused by the fungus Zymoseptoria tritici, is currently the most damaging leaf disease on wheat in Europe  and is one of the most significant diseases on wheat around the world . Z. tritici has a mixed reproductive system, producing airborne ascospores through sexual reproduction that can be disseminated over distances of several kilometers and asexual conidia that are splash-dispersed over spatial scales of only 1-2 meters during the course of a growing season . Z. tritici populations are highly variable within fields as a result of its mixed reproductive system, large effective population sizes, and high levels of gene flow among populations [10, 11]. These properties provide a high evolutionary potential that leads to rapid development of virulence against resistant cultivars [10, 12] as well as resistance to fungicides [10, 13, 14]. STB in Europe is controlled mainly by applying fungicides costing over $1 billion per year , but many European Z. tritici populations have now evolved sufficiently high levels of resistance that fungicides are losing their efficacy [16, 17]. The European Union is planning to ban many fungicides in the near future (EU Regulation 1107/2009). These developments have stimulated new efforts to increase STB resistance through plant breeding.
Many studies have identified strain-specific STB resistance genes that could prove useful in breeding programs (summarized in ). These genes are typically identified and mapped by inoculating a single strain of Z. tritici onto a segregating population derived from a cross between two parents that differ in their resistance to that strain, with disease ratings often made on seedlings under greenhouse conditions. STB resistance in the field is mainly quantitative, but some examples of major gene resistance were identified (e.g., Stb6) that were recently shown to follow the gene-for-gene (GFG) pattern of inheritance [19, 20]. Unfortunately, major STB resistance genes such as Stb6 typically failed within 3-4 years of deployment as a result of pathogen evolution . A different breeding approach that is expected to slow pathogen evolution and be more durable is to make pyramids of quantitative resistance (QR) genes with additive effects [21, 22]. This approach requires the identification and deployment of QR that is effective across a broad cross section of the Z. tritici population as opposed to major gene resistance that works against only a small fraction of the strains found in natural field populations.
Identification of QR is difficult for most pathogens for many reasons including (1) measurement error associated with visual assessments of disease, (2) inherent differences in disease measurements conducted by different people, (3) differences in expression of QR in different environments, (4) the phenotypic effect associated with a QR locus (e.g., 5-10% reduction in disease) which can be too small to detect using disease ratings that operate in larger increments (e.g., 0-9 scales that differentiate disease increments of 10%), and (5) the occurrence of mixed infections by several pathogens under typical field conditions, with overlapping symptoms that often cannot be teased apart (e.g., STB symptoms look very similar to the symptoms associated with tan spot and stagonospora nodorum leaf blotch). These factors combine to create a low heritability for QR that slows progress in accumulating different sources of QR in breeding programs.
The recent development of automated image analysis for STB enabled rapid acquisition of large datasets, including millions of phenotype datapoints that were highly informative under both greenhouse and field conditions, and facilitated the cloning of genes encoding several avirulence effectors, including AvrStb6  and Avr3D1 [23, 24], as well as the Zmr1 gene affecting melanization of Z. tritici colonies and pycnidia [25, 26]. We took advantage of the high levels of fungicide resistance in Swiss populations of Z. tritici by using fungicide treatments to eliminate competing pathogens in a replicated field experiment . The fungicide treatments enabled a pure-pathogen read-out of quantitative resistance to STB caused by a genetically diverse, natural population of Z. tritici in an epidemic that developed under natural field conditions. Here, we use this extensive phenotype dataset in the GWAS to identify 26 chromosome intervals associated with quantitative STB resistance in a broad panel of 335 elite European winter wheat cultivars. Many of these intervals explained 6%-10% of the variance for the associated resistance trait. Several of the intervals contained a relatively small number of annotated genes, including genes known to be associated with disease resistance in wheat or other plants. There was a significant enrichment () for genes encoding putative receptor kinases and kinases within the 17 chromosome segments spanning less than 5 Mbp. Other candidate genes for STB resistance encoded NB-LRR proteins, F-box LRR proteins, sugar transporters, an ABC transporter, superoxide dismutase, and a TCP transcription factor, illustrating how automated image analysis can lead to identification of plausible candidate genes for quantitative disease resistance.
2. Materials and Methods
335 European winter wheat cultivars chosen from the GABI wheat panel  were grown in plots replicated twice as complete blocks at the Field Phenotyping Platform of the ETH research station in Lindau, Switzerland (47.449°N, 8.682°E, 520 masl) . The plots received full agrochemical inputs typically associated with intensive wheat cultivation in Europe, including mineral fertilizers, a stem shortener, and several pesticide applications. Among the pesticides, fungicides comprising five different active ingredients with three modes of action were applied at three time points over the growing season. Additional details associated with the field experiment are given in .
An unusual feature of this experiment is that all STB infection was natural, with the epidemic caused by a highly diverse Z. tritici population that immigrated into the experimental plots via windborne ascospores coming from nearby wheat fields that were treated with fungicides. This local Z. tritici population carried sufficient resistance to all fungicides applied in the experimental plots to enable an STB epidemic to develop despite the intensive fungicide treatments. But other wheat diseases common in this region, including leaf rust, stripe rust, stagonospora nodorum blotch, powdery mildew, and tan spot, were practically absent because the fungicides excluded these pathogens . As a result, we were able to obtain a pure-culture read-out of quantitative STB resistance across all 335 wheat cultivars without confounding effects from other diseases. The local weather during the 2015-2016 growing season was cooler and wetter than usual, providing a highly conducive environment for the development of an STB epidemic. At least six asexual reproduction cycles occurred during the most active period of wheat growth between March and July . Other components of STB epidemiology associated with this experiment were already reported .
All experimental plots were assessed for STB resistance at two time points, (20 May 2016, approximately growth stage GS 41 ) and (4 July 2016, approximately GS 75-85) using automated image analysis of 21,420 scanned leaves infected by Z. tritici [16, 30]. Nearly sixteen infected leaves collected from the same leaf layer in each plot were mounted on A4 paper and scanned at 1,200 dpi using flatbed scanners as described earlier . The scanned images were analyzed using an ImageJ macro script . Automatically generated outputs of the script included percentage of leaf area covered by lesions (PLACL), average pycnidia density within lesions (), and average pycnidia darkness (measured using the 256-point gray scale). To measure pycnidia sizes, we developed a Python (version 3.6.7, https://www.python.org/) program based on the determination of contours of constant brightness in the vicinity of each detected pycnidium with the help of the skimage package (version 0.13, https://scikit-image.org/). Each of these STB-associated phenotypes was analyzed separately in the GWAS. The grand means for each phenotype were calculated based on an average of 60 scanned leaves for each wheat cultivar, including both time points and both replicates for each plot (i.e., four measurements of each trait associated with STB resistance). The actual number of leaves analyzed for each cultivar ranged from 31 to 64, with some leaves omitted because of errors made during the collecting, mounting, or scanning processes. Only 35 cultivars had fewer than 55 leaves included in the analysis. The mean values of PLACL and were 1/x transformed to better fit a normal distribution, yielding a for the Shapiro-Wilk test after transformation. Figure 1 illustrates the steps involved in data acquisition and analysis.
The SNP markers used for the GWAS came from the Illumina 90K SNP array (iSelect, San Diego, USA, ). The majority of the markers on this array were not useful for our experiment because they were not polymorphic in the GABI panel. The remaining markers were positioned on the IWGSC wheat genome  using a BLASTn search with value < 10-30. The position with the lowest value was assigned as the marker position. In the case of ties where it was not possible to unequivocally assign a marker to one of the homeologous chromosomes, the markers were omitted. Additional filtering criteria to choose SNPs for the GWAS were a call rate of >95% per marker, >5% minor allele frequency, and identity by state , using the GenABEL software in the R statistical environment . After filtering, a total of 13,648 high-quality SNP markers were used for the GWAS. Haplotypes were identified using a sliding window of three consecutive SNPs with PLINK  and tested using linear regression models. GWAS Manhattan plots were constructed using R (version 3.5.1, ) with ggplot2 (version 3.1.0. ). Bonferroni thresholds were calculated using (0.05/13,648) yielding a LOD (-log10()) score of 5.44. The fraction of the phenotypic variance associated with the 26 chromosome intervals at or exceeding the Bonferroni threshold was calculated using linear regression models in R (lm function). The adjusted provided a measure of the proportion of the variance explained.
The coordinates of the 26 intervals exceeding the Bonferroni threshold were plotted onto the Chinese Spring reference genome as described earlier and used to compare the positions of the SNPs associated with the STB resistance traits identified in this analysis with the positions of STB resistance identified in earlier studies . The sequence data of the markers associated with STB resistance in earlier studies were retrieved from GrainGenes (https://wheat.pw.usda.gov/GG3/) and then searched using BLASTn against the IWGSC reference  assembly using Unité de Recherche Génomique Info (URGI, https://wheat-urgi.versailles.inra.fr/) for the corresponding chromosome. The position of the best hit was used as the genome position. We then calculated the probability of finding a chance overlap involving at least 16 of the 26 new associations with the previously identified chromosome segments for STB resistance. We treated each new association as an individual point because the new segments were very narrow and calculated the probabilities of overlaps based on a binomial distribution. For each new association, the probability of overlapping with an already known STB resistance interval (which we extended by 5 Mbp at both ends) was 41.2% (i.e., 5.89 Gbp of the 14.30 Gbp genome was included in the previously known intervals, many of which included a large fraction of a chromosome).
Candidate gene identification was based on the gene annotation of the IWGSC v1.0 reference sequence of the wheat landrace Chinese Spring . All high confidence genes in chromosome segments shorter than 5 Mb were identified.
The three fungicide treatments eliminated all competing fungal pathogens, enabling a mono disease readout of the relative degree of quantitative resistance to STB under the field conditions typically used for intensive wheat production in Europe. All STB infection was natural, with an epidemic resulting from at least six cycles of infection by a diverse Z. tritici population that included a high degree of gene and genotype diversity, with infections caused by millions of different Z. tritici strains despite the fungicide applications. As an indicator of the pathogen genetic diversity in these plots, genome sequences of 161 Z. tritici isolates obtained from 21 of the plots revealed 147 unique genome sequences, with all but two of the identified clones found within the same 1.5 m2 plot (Daniel Croll, personal communication). This high level of pathogen diversity was consistent with earlier findings from other naturally infected wheat fields around the world and was expected given that Z. tritici populations experience high levels of recombination [11, 35] that enable different fungicide resistance mutations to segregate and reassort into many different genetic backgrounds in natural populations.
The automated analysis pipeline generated phenotypes for 21,420 leaves, with an average of 30 leaves sampled per plot across the two time points. Nearly 37 m2 of leaf area was analyzed, with approx. 11 m2 scored as damaged by STB. The mean analyzed area for an individual leaf was 17 cm2. The percentage of leaf area covered by STB lesions (PLACL) ranged from 0 to 99 with a mean value of 32. The number of pycnidia found on a leaf ranged from 0 to 4,034, with a mean value of 127. The density of pycnidia within STB lesions () ranged from 0 to 256 per cm2 of lesion, with a mean value of 24. The automatically generated phenotypes included 21,420 measures each of PLACL and and 2.7 million measures each of pycnidia size and pycnidia melanization, yielding a total of >5.44 million automatically measured phenotypes that were not prone to human scoring error.
The quantitative measures of STB severity generated by automated image analysis followed the continuous distribution typically associated with quantitative traits [16, 30]. Earlier analyses of relationships among these traits  showed that resistance that minimizes host damage (PLACL) was largely independent of resistance that minimizes pathogen reproduction (). Hence, the GWAS was conducted independently for each trait. In addition to the traits PLACL and , we measured the average size of pycnidia formed within lesions, which reflects the average size and number of spores contained in each fruiting body  (i.e., pycnidia size is an independent indicator of pathogen reproduction), and the average gray value of pycnidia, which reflects the average amount of melanin accumulated in each fruiting body [25, 30]. Our earlier work indicated that pycnidia melanization is on average greater on wheat cultivars with more resistance to STB [24, 30].
Manhattan plots for PLACL, , pycnidia size, and pycnidia gray value revealed the SNPs with the highest associations for each STB resistance trait (Figure 2). A total of 109 SNPs were at or above the Bonferroni threshold across all traits based on the GWAS. Marker-trait associations were calculated using sliding windows including three consecutive SNPs. Among these, 52 haplotypes were at or above the Bonferroni threshold. Further evaluation of the 52 haplotypes revealed overlaps that were combined to produce a nonredundant set of 26 chromosome segments that explained from 1.9% to 10.6% of the overall phenotypic variance associated with each measured trait (Table 1).
For the PLACL trait that reflects the ability of a wheat cultivar to limit the degree of necrosis caused by an STB infection, 14 SNPs identified 4 different genomic positions distributed across chromosomes 5A, 5B, and 5D with LOD scores exceeding 5.5. Interval 4 on 5D had a LOD score of 9.2 and explained 10.3% of the total variance associated with PLACL (Table 1). For the trait that reflects the ability of a wheat cultivar to restrict Z. tritici reproduction, 51 SNPs identified 13 genomic positions located on chromosomes 2B, 4A, 5D, 6A, 6B, 6D, and 7B, with LOD scores ranging from 5.5 to 7.1. Interval 15 on 6B had the highest LOD score and explained 9.3% of the phenotypic variance associated with . For the pycnidia size trait, three SNPs located on 2B defined a chromosome interval that surpassed the Bonferroni threshold. Interval 18 explained 5.9% of the total variance associated with pycnidia size. For the pycnidia melanization trait, 36 SNPs defined 8 genomic positions located on chromosomes 1A, 2A, 3B, 4D, 5A, 5B, and 7B. Interval 23 on 4D showed the highest LOD score of 8.7 and explained 10.6% of the total variance associated with pycnidia melanization.
The positions of the 26 chromosome segments identified in this experiment were compared to the positions of mapped STB genes reported in earlier publications (summarized in Figure 1 of ). Figure 3 shows that 16 of the 26 chromosome segments identified in our analyses overlapped with or were very close to genomic regions identified in earlier publications. The probability of finding at least 16 associations within or close to the chromosome intervals previously associated with STB resistance was 2.9%. Ten of the chromosome segments were in chromosomal regions that were not previously associated with STB resistance. Among those, two were associated with PLACL (3, 4), five with (7, 8, 9, 16, and 17), and three with pycnidia melanization (20, 23, and 26) (Table 1).
Seventeen of the 26 chromosome segments were smaller than 5 Mb. For these, we identified putative candidate genes responsible for STB resistance based on the wheat reference genome sequence of Chinese Spring . In total, the 17 intervals spanned 26.2 Mb and contained 173 high confidence genes (Supplementary Table 1). There was a significant enrichment () for genes encoding putative receptor kinases and kinases within these 17 chromosome segments. Receptor kinase genes were recently shown to play major roles in disease resistance in cereals [19, 37, 38], including the Stb6 gene encoding resistance to STB . Five of the chromosome segments contained four or fewer genes, with three of these segments (19, 20, and 24) associated with pycnidia gray value and two segments (2, 3) associated with PLACL. The smallest chromosome segment (20) encompassed 28 kb on chromosome 1A and contained a single gene in Chinese Spring (TraesCS1A01G277000) encoding a putative solute carrier family 35 member. The 99.7 kb segment 24 on the long arm of chromosome 5A also had a single gene (TraesCS5A01G524800) encoding a putative 4-hydroxy-tetrahydrodipicolinate reductase, a protein involved in lysine biosynthesis. Intervals 2 (chromosome 5A) and 19 (chromosome 1A) had three candidate genes each, of which a putative kinase gene and a putative nucleotide binding site—leucine-rich repeat (NLR)—represent the most obvious candidates as these categories of genes are known to affect disease resistance . Interval 3 contained four candidate genes, all of which were associated with F-box proteins, a class of proteins often associated with plant defense responses . Other candidate genes known to be involved in disease resistance include sugar transporters (associated with PLACL in intervals 1 and 4), superoxide dismutase (associated with pycnidia size in interval 18), an ABC transporter (associated with PLACL in interval 1), and a TCP transcription factor (associated with in interval 15).
In a year that was highly conducive for the development of an STB epidemic, we combined a novel automated image analysis tool that could differentiate independent components of STB severity with the high level of fungicide resistance existing in a local Swiss population of Z. tritici to make a quantitative comparison of STB resistance across a broad cross section of elite European winter wheat cultivars. GWAS analyses that coupled these quantitative measures of STB resistance with 13,648 SNP markers enabled identification of 109 SNPs on 13 chromosomes that defined 26 chromosome segments highly associated with STB resistance. Because all STB infection in this experiment was natural, including millions of different pathogen genotypes originating from a recombining population, and the growing season was highly conducive to the development of an STB epidemic, we believe that the SNP markers defining the chromosome intervals associated with the highest levels of STB resistance could be especially useful in European breeding programs aimed at increasing overall levels of STB resistance to the Z. tritici populations found in Europe.
The 26 chromosome intervals associated with STB resistance ranged from 28 kbp to 60 Mbp in size and were distributed across 13 chromosomes, with individual intervals explaining 1.9% to 10.6% of the phenotypic variance for each trait. Some of the intervals were clustered in the same chromosomal region (e.g., intervals 1 and 2 associated with PLACL on chromosome 5A; intervals 13, 14, and 15 associated with pycnidia gray value on chromosome 6B), but most of the intervals were genetically distant from each other. Sixteen of the intervals were embedded within or located very close to chromosomal regions previously associated with STB resistance, but 10 intervals were in genomic regions that had not been associated with STB resistance. Though the known STB resistance intervals identified in earlier studies were quite large (covering in total ~40% of the entire wheat genome), the probability that at least 16 new intervals would overlap with these known intervals by chance was less than 3%. This finding increases our confidence that our new phenotyping method can identify simultaneously (i.e., in a single experiment) many previously defined STB resistance intervals that had been identified individually in studies involving only one pathogen strain. Particularly notable novel regions were the intervals 4, 8, and 9 located on 5D, a chromosome which had not previously been associated with STB resistance , though  found some weak associations with STB resistance on this chromosome. There was no overlap between chromosomal segments associated with host damage (PLACL) and pathogen reproduction (pycnidia density or pycnidia size), indicating that these resistance traits were under independent genetic control as hypothesized earlier . Intervals 4, 14, 15, and 23 had LOD scores at or exceeding 7. These are candidate regions for genes encoding broadly based field resistance to STB that may be especially useful against the genetically diverse Z. tritici populations in Europe.
The chromosome segments identified in our GWAS are much smaller than the intervals defined in earlier work as shown in Figure 2. For example, STB15 was previously mapped to a region that includes most of chromosome 6A (~590 Mbp) while we identified two separate chromosome regions (10 and 11) within the STB15 region that encompass only ~8.5 Mbp. Similarly, STB1 was mapped to a region that covered ~69.9 Mbp on chromosome 5B while interval 25 covers only ~13 Mbp within this region. The smaller intervals detected in our GWAS reflect the much higher marker density used in our experiment coupled with more accurate knowledge of marker positions coming from the new wheat genome assembly. Other contributors to the small intervals were the more accurate quantitative phenotypes yielding relatively large effect sizes and the haplotype-based GWAS approach that increased the statistical power compared to standard GWAS pipelines.
Seventeen of the 26 chromosome segments identified in the GWAS were less than 5 Mbp in size and contained between 1 and 28 candidate genes annotated in the Chinese Spring reference genome. The 173 genes located in these intervals were significantly enriched for receptor kinases and kinases, including clusters of 6 and 10 kinases found in intervals 18 and 22, respectively. We consider this enrichment to be notable because Stb6, the only cloned STB resistance gene, is a receptor kinase . Also notable was our finding that genes encoding receptor kinases are strongly upregulated during infection by all tested strains of Z. tritici . Hence, we hypothesize that some of the receptor kinase genes found in these intervals may be responsible for the STB resistance we observed. The interval 4 associated with PLACL explained 10.3% of the overall variance and provided the first report of STB resistance on chromosome 5D. This interval contained 12 genes, including three encoding proteins already shown to affect disease resistance, including an NLR, a S/T protein kinase, and a sugar transporter [39, 41]. We hypothesize that one or more of these genes are responsible for the STB resistance in this chromosome segment. Other interesting candidate genes found in the 17 intervals encode an ABC transporter, a TCP transcription factor, and superoxide dismutase. The Lr34 gene encoding quantitative resistance to leaf rust and other diseases in wheat was shown to be an ABC transporter . Superoxide dismutases are involved in synthesis of hydrogen peroxide, which was already shown to be involved in wheat’s defense response against STB . TCP transcription factors were shown recently to be important components of the signaling pathway involved in systemic acquired resistance . Segment 24, which explained 8% of the variance in pycnidia melanization and lies within the QTL9 region identified in earlier mapping studies, contained a single gene encoding a protein involved in lysine biosynthesis. Recent work on the wheat pathogen Cochliobolus sativus showed that lysine was essential for melanin biosynthesis  and lysine was recently shown to be essential for virulence in Z. tritici . We conclude from this analysis that many of the genes found in the intervals identified in the GWAS are plausible candidates to explain the observed phenotypes associated with STB resistance, but functional validation studies will be needed to confirm whether any of these genes actually play a role in resistance.
Earlier field trials also used association mapping to identify genetic markers associated with STB resistance [27, 47, 48]. In all of these trials, the experimental plots were inoculated with a small number of Z. tritici isolates that were sprayed when all wheat genotypes had fully extended flag leaves (i.e., ) a few weeks before scoring for STB resistance. As a result, the associations identified in those experiments are likely to be strain-specific and represent the outcome of a single cycle of infection based on a high dose of artificially applied blastospore inoculum. Similarly, most experiments that identified STB genes with major effects were based on greenhouse inoculations of seedlings by a single pathogen strain and used disease scores made at a single point in time, leading to identification of genes that encode seedling resistance to the strain used in the experiment. It is now clear that natural field infections of STB are caused by many millions of Z. tritici strains, with a different strain occurring on each infected leaf, on average, and with most leaves infected by more than one strain [10, 49]. The significant STB resistance associations identified in our experiment were based on a natural epidemic that included at least six cycles of pycnidiospore infection by a highly diverse population of the pathogen and included two time points during epidemic development. We believe that the STB resistance identified in our experiment is more likely to be broadly applicable under natural field conditions and hence more useful in breeding programs aiming for stable STB resistance.
An important and novel aspect of our experiment was the use of an automated image analysis pipeline for phenotyping that eliminated human scoring bias while generating millions of accurate phenotype datapoints. As is the case for many plant diseases , the traditional visual assessment of STB typically generates a single number on a 0-9 scale  that tries to integrate the totality of disease in a particular plot, often relative to other plots in the same field or trial. Visual assessments are fast, often requiring less than one minute per plot to produce a measurement, but are prone to variation caused by fatigue, changes in lighting over the course of a day, and differences in opinion among different scorers. The automated image analyses allowed us to simultaneously assess four quantitative phenotypes that could not be accurately measured by the eye. A traditional visual assessment would have generated a total of 4 STB measurements per cultivar to use in the GWAS. Our automated analyses generated an average of over 16,000 STB measurements per cultivar.
The leaf scanning method required significant manual labor (approx. 360 hours in total) to harvest infected leaves, mount them on paper, and scan them on flatbed scanners (Figure 1). A key element of our success was the ability to reproducibly count and measure properties (size and melanization) associated with millions of pycnidia, which measure only 0.5-1.5 mm across. Reliable identification of pycnidia required high-resolution images (1,200 dpi) acquired under uniform, high-intensity lighting, in this case provided by inexpensive flatbed scanners. Further automation of data collection done completely in the field will present new challenges. Flying drones at a height close to the canopy surface may be able to image pycnidia on the top leaf layer, but these drones would need to be coupled with strobe lighting systems and cameras that could produce images with sufficient magnification to reliably identify pycnidia over a significant depth of crop canopy. Ground-based mobile robots would make more sense for data acquisition in row crops such as maize, but wheat fields would be more difficult due to the high plant density. In the near future, we consider it likely that modest amounts of field data (10s-100s of leaves) could be acquired using mobile phone cameras outfitted with clip-on macro lenses or hand-held computers attached to portable single-leaf scanners. We can also envision portable high-throughput conveyor belt leaf scanning systems fit in the bed of a pickup truck that could scan more than 10,000 leaves in a day.
The detailed phenotype data generated by our automated analysis pipeline enabled us to separate different components of STB resistance, in particular allowing us to separate STB resistance that affects host damage (PLACL) from STB resistance that affects pathogen reproduction ( and pycnidia size). We believe that resistance affecting pathogen reproduction is likely to be more effective in the long run for several reasons, including the following: (1) Our earlier analyses  showed that measures of pathogen reproduction () early in the growing season were the best predictors of host damage (PLACL) late in the growing season, showing that resistance that reduces pathogen reproduction is likely to minimize yield losses caused by STB. (2) A decrease in pathogen reproduction diminishes the amount of inoculum available to cause new cycles of infection, which will lower the transmission rate (i.e., decrease the basic reproductive number, ) during each infection cycle and result in less overall infection by the end of the epidemic. (3) A decrease in pathogen inoculum will lead to a decrease in the pathogen population size, which will decrease the overall genetic diversity and provide fewer opportunities for favorable mutations (e.g., for fungicide resistance or gain of virulence) to emerge. This should lower the overall evolutionary potential of the pathogen population .
Recombining the SNP markers associated with the STB resistance intervals identified in this experiment may accelerate breeding efforts aimed at increasing quantitative resistance to STB in European wheat. We showed that resistance affecting leaf damage (PLACL) is genetically distinct from resistance affecting pathogen reproduction (). We consider it likely that these different resistance phenotypes reflect different underlying mechanisms of STB resistance. We hypothesize that PLACL reflects the additive actions of toxin sensitivity genes that interact with host-specific toxins produced by the pathogen, as shown for Parastagonospora nodorum on wheat [53, 54], while pycnidia density reflects the additive actions of quantitative resistance genes that recognize pathogen effectors (e.g., ). Under this scenario, breeders should aim to recombine these two forms of resistance into the same genetic background, bringing together different forms of resistance that may be more durable when deployed together than when either mechanism is deployed in isolation. We anticipate that functional analyses of the most compelling candidate genes identified in this experiment will enable us to identify new genes underlying the different STB resistance traits. Our experiment illustrates how high-throughput automated phenotyping can accelerate breeding for quantitative disease resistance.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
B.A.M. conceived and managed the project. B.A.M., A.M., and B.S. provided funding for the project. A.H., S.F., L.M., D.P., and P.K. performed the experiments. S.Y., A.M., P.K., S.G.K., M.A., and K.Y. analyzed the data. B.A.M., S.Y., and S.G.K. wrote the manuscript. Steven Yates and Alexey Mikaberidze contributed equally to this work.
H. Zellweger managed the wheat trial. Marion Roeder from IPK Gatersleben provided seeds and marker information for the GABI wheat panel. STB research in BAM’s lab was supported by the Swiss National Science Foundation (grants 155955, 134755, 104145, and 56874) and the ETH Zurich Research Commission (grants 12-03 and 15-02). AM and PK were supported by the Swiss National Science Foundation through Ambizione grant PZ00P3_161453.
Table S1: candidate genes in the chromosome intervals defined by the GWAS. (Supplementary Materials)
- S. G. Milner, M. Jost, S. Taketa et al., “Genebank genomics highlights the diversity of a global barley collection,” Nature Genetics, vol. 51, no. 2, pp. 319–326, 2019.
- K. Yano, E. Yamamoto, K. Aya et al., “Genome-wide association study using whole-genome sequencing rapidly identifies new genes influencing agronomic traits in rice,” Nature Genetics, vol. 48, no. 8, pp. 927–934, 2016.
- R. T. Furbank and M. Tester, “Phenomics–technologies to relieve the phenotyping bottleneck,” Trends in Plant Science, vol. 16, no. 12, pp. 635–644, 2011.
- J. L. Araus and J. E. Cairns, “Field high-throughput phenotyping: the new crop breeding frontier,” Trends in Plant Science, vol. 19, no. 1, pp. 52–61, 2014.
- S. Joalland, C. Screpanti, H. Varella et al., “Aerial and ground based sensing of tolerance to beet cyst nematode in sugar beet,” Remote Sensing, vol. 10, no. 5, p. 787, 2018.
- J. Crain, S. Mondal, J. Rutkoski, R. P. Singh, and J. Poland, “Combining high-throughput phenotyping and genomic information to increase prediction and selection accuracy in wheat breeding,” The Plant Genome, vol. 11, no. 1, article 170043, 2018.
- R. Wedeking, A. K. Mahlein, U. Steiner, E. C. Oerke, H. E. Goldbach, and M. A. Wimmer, “Osmotic adjustment of young sugar beets (Beta vulgaris) under progressive drought stress and subsequent rewatering assessed by metabolite analysis and infrared thermography,” Functional Plant Biology, vol. 44, no. 1, pp. 119–133, 2017.
- L. N. Jørgensen, M. S. Hovmøller, J. G. Hansen et al., “IPM strategies and their dilemmas including an introduction to www.eurowheat.org,” Journal of Integrative Agriculture, vol. 13, no. 2, pp. 265–281, 2014.
- S. Savary, L. Willocquet, S. J. Pethybridge, P. Esker, N. McRoberts, and A. Nelson, “The global burden of pathogens and pests on major food crops,” Nature Ecology & Evolution, vol. 3, no. 3, pp. 430–439, 2019.
- B. A. McDonald and C. C. Mundt, “How knowledge of pathogen population biology informs management of Septoria tritici blotch,” Phytopathology, vol. 106, no. 9, pp. 948–955, 2016.
- J. Zhan, R. E. Pettway, and B. A. McDonald, “The global genetic structure of the wheat pathogen Mycosphaerella graminicola is characterized by high nuclear diversity, low mitochondrial diversity, regular recombination, and gene flow,” Fungal Genetics and Biology, vol. 38, no. 3, pp. 286–297, 2003.
- C. Cowger, M. E. Hoffer, and C. C. Mundt, “Specific adaptation by Mycosphaerella graminicola to a resistant wheat cultivar,” Plant Pathology, vol. 49, no. 4, pp. 445–451, 2000.
- L. K. Estep, S. F. F. Torriani, M. Zala et al., “Emergence and early evolution of fungicide resistance in North American populations of Zymoseptoria tritici,” Plant Pathology, vol. 64, no. 4, pp. 961–971, 2015.
- L. E. Hayes, K. E. Sackett, N. P. Anderson, M. D. Flowers, and C. C. Mundt, “Evidence of selection for fungicide resistance in Zymoseptoria tritici populations on wheat in western Oregon,” Plant Disease, vol. 100, no. 2, pp. 483–489, 2016.
- S. F. F. Torriani, J. P. E. Melichar, C. Mills, N. Pain, H. Sierotzki, and M. Courbot, “Zymoseptoria tritici: a major threat to wheat production, integrated approaches to control,” Fungal Genetics and Biology, vol. 79, pp. 8–12, 2015.
- P. Karisto, A. Hund, K. Yu et al., “Ranking quantitative resistance to septoria tritici blotch in elite wheat cultivars using automated image analysis,” Phytopathology, vol. 108, no. 5, pp. 568–581, 2018.
- H. J. Cools and B. A. Fraaije, “Update on mechanisms of azole resistance in Mycosphaerella graminicola and implications for future control,” Pest Management Science, vol. 69, no. 2, pp. 150–155, 2013.
- J. K. M. Brown, L. Chartrain, P. Lasserre-Zuber, and C. Saintenac, “Genetics of resistance to Zymoseptoria tritici and applications to wheat breeding,” Fungal Genetics and Biology, vol. 79, pp. 33–41, 2015.
- C. Saintenac, W. S. Lee, F. Cambon et al., “Wheat receptor-kinase-like protein Stb6 controls gene-for-gene resistance to fungal pathogen Zymoseptoria tritici,” Nature Genetics, vol. 50, no. 3, pp. 368–374, 2018.
- Z. Zhong, T. C. Marcel, F. E. Hartmann et al., “A small secreted protein in Zymoseptoria tritici is responsible for avirulence on wheat cultivars carrying the Stb6 resistance gene,” New Phytologist, vol. 214, no. 2, pp. 619–631, 2017.
- R. P. Singh, P. K. Singh, J. Rutkoski et al., “Disease impact on wheat yield potential and prospects of genetic control,” Annual Review of Phytopathology, vol. 54, no. 1, pp. 303–322, 2016.
- C. C. Mundt, “Pyramiding for resistance durability: theory and practice,” Phytopathology, vol. 108, no. 7, pp. 792–802, 2018.
- L. Meile, D. Croll, P. C. Brunner et al., “A fungal avirulence factor encoded in a highly plastic genomic region triggers partial resistance to septoria tritici blotch,” New Phytologist, vol. 219, no. 3, pp. 1048–1061, 2018.
- E. l. Stewart, D. Croll, M. H. Lendenmann et al., “Quantitative trait locus mapping reveals complex genetic architecture of quantitative virulence in the wheat pathogen Zymoseptoria tritici,” Molecular Plant Pathology, vol. 19, no. 1, pp. 201–216, 2018.
- P. Krishnan, L. Meile, C. Plissonneau et al., “Transposable element insertions shape gene regulation and melanin production in a fungal pathogen of wheat,” BMC Biology, vol. 16, no. 1, p. 78, 2018.
- M. H. Lendenmann, D. Croll, E. L. Stewart, and B. A. McDonald, “Quantitative trait locus mapping of melanization in the plant pathogenic fungus Zymoseptoria tritici,” G3: Genes, Genomes, Genetics, vol. 4, no. 12, pp. 2519–2533, 2014.
- S. Kollers, B. Rodemann, J. Ling et al., “Genetic architecture of resistance to Septoria tritici blotch (Mycosphaerella graminicola) in European winter wheat,” Molecular Breeding, vol. 32, no. 2, pp. 411–423, 2013.
- N. Kirchgessner, F. Liebisch, K. Yu et al., “The ETH field phenotyping platform FIP: a cable-suspended multi-sensor system,” Functional Plant Biology, vol. 44, no. 1, pp. 154–168, 2017.
- J. C. Zadoks, T. T. Chang, and C. F. Konzak, “A decimal code for the growth stages of cereals,” Weed Research, vol. 14, no. 6, pp. 415–421, 1974.
- E. L. Stewart, C. H. Hagerty, A. Mikaberidze, C. C. Mundt, Z. Zhong, and B. A. McDonald, “An improved method for measuring quantitative resistance to the wheat pathogen Zymoseptoria tritici using high-throughput automated image analysis,” Phytopathology, vol. 106, no. 7, pp. 782–788, 2016.
- S. Wang, D. Wong, K. Forrest et al., “Characterization of polyploid wheat genomic diversity using a high-density 90 000 single nucleotide polymorphism array,” Plant Biotechnology Journal, vol. 12, no. 6, pp. 787–796, 2014.
- International Wheat Genome Sequencing Consortium (IWGSC), R. Appels, K. Eversole et al., “Shifting the limits in wheat research and breeding using a fully annotated reference genome,” Science, vol. 361, no. 6403, article eaar 7191, 2018.
- Y. S. Aulchenko, S. Ripke, A. Isaacs, and C. M. Van Duijn, “GenABEL: an R library for genome-wide association analysis,” Bioinformatics, vol. 23, no. 10, pp. 1294–1296, 2007.
- S. Purcell, B. Neale, K. Todd-Brown et al., “PLINK: a tool set for whole-genome association and population-based linkage analyses,” The American Journal of Human Genetics, vol. 81, no. 3, pp. 559–575, 2007.
- R Core Team, R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria, 2018, https://www.R-project.org/.
- H. Wickham, ggplot2: Elegant Graphics for Data Analysis, Springer, 2016.
- B. Keller and S. G. Krattinger, “A new player in race-specific resistance,” Nature Plants, vol. 4, no. 4, pp. 197-198, 2018.
- X. Ma, B. Keller, B. A. McDonald, J. Palma-Guerrero, and T. Wicker, “Comparative transcriptomics reveals how wheat responds to infection by Zymoseptoria tritici,” Molecular Plant-Microbe Interactions, vol. 31, no. 4, pp. 420–431, 2018.
- P. N. Dodds and J. P. Rathjen, “Plant immunity: towards an integrated view of plant-pathogen interactions,” Nature Reviews Genetics, vol. 11, no. 8, pp. 539–548, 2010.
- H. A. van den Burg, D. I. Tsitsigiannis, O. Rowland et al., “The F-box protein ACRE189/ACIF1 regulates cell death and defense responses activated during pathogen recognition in tobacco and tomato,” Plant Cell, vol. 20, no. 3, pp. 697–719, 2008.
- J. W. Moore, S. Herrera-Foessel, C. Lan et al., “A recently evolved hexose transporter variant confers resistance to multiple pathogens in wheat,” Nature Genetics, vol. 47, no. 12, pp. 1494–1498, 2015.
- S. G. Krattinger, E. S. Lagudah, W. Spielmeyer et al., “A putative ABC transporter confers durable resistance to multiple fungal pathogens in wheat,” Science, vol. 323, no. 5919, pp. 1360–1363, 2009.
- N. P. Shetty, R. Mehrabi, H. Lütken et al., “Role of hydrogen peroxide during the interaction between the hemibiotrophic fungal pathogen Septoria tritici and wheat,” New Phytologist, vol. 174, no. 3, pp. 637–647, 2007.
- M. Li, H. Chen, J. Chen et al., “TCP transcription factors interact with NPR1 and contribute redundantly to systemic acquired resistance,” Frontiers in Plant Science, vol. 9, 2018.
- Y. Leng and S. Zhong, “Sfp-type 4-phosphopantetheinyl transferase is required for lysine synthesis, tolerance to oxidative stress and virulence in the plant pathogenic fungus Cochliobolus sativus,” Molecular Plant Pathology, vol. 13, no. 4, pp. 375–387, 2012.
- M. C. Derbyshire, A. M. Gohari, R. Mehrabi et al., “Phosphopantetheinyl transferase (Ppt)-mediated biosynthesis of lysine, but not siderophores or DHN melanin, is required for virulence of Zymoseptoria tritici on wheat,” Scientific Reports, vol. 8, no. 1, article 17069, p. 17069, 2018.
- T. Miedaner, Y. Zhao, M. Gowda et al., “Genetic architecture of resistance to septoria tritici blotch in European wheat,” BMC Genomics, vol. 14, no. 1, p. 858, 2013.
- Q. H. Muqaddasi, Y. Zhao, B. Rodemann, J. Plieske, M. W. Ganal, and M. S. Röder, “Genome-wide association mapping and prediction of adult stage Septoria tritici blotch infection in European winter wheat via high-density marker arrays,” The Plant Genome, vol. 12, no. 1, article 180029, 2019.
- C. C. Linde, J. Zhan, and B. A. McDonald, “Population structure of Mycosphaerella graminicola: from lesions to continents,” Phytopathology, vol. 92, no. 9, pp. 946–955, 2002.
- E. E. Saari and J. M. Prescott, “A scale for appraising the foliar intensity of wheat diseases,” Plant Disease Reporter, vol. 59, pp. 377–380, 1975.
- Z. Eyal, A. L. Scharen, J. M. Prescott, and M. van Ginkel, The Septoria Diseases of Wheat: Concepts and Methods of Disease Management, CIMMYT, 1987.
- B. A. McDonald and C. Linde, “Pathogen population genetics, evolutionary potential, and durable resistance,” Annual Review of Phytopathology, vol. 40, no. 1, pp. 349–379, 2002.
- T. L. Friesen, Z. Zhang, P. S. Solomon, R. P. Oliver, and J. D. Faris, “Characterization of the interaction of a novel Stagonospora nodorum host-selective toxin with a wheat susceptibility gene,” Plant Physiology, vol. 146, no. 2, pp. 682–693, 2008.
- R. P. Oliver, T. L. Friesen, J. D. Faris, and P. S. Solomon, “Stagonospora nodorum: from pathology to genomics and host resistance,” Annual Review of Phytopathology, vol. 50, no. 1, pp. 23–43, 2012.