Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Genome-wide association analysis of the strength of the MAMP-elicited defense response and resistance to target leaf spot in sorghum

## Abstract

Plants have the capacity to respond to conserved molecular features known as microbe-associated molecular patterns (MAMPs). The goal of this work was to assess variation in the MAMP response in sorghum, to map loci associated with this variation, and to investigate possible connections with variation in quantitative disease resistance. Using an assay that measures the production of reactive oxygen species, we assessed variation in the MAMP response in a sorghum association mapping population known as the sorghum conversion population (SCP). We identified consistent variation for the response to chitin and flg22—an epitope of flagellin. We identified two SNP loci associated with variation in the flg22 response and one with the chitin response. We also assessed resistance to Target Leaf Spot (TLS) disease caused by the necrotrophic fungus Bipolaris cookei in the SCP. We identified one strong association on chromosome 5 near a previously characterized disease resistance gene. A moderately significant correlation was observed between stronger flg22 response and lower TLS resistance. Possible reasons for this are discussed.

## Introduction

Sorghum (Sorghum bicolor) has a diploid genome of ~ 730 Mb with 10 chromosomes1,2. It is a widely-grown cereal crop used as feed or silage for animal consumption, for bio-fuel production and as gluten-free grain for human consumption and is better adapted to grow under high heat and drought conditions than other agriculturally important crops like corn and wheat. These agronomically-important traits make the species an attractive crop for the mass production of grains and bio-fuel under challenging growing conditions.

The warm and humid conditions under which much sorghum is grown support the growth of a wide variety of foliar fungi. Among many diseases, target leaf spot (TLS) caused by the necrotrophic fungus Bipolaris cookei, is one of the most economically-important fungal diseases of sorghum in the southeastern US, causing major yield losses3. TLS causes distinctive oval or elliptical reddish-purple spots that eventually coalesce during disease progression.

The genetic basis of TLS resistance in sorghum has been the subject of several studies. A major recessive resistance gene ds1 on chromosome 5 was identified as a loss of function allele of a gene encoding a leucine-rich repeat receptor kinase4. Other work identified QTL for three different fungal diseases; target leaf spot, zonate leaf spot and drechslera leaf blight, co-localized on chromosome 65. A third study identified a TLS resistance QTL on chromosome 3, as well as the previously-reported chromosome 6 QTL6. A recent study identified novel QTL on chromosomes 3, 4 and 9 as well as a strong QTL on chromosome 5 near the ds1 locus7.

Plants possess cell-surface receptors known as pattern recognition receptors (PRRs) that mediate recognition of highly conserved structural molecules associated with microbes known as microbe-associated molecular patterns (MAMPs). The two best-studied MAMPs are bacterial flagellin, especially its flg22 epitope, and chitin, a component of the fungal cell wall8,9. MAMP recognition elicits a basal response at the infection site known as MAMP-triggered immunity (MTI) which often includes phenomena such as callose deposition, changes in membrane ion flux, changes in phytohormone concentrations, induction or repression of plant defense-related genes, and production of reactive oxygen species (ROS) and nitric oxide (NO)10. In some cases, a pathogen adapted to a particular host can overcome MTI by producing so-called effector proteins, which are usually introduced into the cytoplasm and may suppress MTI. Effectors are sometimes recognized by cytoplasmic receptors known as R proteins, eliciting a strong response known as effector-triggered immunity (ETI) which is quantitatively stronger though, qualitatively somewhat similar to MTI11,12.

Non-host resistance can be defined as: “Resistance shown by an entire plant species against all known genetic variants (or isolates) of a specific parasite or pathogen”13,14. It has been hypothesized that MTI is a significant cause of non-host resistance, as most non-adapted pathogens cannot subvert the MTI-based defenses of their non-host plants15.

Host resistance can be subdivided into qualitative and quantitative resistance. Qualitative resistance is typically based on the action of a single, large-effect gene, while quantitative resistance is mediated by large numbers of small-effect genes16,17. There is some evidence that variation in the strength of MTI may underlie some part of quantitative resistance. Genetic variation in the strength of the MTI response has been documented in a number of plant species including brassicas18,19,20,21, maize22, soybean23, and tomato24 and in several of these cases, QTL controlling the variation were identified. In particular, the fact that several genes resembling PRRs confer quantitative resistance in various plant species16 and that the strength of flg22 perception is negatively correlated with susceptibility to Pseudomonas syringae in Arabidopsis20 suggest that there may be some connection between variation in the MTI response and quantitative disease resistance. However, the relationship between these traits is not well understood, especially in crop plants. The objectives of this study were to characterize the genetic control of the MAMP response and TLS resistance in a diverse panel of sorghum germplasm and to determine if there was evidence of shared genetic control of these traits. Specifically, we wanted to determine whether a stronger MAMP response was indicative of stronger quantitative resistance.

## Materials and methods

### Plant and pathogen materials

A sorghum association mapping population known as the sorghum conversion population (SCP) was provided by Dr. Pat Brown at the University of Illinois (now at UC Davis). It has been described previously25 and is a collection of diverse lines converted to photoperiod-insensitivity and smaller stature to facilitate the growth and development of the plants in US environments26. 510 lines from this population were used in this study although due to bad germination and other quality control issues, not all the lines were used in the analysis of all three traits. Ultimately data from 345 lines were used for the analysis of the chitin response, 472 lines for the flg22 response, and 456 for TLS resistance. B. cookei strain LSLP18 was obtained from Dr. Burt Bluhm at the University of Arkansas.

### MAMP response measurement

Two different MAMPs were used in this study flg22, (Genscript catalog# RP19986), and chitin (Sigma catalog # C9752). Sorghum plants were grown in inserts laid on flats filled with soil (33% Sunshine Redi-Earth Pro Growing Mix) in the greenhouse. Plants were watered the day before sample collection to avoid extra leaf moisture on the day of collection.

The lines were randomized and, for logistical reasons, were planted in batches of 60 lines. For each line, three ‘pots’ were planted with two seeds per line. Subsequent batches were planted as soon as the previous batch had been processed until the entire population had been assessed. Two experimental runs were conducted for both MAMPs with genotypes re-randomized in each of the two runs.

ROS assays were carried out as previously described27. Briefly, for each line, six seeds were planted in 3 different pots. From the resulting seedlings, three were selected based on uniformity. Seedlings that looked unusual or were significantly taller or shorter than the majority were not used. Four leaf discs of 3 mm diameter were excised from the broadest part of the 4th leaf of three different 15-day old sorghum plants. One disc per leaf from two plants and two discs from one plant, with the second disc becoming the water control (see below). The discs were individually floated on 50 µl H20 in a black 96-well plate, sealed with an aluminum seal to avoid exposure to light, and kept at room temperature overnight. The next morning a reaction solution was made using 2 mg/ml chemiluminescent probe L-012 (Wako, catalog # 120-04891), 2 mg/ml horseradish peroxidase (Type VI-A, Sigma-Aldrich, catalog # P6782), and 100 mg/ml Chitin or 2 μM of Flg22. 50 µl of this reaction solution was added to three of the four wells. The fourth well was a mock control, to which the reaction solution excluding the MAMP was added. Four blank wells containing only water were also included in each plate.

After adding the reaction solution, the luminescence was measured using SynergyTM 2 multi-detection microplate reader (BioTek) every 2 min for 1 hr. The plate reader takes luminescence measurements every 2 min during this 1 h. The sum of all 31 readings was calculated to give the value for each well. The estimated value for the MAMP response for each genotype was calculated as (average luminescence value of the three experimental wells—the mock well value) -minus the average blank well value. The blank well values were consistently close to zero.

Leaf discs of Nicotiana benthamiana, one high responsive sorghum line (SC0003), and one low responsive sorghum line (PI 6069) were also included as controls in each 96-well plate for quality control purposes.

### B. cookei inoculum preparation and inoculation

B. cookei inoculum was prepared as described previously28. Briefly, sorghum grains were soaked in water for three days, rinsed, scooped into 1L conical flasks and autoclaved for an hour at 15psi and 121 °C. The grains were then inoculated with about 5 ml of macerated mycelia from a fresh culture of B. cookei LSLP18 isolates and left for 2 weeks at room temperature, shaking the flasks every 3 days. After 2 weeks, the fungus infested sorghum grains were air-dried and then stored at 4 °C until field inoculation. The same inoculum was used for the entire trial and made fresh every year. For inoculation, 6–10 infested grains were placed into the whorl of 4–5 week old sorghum plants. The spores produced from these fungi initiated infection in the young sorghum plants within a week.

### Seed preparation

Before planting in the field sorghum seed was treated with a fungicide, insecticide, and safener mixture containing ~ 1% Spirato 480 FS fungicide, 4% Sebring 480 FS fungicide, 3% Sorpro 940 ES seed safener. Then the seeds were air-dried for 3 days which provided a thin coating of this mix around the seeds. The safener allowed the use of the herbicide Dual Magnum as a pre-emergence treatment.

### Evaluation of Target Leaf Spot resistance

The SCP was planted at the Central Crops Research Station in Clayton, NC on June 14–15 2017 and June 20, 2018 in a randomized complete block design with two experimental replications in each case. Experiments were planted in 1.8 m single rows with a 0.9 m row width using 10 seeds per plot. Two border rows were planted around the periphery of each experiment to prevent edge effects. The experiments were inoculated on July 20, 2017 and July 20, 2018 at which point the sorghum plants were at growth stage 3. Ratings were taken on a one to nine scales (Fig. S2), where plants showing no signs of disease were scored as a nine and completely dead plants were scored as one (Fig. S2). Two ratings were taken in 2017 and four readings in 2018 starting two weeks after inoculation each year. sAUDPC (standardized area under disease progression curve) was calculated as described previously29,30.

### Statistical analyses

All statistical analysis of phenotypic data was performed using SAS V9.4 software. The LSmeans of two replicates for each year were calculated and these were used in turn to calculate the overall LS means. Analysis of variance (ANOVA) and least square (LS)means were calculated using the Proc Mix and Proc GLM procedure in SAS respectively. Correlations were calculated using the CORR procedure of SAS31.

### Phenotypic data transformation for association analysis

The phenotypic distribution was right skewed for flg22 and chitin-elicited ROS response traits. From a simple ANOVA we determined that higher predicted score values of these phenotypic traits were moderately associated with higher residuals. Therefore, natural logarithm transformation and root square transformation were performed using the raw scores of chitin and flg22, respectively. After transformation, the phenotypic distribution of each trait was less skewed and the relationship between residual and predicted values was improved. Transformed data were used in the downstream association analysis. Each trial was analyzed separately with SAS mixed model procedure in SAS software version 9.4 (SAS institute. 2019). For the chitin and flg22 response, a best linear unbiased estimator (BLUE) was obtained to estimate each line mean phenotypic value by a mixed model considering inbred lines as fixed effects and replications as random effects. Similarly, for TLS, the data of both years were combined by using a mixed linear model across years considering years and replication within years as random effects and inbred lines as fixed. All the original phenotypic data used for analysis is provided in File S1.

### Genotypic data

All genotypic data used for this study are available upon request from the corresponding author. Genotypic data for the SCP were obtained from Dr. Tiffany Jamann and Dr. Patrick Brown (University of Illinois). The original array consisted of ~ 1.12 million SNPs derived from whole-genome sequencing. We used the genotypes of each set of plants with phenotypic data described above. Each data set was first filtered to exclude SNPs with less than 5% minor allele frequency (MAF) and more than 10% heterozygosity. Linkage disequilibrium-based pruning of genotypic data was performed in software Plink v1.932. After pruning, a set of ~ 58 K SNPs were used to compute the kinship matrix in Tassel 533. Pruned data was also used to perform principal component analysis (PCA) with JMP Genomics 9 (SAS, Institute. 2019). Based on PCA, approximately 20% of the variability was accounted for by the first three principal components. To control for population structure in the association analysis, we excluded 37 inbred lines that explain more than 7% of the variability (Fig. S1). After removing outlier inbred lines, a second filter (< 5% MAF and > 10% heterozygote sites) was performed in each data set. The final sub-set of genotypes used in Genome wide Association Studies (GWAS) contained ~ 755 K SNPs for TLS and flg22 and ~ 750 K SNPs for Chitin. The genotypic datasets used for the analysis of the TLS, flg22 and chitin analyses are available from the corresponding authors.

### Association analysis

Genome-wide association analysis based on a mixed linear model (MLM) was performed in Tassel 533. The MLM model used was y = Xβ + Zu + e where y is the vector of phenotypes (BLUEs), β is a vector of fixed effects, including the SNP marker, tested, u is a vector of random additive effects (inbred lines), X and Z represent matrices, and e is a vector of random residuals. The variance of random line effects was modeled as Var(u) = K $${\sigma }_{a}^{2}$$, where K is the n × n matrix of pairwise kinship coefficient and $${\sigma }_{a}^{2}$$ is the estimated additive genetic variance. A threshold to declare significance of 1/m where m is the number of markers tested was used34,35.

### Candidate gene selection

Genes within 100 Kb of the highly significant markers were considered candidate genes. Identification and annotation of the candidate genes were performed using the maize BTx623 reference genome v3 available on the Ensembl Plants browser. Functional annotation of the candidate gene was based on EnsemblPlants and Gramene annotation.

## Results and discussion

### Evaluation and mapping of TLS resistance

The SCP was assessed in the field for TLS resistance in 2017 and 2018, in randomized complete blocks with two complete replications per year. We observed substantial variation in TLS resistance in the SCP (Fig. 1A). The two replicates in each year were significantly correlated (0.52 and 0.68 in 2017 and 2018 respectively, p < 0.0001) and the LSmean scores were significantly correlated between years (0.45, p < 0.0001) (Table 1). ANOVA analysis indicated that the genotype effect was significant (Table 2).

Association analysis using the LSMeans of the 492 lines that were scored identified a single highly significant association on chromosome 5 (Fig. 1B). Table 3 shows the parameters associated with this locus and details predicted genes located 100 Kb either side in the Btx623 genome. One of these genes is ds1, a leucine-rich repeat serine/threonine protein kinase gene that was previously reported as a major TLS resistance gene4. It seems very likely that this gene underlies the major QTL identified in the SCP.

### Evaluation and mapping of the MAMP response

To assess variation in basal immune response, we measured ROS production in response to flg22 and chitin treatment in the SCP in two full replications. We observed significant variation in response to both MAMPs (Fig. 2). Significant correlations were observed between replicates in both cases (0.5 and 0.38 for flg22 and chitin respectively, p < 0.0001, Table 1). ANOVA indicated that genotype effects were highly significant for both traits (Table 1).

The LSmeans of the flg22 and chitin responses were not significantly correlated, though replicate 2 of flg22 response was somewhat correlated with rep1 and rep2 of the chitin response (0.20 and 0.17 respectively, p < 0.01 and < 0.05). In previous work, we observed significant correlations in flg22 and chitin responses measured using the same ROS plate assay as well as shared QTL, in a maize recombinant inbred line mapping population22. The lack of correlation here is therefore somewhat surprising. It is not clear whether this reflects fundamental differences between the maize and sorghum MAMP responses. Vetter et al21 found a negligible correlation in plant growth responses between the bacterial MAMPs EF-Tu and flagellin in Arabidopsis, we are not aware of other published work comparing variation in the responses to two different MAMPs.

The phenotypic data were transformed as described and used for association analysis. Q-Q plots did not indicate an excess of false positives (Fig. S3). Two associations with the flg22 response were detected on chromosome 4. The significance threshold was calculated using a Bonferroni multiple comparison test correction which is based on the number of markers used. Since we used a relatively high number of markers (more than 750,000) this threshold was consequently relatively conservative. It should be noted that the associations with flg22 were below the threshold for significance we used but we are nevertheless reporting them as they are the highest associations detected and, given the conservative significance threshold used, are nevertheless likely to reflect real associations. One significantly associated locus was detected on chromosome 5 for the chitin response (Fig. 3). Table 3 shows the parameters associated with these associated loci and details predicted genes. While it is premature to assign causation, it is interesting to note that several of the candidate genes associated with flg22- and chitin-induced responses have homology to genes involved in the defense response or disease resistance in other systems. For instance, the durable wheat rust resistance gene LR34 is an ABC transporter36 while genes involved in the auxin response37 and the ubiquitin-mediated protein-degradation38 pathway have been implicated in disease resistance in other systems.

### Comparison of TLS resistance and MAMP response data

To understand whether variation in the response elicited by flg22 and chitin is connected to TLS resistance, we looked for correlation between MAMP response and TLS disease scores. Despite chitin being an integral component of the fungal cell wall and TLS being a fungal disease, we did not observe a significant correlation between the traits. We did observe a small but moderately-significant negative correlation between the flg22 response and TLS scores (− 0.13*, p value < 0.05), indicating that a higher flg22 response was somewhat associated with higher susceptibility. This was unexpected both because flg22 is a bacterial MAMP and TLS is a fungal disease and because we were expecting an association between increased MAMP response and increased resistance. Instead, we observed an opposite relationship, albeit quite weak. Two possible explanations occur to us. Since the correlation is relatively low, this may not be a meaningful correlation. Alternatively, several necrotrophic pathogens of a similar type to Bipolaris cookei have effectors that both induce ETI and facilitate pathogenesis. It appears that in these cases elicitation of HR allows the pathogen to grow on the resulting dead host cells39. It is possible that this correlation is due to a similar subversion of the plant defense machinery.

A recent companion study measuring the flg22 response and TLS resistance in two sorghum recombinant inbred line (RIL) populations did not identify correlations between the 2 traits or any colocalizing QTL7. In the current study, we used the SCP which provided higher genetic and phenotypic diversity than had been available from the two RIL populations but, overall this study also did not produce evidence to support our original hypothesis that a stronger MAMP response is predictive of stronger QDR. However, there are a number of caveats that make it impossible to draw general conclusions.

Perhaps the major caveat is that quantitation of the MTI response is complex. It depends on what MAMP is used and how the response is measured. Low correlations between responses to different MAMPs have been reported previously19,21,24, although, as mentioned above, Zhang et al22 observed a significant correlation between flg22 and chitin responses in maize. Moreover, the MAMP response can be quantified in a number of different ways, including measurement of ROS or NO production, MAP kinase phosphorylation, mRNA accumulation levels, lignin and cell wall-bound phenols, callose deposition, seedling growth inhibition and MAMP-induced pathogen resistance18,19,40. Relative line rankings vary significantly depending on the assay used19,22. Our preliminary data also suggests that the MAMP response varies with the age of the plant and the individual leaf on the plant. Essentially, quantification of the MAMP response is complex and inferences may vary significantly depending on how the response is elicited and how measured7.

The other major caveat is of course resistance to only one disease was assessed. It may be that resistance to certain diseases, perhaps those that are less well adapted to the host and cannot completely suppress basal resistance mechanisms, may be more associated with the MAMP response. As more diseases are assessed on the SCP we may be able to re-evaluate our hypothesis in the light of multiple comparisons.

## References

1. 1.

Price, H. J. et al. Genome evolution in the genus Sorghum (Poaceae). Ann. Bot. 95, 219–227 (2005).

2. 2.

Paterson, A. H. et al. The Sorghum bicolor genome and the diversification of grasses. Nature 457, 551–556 (2009).

3. 3.

Zaccaron, A. Z. & Bluhm, B. H. The genome sequence of Bipolaris cookei reveals mechanisms of pathogenesis underlying target leaf spot of sorghum. Sci. Rep. 7, 17217 (2017).

4. 4.

Kawahigashi, H. et al. Positional cloning of ds1, the target leaf spot resistance gene against Bipolaris sorghicola in sorghum. Theor. Appl. Genet. 123, 131–142 (2011).

5. 5.

Mohan, S. M. et al. Co-localization of quantitative trait loci for foliar disease resistance in sorghum. Plant Breed. 128, 532–535 (2009).

6. 6.

Murali Mohan, S. et al. Identification of quantitative trait loci associated with resistance to foliar diseases in sorghum [Sorghum bicolor (L.) Moench]. Euphytica 176, 199–211 (2010).

7. 7.

Kimball, J. et al. Identification of QTL for Target Leaf Spot resistance in Sorghum bicolor and investigation of relationships between disease resistance and variation in the MAMP response. Sci. Rep. 9, 18285 (2019).

8. 8.

Zipfel, C. et al. Bacterial disease resistance in Arabidopsis through flagellin perception. Nature 428, 764–767 (2004).

9. 9.

Kaku, H. et al. Plant cells recognize chitin fragments for defense signaling through a plasma membrane receptor. Proc. Natl. Acad Sci. U. S. A. 103, 11086–11091 (2006).

10. 10.

Newman, M.-A., Sundelin, T., Nielsen, J. T. & Erbs, G. MAMP (microbe-associated molecular pattern) triggered immunity in plants. Front. Plant Sci. 4, 139 (2013).

11. 11.

Bent, A. F. & Mackey, D. Elicitors, effectors, and r genes: the new paradigm and a lifetime supply of questions. Annu. Rev. Phytopathol. 45, 399–436 (2007).

12. 12.

Bittel, P. & Robatzek, S. Microbe-associated molecular patterns (MAMPs) probe plant immunity. Curr. Opin. Plant Biol. 10, 335–341 (2007).

13. 13.

Heath, M. C. Nonhost resistance and nonspecific plant defenses. Curr. Opin. Plant Biol. 3, 315–319 (2000).

14. 14.

Heath, M. C. A generalized concept of host-parasite specificity. Phytopathology 71, 1121–1123 (1981).

15. 15.

Lipka, U., Fuchs, R. & Lipka, V. Arabidopsis non-host resistance to powdery mildews. Curr. Opin. Plant Biol. 11, 404–411 (2008).

16. 16.

Nelson, R., Wiesner-Hanks, T., Wisser, R. & Balint-Kurti, P. Navigating complexity to breed disease-resistant crops. Nat. Rev. Genet. 19, 21–33 (2018).

17. 17.

Poland, J. A., Balint-Kurti, P. J., Wisser, R. J., Pratt, R. C. & Nelson, R. J. Shades of gray: the world of quantitative disease resistance. Trends Plant Sci. 14, 21–29 (2009).

18. 18.

Lloyd, S. R., Schoonbeek, H.-J., Trick, M., Zipfel, C. & Ridout, C. J. Methods to study PAMP-triggered immunity in brassica species. Mol. Plant Microbe Interact. 27, 286–295 (2014).

19. 19.

Lloyd, S. R., Ridout, C. J. & Schoonbeek, H.-j. Methods to Quantify PAMP-Triggered Oxidative Burst, MAP Kinase Phosphorylation, Gene Expression, and Lignification in Brassicas. In Plant Pattern Recognition Receptors: Methods and Protocols 325–335 (2017).

20. 20.

Vetter, M. M. et al. Flagellin perception varies quantitatively in Arabidopsis thaliana and its relatives. Mol. Biol. Evol. 29, 1655–1667 (2012).

21. 21.

Vetter, M., Karasov, T. L. & Bergelson, J. Differentiation between MAMP triggered defenses in Arabidopsis thaliana. PLoS Genet. 12, e1006068 (2016).

22. 22.

Zhang, X., Valdés-López, O., Arellano, C., Stacey, G. & Balint-Kurti, P. Genetic dissection of the maize (Zea mays L.) MAMP response. Theor. Appl. Genet. 130, 1155–1168 (2017).

23. 23.

Valdes-Lopez, O. et al. Identification of quantitative trait loci controlling gene expression during the innate immunity response of soybean. Plant Physiol. 157, 1975–1986 (2011).

24. 24.

Hind, S. R. et al. Tomato receptor FLAGELLIN-SENSING 3 binds flgII-28 and activates the plant immune system. Nat. Plants 2, 16128 (2016).

25. 25.

Thurber, C. S., Ma, J. M., Higgins, R. H. & Brown, P. J. Retrospective genomic analysis of sorghum adaptation to temperate-zone grain production. Genome Biol. 14, R68 (2013).

26. 26.

Stephens, J. C., Miller, F. R. & Rosenow, D. T. Conversion of alien sorghums to early combine genotypes1. Crop Sci. 7, 396–396 (1967).

27. 27.

Samira, R. et al. Quantifying MAMP-induced production of reactive oxygen species in sorghum and maize. Bio-protocol 9, e3304 (2019).

28. 28.

Sermons, S. M. & Balint-Kurti, P. J. Large scale field inoculation and scoring of maize southern leaf blight and other maize foliar fungal diseases. Bio-protocol 8, e2745 (2018).

29. 29.

Campbell, C. L. & Madden, L. V. Introduction to Plant Disease Epidemiology 192–194 (Wiley, New York, 1990).

30. 30.

Shaner, G. & Finney, P. E. The effect of nitrogen fertilizer on expression of slow mildewing resistance in Knox wheat. Phytopathology 67, 1051–1056 (1977).

31. 31.

SAS Institute, I. SAS 9.2 Help and Documentation,. (SAS, Cary , NC, 2000–2004).

32. 32.

Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575 (2007).

33. 33.

Bradbury, P. et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635 (2007).

34. 34.

Xu, Y. et al. Genome-wide association mapping of starch pasting properties in maize using single-locus and multi-locus models. Front. Plant Sci. 9, 1311 (2018).

35. 35.

Zhang, Y.-M., Jia, Z. & Dunwell, J. M. Editorial: the applications of new multi-locus GWAS methodologies in the genetic dissection of complex traits. Front. Plant Sci. 10 (2019).

36. 36.

Krattinger, S. G. et al. A putative ABC transporter confers durable resistance to multiple fungal pathogens in wheat. Science 323, 1360–1363 (2009).

37. 37.

Ding, X. et al. Activation of the indole-3-acetic acid-amido synthetase GH3-8 suppresses expansin expression and promotes salicylate- and jasmonate-independent basal immunity in rice. Plant Cell 20, 228–240 (2008).

38. 38.

Goritschnig, S., Zhang, Y. & Li, X. The ubiquitin pathway is required for innate immunity in Arabidopsis. Plant J. 49, 540–551 (2007).

39. 39.

Lorang, J. M. Necrotrophic exploitation and subversion of plant defense: a lifestyle or just a phase, and implications in breeding resistance. Phytopathology 109, 332–346 (2019).

40. 40.

Zeidler, D. et al. Innate immunity in Arabidopsis thaliana: Lipopolysaccharides activate nitric oxide synthase (NOS) and induce defense genes. Proc. Natl. Acad. Sci. U.S.A. 101, 15811–15816 (2004).

## Acknowledgements

We thank Dr. Steve Kresovich for advice on seed handling and agronomic practice, and Dr. Burt Bluhm for providing isolates of B. cookei. We thank Cathy Herring and the staff at Central Crops Research Station for their work facilitating the field trials. Dr. Shannon Sermons and Greg Marshall assisted with several aspects of the field research. This work was funded by the DOE Plant Feedstock Genomics for Bioenergy program grants # DE-SC0014116 and DE-SC0019189.

## Author information

Authors

### Contributions

R.S., J.A.K., G.S. and P.B.K. planned the experiments, RS, JAK and PJBK performed the experiments and wrote the manuscript, R.S., J.A.K., L.S.F.L. and J.H. conducted the analyses. L.S.F.L. and R.S. prepared the figures. T.M.J. and P.J.B. provided seeds and unpublished genotypic datasets. All authors reviewed and edited the manuscript.

### Corresponding authors

Correspondence to Jennifer A. Kimball or Peter J. Balint-Kurti.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

### Publisher's note

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

## Rights and permissions

Reprints and Permissions

Samira, R., Kimball, J.A., Samayoa, L.F. et al. Genome-wide association analysis of the strength of the MAMP-elicited defense response and resistance to target leaf spot in sorghum. Sci Rep 10, 20817 (2020). https://doi.org/10.1038/s41598-020-77684-w

• Accepted:

• Published: