# Genome-Wide Association Study for Resistance to Rhynchosporium in a Diverse Collection of Spring Barley Germplasm

Jean-Noël Thauvin, Joanne Russell, Dominique Vequaud, Mark Looseley, Micha Bayer, Pierre-Marie Le Roux, Pierre Pin, Robbie Waugh, Anna Avrova

► To cite this version:

Jean-Noël Thauvin, Joanne Russell, Dominique Vequaud, Mark Looseley, Micha Bayer, et al. Genome-Wide Association Study for Resistance to Rhynchosporium in a Diverse Collection of Spring Barley Germplasm. Agronomy, 2022, 12 (4), 10.3390/agronomy12040782. [https://doi.org/10.3390/agronomy12040782](https://doi.org/10.3390/agronomy12040782)

HAL Id: hal-03712623

https://hal.inrae.fr/hal-03712623

Submitted on 4 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Genome-Wide Association Study for Resistance to Rhynchosporium in a Diverse Collection of Spring Barley Germplasm

Jean-Noël Thauvin <fnref n="1" />,<fnref n="2" />,*,†, Joanne Russell <fnref n="1" />, Dominique Vequaud <fnref n="3" />, Mark Looseley <fnref n="1" />,‡ ID, Micha Bayer <fnref n="1" /> ID,
Pierre-Marie Le Roux <fnref n="3" />, Pierre Pin <fnref n="3" />, Robbie Waugh <fnref n="1" />,<fnref n="2" /> and Anna Avrova <fnref n="1" />,* ID

<figure><img src="image_2.png" /></figure>

### check for updates

**Citation:** Thauvin, J.-N.; Russell, J.; Vequaud, D.; Looseley, M.; Bayer, M.; Le Roux, P.-M.; Pin, P.; Waugh, R.; Avrova, A. Genome-Wide Association Study for Resistance to Rhynchosporium in a Diverse Collection of Spring Barley Germplasm. *Agronomy* **2022**, *12*, 782. https://doi.org/10.3390/agronomy12040782

Academic Editors: Jerzy Henryk Czembor, Alan H. Schulman and Guoping Zhang

Received: 17 February 2022

Accepted: 21 March 2022

Published: 24 March 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Rhynchosporium is one of the main biotic stresses on barley production worldwide. A set of 312 spring barley accessions was tested in four different locations over 3 years, to identify novel genetic resistances to rhynchosporium and to explore the allelic diversity for resistance genes present in this global germplasm collection. High-density genotypes from exome capture and RNA-seq were used to conduct high-resolution association mapping. Seven quantitative trait loci (QTL) were detected, including one in the *Rrs2* region, amongst five containing known resistances. Relatively short physical intervals harbouring these resistances were proposed, providing a platform for the identification of underlying genes and tightly linked genetic markers for use in marker assisted selection. Genes encoding kinases were present in four of the QTL, in addition to *Rrs1* and *Rrs18*, two loci known to contribute to rhynchosporium resistance. The frequencies and distributions of these novel and known QTL were superimposed on the regional origin of the landrace genotypes comprising the genome-wide association studies (GWAS) panel, highlighting the value of genetic resources as a source of diverse genetically controlled resistance to rhynchosporium. The detected QTL along with their linked genetic markers, could be exploited either directly for breeding purposes or for candidate gene identification in future studies.

**Keywords:** barley; landraces; rhynchosporium; resistance; QTL; GWAS; SNP; genetic markers

## 1. Introduction

Barley is the fourth most important cereal crop worldwide, by cultivated area, with 52 million hectares grown in 2020 [1]. One of its major diseases is rhynchosporium or scald, caused by the hemibiotrophic fungal pathogen *Rhynchosporium commune* [2]. Rhynchosporium develops in cold and wet conditions, with primary inoculum originating from infected seeds or straw and further spore dispersal through water splashes [3]. The pathogen is highly genetically variable and 58% of the global genetic diversity has been reported to be represented within individual fields [4]. In the United Kingdom, rhynchosporium is the most damaging disease of barley and is mainly managed by fungicides, with two sprays recommended for spring barley [5]. The number of effective fungicides available is decreasing, due to the withdrawal of products or loss of efficacy, meaning that resistant barley cultivars will be key for sustainable barley production [6,7].

<footnote n="1">¹ The James Hutton Institute, Invergowrie, Dundee DD2 5DA, UK; joanne.russell@hutton.ac.uk (J.R.); mark.looseley@xelect.co.uk (M.L.); micha.bayer@hutton.ac.uk (M.B.); robbie.waugh@hutton.ac.uk (R.W.)</footnote>

<footnote n="2">² Division of Plant Sciences, University of Dundee (at JHI), Invergowrie, Dundee DD2 5WE, UK</footnote>

<footnote n="3">³ Secobra Recherches SAS, Centre de Bois-Henry, 78580 Maule, France; dominique.vequaud@secobra.com (D.V.); pierremarie.leroux@secobra.com (P.-M.L.R.); pierre.pin@secobra.com (P.P.)</footnote>

<footnote>* Correspondence: jeannoel.thauvin@gmail.com (J.-N.T.); anna.avrova@hutton.ac.uk (A.A.);</footnote>

<footnote>† Current address: INRAE, UMR 1095 Genetics, Diversity and Ecophysiology of Cereals (GDEC), 5 Chemin de Beaulieu, 63000 Clermont-Ferrand, France.</footnote>

<footnote>‡ Current address: Xelect Ltd., Horizon House, Abbey Walk, St Andrews KY16 9LB, UK.</footnote>
A better understanding of genetic resistance to rhynchosporium is needed to improve
barley cultivars and the sustainability of the crop. Indeed, very little is known about
the actual genes responsible for resistance to rhynchosporium. The major resistance loci
are collectively known as *Rrs* (for Reaction to *Rhynchosporium secalis*, the former name of
the species) [8], but so far candidate genes have been published for only three *Rrs* genes.
Receptor-like kinases were suggested as causative genes for *Rrs1*<sub>Rh4</sub> [9] and *Rrs18* [10].
In contrast, *Rrs2* has been fine mapped to a 0.08 cm interval containing a cluster of Pectin
Esterase Inhibitor genes, but none of them individually trigger the *Rrs2* resistance re-
sponse [11,12]. A receptor kinase-like protein has also been suggested as a candidate gene
for a quantitative trait locus (QTL) QSc.VR4 [13]. In total, 11 major resistance genes (*Rrs*
genes) have been described, but numerous other QTL have been identified and a full review
was conducted in 2020 [14]. The authors listed 148 QTL, many of which are mapped to
genetically large or overlapping intervals. So far, truly diagnostic markers based on single
nucleotide polymorphisms (SNPs) have only been published for *Rrs1*<sub>Rh4</sub> [9].

European spring barley germplasm is relatively susceptible to rhynchosporium, with
a prior association mapping study identifying only one effective major resistance gene
(*Rrs1*<sub>Rh4</sub>) segregating in elite spring barley cultivars [15]. As the pathogen is highly diverse
and can adapt very quickly to defeat a single resistance gene [16], the identification and
characterization of novel sources of resistance is necessary to provide long-term sustainable
protection. Therefore, the use of diverse germplasm, including landraces, will be a valuable
source of untapped genetic variability for rhynchosporium resistance for use in barley
breeding [17].

Genome-wide association studies (GWAS) are a powerful tool to map regions of the
genome associated with traits in germplasm collections. Indeed, GWAS can access con-
siderable diversity and numerous historic recombination events, which results in reduced
linkage disequilibrium (LD), providing greater mapping resolution [18]. Recently this ap-
proach has been used in a number of studies to identify resistance loci in barley [15,19,20].
In European spring barley germplasm, the *Rrs1*<sub>Rh4</sub> locus was shown to be the main con-
tributor to resistance, with a further 15 resistance QTL showing minor effects [9]. In a
collection of Ethiopian, North American, and the International Center for Agricultural
Research in the Dry Areas (ICARDA) barley lines, multiple resistance QTL were found with
17 marker-trait associations (MTA) in 16 genomic regions [19]. The value of Scottish bere
barley landraces to mine for rhynchosporium resistance genes was also pointed out and
eight genetic regions associated with resistance were identified [20]. These studies used
three, five, and two field trials, respectively, for their GWAS analysis. A further eight QTL
originating from wild barley accessions were identified by GWAS in the Nested Association
Mapping population HEB-25 [21]. GWAS was also used to map nine rhynchosporium
resistance regions, including three new regions, in four multi-parent advanced generation
inter-cross (MAGIC) populations tested in northern Europe [22].

The aim of this study was to determine whether using a highly diverse legacy collection of barley, from across the geographical range of the species, combined with over 100,000 physically located SNPs, variable field trials, and state of the art genome-wide analyses, could lead to detection of novel resistance QTL and narrow down the physical intervals for the previously published resistance QTL. A further important aim was to identify resistant accessions that could provide novel sources of pre-breeding material for the breeding community. SNP markers associated with resistance QTL detected here will help to improve resistance against this dangerous barley pathogen.

**2. Materials and Methods**

2.1. *Plant Material*

A diverse collection of 312 accessions (Supplementary Table S1) chosen to represent
the maximum diversity among the worldwide spring barley germplasm was assembled
from three different sources. The majority (229 accessions) came from the Wheat and Barley
Legacy (WHEALBI) collection [23]. The other accessions were landraces from Syria and
Jordan (60 accessions), previously identified as highly resistant to rhynchosporium [9], and lastly a set of landraces from the Vavilov Institute collection (15 accessions) [24]. Eight cultivars were added as resistant and susceptible controls for known resistances or susceptible cultivars: Armelle (Rrs1<sub>Brier</sub>), Atlas, Atlas46 (Rrs1<sub>Rh3</sub>, Rrs2), Belgravia (Rrs1<sub>Rh4</sub>), Digger (Rrs2), Doyen (relatively susceptible control), Optic (susceptible control) and Retriever (Rrs1<sub>Rh4</sub>).

## 2.2. Field Trials

Four sites were used for the field trials between the 2018 and 2020 growing seasons. Table 1 provides details of the locations and the number of field trials carried out at each site. All the sites were autumn-sown (end of November), except for Maule in 2018, where a spring sown (beginning of March) field trial was also conducted. Each field trial had a single replicate block design with 1 m² mini-plots, except for Dundee where the field had a double replicate randomised block design. Rhynchosporium spread from the naturally present inoculum from residual barley crop debris built up over several years of field trials [9,25] and, therefore, the pathogen population was not controlled and could be different in each field [3,4]. In order to increase the quantity of primary inoculum, the previous crop was always barley, and the Dundee site had overhead irrigation and spreader rows consisting of a mixture of susceptible winter barley cultivars.

<table><caption>**Table 1.** Field trials of barley collection for resistance to rhynchosporium.</caption><thead><tr><th>Site</th><th>Soil Type</th><th>GPS Coordinates</th><th>Number of Field Trials</th><th>Years and Sowing Season (A = Autumn, S = Spring)</th></tr></thead><tbody><tr><td>Dundee (Scotland)</td><td>Sandy loam</td><td>56.45 N, -3.07 E</td><td>2</td><td>2019 (A), 2020 (A)</td></tr><tr><td>Dunmow (England)</td><td>Loam</td><td>51.88 N, 0.41 E</td><td>1</td><td>2019 (A)</td></tr><tr><td>Maule (France)</td><td>Clay loam</td><td>48.92 N, 1.82 E</td><td>3</td><td>2018 (A), 2018 (S), 2019 (A)</td></tr><tr><td>Pithiviers (France)</td><td>Sandy loam</td><td>48.12 N, 2.14 E</td><td>2</td><td>2018 (A), 2019 (A)</td></tr></tbody></table>

## 2.3. Field Phenotyping and Data Preparation

The rhynchosporium phenotyping was carried out toward the end of the growing season, in May and June (post anthesis). Rhynchosporium levels were scored on a 1 to 9 scale provided in Table 2, where 1 represents no visible symptoms in the entire plot and 9, where all the leaves were dead due to rhynchosporium [25,26]. Phenotypes provided in Supplementary Table S1 represent the disease level across the canopy. Most of the field trials were scored once during the season, although Dundee 2020, Maule 2018 and Pithiviers 2018 were scored on two separate occasions to take the progress of the infection into account. In these cases, the average of the two scores was used for further analysis. For the fields with two replicates of each accession, the average of the two scores was used for the analysis.

<table><caption>**Table 2.** Scoring scale for rhynchosporium.</caption><thead><tr><th>Score</th><th>Description</th></tr></thead><tbody><tr><td>1</td><td>No infection observed</td></tr><tr><td>2</td><td>1% infection on lower leaves</td></tr><tr><td>3</td><td>5% infection on lower leaves</td></tr><tr><td>4</td><td>25% infection on lower leaves</td></tr><tr><td>5</td><td>50% infection on lower leaves</td></tr><tr><td>6</td><td>Leaves appear 1/2 infected 1/2 green</td></tr><tr><td>7</td><td>Leaves appear more infected than green</td></tr><tr><td>8</td><td>Very little green leaf tissue left</td></tr><tr><td>9</td><td>Leaves dead no green leaf left</td></tr></tbody></table>
In addition to rhynchosporium, heading date and plant height were scored in subsets of five and six field trials respectively. Indeed, plant height, and earliness are known to be involved in disease escape from rhynchosporium [3] and information on these traits can help the identification of rhynchosporium resistance genes per se. The heading date was scored as the day when half of the ears in a plot were becoming visible. Heading date is a relatively reliable and simple method to score earliness, it corresponds to Zadoks stage Z55 [27]. Plant height was measured in cm from the soil to the bottom of the spike at the beginning of ripening (Zadoks stage Z91). The best linear unbiased estimators (BLUE) for plant height and heading date were calculated for all the accessions, with the lmer function from the lme4 R package [28].

A multiple regression was then applied to explain the rhynchosporium phenotypes for the 8 field trials depending on the BLUE-height and BLUE-heading date, following the model below:

$$ \text{Rhynchosporium}_{ij} = a \times \text{BLUE-height}_{j} + b \times \text{BLUE-heading date}_{j} + \varepsilon_{ij} $$

where i = field trial, j = genotype, a and b are coefficients, and εij are the residuals for the field trial i and the genotype j.

The residuals of this model were then used as the rhynchosporium phenotypes to analyse in the association studies of the eight field trials. The use of a regression to prepare a dataset to analyse is a common practice when two quantitative variables are linked and is commonly used in the case of grain protein content in wheat, with the residuals of the regression protein content depending on yield called the grain protein deviation [29].

## 2.4. Genotyping

The accessions from the WHEALBI collection have legacy exome capture sequence data available. The exome capture data were generated using Illumina short read sequencing of predicted coding areas of the genome. A description of this and the population structure and LD pattern is available from a previous study [23]. The WHEALBI exome capture reads were mapped onto the Morex 2017 assembly [30] using BWA-MEM [31]. A SNP calling was performed with the GATK [32] following GATK's best practices pipeline [33]. The raw variant calls were filtered with custom Java version 6 USA. code using the following criteria:

1. More than 8× coverage for at least 50% of the samples (to ensure robust SNP and genotype calls);

2. More than 95% of samples represented at SNP locus (for maximum sample representation);

3. Minor allele frequency of 1% (to exclude SNPs based on very rare alleles);

4. Quality score: less than 1/1000 chance of the SNP having been called by error (for SNP robustness).

All the other accessions were RNA-sequenced (RNA-seq) except for Atlas46, SLB-22-014 and SLB-34-076 accessions for which whole genome shotgun (WGS) sequencing was carried out (>30× coverage). The RNA and genomic DNA extraction was done using the second leaf from 3–4 weeks old seedling using the Qiagen RNeasy Plant mini kit and DNeasy kit, following the manufacturer's instructions. DNA and RNA integrity was tested by gel electrophoresis. DNA and RNA yield was measured using PicoGreen (Turner BioSystems Inc., Sunnyvale, CA, USA). RNA samples were DNase I treated using the Ambion Turbo DNA-free kit (Fisher Scientific UK Ltd., Loughborough, UK), following the manufacturer's protocol. All the sequencing was done using Illumina paired-end sequencing. The RNA-seq gave on average 63 million reads (from 1.7 to 319 million reads depending on accession). The reads were mapped onto the Morex 2017 genome assembly [30] using STAR [34]. Both the RNA-seq and WGS data were used together for the SNP calling with the GATK. SNPs with less than 50% missing data in this dataset were selected to be merged with the exome capture dataset.

Hereafter, the data analysis was carried out using R 3.6.1 'Action of the toes' [35].
The diagnostic SNP chr3H_490244130 identified by Looseley et al. [9] was used to determine the list of accessions containing *Rrs1*<sub>Rh4</sub>.

The merging of the different SNP datasets resulted in 115,172 SNPs (Supplementary Table S2). An imputation of the missing data was necessary to avoid the clustering of the accessions depending on the genotyping method used. The imputation was done with the knnimpute function from the R package scrime [36].

## 2.5. GWAS Using Single-Locus and Multi-Locus Models

The association mapping was carried out for each phenotyping dataset using three different methods: one single-locus GWAS method (EMMAX) and two multi-locus methods, where the SNPs were not tested alone, but in combinations to take into account epistasis. First, the GWAS was performed using an efficient mixed model association expedited, EMMAX, implemented by the rrBLUP R package [37]. In this model, the population structure is estimated in a first step and used as fixed effect in the mixed linear model [38]. The two multi-locus GWAS methods used were FASTmrMLM and BLINK. The multi-locus GWAS methods were designed to have higher statistical power and lower false positive results for multigenic traits [39]. FASTmrMLM was implemented with the mrMLM R package [40]. This method works in two stages, first the whole genome is scanned to detect the most significant markers, and in a second step, a few SNPs are selected and tested together in the multi locus model. The second multi-locus model used in this study was BLINK: Bayesian-information and Linkage Disequilibrium Iteratively Nested Keyway [41]. BLINK was implemented using the BLINK-R package [42]. Single-locus and multi-locus GWAS models are considered to be complementary depending on genetic architecture and population structure within the germplasm collection [43,44].

For all the GWAS methods, the SNPs detected were considered as having a significant impact on rhynchosporium scores when $-\log_{10}P > 5$. This significance threshold is relatively high and was set to avoid false positives in a context where GWAS for the same trait was run with three different methods potentially generating a high number of MTAs. The population structure was assessed with a cluster analysis based on Principal Component Analysis (PCA), and sub-populations were determined using the R 'kmeans' function. The LD decay was then calculated around each SNP detected in a GWAS, taking into account the population structure, using the LDcorSV package [45], as done by Bustos-Korts et al. [23]. The LD was averaged using a sliding window of 10 markers to estimate the local LD decay. All the SNPs detected with the different methods have been gathered in a single dataset. The SNPs with overlapping intervals for a LD decay of $R^2 = 0.2$ were considered as detecting the same QTL, as recommended by Alqudah et al. [18]. QTL were considered as consistent and kept in the final results if they were detected within at least two field trials.

# 3. Results and Discussion

## 3.1. Genotyping

The cluster analysis of the genotyping datasets composed of 115 172 SNPs per accession identified six groups within the landrace collection (Figure 1). Cluster 1 corresponds mainly to six rowed landraces of European origin. Cluster 2 is composed of Ethiopian landraces. Cluster 3 is mainly Asian landraces, from the Middle East to India and Korea. Cluster 4 is mainly European two-rowed landraces. Cluster 5 corresponds to six-rowed landraces from the Mediterranean region, and finally cluster 6 is mainly two-rowed landraces from the Middle East. The fact that barley accessions were clustered by their region of origin independent of the genotyping method used confirms the successful merger of the different SNP calling datasets.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Principal Component Analysis (PCA) graph of individuals for the genotyping dataset composed of 115,172 single nucleotide polymorphisms (SNPs) on the collection of 312 barley accessions, showing the results of the cluster analysis on the first two principal components (PC1 and PC2). Figures in brackets following the axis labels indicate the percentage of the total genotypic variance accounted for by the corresponding principal component. Cluster 1: mainly six rowed landraces of European origin. Cluster 2: Ethiopian landraces. Cluster 3: mainly Asian landraces, from the Middle East to India and Korea. Cluster 4: mainly European two-rowed landraces. Cluster 5: six-rowed landraces from the Mediterranean region. Cluster 6: mainly two-rowed landraces from the Middle East.</figcaption></figure>

The diagnostic SNP for *Rrs1*<sub>Rh4</sub> [9] was used to identify a list of 32 accessions containing this major resistance gene to rhynchosporium (Supplementary Table S1). This list includes Belgravia, Retriever, and 13 Syrian/Jordanian landraces previously shown to contain *Rrs1*<sub>Rh4</sub> [9] as well as 1 more Syrian landrace, 3 landraces from Vavilov Institute collection, and 13 landraces from the WHEALBI collection. The majority of barley landraces containing *Rrs1*<sub>Rh4</sub> originated in the Mediterranean region (Supplementary Table S1).

## 3.2. *Phenotyping*

The rhynchosporium scores showed a high level of variability between genotypes. To ensure rhynchosporium infection most field trials were autumn sown as barley sown in autumn in the locations used for field trials was consistently infected with *R. commune*. This is mainly due to relatively mild and wet conditions in winter and early spring in the UK and France which favour *R. commune* development. Rhynchosporium infection
occurred in each field trial, although the levels of disease severity and distribution patterns
were different (Figure 2). Indeed, Maule 2018, Pithiviers 2018, Dundee 2019 and Pithiviers
2019 reached higher infection levels than the other field trials. The only spring sown trial
used in this study, Maule 2018 spring sowing, showed one of the lowest levels of infection,
comparable only to Maule 2019 (Figure 2 and Supplementary Table S1). The shape of
the distribution of scores varied for each field trial, reflecting different weather patterns,
inoculum or other environmental conditions. The collection used in this study was also
highly diverse in relation to plant height, as it contains both tall and dwarf accessions, and
earliness, reflected by the heading date scores. The BLUE heading date and plant height
showed Gaussian distribution of the landraces, with an average heading date on 20 May
within an 18-day range, and an average plant height of 64 cm within a range of 31 to 94 cm
(Figure 3, Supplementary Table S1).

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Histograms representing the distribution of the rhynchosporium scores within the eight field trials, on a 1–9 scale, where 1 = no symptoms and 9 = all leaves infected. The field trials are (a) Dundee 2020, (b) Dundee 2019, (c) Maule 2019, (d) Pithiviers 2019, (e) Dunmow 2019, (f) Maule 2018, (g) Pithiviers 2018, (h) Maule spring sowing 2018.</figcaption></figure>

The highest correlation coefficients were between the rhynchosporium phenotypes in
two Dundee trials (Pearson correlation = 0.8) and between Maule 2018 and Pithiviers 2018
trials (0.78), followed by correlation between Maule 2018 and both Dundee trials (0.74), and
Pithiviers 2018 and Dundee 2019 (0.73) (Figure 4). By comparison, the correlation between
replicates for the fields with two replicates, Dundee 2019 and Dundee 2020, was 0.79 and
0.62, respectively (this difference between years can be explained by a higher level of infec-
tion at the end of the season in Dundee in 2019 compared with 2020, resulting in a more
reliable phenotyping). The correlation coefficients between rhynchosporium scores for the
remaining trials were moderate with the lowest correlation between Maule 2018 spring
sowing and all the autumn sown trials. Heading date and plant height were relatively
highly correlated (0.5) and significantly negatively correlated with the rhynchosporium
scores for all the field trials, except Maule 2018 spring sowing (Figure 4). This is in agree-
ment with plant height and earliness being linked to disease escape to rhynchosporium [3].
One of the resistance QTL identified by Looseley et al. [15] co-localised with the well-
characterised semi-dwarfing gene *sdw1* [46], which was previously shown to contribute to disease escape [25].

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Histograms representing the distribution of the best linear unbiased estimator (BLUE) for (a) the heading date (measured in days from 1 May) and (b) plant height (measured in cm between the soil and the bottom of the spike).</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Pearson correlation table between the rhynchosporium phenotypes of the eight field trials, and the best linear unbiased estimator for heading date and height (upper triangle). *p*-values are indicated with red stars where: \* = *p*-value < 0.05, \*\* = *p*-value < 0.01, and \*\*\* = *p*-value < 0.001. The names of the phenotypes are written diagonally, and the phenotypes are plotted depending on each other in the lower triangle.</figcaption></figure>

### 3.3. *Genome-Wide Association Studies*

The GWAS with the three methods EMMAX, BLINK, and FASTmrMLM resulted in seven QTL detected in at least two field trials. These QTL are located on chromosomes 3H, 5H, and 7H, with one, two, and four QTL detected on these chromosomes, respectively (Table 2 and Supplementary Figures S1–S16). All the QTL detected on chromosomes 3H and
7H correspond to genomic regions with published QTL for resistances to rhynchosporium. The QTL-7H-1 was detected in three field trials (Table 2) and was located within the physical interval for *Rrs2* [11]. The QTL-7H-1 allele associated with resistance was the minor allele in the collection, with a frequency of 0.22. The QTL-7H-2 is a QTL with the highest P-value and colocalises with QTL mapped in an Australian cultivar [47]. QTL-3H-1, QTL-7H-3, and QTL-7H-4 are also in areas where QTL were previously identified based on bi-parental population studies [10,19,22,48–51]. However, on chromosome 5H, QTL-5H-1 and QTL-5H-2 were detected in regions with no known resistances to rhynchosporium. Unlike chromosomes 2H, 3H, 4H, 6H, and 7H, chromosome 5H contain only five other rhynchosporium resistance QTL identified so far: Qsc5H-Shyri (0.9 Mb) [52], qSUK7_5 (456–526 Mb) [10], QTLCW5H.1 (555.7 Mb) [26], Qsc5H.3 (587.4 Mb) [19], and Qsc_5H_1 (638.2–659.5 Mb) [22]. The first two QTL on chromosome 7H are the only detected loci where the allele associated with resistance is the minor allele, the other putative resistances appear more widespread with frequencies of 0.64 to 0.91 in the landrace collection (Table 3).

<table><caption>**Table 3.** QTL for rhynchosporium scores detected by genome-wide association study (GWAS) in at least two of the eight field trials of the study.</caption><thead><tr><th>QTL Name</th><th>Chr<sup>a</sup></th><th>Physical Interval<sup>b</sup></th><th>Top SNP<sup>c</sup></th><th>Max<br/>-log<sub>10</sub>P<sup>d</sup></th><th>RAF<sup>e</sup></th><th>Nb<br/>Genes<sup>f</sup></th><th>Dataset<sup>g</sup></th><th>Methods<sup>h</sup></th><th>References<sup>i</sup></th></tr></thead><tbody><tr><td>QTL-3H-1</td><td>chr3H</td><td>673759655-674277867</td><td>chr3H_674155754</td><td>6.08</td><td>0.84</td><td>16</td><td>Dundee19; Pithiviers18</td><td>BLINK; FASTM-rMLM</td><td>QRh.S42-3H.a [49], QTLRS6 and QTLRS7b [48]</td></tr><tr><td>QTL-5H-1</td><td>chr5H</td><td>5560028-5561932</td><td>chr5H_5561636</td><td>5.01</td><td>0.85</td><td>1</td><td>Dundee19; Maule19</td><td>BLINK; FASTM-rMLM</td><td>reported in the current study only</td></tr><tr><td>QTL-5H-2</td><td>chr5H</td><td>394034355-423047705</td><td>chr5H_400584464</td><td>5.52</td><td>0.91</td><td>178</td><td>Dunmow19; Dundee19</td><td>BLINK; EMMAX</td><td>reported in the current study only</td></tr><tr><td>QTL-7H-1</td><td>chr7H</td><td>5389278-5394956</td><td>chr7H_5389278</td><td>5.31</td><td>0.22</td><td>2</td><td>Maule18; Dundee20; Dundee19</td><td>EMMAX; FASTM-rMLM</td><td>Rrs2 [11], Qsc_7H_1 [22]</td></tr><tr><td>QTL-7H-2</td><td>chr7H</td><td>68708954-68768630</td><td>chr7H_68766673</td><td>7.28</td><td>0.25</td><td>3</td><td>Dundee19; Pithiviers19</td><td>BLINK; FASTM-rMLM</td><td>QTLSR-7H-2017 [47]</td></tr><tr><td>QTL-7H-3</td><td>chr7H</td><td>113870715-193866189</td><td>chr7H_151131820</td><td>5.59</td><td>0.85</td><td>512</td><td>Dundee19; Dundee20</td><td>BLINK; FASTM-rMLM</td><td>QTLVixenRrs7H271 [50], Qsc-7H.3 [19]</td></tr><tr><td>QTL-7H-4</td><td>chr7H</td><td>638442686-638453794</td><td>chr7H_638446781</td><td>6.18</td><td>0.64</td><td>1</td><td>Dundee19; Maule19</td><td>BLINK; FASTM-rMLM</td><td>qS271_7 [10], Rrs15 [51], Qsc_7H_2 [22]</td></tr></tbody></table>

(a) Chromosome, (b) physical interval based on Morex V1 [30], (c) most significant single nucleotide polymorphism (SNP), (d) decimal logarithm of the *p*-value for the most significant SNP, (e) resistance allele frequency, (f) number of genes in the interval, (g) field trial dataset in which the QTL has been detected, (h) GWAS method where the QTL has a significant effect, (i) published QTL at the same loci.

This study used the highest number of markers for GWAS for resistance to rhynchosporium so far, 115.172 SNPs compared to under 9000 markers used by Looseley et al. [15], 23.549 SNPs used by Daba et al. [19], and 37.242 makes used by Cope et al. [20]. This resulted in identification of relatively narrow QTL intervals with the majority of them being under 1 Mb (Table 2). The smallest QTL interval is less than 0.01 Mb for QTL-7H-1 within the *Rrs2* interval. Only QTL-7H-3 interval stretches to 80 Mb, while QTL-5H-2 is 29 Mb. The intervals contain on average 102 genes, ranging from 1 gene for QTL-5H-1 and QTL-7H-4 to 512 genes for QTL-7H-3 (Table 2). The list of the genes present in each interval in the Morex 2017 genome annotation [30] is provided in Supplementary Table S3.

Four of the QTL (QTL-3H-1, QTL-5H-1, QTL-5H-2 and QTL-7H-3) contain kinases and a receptor-like kinase (RLK) is the only gene within the QTL-5H-1 interval. As *R. commune* is known to colonise the plant extracellular space, R genes recognising its effectors, or their effect on the plant, are likely to be cell surface-localised RLKs or receptor-like proteins
(RLPs) [2,53,54]. Therefore RLKs have been proposed as candidate genes for *Rrs1*, *Rrs18* and QSc.VR4 [9,10,13]. Previously, RLPs and RLKs were shown to be responsible for tomato *Cf-2*, *Cf-4* and *Cf-9* resistance against *C. fulvum* [55], oilseed rape resistance against *L. maculans* [56], wheat *Stb6* and *Stb16q* resistance against *Zymoseptoria tritici* [54,57], and rice *OsWAK1* resistance against *Magnaporthe oryzae* [58].

Another putative candidate gene for rhynchosporium resistance within QTL-3H-1 is a
germin-like protein (GLP). GLPs were previously found within another rhynchosporium
resistance QTL [13]. GLPs were shown to be involved in basal host resistance against
several fungal pathogens of cereals in barley [59,60], wheat [60] and rice [61].

It is not clear which genes could be putative candidate genes for QTL-7H-1 overlapping
with *Rrs2* and QTL-7H-4 overlapping with *Rrs15* as the genome annotation for these regions
does not contain any kinases. QTL-7H-2 overlapping with QTLSR-7H-2017 [47] is likely to
contain a quantitative resistance gene that does not have to be a RLK or RLP. While there
are only three annotated genes within this interval, none of them represent an obvious
candidate for a quantitative resistance gene. Considering that the genome annotation is
based on the North-American cultivar Morex genome assembly, which is not expected to
contain many of the quantitative resistance genes or any R genes to rhynchosporium, and
might not even have alleles of these genes, the gene content of these QTL would have to be
examined in the genotypes containing these resistances. GWAS ability to detect MTAs relies
on the allele frequency at QTL and, therefore, it is likely that rare resistance genes were
missed during this analysis. Despite the presence of diagnostic SNP for *Rrs1*<sub>Rh4</sub> in 10% of the
accessions used in this study, only a weak effect of this resistance, with −log<sub>10</sub>P reaching a
value above the significance threshold in only one of the field trials (*p*-value = 3.15 × 10<sup>-10</sup>
for Pithiviers 2018 with BLINK), was observed (Supplementary Figure S5). Previously
*Rrs1*<sub>Rh4</sub> locus was shown to be the main contributor to resistance in European spring barley
germplasm [15]. Interestingly, the disease nursery trials analysed in that study took place
at the same field in Dundee, UK between 2013 and 2015 but using a different germplasm
collection. This may indicate a loss of *Rrs1*<sub>Rh4</sub> efficiency in this field following the growth
of barley accessions containing this resistance gene over several years. The relatively low
effect of the *Rrs1*<sub>Rh4</sub> in field trials might also be a consequence of the wide use of this
resistance gene in commercially grown cultivars [15].

3.4. *Geographical Distribution of Resistance Alleles*

The countries of origin of the landraces were combined into wider geographical areas to study the regional frequency of the resistance QTL. The regions attributed to each accession are presented in the Supplementary Table S1. The Middle East region provided the highest number of genotypes, with 108 accessions, followed by Europe (87 accessions), Asia (52 accessions), North Africa (29 accessions), Ethiopia (18 accessions), North America (10 accessions), and South America (8 accessions). The frequencies of the alleles associated with resistance for each region are presented in Table 4.

The Ethiopian landraces carried the highest number of resistance QTL (the allele associated with resistance is the major allele for all QTL except in the case of QTL-7H-1, for which all the Ethiopian landraces possess an allele associated with susceptibility) while the Asian landraces showed the lowest number. The landraces from Asia and Europe had frequencies for alleles associated with resistance below 50% for 4 and 3 QTL located on chromosome 7H respectively. The landraces from the Middle East carried all the alleles associated with resistance at relatively high proportions (39% minimum). The allele frequencies for the landraces from North Africa, North America, and South America are based on very small numbers of accessions but still reflect the global trend of QTL-7H-1 and QTL-7H-2 having the lowest frequencies for alleles associated with resistance.
<table><caption>**Table 4.** Frequencies of the alleles associated with resistance, for the QTL most significant SNPs, across geographical regions of origin of the landraces.</caption><thead><tr><th>QTL Most Significant SNP</th><th>Asia</th><th>Ethiopia</th><th>Europe</th><th>Middle East</th><th>North Africa</th><th>North America</th><th>South America</th><th>Total</th></tr></thead><tbody><tr><td>chr3H_674155754</td><td>0.88</td><td>1.00</td><td>0.85</td><td>0.77</td><td>0.83</td><td>0.90</td><td>1.00</td><td>0.84</td></tr><tr><td>chr5H_5561636</td><td>0.50</td><td>1.00</td><td>0.98</td><td>0.91</td><td>0.72</td><td>1.00</td><td>1.00</td><td>0.85</td></tr><tr><td>chr5H_400584464</td><td>0.77</td><td>1.00</td><td>0.90</td><td>0.96</td><td>0.97</td><td>1.00</td><td>0.88</td><td>0.91</td></tr><tr><td>chr7H_5389278</td><td>0.17</td><td>0.00</td><td>0.06</td><td>0.45</td><td>0.07</td><td>0.30</td><td>0.13</td><td>0.22</td></tr><tr><td>chr7H_68766673</td><td>0.15</td><td>0.83</td><td>0.03</td><td>0.39</td><td>0.14</td><td>0.30</td><td>0.38</td><td>0.25</td></tr><tr><td>chr7H_151131820</td><td>0.48</td><td>0.94</td><td>0.93</td><td>0.93</td><td>0.79</td><td>1.00</td><td>1.00</td><td>0.85</td></tr><tr><td>chr7H_638446781</td><td>0.33</td><td>0.83</td><td>0.41</td><td>0.88</td><td>0.90</td><td>0.60</td><td>0.75</td><td>0.64</td></tr><tr><td>Average frequency</td><td>0.47</td><td>0.80</td><td>0.59</td><td>0.76</td><td>0.63</td><td>0.73</td><td>0.73</td><td>0.65</td></tr></tbody></table>

Geographical distribution of alleles associated with resistance in the global barley landrace collection used in this study clearly demonstrated the presence of resistance to rhynchosporium in barley accessions from the primary (Fertile Crescent) and secondary (Ethiopia) centres of barley diversity. This further supports the suggestion that at least some resistance to rhynchosporium has evolved in the Fertile Crescent [9] despite the theory of *R. commune* originating in Northern Europe [4]. Therefore, global barley accessions, including those from the Fertile Crescent and Ethiopia, represent a valuable source of resistance to rhynchosporium yet to be fully utilised by European breeders.

### 4. Conclusions

Association mapping is a powerful tool for detection of genomic regions associated with a trait of interest. Seven QTL were detected in at least two field trials in this study. The use of historical recombination events in germplasm collection allows mapping of traits to relatively small intervals and exploring the diversity of QTL available in a wide germplasm collection. However, the QTL detected by GWAS need to be confirmed by other methods. In our case, five of the QTL detected were previously described in bi-parental populations, and the importance of our study is in identification of relatively narrow intervals for these QTL. Our results provide insight into potential candidate genes annotated in these intervals, though the genetic content of these QTL should be investigated in cultivars where their presence have been confirmed. Further work is still needed to clone the resistance genes underpinning these QTL, but our study provides insight into the location of rhynchosporium resistance QTL in a spring barley collection representing global diversity and into the geographical distribution of the resistances to rhynchosporium. Achieving long-lasting protection against rhynchosporium would benefit from combining a wide range of complementary resistances. This requires further efforts in identification of diagnostic markers to make marker-assisted breeding for resistance to rhynchosporium a reality. The QTL identified in this study and the accessions containing alleles associated with resistance further expand the resources available for breeding varieties required for sustainable barley production.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agronomy12040782/s1, Figures S1-S16: Manhattan plots of the GWAS with the BLINK and EMMAX methods (FASTmrMLM do not produce Manhattan plots for its final results); Table S1: Details of the barley accessions used in the study, and phenotypes scored in the field trials; Table S2: Genotyping matrix containing alleles at the 115,172 SNPs used in the GWAS; Table S3: Gene content of the QTL detected in the GWAS, based on the Morex V1 genome annotation.

**Author Contributions:** Conceptualization, A.A., D.V., J.R. and M.L.; methodology, J.-N.T., M.L., M.B., A.A. and J.R.; software, M.B., J.-N.T. and M.L.; validation, J.-N.T. and M.B.; formal analysis, J.-N.T. and M.B.; investigation, J.-N.T. and A.A.; resources, A.A., J.R., D.V., M.L. and M.B.; data acquisition and interpretation, J.-N.T., A.A., P.-M.L.R., P.P., J.R., M.B., M.L. and R.W.; data curation, J.-N.T. and M.B.; writing—original draft preparation, J.-N.T. and A.A.; writing—review and editing, J.-N.T.,
A.A., J.R., M.B., M.L. and R.W.; visualization, J.-N.T.; supervision, A.A., J.R., M.L. and R.W.; project administration, A.A. and D.V.; funding acquisition, A.A., J.R., M.L. and D.V. All authors have read and agreed to the published version of the manuscript.

**Funding:** This project was funded by the Agriculture and Horticulture Development Board (AHDB). J.R., M.L., M.B., R.W., and A.A. were supported by the Scottish Government: Rural and Environment Science and Analytical Services (RESAS).

**Data Availability Statement:** The data presented in this study are available within the article and Supplementary Materials.

**Acknowledgments:** We would like to acknowledge Secobra Recherches for conducting the seed multiplication and the field trials at Maule, Pithiviers, and Dunmow. We are also thankful to Richard Keith, Christopher Warden, and Derek Matthews for the management of the field at the James Hutton Institute, Dundee.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. FAO. FAOSTAT Database. 2020. Available online: https://www.fao.org/faostat (accessed on 15 February 2022).

2. Avrova, A.; Knogge, W. *Rhynchosporium commune*: A persistent threat to barley cultivation. *Mol. Plant Pathol.* **2012**, *13*, 986–997.
[PubMed]

3. Zhan, J.; Fitt, B.D.L.; Pinnschmidt, H.O.; Oxley, S.J.P.; Newton, A.C. Resistance, epidemiology and sustainable management of *Rhynchosporium secalis* populations on barley. *Plant Pathol.* **2008**, *57*, 1–14. [CrossRef]

4. Zaffarano, P.L.; McDonald, B.A.; Zala, M.; Linde, C.C. Global Hierarchical Gene Diversity Analysis Suggests the Fertile Crescent Is Not the Center of Origin of the Barley Scald Pathogen *Rhynchosporium secalis*. *Phytopathology* **2006**, *96*, 941–950. [CrossRef]
[PubMed]

5. AHDB. *Barley Disease Management Guide*; Agriculture and Horticulture Development Board; AHDB: Warwickshire, UK, 2016.

6. Cooke, L.R.; Locke, T.; Lockley, K.D.; Phillips, A.N.; Sadiq, M.D.S.; Coll, R.; Black, L.; Taggart, P.J.; Mercer, P.C. The effect of fungicide programmes based on epoxiconazole on the control and DMI sensitivity of *Rhynchosporium secalis* in winter barley. *Crop Prot.* **2004**, *23*, 393–406. [CrossRef]

7. Lamichhane, J.R.; Dachbrodt-Saaydeh, S.; Kudsk, P.; Messéan, A. Toward a Reduced Reliance on Conventional Pesticides in European Agriculture. *Plant Dis.* **2016**, *100*, 10–24. [CrossRef]

8. Zaffarano, P.L.; McDonald, B.A.; Linde, C.C. Two new species of *Rhynchosporium*. *Mycologia* **2011**, *103*, 195–202. [CrossRef]

9. Looseley, M.; Griffe, L.; Buttner, B.; Wright, K.; Bayer, M.; Coulter, M.; Thauvin, J.-N.; Middlefell-Williams, J.; Maluk, M.; Okpo, A.; et al. Characterisation of barley landraces from Syria and Jordan for resistance to rhynchosporium and identification of diagnostic markers for *Rrs1Rhl4*. *Theor. Appl. Genet.* **2020**, *133*, 1243–1264. [CrossRef]

10. Coulter, M.; Buttner, B.; Hofmann, K.; Bayer, M.; Ramsay, L.; Schweizer, G.; Waugh, R.; Looseley, M.; Avrova, A. Characterisation of barley resistance to rhynchosporium on chromosome 6HS. *Theor. Appl. Genet.* **2019**, *132*, 1089–1107. [CrossRef]

11. Hanemann, A.; Schweizer, G.F.; Cossu, R.; Wicker, T.; Röder, M.S. Fine mapping, physical mapping and development of diagnostic markers for the *Rrs2* scald resistance gene in barley. *Theor. Appl. Genet.* **2009**, *119*, 1507–1522. [CrossRef]

12. Marzin, S.; Hanemann, A.; Sharma, S.; Hensel, G.; Kumlehn, J.; Schweizer, G.; Röder, M.S. Are pectin esterase inhibitor genes involved in mediating resistance to *Rhynchosporium commune* in barley? *PLoS ONE* **2016**, *11*, e0150485. [CrossRef]

13. Wang, Y.; Xu, Y.; Gupta, S.; Zhou, Y.; Wallwork, H.; Zhou, G.; Broughton, S.; Zhang, X.-Q.; Tan, C.; Westcott, S.; et al. Fine mapping QSc.VR4, an effective and stable scald resistance locus in barley (*Hordeum vulgare* L.), to a 0.38-Mb region enriched with LRR-RLK and GLP genes. *Theor. Appl. Genet.* **2020**, *133*, 2307–2321. [CrossRef] [PubMed]

14. Zhang, X.; Ovenden, B.; Milgate, A. Recent insights into barley and *Rhynchosporium commune* interactions. *Mol. Plant Pathol.* **2020**, *21*, 1111–1128. [CrossRef] [PubMed]

15. Looseley, M.E.; Griffe, L.L.; Buttner, B.; Wright, K.M.; Middlefell-Williams, J.; Bull, H.; Shaw, P.; Macaulay, M.; Booth, A.; Schweizer, G.; et al. Resistance to *Rhynchosporium commune* in a collection of European spring barley germplasm. *Theor. Appl. Genet.* **2018**, *131*, 2513–2528. [CrossRef] [PubMed]

16. Zhan, J.; Yang, L.; Shang, L.; Newton, A.C. Pathogen Populations Evolve to Greater Race Complexity in Agricultural Systems— Evidence from Analysis of *Rhynchosporium secalis* Virulence Data. *PLoS ONE* **2012**, *7*, e38611. [CrossRef]

17. Milner, S.G.; Jost, M.; Taketa, S.; Mazón, E.R.; Himmelbach, A.; Oppermann, M.; Weise, S.; Knüpffer, H.; Basterrechea, M.; König, P.; et al. Genebank genomics highlights the diversity of a global barley collection. *Nat. Genet.* **2019**, *51*, 319–326. [CrossRef]

18. Alqudah, A.M.; Sallam, A.H.; Baenziger, P.S.; Borner, A. GWAS: Fast-forwarding gene identification and characterization in temperate Cereals: Lessons from Barley—A review. *J. Adv. Res.* **2020**, *22*, 119–135. [CrossRef] [PubMed]

19. Daba, S.D.; Horsley, R.; Brueggeman, R.; Chao, S.; Mohammadi, M. Genome-wide Association Studies and Candidate Gene Identification for Leaf Scald and Net Blotch in Barley (*Hordeum vulgare* L.). *Plant Dis.* **2019**, *103*, 880–889. [CrossRef]
20. Cope, J.E.; Norton, G.J.; George, T.S.; Newton, A.C. Identifying potential novel resistance to the foliar disease 'Scald' (*Rhynchosporium commune*) in a population of Scottish Bere barley landrace (*Hordeum vulgare* L.). *J. Plant Dis. Prot.* **2021**, *128*, 999–1012.
[CrossRef]

21. Buttner, B.; Draba, V.; Pillen, K.; Schweizer, G.; Maurer, A. Identification of QTLs conferring resistance to scald (*Rhynchosporium commune*) in the barley nested association mapping population HEB-25. *BMC Genom.* **2020**, *21*, 837. [CrossRef]

22. Hautsalo, J.; Novakazi, F.; Jalli, M.; Goransson, M.; Mannisen, O.; Isolahti, M.; Reitan, L.; Bergersen, S.; Krusell, L.; Damsgard Robertsen, C.; et al. Pyramiding of scald resistance genes in four spring barley MAGIC populations. *Theor. Appl. Genet.* **2021**, *134*, 3829–3843. [CrossRef]

23. Bustos-Korts, D.; Dawson, I.K.; Russell, J.; Tondelli, A.; Guerra, D.; Ferrandi, C.; Nicolazzi, E.L.; Molnar-Lang, M.; Ozkan, H.; Megyeri, M.; et al. Exome sequences and multi-environment field trials elucidate the genetic basis of adaptation in barley. *Plant J.* **2019**, *99*, 1172–1191. [CrossRef] [PubMed]

24. Soboleva, O.N.; Konovalova, G.S.; Radchenko, E.E. Genetic control of scald resistance in barley landraces. *Vavilov J. Genet. Breed.* **2016**, *20*, 616–622. [CrossRef]

25. Looseley, M.E.; Newton, A.C.; Atkins, S.D.; Fitt, B.D.; Fraaije, B.A.; Thomas, W.T.B.; Keith, R.; Macaulay, M.; Lynott, J.; Harrap, D. Genetic basis of control of *Rhynchosporium secalis* infection and symptom expression in barley. *Euphytica* **2012**, *184*, 47–56. [CrossRef]

26. Looseley, M.E.; Keith, R.; Guy, D.; Barral-Baron, G.; Thirugnanasambandam, A.; Harrap, D.; Werner, P.; Newton, A.C. Genetic mapping of resistance to *Rhynchosporium commune* and characterisation of early infection in a winter barley mapping population. *Euphytica* **2015**, *203*, 337–347. [CrossRef]

27. Zadoks, J.C.; Chang, T.T.; Konzak, C.F. A decimal code for the growth stages of cereals. *Weed Res.* **1974**, *14*, 415–421.

28. Bates, D.; Machler, M.; Bolker, B.; Walker, S. Fitting Linear Mixed-Effects Models Using lme4. *J. Stat. Softw.* **2015**, *7*, 48. [CrossRef]

29. Mosleth, E.F.; Lillehammer, M.; Pellny, T.K.; Wood, A.J.; Riche, A.B.; Hussain, A.; Griffiths, S.; Hawkesford, M.J.; Shewry, P.R. Genetic variation and heritability of grain protein deviation in European wheat genotypes. *Field Crops Res.* **2020**, *255*, 107896. [CrossRef]

30. Mascher, M.; Gundlach, H.; Himmelbach, A.; Beier, S.; Twardziok, S.O.; Wicker, T.; Radchuk, V.; Dockter, C.; Hedley, P.E.; Russell, J.; et al. A chromosome conformation capture ordered sequence of the barley genome. *Nature* **2017**, *544*, 427–433. [CrossRef]

31. Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. *arXiv* **2013**, arXiv:1303.3997.

32. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabrie, S.; Daly, M.; et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. *Genome Res.* **2010**, *20*, 1297–1303. [CrossRef]

33. Van der Auwera, G.A.; Carneiro, M.O.; Hartl, C.; Poplin, R.; Del Angel, G.; Levy-Moonshine, A.; Jordan, T.; Shakir, K.; Roazen, D.; Thibault, J.; et al. From FastQ data to high confidence variant calls: The Genome Analysis Toolkit best practices pipeline. *Curr. Protoc. Bioinform.* **2013**, *43*, 11.10.1–11.10.33. [CrossRef]

34. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. *Bioinformatics* **2013**, *29*, 15–21. [CrossRef] [PubMed]

35. R Core Team. R: A Language and Environment for Statistical Computing. 2019. Available online: https://www.R-project.org/ (accessed on 22 February 2022).

36. Schwender, H.; Fritsch, A. Package Scrime: Analysis of High-Dimensional Categorical Data Such as SNP Data. 2008. Available online: https://cran.r-project.org/web/packages/scrime/ (accessed on 22 March 2022).

37. Endelman, J.B. Ridge regression and other kernels for genomic selection with R package rrBLUP. *Plant Genome* **2011**, *4*, 255–258. [CrossRef]

38. Kang, H.M.; Sul, J.H.; Service, S.K.; Zaitlen, N.A.; Kong, S.-Y.; Freimer, N.B.; Sabatti, C.; Eskin, E. Variance component model to account for sample structure in genome-wide association studies. *Nat. Genet.* **2010**, *42*, 348–354. [CrossRef] [PubMed]

39. Wang, S.-B.; Feng, J.-Y.; Ren, W.-L.; Huang, B.; Zhou, L.; Wen, Y.-J.; Zhang, J.; Dunwell, J.M.; Xu, S.; Zhang, Y.-M. Improving power and accuracy of genome-wide association studies via a multi-locus mixed linear model methodology. *Nat. Sci. Rep.* **2016**, *6*, 19444. [CrossRef]

40. Tamba, C.L.; Zhang, Y.-M. A fast mrMLM algorithm for multi-locus genome-wide association studies. *BioRxiv* **2018**. [CrossRef]

41. Huang, M.; Liu, X.; Zhou, Y.; Summers, R.M.; Zhang, Z. BLINK: A package for the next level of genome-wide association studies with both individuals and markers in the millions. *Gigascience* **2019**, *8*, 154. [CrossRef]

42. Wang, J.; Zhang, Z. GAPIT Version 3: Boosting Power and Accuracy for Genomic Association and Prediction. *Genom. Proteom. Bioinform.* **2021**. [CrossRef]

43. Tibbs Cortes, L.; Zhang, Z.; Yu, J. Status and prospects of genome-wide association studies in plants. *Plant Genome* **2021**, *14*, 20077. [CrossRef]

44. Merrick, L.F.; Burke, A.B.; Zhang, Z.; Carter, A.H. Comparison of Single-Trait and Multi-Trait Genome-Wide Association Models and Inclusion of Correlated Traits in the Dissection of the Genetic Architecture of a Complex Trait in a Breeding Program. *BioRxiv* **2021**. Available online: https://www.biorxiv.org/content/10.1101/2021.08.23.457367v2 (accessed on 22 February 2022).

45. Mangin, B.; Siberchicot, A.; Nicolas, S.; Doligez, A.; This, P.; Cierco-Ayrolles, C. Novel measures of linkage disequilibrium that correct the bias due to population structure and relatedness. *Heredity* **2012**, *108*, 285–291. [CrossRef]
46. Xu, Y.; Jia, Q.; Zhou, G.; Zhang, X.-Q.; Angessa, T.; Broughton, S.; Yan, G.; Zhang, W.; Li, C. Characterization of the *sdw1* semi-dwarf gene in barley. *BMC Plant Biol.* **2017**, *17*, 11. [CrossRef]

47. Zhang, X.; Ovenden, B.; Orchard, B.A.; Zhou, M.; Park, R.F.; Singh, D.; Milgate, A. Bivariate analysis of barley scald resistance with relative maturity reveals a new major QTL on chromosome 3H. *Nat. Sci. Rep.* **2019**, *9*, 19258. [CrossRef]

48. Sayed, H.; Backes, G.; Kayyal, H.; Yahyaoui, A.; Ceccarelli, S.; Grando, S.; Jahoor, A.; Baum, M. New molecular markers linked to qualitative and quantitative powdery mildew and scald resistance genes in barley for dry areas. *Euphytica* **2004**, *135*, 225–228. [CrossRef]

49. von Korff, M.; Wang, H.; Léon, J.; Pillen, K. AB-QTL analysis in spring barley. I. Detection of resistance genes against powdery mildew, leaf rust and scald introgressed from wild barley. *Theor. Appl. Genet.* **2005**, *111*, 583–590. [CrossRef] [PubMed]

50. Wagner, C.; Schweizer, G.; Kramer, M.; Dehmer-Badani, A.G.; Ordon, F.; Friedt, W. The complex quantitative barley-Rhynchosporium secalis interaction: Newly identified QTL may represent already known resistance genes. *Theor. Appl. Genet.* **2008**, *118*, 113–122. [CrossRef]

51. Genger, R.K.; Nesbitt, K.; Brown, A.H.D.; Abbott, D.C.; Burdon, J.J. A novel barley scald resistance gene: Genetic mapping of the *Rrs15* scald resistance gene derived from wild barley, *Hordeum vulgare* ssp. *spontaneum*. *Plant Breed.* **2005**, *124*, 137–141. [CrossRef]

52. Zantinge, J.; Xue, S.; Holtz, M.; Xi, K.; Juskiw, P. The identification of multiple SNP markers for scald resistance in spring barley through restriction-site associated sequencing. *Euphytica* **2019**, *215*, 8. [CrossRef]

53. Stotz, H.; Mitrousia, G.; de Wit, P.; Fitt, B.D.L. Effector-triggered defence against apoplastic fungal pathogens. *Trends Plant Sci.* **2014**, *19*, 491–500. [CrossRef]

54. Saintenac, C.; Lee, W.S.; Cambon, F.; Rudd, J.J.; King, R.C.; Marande, W.; Powers, S.J.; Bergès, H.; Phillips, A.L.; Uauy, C.; et al. Wheat receptor kinase-like protein *Stb6* controls gene-for-gene resistance to fungal pathogen *Zymoseptoria tritici*. *Nat. Genet.* **2018**, *50*, 368–374. [CrossRef]

55. Liebrand, T.W.H.; van den Berg, G.C.M.; Zhang, Z.; Smit, P.; Cordewener, J.H.G.; America, A.H.P.; Sklenar, J.; Jones, A.M.E.; Tameling, W.I.L.; Robatzek, S.; et al. Receptor like kinase SOBIR1/EVR interacts with receptor-like proteins in plant immunity against fungal infection. *Proc. Natl. Acad. Sci. USA* **2013**, *110*, 10010–10015. [CrossRef]

56. Larkan, N.J.; Lydiate, D.J.; Parkin, I.A.P.; Nelson, M.N.; Epp, D.J.; Cowling, W.A.; Rimmer, S.R.; Borhan, M.H. The *Brassica napus* blackleg resistance gene *LepR3* encodes a receptor-like protein triggered by the *Leptosphaeria maculans* effector *AVRLM1*. *New. Phytol.* **2013**, *197*, 595–599. [CrossRef] [PubMed]

57. Saintenac, C.; Cambon, F.; Aouini, L.; Verstappen, E.; Ghaffary, S.M.T.; Poucet, T.; Marande, W.; Berges, H.; Xu, S.; Jaouannet, M.; et al. A wheat cysteine-rich receptor-like kinase confers broad-spectrum resistance against *Septoria tritici* blotch. *Nat. Commun.* **2021**, *12*, 433. [CrossRef] [PubMed]

58. Li, H.; Zhou, S.-Y.; Zhao, W.-S.; Su, S.-C.; Peng, Y.-L. A novel wall-associated receptor-like protein kinase gene, OsWAK1, plays important roles in rice blast disease resistance. *Plant Mol. Biol.* **2009**, *69*, 337–346. [CrossRef] [PubMed]

59. Zimmermann, G.; Baumlein, H.; Mock, H.-P.; Himmelbach, A.; Schweizer, P. The multigene family encoding germin-like proteins of barley. Regulation and function in Basal host resistance. *Plant Physiol.* **2006**, *142*, 181–192. [CrossRef]

60. Christensen, A.B.; Thordal-Christensen, H.; Zimmermann, G.; Gjetting, T.; Lyngkjær, M.F.; Dudler, R.; Schweizer, P. The Germinlike Protein GLP4 Exhibits Superoxide Dismutase Activity and Is an Important Component of Quantitative Resistance in Wheat and Barley. *MPMI* **2007**, *17*, 109–117. [CrossRef]

61. Manosalva, P.M.; Davidson, R.M.; Liu, B.; Zhu, X.; Hulbert, S.H.; Leung, H.; Leach, J.E. A Germin-Like Protein Gene Family Functions as a Complex Quantitative Trait Locus Conferring Broad-Spectrum Disease Resistance in Rice. *Plant Physiol.* **2009**, *149*, 286–296. [CrossRef]