High Resolution Mapping of RphMBR1012 Conferring Resistance to Puccinia hordei in Barley (Hordeum vulgare L.)

Isolation of disease resistance genes in barley was hampered by the large genome size, but has become easy due to the availability of the reference genome sequence. During the last years, many genomic resources, e.g., the Illumina 9K iSelect, the 50K Infinium arrays, the Barley Genome Zipper, POPSEQ, and genotyping by sequencing (GBS), were developed that enable enhanced gene isolation in combination with the barley genome sequence. In the present study, we developed a fine map of the barley leaf rust resistance gene RphMBR1012. 537 segmental homozygous recombinant inbred lines (RILs) derived from 4775 F2-plants were used to construct a high-resolution mapping population (HRMP). The Barley Genome Zipper, the 9K iSelect chip, the 50K Infinium chip and GBS were used to develop 56 molecular markers located in the target interval of 8 cM. This interval was narrowed down to about 0.07 cM corresponding to 0.44 Mb of the barley reference genome. Eleven low-confidence and 18 high-confidence genes were identified in this interval. Five of these are putative disease resistance genes and were subjected to allele-specific sequencing. In addition, comparison of the genetic map and the reference genome revealed an inversion of 1.34 Mb located distally to the resistance locus. In conclusion, the barley reference sequence and the respective gene annotation delivered detailed information about the physical size of the target interval, the genes located in the target interval and facilitated the efficient development of molecular markers for marker-assisted selection for RphMBR1012.


INTRODUCTION
Leaf rust of barley is a serious disease caused by the biotrophic fungus Puccinia hordei Otth., which, under favorable conditions, may cause yield losses of up to 62% , while in general loses are about 15-25% (Whelan et al., 1997). Symptoms of leaf rust vary from small chlorotic flecks to large orange-brown pustules of up to 0.5 mm in size, often surrounded by green islands (Clifford, 1985). Although several resistance genes in barley have been identified, the major challenge in control of barley leaf rust is the breakdown of resistance caused by mutations in effector (avirulence) genes of the pathogen, leading to occurrence of new virulent races on previously resistant plant cultivars in a short period of time (Park, 2003). Therefore, to combat leaf rust epidemics caused by newly occurring/generated virulent races and to achieve a sustainable disease control, the employment of new resistance genes using functional molecular markers in breeding schemes as well as the isolation of known ones in order to get detailed information on the structure and function is of prime importance. Furthermore, isolation of known resistance genes is a prerequisite that allow an efficient allele mining of genetic resources (Li et al., 2016) as well as allele editing, e.g., by CRISPR/Cas9 (Wang et al., 2014).
Molecular markers have been widely used in barley breeding for mapping of genes, marker-assisted selection, as well as in positional isolation of genes (Stein and Graner, 2005;Perovic et al., 2018). The most abundant molecular markers are single nucleotide polymorphism (SNP). Employing next generation sequencing (Ganal et al., 2018), SNPs are easily detectable in a high throughput manner and are therefore currently the markers of choice. The number of available SNP markers rapidly increased from about 180 EST markers to about 6,800 SNPs on the 9K Illumina iSelect chip up to 44,040 SNPs on the 50K Illumina Infinium array (Kota et al., 2003;Rostoks et al., 2005;Stein et al., 2007;Close et al., 2009;Muñoz-Amatriaín et al., 2011;Comadran et al., 2012;Bayer et al., 2017). The barley Genome Zipper (GZ) assembled 86% of the barley genes in a putative linear order (Mayer et al., 2011). Population sequencing methodology (POPSEQ) was developed as an integrated method to create a linear order of contigs using whole-genome-shotgun sequencing (WGS) data that resulted in the first ultra-high density map of the barley genome (Mayer et al., 2011;Mascher et al., 2013a). Assessment of the GZ and POPSEQ by Silvar et al. (2015) at seven loci mapped with higher genetic resolution revealed an accuracy of 97.8% with respect to the GZ and 99.3% to POPSEQ in comparison to consensus genetic maps. In addition to the above mentioned resources, advances in target capture/enrichment and next-generation sequencing, like GBS (Poland et al., 2012), exome capture (Mascher et al., 2013b), and barley reference genome sequence (Mascher et al., 2017) are available for marker development.
Although high resolution mapping allows precise zooming into targeted loci, the un-even distribution of crossovers along chromosomes (International Barley Genome Sequencing Consortium [IBSC], 2012) and the large variation in the genetic/physical ratio across the genome  often hampers high-resolution genetic dissection. In barley, pericentromeric regions (pCENR) comprise at least 48% of the physical genome but harbor only 14-22% of the total barley gene content (Mascher et al., 2017). The other extreme are hotspots of high recombination rates in telomeric regions (Bhakta et al., 2015). In case of the locus of Ryd3, which is located in a centromeric region, the physical/genetic ratio has been estimated at 14-60 Mb/cM, while the genome-wide average is 4.4 Mb/cM (Lüpken et al., 2014). At the rym4/rym5 locus, the ratio of physical to genetic distances was in the range between 0.8 and 2.3 Mb cM and have increased to over 30 Mb cM, although the gene has been mapped on the telomeric region of chromosome 3H (Stein and Graner, 2005). This indicates that a large number of meiotic events is essential for a sufficient genetic resolution to detect recombination events in close vicinity to the targeted genes, and highlights the need for development of HRMPs.
Diverse collections of barley germplasm were evaluated for detecting new sources of leaf rust resistance (Perovic et al., 2003). In this respect, the Rph MBR1012 gene was mapped on the short arm of chromosome 1H (König et al., 2012), where only Rph4 has previously been localized (McDaniel and Hathcock, 1969;König et al., 2012). Prior to the recently cloned gene Rph1 by Dracatos et al. (2018), all efforts to isolate leaf rust resistance genes in barley were unsuccessful. An example of unsuccessful isolation is the case of Rph7 (Brunner et al., 2000;Scherrer et al., 2005). Hence, positional cloning is still one of the most efficient and reliable approaches to isolate a resistance gene in crop species with large genomes, such as wheat and barley (Krattinger et al., 2009). In barley, up to now five genes conferring resistance to fungal and viral pathogens have been isolated through map-based cloning, comprising mlo (Büschges et al., 1997;Simons et al., 1997), Mla6 (Halterman et al., 2001), Rpg1 (Brueggeman et al., 2002), rym4/rym5 (Pellio et al., 2005) and rym11 (Yang et al., 2014).
The aims of this study were to: (i) develop at HRMP for the Rph MBR1012 resistance gene, (ii) saturate the locus using all available state-of-the-art genomic resources i.e., GBS, 50K Infinium and the barley reference genome, (iii) anchor the genetic map to the barley reference sequence (iv) characterize the putative candidate rust resistance genes by allele specific re-sequencing and (v) test the developed markers for their diagnostic value.

Plant Material and Construction of a High-Resolution Mapping Population
For high resolution mapping of Rph MBR1012 , a segregating population comprising of 4,775 F 2 plants was constructed based on crosses between five DH-lines namely, the resistant (R) DH3/6 and DH3/127 and the susceptible (S) DH3/9, DH3/62 and DH3/74, which were derived from the original cross between the parental line MBR1012 (resistant) and Scarlett (susceptible). Based on these five DH-lines four crosses were conducted, namely DH3/74 (S) × DH3/6 (R), DH3/74 (S) × DH3/127 (R), DH3/6 (R) × DH3/9 (S) and DH3/62 (S) × DH3/127 (R) ( Table 1). In order to identify recombinants, F 2 plants were analyzed using two flanking co-dominant SSRs, i.e., QBS94 (distal) and QBS113 (proximal) . Respective markers were analyzed by capillary electrophoresis at the genetic analyzer ABI PRISM R 3100 (Applied Biosystems, Darmstadt, Germany). From identified heterozygous recombinant F 2 plants in target interval, 12 progeny plants, representing F 3 families were sown in 96 Quick pot plates. Genomic DNA of 10 days old plantlets was extracted in F 2 and F 3 according to Dorokhov and Klocke (1997). The quality of the extracted genomic DNA was checked by electrophoresis on 1% agarose gel and latter quantified by using the NanoDrop ND-100 spectrophotometer (PeQLab, Erlangen, Germany). By this approach, a HRMP of 537 recombinant inbred lines (RILs) was developed and subsequently used for marker saturation and resistance testing. Genomic DNA of the selected segmental homozygous RILs was extracted using the Miniprep method according to Stein et al. (2001). DNA of all samples was adjusted to a final concentration of 20 ng/µl. Furthermore, F 3 recombinant plants were self-fertilized and as F 4 segmental RILs used for phenotyping and genotyping with newly developed PCR based markers.

Inoculum Preparation
Fresh urediniospores of leaf rust isolate I-80 were prepared by artificial inoculation at the two-leaf stage of Hordeum vulgare cultivar Grossklappige, which is highly susceptible to the majority of P. hordei isolates. Inoculated plants were covered with plastic for 24 h at 18 • C to ensure a moist environment. After 15 days, rust urediniospores were harvested and used for inoculation of RILs seedlings.

Resistance Tests
Resistance tests were carried out in the greenhouse by inoculation of RILs along with the two H. vulgare parental lines, i.e., MBR1012 (resistant), Scarlett (susceptible) and susceptible (DH3/62) and resistant (DH3/127) DH-lines as well as the cv. Grossklappige as a control. Three plants per segmental RILs were sown in 96 Quick pot trays and 10 days old plantlets were inoculated with fresh I-80 urediniospores according to Ivandic et al. (1998). Briefly, 10 mg of fresh spores were used per 100 plants and mixed with white clay (Laborchemie Apolda, Germany), (1:3). The inoculated plants were kept at 18 • C and covered with plastic for 24 h, providing a moist environment for successful infection. All plants were scored at two time points, i.e., 10 and 13 days post-inoculation (dpi) according to Levine and Cherewick (1952). Segregation of resistant and susceptible plants was analyzed using the Chi-square (χ 2 ) tests for goodness-of-fit to the expected Mendelian segregation ratios.

Marker Development
For marker saturation, initially 6 Simple Sequence Repeats (SSRs), 7 size polymorphism and 24 SNPs markers derived from the barley GZ and 9K iSelect high-density custom genotyping bead chip were used for random saturation of the large interval of about 8 cM (Perovic et al., personal communication), while the Illumina 50K Infinium array and Genotyping By Sequencing (GBS) were used in combination with the barley reference sequence (Mascher et al., 2017) for very precise marker saturation within an interval of 0.1 cM of the locus in this study (Supplementary Table S1).

50K iSelect Illumina SNP Array
The genomic DNA of parental lines, two DH-lines and two RILs from HRMP (carrying critical recombination within the resistance locus region) were used for the identification of polymorphic SNPs derived from the 50K Infinium array (TraitGenetics Gatersleben, Germany). The polymorphic SNPs located in the target interval were converted into Kompetitive Allele Specific PCR (KASP) assays by designing the two allele-specific forward primers, and one common reverse primer spanning the sequence of interest carrying the SNP position using Primer3 v. 0.4.0 1 (Koressaar and Remm, 2007;Untergasser et al., 2012). KASP markers were then used for genotyping of the HRMP.

Genotyping-by-Sequencing (GBS)
The same lines as for the 50K array were used for GBS screening. A 20 ng/µl of genomic DNA of each line was used for GBS according to Wendler et al. (2014). Sequencing of selected lines was done on Illumina R MiSeq TM (Illumina, San Diego, United States). Sequencing data were analyzed using the Galaxy platform (Blankenberg et al., 2001;Giardine et al., 2005;Goecks et al., 2010) implemented at the JKI. After adapter and quality trimming (trim galore version 0.2.8.1; quality < 30, read length > 50), read mapping of the GBS data was executed using BWA version 0.7.15-r1140 (Li and Durbin, 2009) with standard settings to map the reads to the pseudomolecules of barley (Mascher et al., 2017). SNP calling was performed using mpileup version 1.2 (Li and Durbin, 2009), with genotype likelihood computation. Missing data was imputed with Beagle v4.1 (Browning and Browning, 2016). Biallelic SNPs were detected and subsequently filtered for differences between the resistant and susceptible parental lines and a minimum coverage of five reads per SNP using SnpSift version 4.2 (Cingolani et al., 2012). KASP markers were designed for polymorphic SNPs positioned in the target region 2 .
The following PCR conditions were used for all SSRs, size polymorphism and CAPS markers: denaturation at 94 • C for 5 min followed by 12 cycles at 94 • C for 30 s, annealing at 62 • C to 56 • C (-0.5 • C/cycle) for 30 s, extension 30 s at 72 • C, 94 • C for 30 s, 56 • C for 30 s, 72 • C 30 s, 35 cycles, final extension at 72 • C for 10 min.
The PCR amplification condition for KASP markers were: 10 min at 94 • C, followed by 10 cycles: 94 • C for 20 s, annealing at 61 • C to 55 • C (-0.6 • C/cycle) for 60 s, followed by 26 cycles: 94 • C for 20 s, 55 • C for 60 s, 30 • C 60 s. The real-time PCR machine was used to detect the fluorescence from HEX and FAM on plate reads. After thermal cycling was completed, the fluorescent signal was detected by reading the plate in the qPCR machine at 37 • C. At the end of the run the results were shown in the data analysis software under "Allelic Discrimination." The software automatically showed the clusters for the alleles for samples based on their position in the allelic discrimination plot (LGC, Guide to running KASP genotyping on the BIO-RAD CFX-series instruments').

Linkage Analysis
Linkage analysis was performed by dividing the number of the recombination events with the number of analyzed gametes, multiplied with 100. The recombination frequency was used for the genetic linkage map construction and visualized using MapChart (Voorrips, 2002)

Anchoring the Genetic Map to the Barley Reference Sequence
All 56 markers used for construction of the HRMP were anchored to the barley Reference genome sequence (Mascher et al., 2017). All sequences including forward and revers primers were blasted against the barley reference genome sequence 4 using BLASTN algorithm applying default parameters. Obtained physical positions of mapped markers were visualized using software MapChart (Voorrips, 2002).

Use of the Barley Reference Sequence for the Identification of Candidate Genes
Marker positions in the barley reference sequence were used to determine the target interval of the resistance gene locus and to extract putative candidate genes 5 . After defining the genomic region of the resistance locus at the barley reference sequence, High-Confidence (HC) and Low-Confidence (LC) genes including Exon-intron boundaries were extracted from the available annotation (Mascher et al., 2017). The reconstruction of the gene intron-exon-structure was performed using the internet platform "Splign" 6 from NCBI, which allows alignment of mRNA to genomic sequence (Kapustin et al., 2008).

Allele Specific Re-sequencing of Candidate Genes
The allele specific re-sequencing of candidate genes was conducted for 18 high and 11 low confidence genes positioned in the candidate interval. Online software Primer3 v. 0.4.0 (see footnote 1) (Koressaar and Remm, 2007;Untergasser et al., 2012) setting the parameters at 20-22 bp, temperature 58-62 • C and product size of 350 bp was used for primer design, which subsequently were then tested for their specificity for chromosome 1H using the barley blast server(see footnote 4) against the barley pseudomolecules according to Mascher et al. (2017). In the first round of low pass resequencing, a set of 36 primer pairs were designed covering all 29 high and low confidence genes. In the second round of the experiment, 25 primer pairs were designed in order to sequence the full length of five disease resistance genes. To sequence the entire gene, Morex contigs including the gene sequence of each disease resistance gene were identified using (see footnote 4) allowing to design primers at least 20 bases upstream of the start codon and 20 bases downstream of the stop codon. Moreover, the primers should overlap to ensure that there are no gaps between the fragments after sequence analysis. Next, obtained PCR products of the same size were subjected for sequencing. PCR fragments were separated by agarose gel electrophoresis and analyzed using the imaging system Gel Doce TM XR and the Quantity One R 1-D analysis software (4.6.2) (Bio-Rad, Hercules, United States) and subsequently sequenced by the company Microsynth AG (Balgach, Switzerland) using the Sanger sequencing method (Sanger et al., 1977). Obtained sequences were edited and analyzed using Sequencher 5.1 software (Gene Codes, Ann Arbor, MI, United States) using default parameters.
Functional analysis of identified polymorphisms between parental lines (MBR1012 and Scarlett) was done using the multiple sequence alignment program, MAFFT by default parameters (Katoh and Standley, 2013).
Phenotypic analysis of resistance to Rph MBR1012 showed a segregation of 261 resistant and 276 susceptible RILs and revealed TABLE 3 | Selected Bowman lines and parental lines carrying 25 known Rph genes for diagnostic value evaluation of developed markers linked to the resistance locus, size of alleles and restriction patterns.  the expected 1r:1s segregation ratio among these RILs. Chisquare test (χ 2 1:1 = 0.4189, df = 1, p < 0.05) for goodness of fit indicated that the resistance in MBR1012 is monogenically controlled (Figure 1 and Table 1).

Locus and Anchoring to the Barley Reference Sequence
A fine map of the Rph MBR1012 was constructed using the set of 537 segmental homozygous RILs (Figure 2). Marker saturation of the HRMP resulted in reducing the target interval to 0.1 cM. After screening parental lines using the 50K chip and GBS, 19 new polymorphisms were identified in the target region of 0.1 cM. The 50K screen revealed in total, a set of 40,777 scoreable SNPs at the barley genome (Figure 3). Out of these, 14,616 SNPs showed homozygous polymorphisms between resistant and susceptible genotypes. Thirty-nine SNPs were located at the large interval of 8.0 cM on chromosome 1HS, and four SNPs were located at the closest target interval comprising 0.1 cM. These SNPs were converted into KASP markers and mapped on the whole HRMP population (Supplementary Table S1).
Genotyping by sequencing analysis yielded 48,226 SNPs distributed over all seven barley chromosome, of which 37,287 showed homozygous polymorphisms between resistant and susceptible lines (Figure 3). Out of these, 80 polymorphic markers were located in the larger interval, flanked by QBS94 and QBS113 (8.0 cM) and 15 SNPs were identified in the shortened interval of 0.1 cM. KASP markers were designed for all 15 SNPs and used for genotyping of the 537 RILs (Supplementary Table S1).
Mapping of all mentioned markers showed that the Rph MBR1012 locus is located in a region of 0.07 cM between tightly linked markers QBS127 (SNP) and QBS98 (size polymorphism) at 0.020% (distal) and 0.050% (proximal) recombination of the Rph MBR1012 locus. Thus, the target interval was shortened from 0.1% recombination to 0.07% recombination (Figure 4). A high-density genetic map revealed ten markers co-segregating (QBS128, QBS129, QBS130, QBS131, QBS132, GBS626, GBR534, GBS546, QBS116 and QBS117) within the Rph MBR1012 locus (Figure 4). Moreover, recombination distribution in the target interval was uneven varying from 0.58 to 0.60 Mb/cM proximally and distally, respectively to the resistance locus, to 7.26 Mb/cM at the Rph MBR1012 locus (Figure 2). Marker saturation also revealed a high number of recombination between markers QBS96 and QBS71 in the distal region of the interval, i.e., 177 recombination events and 112 recombinations between markers QBS112 and QBS113 located proximally (Figure 5). However, the analysis allowed narrowing the Rph MBR1012 locus to a region comprising a limited number of candidate genes.
BLAST searches against the barley reference sequence revealed that the mapped markers were in a nearly perfect co-linear order. However, 15 markers within 1.34 Mb in the distal part of chromosome 1HS showed a marker inversion (Figure 4). The BLAST searches also indicated only one hit for 20 markers (12 SNPs, 4 SSRs and 4 size polymorphism) on chromosome 1H and two or more hits for 17 markers (12 SNPs, 2 SSRs and 3 size polymorphism). The physical size of the large target interval of 8.0 cM between the flanking markers QBS94 and QBS113 encompassed 6.24 Mb. This region harbors 299 genes of which 183 are high confidence (HC) genes and 116 are low confidence genes. Based on the sequence annotation of HC and LC genes, 23 genes were disease resistance proteins and three were annotated as powdery mildew resistance proteins (Supplementary Table S2). Likewise, physical size of the shortened interval carrying Rph MBR1012 flanked between QBS127 and QBS98 was estimated to 0.44 Mb (Figure 4). In this interval 11 low confidence and 18 high confidence (HC) genes were detected (Supplementary Table S2). Fifteen of these genes are functionally annotated and five of them are related to pathogen resistance, i.e., HORVU1Hr1G000830 (disease resistance protein), HORVU1Hr1G000840 (powdery mildew resistance protein PM3 variant), HORVU1Hr1G000860 (disease resistance protein), HORVU1Hr1G000900 (disease resistance protein) and HORVU1Hr1G000910 (disease resistance protein) (Supplementary Table S2). The markers QBS128 and QBS130 are exactly located at two disease resistance genes, namely HORVU1Hr1G000830 and HORVU1Hr1G000910.
Furthermore, the available barley annotation (Mascher et al., 2017) revealed a mosaic structure of exon and intron fragments only for two disease resistance genes, namely HORVU1Hr1G000830 and HORVU1Hr1G000860, while the three other disease resistance genes (HORVU1Hr1G000840, HORVU1Hr1G000900 and HORVU1Hr1G000910) only have one coding exon (Figure 6).

Testing Diagnostic Value of Developed Markers
Diagnostic assessment of markers co-segregating markers with the Rph MBR1012 was conducted. However, out of ten tested markers only six showed clear a allele differentiation, whereas for four markers, i.e., QBS129, QBS130, QBS131 and QBS132 had to be excluded. The number of alleles detected varied from two alleles for markers QBS116, QBS117, QBS128, QBS130, GBS546, GBS626 and seven alleles for GBR534. For the two markers GBS546 and GBS626 most of the cultivars/lines showed the same allele as the susceptible parental line Scarlett with 80.32 and 83.60% accuracy, respectively. Marker QBS117 with 9.8% accuracy for Rph MBR1012 has no diagnostic value to trigger this gene. Other tested markers were also of limited value for markerassisted selection ( Table 3).

Allele Specific Re-sequencing of Candidate Genes
Allele specific re-sequencing for all 29 putative genes located on the pseudomolecule of chromosome 1H from 2,206,515 to 2,763,382 bp located in a narrowed interval comprising 0.44 Mb was conducted twice. In the first round of low pass resequencing, a set of 36 primer pairs were designed, 33 primer pairs amplified products in both parental lines, one was dominant by amplifying products in Scarlett and two were dominant for MBR1012 and did not produce any fragment on Scarlett. For two genes no specific primer on chromosome 1HS could be designed due to the high similarity of the sequences of these genes (e.g., gene HORVU1Hr1G000820.1: on chromosome 4H, 1863 bp of 1866 bp identical to chromosome 1H). Out of 36 primer pairs, 24 primer pairs were functional, while 12 primer pairs were not functional, since PCR products gave multiple bands, smear or present/absent patterns. Finally, 24 PCR amplicons of the functional primer pairs were sequenced. Moreover, markers for which polymorphisms were based on size polymorphism of polymerase chain reaction (PCR) fragments between parental lines (HORVU1Hr1G000910.9_s3958_as4143 and HORVU1Hr1G001060.1_s173_as480) were directly mapped into the HRMP population. By editing the sequence data, sequence of 18 amplicons could be aligned in both parental lines while for six fragments no alignments were achieved due to the low quality of the sequence data or obtained heterozygous signals (Supplementary Table S3).
Next, for whole length amplification and re-sequencing of five disease resistance genes in the target interval, 25 new primer pairs were designed (Supplementary Table S4). Out of 25 designed primers, 23 amplified products in both parental lines. From this experiment, 12 PCR products were sequenced (Supplementary Table S3). Finally, for 31,204 bp of all 29 candidate genes 61 primer pairs were designed, yielding DNA sequence information for 17,107 bp in MBR1012 and 16,963 bp in Scarlett. Using this sequence data, 259 SNPs were identified for disease resistance genes from the target interval. Moreover, from gene HORVU1Hr1G000900.5 (Disease resistance protein) a large deletion (InDel) was identified in Scarlett ranging from 26 to 222 bp. Seven SNPs for HORVU1Hr1G000830.3, nine for HORVU1Hr1G000860.7 and 243 for HORVU1Hr1G000900.5 were identified (Supplementary Table S3). For two resistance genes i.e., HORVU1Hr1G000840.1 and HORVU1Hr1G000910.9 no SNP/InDel were identified. Functional annotation of defined  SNPs between parental lines, MBR1012 and Scarlett, revealed synonymous mutations for 11 SNPs whereas for 17 SNPs amino acid substitutions were detected. For two SNPs the arginine amino acid changed to a stop codon (TGA) ( Table 4). Multiple alignment also revealed polymorphisms between the parents and barley reference sequence (Supplementary Table S5).

DISCUSSION
Leaf rust is an important fungal disease affecting barley production (Park, 2003). Fungicide application is an option to reduce yield losses but is not always efficient and cannot be considered as a sustainable disease management . Thus, growing of resistant cultivars is the most economical and environmental friendly way to reduce yield losses caused by leaf rust (Kolmer, 1996). However, disease resistance provided by major Rph genes is often overcome due to the emergence of new P. hordei pathotypes (Niks, 1982;Steffenson et al., 1993;Park, 2003) indicating the need for introducing new sources of resistance into barley breeding as well as the need for isolating known ones toward deciphering the structure and function offering the possibility of developing functional markers for breeding and create new alleles by e.g., CRISPR/Cas9 (Kumar et al., 2018).
In this study we have shown the efficient use of the barley reference sequence in physical mapping and especially in marker saturation. Previously, Perovic et al. (2003) demonstrated that the barley landrace MBR1012 is resistant to the barley leaf rust isolate I-80, which later was mapped using 14 SSRs and three SNPs markers on barley chromosome 1HS (König et al., 2012). A null allele of the SSR marker GBMS187 was identified as the closest linked marker at 0.8 cM proximal to the resistance gene. The allelic status of Rph MBR1012 andRph4 (McDaniel andHathcock, 1969), two genes mapped on the short arm on barley chromosome 1HS, is part of an ongoing experiment (Perovic et al., in preparation). The phenotypic evaluation conducted here revealed a hypersensitive reaction of the Rph MBR1012 resistance gene (Figure 1), while the genetic analysis demonstrates that by using genetically mapped markers in combination with the genome sequence information (Mascher et al., 2017) the physical position of this locus can be determined easily. An initial size of the locus of 6.25 Mb that was estimated based on the published map was further downsized by the use of new marker resources and by increasing the genetic resolution.
For many years, mapping of resistance genes relied on the use of various molecular markers i.e., restriction fragment length polymorphism (RFLP) (Graner et al., 1991;Kleinhofs et al., 1993), random amplified polymorphic DNAs (RAPDs) (Williams et al., 1990;Chalmers et al., 1993), amplified fragment length polymorphism (AFLPs) (Vos et al., 1995;Qi et al., 1998) and SSRs (Ramsay et al., 2000;Varshney et al., 2007). For instance, the powdery mildew resistance gene mlo was identified by a combined use of RFLP and AFLP markers which was the first gene, isolated by map-based cloning in barley (Büschges et al., 1997). AFLP, RAPD and RFLP-derived markers were also used to saturate the Mla region (Wei et al., 1999). However, using these marker systems, gene isolation was a laborious and time consuming effort. Advances in molecular marker technologies as well as the previous version of the barley genome sequence already facilitated an accelerated fine mapping of disease resistance genes (Lüpken et al., 2013(Lüpken et al., , 2014Yang et al., 2014). New Illumina SNP genotyping assays, namely 9K and 50K (Comadran et al., 2012;Bayer et al., 2017), together with GBS (Poland et al., 2012) opened a new way for a more efficient and faster marker saturation of target loci in barley. In our study, above mentioned marker resources were used for a first marker saturation of Rph MBR1012 . During the simultaneous construction of a fine map and an initial marker saturation a set of 37 GZ and 9K iSelect SNP markers were randomly selected and mapped to our target interval of 8.0 cM reducing the target interval to 0.1 cM flanked by QBS97 and QBS98. Subsequently, the newly developed high-density barley 50K Infinium SNP markers (Bayer et al., 2017) and GBS markers, which were selected using the reference sequence in the shortened candidate interval (0.1 cM), resulted in the identification of nineteen additional polymorphic SNPs. These markers were converted into KASP markers and the Rph MBR1012 locus was genetically further narrowed into an interval of 0.07 cM between the markers QBS127 and QBS98. In the target interval, ten markers i.e., QBS128, QBS129, QBS130, GBS626, GBR534, GBS546, QBS116, QBS117, QBS131 and QBS132 spanning 0.07 cM genetic distance between QBS127 (at 0.02 cM) and QB98 (at 0.05 cM) were cosegregating. Seven out of the ten co-segregating markers, namely QBS116 (50K), QBS117 (50K), QBS128 (GBS), QBS130 (GBS), QBS131 (GBS), GBS546 and GBR534, were located in five genes in the target interval. Fine mapping of resistance genes is a bottleneck in gene isolation due to the presence of many genes within target intervals, an uneven recombination frequency and a lack of molecular markers, (Stein and Graner, 2005). The fine map for the Rph MBR1012 region constructed in this study was based on a set of 56 molecular markers including four InDel, three present/absent, six SSRs, and 43 SNPs markers. Even though Rph MBR1012 is located in the telomeric region, it turned out that recombination events are not distributed continuously along this region. Although Rph MBR1012 is surrounded with two highly recombining regions at the telomere of chromosome 1HS, 0.58 and 0.6 Mb/cM, the locus is in very unfavorable region of 7.28 Mb/cM with a high number of co-segregating markers, again elucidating that the potential of map-based cloning still depends on the genomic context around the gene of interest. Uneven distribution of recombination frequencies along the genome Akhunov et al., 2003) and differences in local recombination rates, may cause regions even without any recombination over large physical distances which are not suited to map based cloning (Qi and Gill, 2001;Neu et al., 2002). Consequently, the efficiency of an effort of increasing the population size has always to be considered.
Genome-wide studies and multiple gene surveys recorded variation of the SNP frequencies in barley from one SNP per 240 bp, per 200 bp, and per 189 bp (Shavrukov, 2016). In contrast, one SNP per 7 bp in the leaf rust resistance Rph7 gene region evidently showed the usefulness of high-density SNP markers for the purpose of gene isolation in barley (Scherrer et al., 2005).
In addition to the resources used in our study the following genomic resources for marker saturation nowadays may be used: exome sequencing (Mascher et al., 2013b), RNA sequencing (RNAseq) (Wang et al., 2009) which is based on transcriptome profiling, resistance gene enrichment sequencing (RenSeq) (Andolfo et al., 2014) and WGS. These methods may serve to enhance the detection of polymorphism in the genome and to develop markers toward gene isolation in a short period of time. More recently, MutRenSeq that combines the complexity reduction of R gene targeted enrichment sequencing and computational analysis based on comparative genomics provides a tool for the rapid cloning of disease resistance (R) genes in plants (Steuernagel et al., 2016;Dracatos et al., 2018).
Anchoring of the markers against the barley reference sequence elucidated the physical size of 0.44 Mb for the interval harboring Rph MBR1012 . The order of all mapped markers were inconsistent with the order in the barley physical map . However, a large rearrangement of 15 markers within 1.34 Mb in the distal part of chromosome 1H was observed. This inversion is based on non-fixed orientation of the BAC-based sequence contig within the small scaffold having only one the anchor point (personal communication Martin Mascher).
Twenty-nine annotated genes were identified within the narrowed down interval between markers QBS127 and QBS98 comprising five disease associated resistance genes (R genes) which support the prior observation that many barley resistance genes are located distally in regions with high recombination frequency (International Barley Genome Sequencing Consortium [IBSC], 2012). It has been indicated, that more than 80% of all known R genes are of the NBS-LRR type (nucleotide-binding leucine rich repeat) (Shao et al., 2016). LRR domains have particular function in plant-pathogen recognition (Hong and Zhang, 2016). The annotation using Blastx against the non-redundant protein database of NCBI also indicates the presence of the NBS-LRR domain in all five disease resistance genes in the target interval. Disease resistance genes located in the target interval tend to cluster which is typical for NBS-LRR based resistance gene analogs (DeYoung and Innes, 2006). Since P. hordei is a biotrophic fungi and the fact that NBS-LRR resistance genes are only effective in conferring resistance to biotrophic or hemibiotrophic pathogens, but not against necrotrophic pathogens (Belkhadir et al., 2004) provides evidence that resistance is due to a gene carrying the NBS-LRR motif. Hence, full length re-sequencing of five disease resistance genes in parental lines was conducted. However, more than 80% similarity in the sequences of R genes considerably hampered sequencing, therefore in order to obtain a complete sequence of the disease resistance genes, new primer design will be conducted.
Marker validation of seven co-segregating markers in 51 already tested barley lines (König et al., 2012), as well as 12 other barley cultivars/lines, gave hint that new markers identified in this study are not all diagnostic for Rph MBR1012 . Based on our study, the markers GBS546 and GBS626 with 80.32 and 83.60% accuracies in prediction of Rph MBR1012 are the best diagnostic markers and facilitate faster and easier detection of Rph MBR1012 (and putative alleles) in barley breeding lines. Out of the selected markers QBS128 (HORVU1Hr1G000830/Disease resistance protein), QBS130 (HORVU1Hr1G000910/Disease resistance protein), QBS116, QBS117 and GBR534 (HORVU1Hr1G000940/copper ion binding), and marker GBS546 (HORVU1Hr1G000930/Low molecular weight glutenin subunit) were directly derived from putative candidate genes in the target interval but revealed a less diagnostic character. However, the diagnostic Rph MBR1012 markers identified in this study could be very useful not only for discriminating between resistant and susceptible cultivars but also for pyramiding of Rph MBR1012 with other resistance genes to aim a durable resistance in barley cultivars (Sharma Poudel et al., 2018).

CONCLUSION
In summary, by using high-throughput genotyping and sequencing techniques together with the barley reference sequence we succeeded to downsize the Rph MBR1012 target interval to 0.44 Mb between markers QBS127 and QBS98 in comparison to 6.24 Mb in a previous study. This is an indispensable step toward isolation of this gene. Four strategies might be then considered in next step in order to define the loci underlying the resistance gene Rph MBR1012 ; enhancing the map resolution via screening a new set of F 2 plants and using the new SNPs and InDel defined from candidate genes at target interval to develop the new markers for further marker saturation, screening a non-gridded BAC library from donor line MBR1012, overexpression of five detected disease resistance genes in the target interval in a susceptible barley cultivar, e.g., Scarlett and knock out the genes in resistant lines using CRISPR/cas9. The co-segregating and closely linked markers detected in this study, may be useful as probes for BAC library screening and construction of the physical map in MBR1012.

AUTHOR CONTRIBUTIONS
DP, DK, and FO conceived and designed the experiments, provided the experimental material and contributed to study design, subject recruitment and sample preparation. LF, DK, and DP performed the experiments. LF, JK, and DP analyzed the data. LF, JK, HD, FO, and DP interpreted the data. All authors wrote the manuscript, and read and approved the final manuscript.

FUNDING
This research financially was supported by grants from the Federal Ministry of Education and Research, Bundesministerium für Bildung und Forschung (BMBF) and the Deutscher Akademischer Austauschdienst e. V. (DAAD) to LF.