ORIGINAL RESEARCH article
Sec. Plant Breeding
High Resolution Mapping of RphMBR1012 Conferring Resistance to Puccinia hordei in Barley (Hordeum vulgare L.)
- 1Institute for Resistance Research and Stress Tolerance, Federal Research Centre for Cultivated Plants, Julius Kühn-Institute (JKI), Quedlinburg, Germany
- 2Department of Phytopathology and Plant Protection, Institute of Agricultural and Nutrition Sciences, Faculty of Natural Sciences III, Martin Luther University of Halle-Wittenberg, Halle, Germany
- 3Institute for Biosafety in Plant Biotechnology, Federal Research Centre for Cultivated Plants, Julius Kühn-Institute (JKI), Quedlinburg, Germany
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.
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% (Park et al., 2015), 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).
Since the first genetic study on leaf rust resistance (Waterhouse, 1927), 25 Rph (Resistance to P. hordei) genes have been mapped in barley (Kavanagh et al., 2017). Among them, two genes, namely Rph20 and Rph23, mediate an adult plant resistance (APR) (Hickey et al., 2011; Singh et al., 2015), while the remaining 23 (Rph1 to Rph19, Rph21, Rph22, Rph24, and Rph25) establish seedling resistance (Kavanagh et al., 2017). Rph5 and Rph6 on chromosome 3H (Zhong et al., 2003), Rph9 and Rph12 on chromosome 5H (Borovkova et al., 1998) and Rph15 and Rph16 on chromosome 2H have been described as alleles of the same gene (Weerasena et al., 2004). Only Rph7, Rph15, and Rph16 are still effective in Europe (Niks et al., 2000; Perovic et al., 2004) and the number of effective Rph genes available to breeders is decreasing rapidly (Kavanagh et al., 2017). Among all known Rph genes, only Rph1 has been isolated recently, using the newly developed cloning approach called Mutant Chromosome Sequencing (MutChromSeq) (Steuernagel et al., 2016) in combination with genetic mapping (Dracatos et al., 2018).
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 (Künzel et al., 2000) often hampers high-resolution genetic dissection. In barley, peri-centromeric 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 RphMBR1012 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 RphMBR1012 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.
Materials and Methods
Plant Material and Construction of a High-Resolution Mapping Population
For high resolution mapping of RphMBR1012, a segregating population comprising of 4,775 F2 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, F2 plants were analyzed using two flanking co-dominant SSRs, i.e., QBS94 (distal) and QBS113 (proximal) (Perovic et al., 2013). Respective markers were analyzed by capillary electrophoresis at the genetic analyzer ABI PRISM® 3100 (Applied Biosystems, Darmstadt, Germany). From identified heterozygous recombinant F2 plants in target interval, 12 progeny plants, representing F3 families were sown in 96 Quick pot plates. Genomic DNA of 10 days old plantlets was extracted in F2 and F3 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, F3 recombinant plants were self-fertilized and as F4 segmental RILs used for phenotyping and genotyping with newly developed PCR based markers.
Table 1. DH lines and crosses used for the construction of the high resolution mapping population for RphMBR1012.
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 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.
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.01 (Koressaar and Remm, 2007; Untergasser et al., 2012). KASP markers were then used for genotyping of the HRMP.
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® MiSeqTM (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 region2.
The HRMP was genotyped using in total 56 molecular markers derived from the procedures described above. Molecular markers used may be divided in five types as follows: six SSRs based markers from the pyrosequencing assay (Silvar et al., 2011), three dominant present/absent markers, four size polymorphism markers [insertion/deletion polymorphisms (InDels)], 19 KASP markers and 24 Cleaved Amplified Polymorphic Sequences (CAPS) markers. Size polymorphisms markers and SSRs were amplified in a total volume of 10 μl, according to Perovic et al. (2013) and detected either using fluorescently labeled primers (M13) by capillary electrophoresis on the ABI Genetic Analyzer (ABI sequencer, ABI Perkin Elmer, Weiterstadt, Germany), or directly separated on a 1.5% agarose gel. For ABI analysis, 0.1 μl of M13 primer (10.0 pmol/μl)(5′-CACGACGTTGTAAAACGAC-3′) labeled with fluorescent dye was added to the reaction mix. One microliter of diluted PCR product was added to 14 μl of HiDi-Rox mastermix (1.4 ml Hidi and 6 μl Rox) in a total volume of 15 μl. Results were analyzed using the software package GeneMapper v4.0 (Applied Biosystems, Darmstadt, Germany). For 43 sequences, detected SNPs were converted either to KASP markers (see footnote 2) or CAPS markers using NEB cutter v.2.03. KASP reaction was performed in total volume of 5 μl containing 2.5 μl KASP mix (LGC Genomics GmbH, Germany), 0.08 μl forward primer, allele 1 (10.0 pmol/μl, labeled with FAM M13 tail), 0.08 μl forward primer allele 2 (10.0 pmol/μl, labeled with HEX M13 tail), 0.2 μl reverse common primer (10.0 pmol/μl), and 2.2 μl template DNA (20 ng/μl). For CAPS analysis, DNA amplicons were cleaved with the respective restriction endonuclease (Table 2) in a volume of 20 μl, containing 2 μl corresponding 10× buffer, 0.1 μl appropriate enzyme, 7.9 μl HPLC gradient grade water (Carl Roth, Karlsruhe, Germany) and 8/10 μl of the PCR product. Proper temperature was applied according to manufacturer’s instructions for each restriction endonuclease and digestion was done for 3 h.
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 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) software package.
Testing the Diagnostic Value of Co-segregating and Closely Linked Markers
Co-segregating markers in RphMBR1012 locus were tested for their diagnostic value on a set of 63 genotypes comprising 25 selected barley genotypes/lines carrying Rph1 to Rph25, 23 parental lines and 15 Bowman introgression lines carrying Rph1 to Rph15 (Table 3). The diagnostic value of tested co-segregating markers (%) was calculated using the following equation:
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 sequence4 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 genes5. 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. A fragment size of 400 to 1,200 bp was chosen because of the maximum sequencing length. Amplification was done on the parental genotypes MBR1012 and Scarlett, as well as on two DH-lines [DH3/62 (S), DH3/127 (R)]. Amplification reaction was prepared in a total volume of 20 μl containing 2 μl of 10× PCR buffer (Qiagen, Hilden, Germany), 2 μl of 25 mM MgCl2, 0.4 μl of 10 mM dNTPs (Fermentas, Schwerte, Germany), 0.5 μl of each forward and reverse primer (10.0 pmol/μl), 0.16 unit of fire DNA polymerase (5 U/μl), (Qiagen, Hilden, Germany), 12.44 μl HPLC gradient grade water (Carl Roth, Karlsruhe, Germany) and 2 μl of template DNA (20 ng/μl). 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 DoceTM XR and the Quantity One® 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).
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.
Construction of the High-Resolution Mapping Population
Four crosses i.e., DH 3/74 (S) × DH3/6 (R), DH3/74 (S) × DH3/127 (R), DH3/6 (R) × DH3/9 (S) and DH3/62 (S) × DH3/127 (R), were used for the construction of the HRMP (Table 1). In total, of 5,237 F2 plants, 4,775 survived, and from corresponding F3 families 537 recombinant F4 RILs were developed, resulting in an interval harboring the resistance locus of 0.07% recombination. Finally, a genetic resolution of 0.010% recombination was achieved.
Phenotypic analysis of resistance to RphMBR1012 showed a segregation of 261 resistant and 276 susceptible RILs and revealed the expected 1r:1s segregation ratio among these RILs. Chi-square 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).
Figure 1. (A) Macroscopic symptoms of Puccinia hordei isolate I-80 on the resistant parent (MBR1012), the susceptible parent (Scarlett), two DH lines and the F1 7 days post-inoculation. (B) CAPS marker GBS546 originated from a high-confidence gene in the target interval.
Marker Saturation of the RphMBR1012 Locus and Anchoring to the Barley Reference Sequence
A fine map of the RphMBR1012 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.
Figure 2. High-resolution genetic map of RphMBR1012. (A) The genomic region harboring the RphMBR1012. (B) The identification of 537 recombinants and mapping of RphMBR1012 locus based on markers derived from the Genome Zipper and the 9K iSelect chip flanked by QBS97 and QBS98. The blue boxes indicate the physical size. (C) Target region used for marker saturation based on the 50K Infinium chip and GBS markers. Markers derived from the 50K Infinium chip are highlighted with orange and markers derived from GBS are shown in green. Co-segregating markers are indicated by italics. (D) Genetic map of RphMBR1012 locus.
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).
Figure 3. Landscape of the 50K and GBS marker distribution. Track A: gives the seven barley chromosomes. Track B: grey color depicts GBS: (all, 48.226), blue, position of 37,287 polymorphic SNPs between MBR1012 and Scarlett. Track C: distribution of SNP Chip (50K) markers: grey (all, 40,777), green, position of 14,616 polymorphic SNPs between MBR1012 and Scarlett.
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 RphMBR1012 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 RphMBR1012 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 RphMBR1012 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 RphMBR1012 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 RphMBR1012 locus to a region comprising a limited number of candidate genes.
Figure 4. High-density genetic and physical map of the RphMBR1012 region on barley chromosome 1HS based on 56 molecular markers and 537 recombinant inbred lines derived from the cross MBR1012 × Scarlett.
Figure 5. Graphical genotypes of F4 RILs for all 537 recombinant lines carrying cross-over events between QBS94 and QBS113 (8.0 cM). A – (Susceptible genotype = white) and B – (Resistant genotype = hatched) in the target locus indicate the result of the resistance test of recombinant lines. Border of hatched to white shows the recombination position between the MBR1012 allele to the Scarlett allele and white to hatched shows the recombination position between the Scarlett allele to the MBR1012 allele.
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 RphMBR1012 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).
Figure 6. Gene structure of five disease resistance genes positioned in target interval. Colored boxes in genes indicate CDSs (exon) and thin lines indicate introns.
Testing Diagnostic Value of Developed Markers
Diagnostic assessment of markers co-segregating markers with the RphMBR1012 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 RphMBR1012 has no diagnostic value to trigger this gene. Other tested markers were also of limited value for marker-assisted 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).
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 (Park et al., 2015). 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 RphMBR1012 and Rph4 (McDaniel and Hathcock, 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 RphMBR1012 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.
Table 4. Functional annotation of SNPs between parental lines (MBR1012 and Scarlett) originated from candidate genes located within the 0.44 Mb of target interval.
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, 2014; Yang 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 RphMBR1012. 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 RphMBR1012 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 co-segregating. 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 RphMBR1012 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 RphMBR1012 is located in the telomeric region, it turned out that recombination events are not distributed continuously along this region. Although RphMBR1012 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 (Künzel et al., 2000; 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 RphMBR1012. The order of all mapped markers were inconsistent with the order in the barley physical map (Stein et al., 2007). 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 RphMBR1012. Based on our study, the markers GBS546 and GBS626 with 80.32 and 83.60% accuracies in prediction of RphMBR1012 are the best diagnostic markers and facilitate faster and easier detection of RphMBR1012 (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 RphMBR1012 markers identified in this study could be very useful not only for discriminating between resistant and susceptible cultivars but also for pyramiding of RphMBR1012 with other resistance genes to aim a durable resistance in barley cultivars (Sharma Poudel et al., 2018).
In summary, by using high-throughput genotyping and sequencing techniques together with the barley reference sequence we succeeded to downsize the RphMBR1012 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 RphMBR1012; enhancing the map resolution via screening a new set of F2 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.
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.
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.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors acknowledge the excellent technical assistance of Katy Niedung and Marlis Weilepp.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.00640/full#supplementary-material
- ^ http://bioinfo.ut.ee/primer3-0.4.0/
- ^ http://www.lgcgroup.com/
- ^ http://tools.neb.com/NEBcutter2
- ^ http://webblast.ipk-gatersleben.de/barley_ibsc/
- ^ https://plants.ensembl.org/Hordeum_vulgare/Info/Index
- ^ https://www.ncbi.nlm.nih.gov/sutils/splign/splign.cgi?textpage=online&level=form
Akhunov, E. D., Goodyear, A. W., Geng, S., Qi, L. L., Echalier, B., and Gill, B. S. (2003). The organization and rate of evolution of wheat genomes are correlated with recombination rates along chromosome arms. Genome Res. 13, 753–763.
Andolfo, G., Jupe, F., Witek, K., Etherington, G. J., Ercolano, M. R., and Jones, J. D. G. (2014). Defining the full tomato NB-LRR resistance gene repertoire using genomic and cDNA RenSeq. BMC Plant Biol. 14:120. doi: 10.1186/1471-2229-14-120
Bayer, M. M., Rapazote-Flores, P., Ganal, M., Hedley, P. E., Macaulay, M., Plieske, J., et al. (2017). Development and evaluation of a Barley 50k iSelect SNP Array. Front. Plant Sci. 8:1792. doi: 10.3389/fpls.2017.01792
Bhakta, M. S., Jones, V. A., and Vallejos, C. E. (2015). Punctuated distribution of recombination hotspots and demarcation of pericentromeric regions in Phaseolus vulgaris L. PLoS One 10:e0116822. doi: 10.1371/journal.pone.0116822
Blankenberg, D., Von Kuster, G., Coraor, N., Ananda, G., Lazarus, R., and Mangan, M. (2001). Galaxy: aweb-based genome analysis UNIT tool for experimentalists. Curr. Protoc. Mol. Biol. 89, 19.10.1–19.10.21. doi: 10.1002/0471142727.mb1910s89
Borovkova, I. G., Jin, Y., and Steffenson, B. J. (1998). Chromosomal location and genetic relationship of leaf rust resistance genes Rph9 and Rph12 in barley. Phytopathology 88, 76–80. doi: 10.1094/PHYTO.19126.96.36.199
Brueggeman, R., Rostoks, N., Kudrna, D., Kilian, A., Han, F., and Chen, J. (2002). The barley stem rust-resistance gene Rpg1 is a novel disease-resistance gene with homology to receptor kinases. Proc. Natl. Acad. Sci. U.S.A. 99, 9328–9333.
Chalmers, K. J., Barua, U. M., Hackett, C. A., Thomas, R., and Waugh, W. P. (1993). Identification of RAPD markers linked to genetic factors controlling the milling energy requirement of barley. Theor. Appl. Genet. 87, 314–320. doi: 10.1007/BF01184917
Cingolani, P., Platts, A., Wang le, L., Coon, M., Nguyen, T., Wang, L., et al. (2012). A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 6, 80–92. doi: 10.4161/fly.19695
Close, T. J., Bhat, P. R., Lonardi, S., Wu, Y. H., Rostoks, N., and Ramsay, L. (2009). Development and implementation of high-throughput SNP genotyping in barley. BMC Genom. 10:582. doi: 10.1186/1471-2164-10-582
Comadran, J., Kilian, B., Russell, J., Ramsay, L., Stein, N., and Ganal, M. (2012). Natural variation in a homolog of Antirrhinum CENTRORADIALIS contributed to spring growth habit and environmental adaptation in cultivated barley. Nat. Genet. 44, 1388–1392. doi: 10.1038/ng.2447
Dracatos, P. M., Bartos, J., Elmansour, H., Singh, D., Karafiatova, M., and Zhang, P. (2018). The coiled-coil NLR Rph1, confers leaf rust resistance in barley cultivar Sudan. Plant Physiol. 179, 1362–1372. doi: 10.1104/pp.18.01052
Ganal, M. W., Plieske, J., Hohmeyer, A., Polley, A., and Röder, M. S. (2018). “High throughput genotyping for cereal research and breeding,” in Applications of Genetic and Genomic Research in Cereals, eds T. Miedaner and V. Korzun (Sawston: Woodhead Publishing).
Goecks, J., Nekrutenko, A., and Taylor, J. (2010). Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 11:R86. doi: 10.1186/gb-2010-11-8-r86
Halterman, D., Zhou, F., Wei, F., Wise, R. P., and Schulze-Lefert, P. (2001). The MLA6 coiled-coil, NBS-LRR protein confers AvrMla6-dependent resistance specificity to Blumeria graminis f. sp. hordei in barley and wheat. Plant J. 25, 335–348.
Hickey, L. T., Lawson, W., Platz, G. J., Dieters, M., Arief, V. N., German, S., et al. (2011). Mapping Rph20: a gene conferring adult plant resistance to Puccinia hordei in barley. Theor. Appl. Genet. 123, 55–68. doi: 10.1007/s00122-011-1566-z
Jin, Y., Statler, G. D., Franckowiak, J. D., and Steffenson, B. J. (1993). Linkage between leaf rust resistance genes and morphological markers in barley. Phytopathology 83, 230–233. doi: 10.1094/PHYTO.19188.8.131.52
Johnston, P. A., Niks, R. E., Meiyalaghan, V., Blanchet, E., and Pickering, R. (2013). Rph22: mapping of a novel leaf rust resistance gene introgressed from the non-host Hordeum bulbosum L. into cultivated barley (Hordeum vulgare L.). Theor. Appl. Genet. 126, 1613–1625. doi: 10.1007/s00122-013-2078-9
Kavanagh, P. J., Singh, D., Bansal, U. K., and Robert, P. K. (2017). Inheritance and characterization of the new and rare gene Rph25 conferring seedling resistance in Hordeum vulgare against Puccinia hordei. Plant Breed. 136, 908–912.
Kleinhofs, A., Kilian, A., Saghai, Maroof MA, Biyashev, R. M., Hayes, P., Chen, F. Q., et al. (1993). A molecular, isozyme and morphological map of the barley (Hordeum vulgare) genome. Theor. Appl. Genet. 86, 705–712. doi: 10.1007/BF00222660
König, J., Kopahnke, D., Steffenson, B. J., Przulj, N., Romeis, T., Roder, M. S., et al. (2012). Genetic mapping of a leaf rust resistance gene in the former Yugoslavian barley landrace MBR1012. Mol. Breed. 30, 1253–1264.
Kota, R., Varshney, R. K., Prasad, M., Zhang, H., Stein, N., and Graner, A. (2008). EST-derived single nucleotide polymorphism markers for assembling genetic and physical maps of the barley genome. Funct. Integr. Genom. 8, 223–233.
Krattinger, S., Wicker, T., and Keller, B. (2009). “Map-based cloning of genes in triticeae (Wheat and Barley),” in Genetics and Genomics of the Triticeae, eds G. J. Muehlbauer and C. Feuillet (New York, NY: Springer), 337–357.
Kumar, N., Galli, M., Ordon, J., Stuttmann, J., Kogel, K. H., and Imani, J. (2018). Further analysis of barley MORC1 using a highly efficient RNA-guided Cas9 gene-editing system. Plant Biotechnol. J. 16, 1892–1903. doi: 10.1111/pbi.12924
Künzel, G., Korzun, L., and Meister, A. (2000). Cytologically integrated physical restriction fragment length polymorphism maps for the barley genome based on translocation breakpoints. Genetics 154, 397–412.
Li, Y. I., van de Geijn, B., Raj, A., Knowles, D. A., Petti, A. A., Golan, D., et al. (2016). RNA splicing is a primary link between genetic variation and disease. Science 352, 600–604. doi: 10.1126/science.aad9417
Lüpken, T., Stein, N., Perovic, D., Habekuss, A., Kramer, I., and Hahnel, U. (2013). Genomics-based high-resolution mapping of the BaMMV/BaYMV resistance gene rym11 in barley (Hordeum vulgare L.). Theor. Appl. Genet. 126, 1201–1212. doi: 10.1007/s00122-013-2047-3
Mammadov, J. A., Zwonitzer, J. C., Biyashev, R. M., Griffey, C. A., Jin, Y., Steffenson, B. J., et al. (2003). Molecular mapping of leaf rust resistance gene Rph5 in barley. Crop Sci. 43, 388–393. doi: 10.1094/PHYTO.2003.93.5.604
Mascher, M., Gundlach, H., Himmelbach, A., Beier, S., Twardziok, S. O., and Wicker, T. (2017). A chromosome conformation capture ordered sequence of the barley genome. Nature 544, 427–433. doi: 10.1038/nature22043
Mascher, M., Muehlbauer, G. J., Rokhsar, D. S., Chapman, J., Schmutz, J., Barry, K., et al. (2013a). Anchoring and ordering NGS contig assemblies by population sequencing (POPSEQ). Plant J. 76, 718–727. doi: 10.1111/tpj.12319
Mascher, M., Richmond, T. A., Gerhardt, D. J., Himmelbach, A., Clissold, L., and Sampath, D. (2013b). Barley whole exome capture: a tool for genomic research in the genus Hordeum and beyond. Plant J. 76, 494–505. doi: 10.1111/tpj.12294
Mayer, K. F., Martis, M., Hedley, P. E., Simkova, H., Liu, H., and Morris, J. A. (2011). Unlocking the barley genome by chromosomal and comparative genomics. Plant Cell 23, 1249–1263. doi: 10.1105/tpc.110.082537
Muñoz-Amatriaín, M., Moscou, M. J., Bhat, P. R., Svensson, J. T., Bartoš, J., and Suchánková, P. (2011). An improved consensus linkage map of barley based on flow-sorted chromosomes and single nucleotide polymorphism markers. Plant Genome 4, 238–249.
Niks, R. E., Walther, U., Jaiser, H., Martinez, F., Rubiales, D., and Andersen, O. (2000). Resistance against barley leaf rust (Puccinia hordei) in West-European spring barley germplasm. Agronomie 20, 769–782.
Park, R. F., Golegaonkar, P. G., Derevnina, L., Sandhu, K. S., Karaoglu, H., and Elmansour, H. M. (2015). Leaf rust of cultivated barley: pathology and control. Annu. Rev. Phytopathol. 53, 565–589. doi: 10.1146/annurev-phyto-080614-120324
Park, R. F., and Karakousis, A. (2002). Characterization and mapping of gene Rph19 conferring resistance to Puccinia hordei in the cultivar ’Reka 1’ and several Australian barleys. Plant Breed. 121, 232–236.
Pellio, B., Streng, S., Bauer, E., Stein, N., Perovic, D., and Schiemann, A. (2005). High-resolution mapping of the Rym4/Rym5 locus conferring resistance to the barley yellow mosaic virus complex (BaMMV, BaYMV, BaYMV-2) in barley (Hordeum vulgare ssp. vulgare L.). Theor. Appl. Genet. 110, 283–293.
Perovic, D., Doris, K., Brian, S. J., Jutta, F., Janine, K., and Benjamin, K. (2012). “Genetic fine mapping of a novel leaf rust resistance gene and a barley yellow dwarf virus tolerance (BYDV) introgressed from hordeum bulbosum by the use of the 9K iSelect Chip,” in Advance in Barley Sciences, eds G. Zhang, C. Li, and X. Liu (Dordrecht: Springer).
Perovic, D., Kopahnke, D., Habeku, A. ß, Ordon, F., and Serfling, A. (2018). “Marker-based harnessing of genetic diversity to improve resistance of barley to fungal and viral diseases,” in Applications of genetic and genomic research in cereals, eds T. Miedaner and V. Korzun (Amsterdam: Elsevier), 137–164.
Perovic, D., Przulj, N., Milovanovic, M., Prodanovic, S., Perovic, J., Kophanke, D., et al. (2003). “Characterisation of spring barley genetic resources in Yugoslavia,” in Proceedings of the Symposium Dedicated to the 100th Birthday of Rudolf Mansfeld, Vol. 22, (Gatersleben), 301–306.
Perovic, D., Stein, N., Zhang, H., Drescher, A., Prasad, M., Kota, R., et al. (2004). An integrated approach for comparative mapping in rice and barley with special reference to the Rph16 resistance locus. Funct. Integr. Genomics 4, 74–83.
Perovic, J., Silvar, C., Koenig, J., Stein, N., Perovic, D., and Ordon, F. (2013). A versatile fluorescence-based multiplexing assay for CAPS genotyping on capillary electrophoresis systems. Mol. Breed. 32, 61–69.
Pickering, R. A., Malyshev, S., Künzel, G., Johnston, P. A., Korzun, V., Menke, M., et al. (2000). Locating introgressions of Hordeum bulbosum chromatin within the H. vulgare genome. Theor. Appl. Genet. 100, 27–31.
Pickering, R. A., SteVenson, B. J., Hill, A. M., and Borovkova, I. (1998). Association of leaf rust and powdery mildew resistance in a recombinant derived from a Hordeum vulgare £ Hordeum bulbosum hybrid. Plant Breed. 117, 83–84.
Poland, J. A., Brown, P. J., Sorrells, M. E., and Jannink, J. L. (2012). Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach. PLoS One 7:e32253. doi: 10.1371/journal.pone.0032253
Qi, L. L., and Gill, B. S. (2001). High-density physical maps reveal that the dominant male-sterile gene Ms3 is located in a genomic region of low recombination in wheat and is not amenable to map-based cloning. Theor. Appl. Genet. 103, 998–1006.
Roane, C. W., and Starling, T. M. (1967). Inheritance of reaction to Puccinia hordei in barley. II. Gene symbols for loci in differential cultivars. Phytopathology 57, 66–68. doi: 10.1094/PHYTO-95-0884
Rostoks, N., Mudie, S., Cardle, L., Russell, J., Ramsay, L., Booth, A., et al. (2005). Genome-wide SNP discovery and linkage analysis in barley based on genes responsive to abiotic stress. Mol. Genet. Genomics 274, 515–527.
Sandhu, K. S., Forrest, K. L., Kong, S., Bansal, U. K., Singh, D., Hayden, M. J., et al. (2012). Inheritance and molecular mapping of a gene conferring seedling resistance against Puccinia hordei in the barley cultivar Ricardo. Theor. Appl. Genet. 125, 1403–1411. doi: 10.1007/s00122-012-1921-8
Scherrer, B., Isidore, E., Klein, P., Kim, J. S., Bellec, A., Chalhoub, B., et al. (2005). Large intraspecific haplotype variability at the Rph7 locus results from rapid and recent divergence in the barley genome. Plant Cell 17, 361–374.
Sharma Poudel, R., Al-Hashel, A. F., Gross, T., Gross, P., and Brueggeman, R. (2018). Pyramiding rpg4- and Rpg1-mediated stem rust resistance in barley requires the Rrr1 gene for both to function. Front. Plant Sci. 9:1789. doi: 10.3389/fpls.2018.01789
Simons, G., van der Lee, T., Diergaarde, P., van Daelen, R., Groenendijk, J., Frijters, A., et al. (1997). AFLP-based fine mapping of the Mlo gene to a 30-kb DNA segment of the barley genome. Genomics 44, 61–70.
Singh, D., Dracatos, P., Derevnina, L., Zhou, M., and Park, R. F. (2015). Rph23: a new designated additive adult plant resistance gene to leaf rust in barley on chromosome 7H. Plant Breed. 134, 62–69.
Steffenson, B. J., Jin, Y., and Griffey, C. A. (1993). Pathotypes of Puccinia-Hordei with virulence for the barley leaf rust resistance gene Rph7 in the United-States. Plant Dis. 77, 867–869. doi: 10.1094/PDIS-08-12-0785-PDN
Stein, N., Prasad, M., Scholz, U., Thiel, T., Zhang, H., Wolf, M., et al. (2007). A 1,000-loci transcript map of the barley genome: new anchoring points for integrative grass genomics. Theor. Appl. Genet. 114, 823–839.
Steuernagel, B., Periyannan, S. K., Hernández-Pinzón, I., Witek, K., Rouse, M. N., and Yu, G. (2016). Rapid cloning of disease-resistance genes in plants using mutagenesis and sequence capture. Nat. Biotechnol. 34, 652–655. doi: 10.1038/nbt.3543
Wang, Y., Cheng, X., Shan, Q., Zhang, Y., Liu, J., Gao, C., et al. (2014). Simultaneous editing of three homoeoalleles in hexaploid bread wheat confers heritable resistance to powdery mildew. Nat. Biotechnol. 32, 947–951. doi: 10.1038/nbt.2969
Waterhouse, W. L. (1927). “Studies in the inheritance of resistance to leaf rust Puccinia anomal Rostr. In crosses of barley,” in Proceedings of the Royal society of New South Wale, (South Wale), 218–247.
Weerasena, J. S., Steffenson, B. J., and Falk, A. B. (2004). Conversion of an amplified fragment length polymorphism marker into a co-dominant marker in the mapping of the Rph15 gene conferring resistance to barley leaf rust, Puccinia hordei Otth. Theor. Appl. Genet. 108, 712–719.
Wei, F., Gobelman-Werner, K., Morroll, S. M., Kurth, J., Mao, L., Wing, R., et al. (1999). The Mla (powdery mildew) resistance cluster is associated with three NBS-LRR gene families and suppressed recombination within a 240-kb DNA interval on chromosome 5S (1HS) of barley. Genetics 153, 1929–1948.
Wendler, N., Mascher, M., Noh, C., Himmelbach, A., Scholz, U., Ruge-Wehling, B., et al. (2014). Unlocking the secondary gene-pool of barley with next-generation sequencing. Plant Biotechnol. J. 12, 1122–1131. doi: 10.1111/pbi.12219
Whelan, H. G., Gaunt, R. E., and Scott, W. R. (1997). The effect of leaf rust (Puccinia hordei) on yield response in barley (Hordeum vulgareL.) crops with different yield potentials. Plant Pathol. 46, 397–406.
Williams, J. G., Kubelik, A. R., Livak, K. J., Rafalski, J. A., and Tingey, S. V. (1990). DNA polymorphisms amplified by arbitrary primers are useful as genetic markers. Nucleic Acids Res. 18, 6531–6535.
Yang, P., Lupken, T., Habekuss, A., Hensel, G., Steuernagel, B., and Kilian, B. (2014). Protein disulfide isomerase like 5-1 is a susceptibility factor to plant viruses. Proc. Natl. Acad. Sci. U.S.A. 111, 2104–2109. doi: 10.1073/pnas.1320362111
Zhong, S. B., Effertz, R. J., Jin, Y., Franckowiak, J. D., and Steffenson, B. J. (2003). Molecular mapping of the leaf rust resistance gene Rph6 in barley and its linkage relationships with Rph5 and Rph7. Phytopathology 93, 604–609. doi: 10.1094/PHYTO.2003.93.5.604
Ziems, L. A., Hickey, L. T., Platz, G. J., Franckowiak, J. D., Dracatos, P. M., Singh, D., et al. (2017). Characterization of Rph24: a gene conferring adult plant resistance to Puccinia hordei in Barley. Phytopathology 107, 834–841. doi: 10.1094/PHYTO-08-16-0295-R
Keywords: barley, leaf rust resistance gene RphMBR1012, positional isolation, GBS, Infinium 50K
Citation: Fazlikhani L, Keilwagen J, Kopahnke D, Deising H, Ordon F and Perovic D (2019) High Resolution Mapping of RphMBR1012 Conferring Resistance to Puccinia hordei in Barley (Hordeum vulgare L.). Front. Plant Sci. 10:640. doi: 10.3389/fpls.2019.00640
Received: 30 January 2019; Accepted: 29 April 2019;
Published: 22 May 2019.
Edited by:Takaki Yamauchi, PRESTO, Japan Science and Technology Agency, Japan
Reviewed by:Kenji Yano, RIKEN Center for Advanced Intelligence Project (AIP), Japan
Urmil Bansal, University of Sydney, Australia
Copyright © 2019 Fazlikhani, Keilwagen, Kopahnke, Deising, Ordon and Perovic. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dragan Perovic, email@example.com