The Barley HvWRKY6 Transcription Factor Is Required for Resistance Against Pyrenophora teres f. teres

Barley is an important cereal crop worldwide because of its use in the brewing and distilling industry. However, adequate supplies of quality malting barley are threatened by global climate change due to drought in some regions and excess precipitation in others, which facilitates epidemics caused by fungal pathogens. The disease net form net blotch caused by the necrotrophic fungal pathogen Pyrenophora teres f. teres (Ptt) has emerged as a global threat to barley production and diverse populations of Ptt have shown a capacity to overcome deployed genetic resistances. The barley line CI5791 exhibits remarkably effective resistance to diverse Ptt isolates from around the world that maps to two major QTL on chromosomes 3H and 6H. To identify genes involved in this effective resistance, CI5791 seed were γ-irradiated and two mutants, designated CI5791-γ3 and CI5791-γ8, with compromised Ptt resistance were identified from an M2 population. Phenotyping of CI5791-γ3 and -γ8 × Heartland F2 populations showed three resistant to one susceptible segregation ratios and CI5791-γ3 × -γ8 F1 individuals were susceptible, thus these independent mutants are in a single allelic gene. Thirty-four homozygous mutant (susceptible) CI5791-γ3 × Heartland F2 individuals, representing 68 recombinant gametes, were genotyped via PCR genotype by sequencing. The data were used for single marker regression mapping placing the mutation on chromosome 3H within an approximate 75 cM interval encompassing the 3H CI5791 resistance QTL. Sequencing of the mutants and wild-type (WT) CI5791 genomic DNA following exome capture identified independent mutations of the HvWRKY6 transcription factor located on chromosome 3H at ∼50.7 cM, within the genetically delimited region. Post transcriptional gene silencing of HvWRKY6 in barley line CI5791 resulted in Ptt susceptibility, confirming that it functions in NFNB resistance, validating it as the gene underlying the mutant phenotypes. Allele analysis and transcript regulation of HvWRKY6 from resistant and susceptible lines revealed sequence identity and upregulation upon pathogen challenge in all genotypes analyzed, suggesting a conserved transcription factor is involved in the defense against the necrotrophic pathogen. We hypothesize that HvWRKY6 functions as a conserved signaling component of defense mechanisms that restricts Ptt growth in barley.

Barley is an important cereal crop worldwide because of its use in the brewing and distilling industry. However, adequate supplies of quality malting barley are threatened by global climate change due to drought in some regions and excess precipitation in others, which facilitates epidemics caused by fungal pathogens. The disease net form net blotch caused by the necrotrophic fungal pathogen Pyrenophora teres f. teres (Ptt) has emerged as a global threat to barley production and diverse populations of Ptt have shown a capacity to overcome deployed genetic resistances. The barley line CI5791 exhibits remarkably effective resistance to diverse Ptt isolates from around the world that maps to two major QTL on chromosomes 3H and 6H. To identify genes involved in this effective resistance, CI5791 seed were γ-irradiated and two mutants, designated CI5791-γ3 and CI5791-γ8, with compromised Ptt resistance were identified from an M 2 population. Phenotyping of CI5791-γ3 and -γ8 × Heartland F 2 populations showed three resistant to one susceptible segregation ratios and CI5791-γ3 × -γ8 F 1 individuals were susceptible, thus these independent mutants are in a single allelic gene. Thirty-four homozygous mutant (susceptible) CI5791-γ3 × Heartland F 2 individuals, representing 68 recombinant gametes, were genotyped via PCR genotype by sequencing. The data were used for single marker regression mapping placing the mutation on chromosome 3H within an approximate 75 cM interval encompassing the 3H CI5791 resistance QTL. Sequencing of the mutants and wild-type (WT) CI5791 genomic DNA following exome capture identified independent mutations of the HvWRKY6 transcription factor located on chromosome 3H at ∼50.7 cM, within the genetically delimited region. Post transcriptional gene silencing of HvWRKY6 in barley line CI5791 resulted in Ptt susceptibility, confirming that it functions in NFNB resistance, validating it as the gene underlying the mutant phenotypes. Allele analysis and transcript regulation of HvWRKY6 from resistant and susceptible lines revealed sequence identity and upregulation upon

INTRODUCTION
A recent study determined that climate change is a major threat to malt barley production as yield loss is projected due to high temperatures and water deficiency in some growing regions (Xie et al., 2018). However, as climate change brings drought to some regions, others will experience excess precipitation and combined with elevated temperatures, will provide environments that are more conducive to fungal disease epidemics. Thus, without adequate management tools, disease problems can be expected to be exacerbated in these regions. The threat of greater disease epidemics due to rising temperatures in regions with excess precipitation was not accounted for in the predictions of world barley shortages (Xie et al., 2018), thus production shortfalls could be even greater than predicted. For sustainable barley production, the intelligent deployment of durable genetic resistance to important fungal pathogens is critical and is a major focus of breeding programs with the primary goal of releasing high yielding and broadly adapted varieties that produce quality grain across dynamic environments (Horsley et al., 2009).
The net blotch of barley (Hordeum vulgare L.), caused by Pyrenophora teres Drechs. is an economically important foliar disease in barley growing regions worldwide with epidemics causing 10-40% of yield loss when susceptible varieties are grown, but under environmental conditions conducive to disease epidemics losses can reach 100% (Mathre, 1997;Murray and Brennan, 2010). Because the brewing and distilling industries demand quality malting barley, which brings premium prices, producers are concerned by biotic or abiotic stresses that negatively affect yield and quality (Grewal et al., 2008). Foliar infection by P. teres is a major concern as it has a large impact on yield, but foliar and kernel infection can also impact quality (Liu et al., 2011). The most sustainable and environmentally friendly way to manage net blotch is deploying effective genetic resistance into varieties, yet, a better understanding of this complex pathosystem and the quantitative nature of the host pathogen genetic interactions is needed. Practicing the stewardship of effective resistance sources is important when deploying resistance so that resistance is durable and genes are conserved. One way of accomplishing this goal is through gene discovery and subsequent functional analysis which is more practical with the new array of genomic tools available to the barley research community.
Net blotch exists in two forms; net form net blotch (NFNB) caused by Pyrenophora teres f. teres (Ptt) and spot form net blotch (SFNB) caused by Pyrenophora teres f. maculata (Ptm) (Smedegård-Petersen, 1971; Steffenson and Webster, 1992). The symptoms of NFNB first appear as small dark brown necrotic lesions that expand over time forming net like longitudinal and transverse striated necrotic lesions commonly surrounded by chlorosis on susceptible host genotypes ( Figure 1A). The symptoms of SFNB also initially appear as small dark brown necrotic lesions that expand over time producing elliptical necrotic lesions typically surrounded by chlorosis on susceptible host genotypes ( Figure 1B). Although, these two pathogens are morphologically identical (conidia and mycelium), their genetics as well as host-pathogen interactions are considered to be relatively distinct (Liu et al., 2011), thus, are considered to be different diseases. However, contrary to this statement the recent review of barley-Pyrenophora teres genetic interactions showed that 17 of the 19 Ptm resistance/susceptibility loci known to date overlap with Ptt resistance/susceptibility loci (Clare et al., 2020). This may include the 3H QTL reported by Koladia et al. (2017), where they identified the HvWRKY6 transcription factor required for NFNB resistance as there were also SFNB resistance QTL identified by biparental and association mapping within the region. Interestingly, we also generated data showing that the HvWRKY6 gene described here is also required for Ptm resistance (data not presented).
The barley line CIho 5791 (hereafter referred to as CI5791), is an Ethiopian breeding line, that is highly resistant to most Ptt isolates collected from barley growing regions worldwide (Mode and Schaller, 1958;Steffenson and Webster, 1992;Wu et al., 2003;Richards et al., 2016;Koladia et al., 2017). CI5791 consistently exhibits a highly resistant reaction at the seedling stage (secondary leaf) mainly consisting of small pinpoint necrotic lesions ( Figure 1A) and the resistance translates to the field at the adult plant stages. However, CI5791 resistance has been compromised by Canadian, French (Arabi et al., 1992;Akhavan et al., 2016), and Moroccan isolates that are moderately virulent on CI5791 ( Figure 1C). Although, the dominant CI5791 resistance located on chromosome 6H is remarkably broad and effective, it is apparent that pathogen populations have the diversity to overcome this resistance.
distinct Ptt isolates. The major CI5791 6H QTL was shown to be dominantly resistant and was effective against all isolates used in the study. The gene underlying the CI5791 3H QTL was also dominant in nature, but was only effective against two Japanese Ptt isolates, JPT0101 and JPT9901. Interestingly, a similar chromosome 3H QTL was also contributed by Tifang and was shown to confer dominant resistance against Ptt isolates Br. Pteres (Brazil), BB06 (Denmark), 6A (California, United States), and 15A (California, United States).
Exome capture is a cost-effective yet powerful genomics tool that allows for targeted sequencing of the coding regions (exons) of specific specie genomes. This tool has been used for the identification of polymorphism within coding regions that contribute to disease in humans, other animals, and plants (Choi et al., 2009;Raca et al., 2010;Wang et al., 2010;Bamshad et al., 2011;Cosart et al., 2011;Mascher et al., 2013bMascher et al., , 2014Mascher et al., , 2016Warr et al., 2015;Russell et al., 2016). The exome capture array specific to barley (Hordeum vulgare), representing 61.6 mega bases of coding region from the complex ∼5.1 Gb barley genome, is available (Mascher et al., 2013b). This array was used to study the domestication and evolution of barley by resequencing and identifying variants in wild barley (Hordeum sponteneum) and land race (Hordeum vulgare) coding regions, including ancient barley germplasm Russell et al., 2016). Mascher et al. (2014) also identified the barley HvMND gene that governs increased tiller numbers utilizing this array. Schreiber et al. (2019) utilized the array on a highly mutagenized TILLING population of the barley variety Golden Promise to identify and evaluate mutation density. They also assembled a collection of semi-sterile mutants from the population and developed a custom exome capture array of 46 candidate genes to identify potential mutations causing the sterility phenotype (Schreiber et al., 2019). The exome capture methodology has also been applied to resistance gene enrichment sequencing (RenSeq), single-molecule real-time RenSeq (SMRT RenSeq), mutagenesis RenSeq (MutRenSeq), and association genetics (AgRenSeq) technology designed to identify nucleotide binding site-leucine rich repeat (NLR) disease resistance genes (R genes) in plants (Jupe et al., 2013;Steuernagel et al., 2016;Witek et al., 2016;Arora et al., 2019). A cDNA RenSeq method was also utilized to accelerate the identification of R genes in tomato (Andolfo et al., 2014).
The exome capture method coupled with forward genetics screens was considered an efficient genomics tool for the identification of resistance/susceptibility genes. Thus, we utilized this methodology to efficiently identify a gene required for broad and effective CI5791 NFNB resistance underlying the major QTL located on chromosomes 3H and 6H previously identified by Koladia et al. (2017). A CI5791 γ-irradiated mutant population was created and screened for individuals with compromised Ptt resistance. Utilizing a forward genetics screen, exome capture, and comparative sequence analysis, the HvWRKY6 transcription factor (TF) gene underlying the CI5791 chromosome 3H Ptt resistance QTL was identified as being required for broad Ptt resistance. To the best of our knowledge this is the first gene identified that contributes to NFNB resistance in barley.

Mutants Development
The Ethiopian barley line CI5791 is highly resistant to most NFNB isolates collected worldwide. A g-irradiation approach was used to develop a CI5791 mutant population. Briefly, ∼500 g of seed was hydrated in an airtight container with 60% glycerol for about 7-10 days. The hydrated seeds were irradiated with 35 kilorads (350 Gy) of g rays in a Gammator (M38-4, Radiation Machinery Corporation) prior to planting. Approximately, 1,400 M 1 seeds were planted in trays and allowed to self-generate the M 2 generation. Approximately 10,000 M 2 seedlings, derived from the original 1,400 M 1 individuals were screened by inoculation with the Ptt isolate LDNH04Ptt19 (hereafter referred to as LDN) collected from Langdon, North Dakota. Planting, inoculum preparation, inoculation, and disease evaluations were performed as described in Friesen et al. (2006). After identifying putative mutants, these seedlings were transplanted to 15.24 cm (6 inch) pots and allowed to self-generate M 3 generation seeds. The M 3 generation seeds were planted in cone-tainers and screened with Ptt isolate 0-1 using three replicates with WT CI5791 and Heartland as the resistant checks and the susceptible check Robust. Ptt isolate 0-1 is a Canadian isolate collected from Ontario (Weiland et al., 1999;Wyatt et al., 2018) and has similar virulence as the isolate LDN which was used to identify the original CI5791-γ3 and CI5791-γ8 mutants from the M 2 generation. The planting, inoculum preparation, inoculation, and disease reading were performed as described in Friesen et al. (2006).

Mapping Populations and Phenotyping
Two F 2 mapping populations were developed by crossing CI5791-γ3 and CI5791-γ8 homozygous mutant M 3 individuals with the NFNB resistant barley line Heartland. Heartland is a sixrowed spring feed barley that was developed at the Agriculture Canada Research Station, Brandon, Manitoba and registered and released in 1984 (Therrien and Wolfe, 1985). Heartland was shown to be resistant to three major Canadian races of Ptt before its release. Heartland is hypothesized to contain a similar dominant resistance gene as CI5791 at the chromosome 6H locus designated as Rpt5. The planting, inoculum preparation, inoculation, and disease reading for the F 2 individuals from each of the CI5791-γ3 and CI5791-γ8 × Heartland populations were performed as described in Friesen et al. (2006). The disease reading was performed 7 days after inoculation (DAI) using a 1-5 rating scale (Neupane et al., 2015) for CI5791-γ3 × Heartland due to the phenotypic resemblance to Ptm infection and a 1-10 rating scale developed by Tekauz (1985) for CI5791-γ8 × Heartland. The CI5791-γ3 × Heartland F 2 susceptible individuals using a rating cutoff of > 2 representing the homozygous CI5791 genotype at the mutant region were used for mapping the gene. We utilized PCR genotyping by sequencing (PCR-GBS) to genotype CI5791-γ3 × Heartland F 2 homozygous susceptible lines (a total of 34 lines representing 68 recombinant gametes) on an Ion Torrent TM PGM: a PCR-GBS marker panel designed for polymorphism between CV. Tradition and barley line PI67381 consisting of 365 markers (Supplementary Table 1) was used to genotype all 34 susceptible F 2 lines. Primer development, DNA extraction, PCR cycle parameters, library preparation, and sequencing on the Ion Torrent TM PGM were performed as previously described in Richards et al. (2016).
The disease severity of the 34 CI5791-γ3 × Heartland F 2 homozygous susceptible lines along with the genotypic data were used for QTL mapping using MapDisto 2.0 (Heffelfinger et al., 2017) and QGene 4.0 (Joehanes and Nelson, 2008). The individual SNP calls were filtered for a minimum genotype quality of 10, and a minimum read depth of 3. The markers with more than 30% missing data and MAF < 25% were removed from further analysis. Single marker regression was used to identify the susceptible QTL in the γ3 × Heartland F 2 population. CI5791 and the two mutants were also phenotyped with the two Moroccan Ptt isolates SM36-2 and SM36-3 that were shown to be moderately virulent on CI5791. The spot form net blotch susceptible barley line Tradition was inoculated with Ptm isolate FGOB10Ptm-1 as previously described in Neupane et al. (2015) to compare a typical SFNB phenotype with the NFNB phenotypes on the CI5791-γ3 and CI5791-γ8 mutants (Figure 1).

Allelism Tests
Reciprocal crosses were made between CI5791-γ3 and CI5791-γ8 to determine if the two putative independent mutants were allelic. Six F 1 s of CI5791-γ3 × CI5791-γ8 and ten F 1 s of CI5791-γ8 × CI5791-γ3 were phenotyped as previously described using the 1-10 rating scale developed by Tekauz (1985) with Ptt isolate 0-1. The F 1 individuals from both crosses were also genotyped utilizing the STS markers targeting the putative mutant gene underlying the region delimited by genetic mapping. The primers FIGURE 2 | Linkage map developed using 34 CI5791-γ3 × Heartland F 2 susceptible individuals (representing 68 recombinant gametes) developed with 123 polymorphic SNP markers showing the seven barley chromosomes. For chromosome 3, the QTL map of resistance/susceptibility to Pyrenophora teres f. teres isolate 0-1 using single marker regression analysis showing the only significant peak (red dashed line) with a LOD score of 71. The X-axis represents LOD values and the Y -axis represents the PCR-GBS SNP markers. The most significant marker was 11_10444 (red boxes) positioned at 78.9 cM with iSelect positions used to develop the map and 74.99 cM with POPSEQ positions with a LOD value of 71. The red region filled in on the chromosome 3 map flanked by the markers 11_20742 and 11_21493 (white boxes) shows the high confidence interval region containing the CI5791-γ3 mutation. The gray box shows the comparative region mapped by Koladia et al. (2017) with the flanking (white boxes) and most significant (red box) markers from our mapping shown relative to their map. The HvWRKY6 gene is shown within the QTL detected in their CI5791 × Tifang biparental mapping population. Markers and the HvWRKY6 gene marked with an asterisk were not mapped by Koladia et al. (2017) but were placed on the map based on their POPSEQ positions for comparison. The markers without an asterisk were mapped by Koladia et al. (2017) but were given POPSEQ positions.
Frontiers in Genetics | www.frontiersin.org used for STS marker development are described below under the STS marker development and mutation validation section.

Exome Capture Sequencing and Analysis
The exome capture, sequencing, and analysis of the CI5791-γ3 and CI5791-γ8 mutants and CI5791 WT were performed as thoroughly described in Solanki et al. (2019). The POPSEQ positions of the markers flanking the QTL identified in the segregating F 2 population (Figure 2), described above, were obtained and used to identify exome capture targets within the mapped region. BAM files from the analysis were imported into CLC Genomics Workbench version 8.0.3 (Qiagen) for the visualization of sequence alignments (Figure 3).

STS Marker Development and Mutation Validation
Based on the identified nucleotide deletion in the CI5791-γ8 mutant specific sequence tagged site (STS), markers was developed. Oligonucleotide primer pairs were developed specific to CI5791 WT and the CI5791-γ8 mutants. Primers WRKY6-F1 (5 -GCCGCTGGTTCTCGTCG TTCATGCG-3 ) and WRKY6-Wt-R1 (5 -TAGTCGACGACGACGGGGCGTCCC-3 ) only produced an amplicon from CI5791 WT (Figure 3) whereas the primer combination of WRKY6-F1 and WRKY6-Mt-R2 (5 -TAGTCGACGACGACGGGGCGT CCG-3 ) only produced an amplicon from the CI5791-γ8 mutant due to designing specificity in the 3 bases at the 3 terminus of the primer that are specific to the 1 bp deletion discovered in the CI5791-γ8 mutant from the exome capture data. The PCR was optimized so the discriminant amplicons were specific to the WT or mutant genotypes. The PCR amplification program was set as: denaturation at 95 o C for 5 min, 25 cycles of 95 o C for 30 s, 76 o C for 1 min, and 76 o C for 30 s, and a final extension of 72 o C for 5 min. Wild-type CI5791, Heartland, CI5791-γ3, CI5791-γ8, homozygous susceptible F 2 individuals from two populations, and 15 randomly selected resistance F 2 lines (CI5791-γ8 × Heartland) were genotyped with the WT and mutant specific primers. The F 1 reciprocal cross between CI5791-γ3 and CI5791-γ8 were also genotyped with these primer sets. All PCR amplicons were visualized on a 1% agarose gel with GelRed R (Biotium).

HvWRKY6 Allele Sequencing and Analysis
To compare if there was the presence of any allelic variation between resistant and susceptible barley cultivars, we sequenced FIGURE 3 | Alignment of exome capture reads from wild-type (WT) CI5791, and the CI5791-γ3 and CI5791-γ8 mutants to the HORVU3Hr1G033740.2 (HvWRKY6) gene model from the cultivar Morex genome sequence. The left shows the reaction of WT CI5791, CI5791-γ3, and CI5791-γ8 to Pyrenohpora teres f. teres isolate 0-1. The center shows the pile up of reads from each genotype to the HORVU3Hr1G033740.2 (HvWRKY6) gene model. The alignment shows only three reads in CI5791-γ3, and ∼2,000 reads in both CI5791-γ8 and WT CI5791 confirming the complete gene deletion in CI5791-γ3. The CI5791-γ8 pileup shows a single bp deletion in the 2nd exon. Reverse PCR primers were developed from the deletion region that produced WT specific amplicons and CI5791-γ8 specific amplicons as shown in the bottom right. The WT specific primers amplified only from Heartland and WT CI5791. The mutant specific primers amplified only from CI5791-γ8 and no amplification in CI5791-γ3 as it is missing the entire gene sequence. The right side shows the predicted proteins encoded by the HvWRKY6 gene in each genotype.
Frontiers in Genetics | www.frontiersin.org HvWRKY6 from CI5791 (resistant), cvs. Tifang (susceptible), and Morex (moderately susceptible). We designed four primer pairs at 1 kb intervals to sequence the entire gene including the promoter region (∼3,544 bp) (Supplementary Table 2). The gDNA extractions were performed as previously described in the "DNA Extraction, Exome Capture, and Sequencing" section above and were quantified using the Qubit TM 2.0 Fluorometer with a Qubit TM dsDNA Broad Sensitivity Kit (Thermo Fisher Scientific). PCR parameters were initial denaturation at 95 • C for 5 min, 35 cycles of 94 • C for 30 s, 60 • C for 60 s, and 72 • C for 60 s, followed by a final extension at 72 • C for 5 min. PCR amplicons were visualized on a 1% agarose gel containing GelRed R (Biotium) and purified using an E.Z.N.A. R Cycle Pure centrifugation columns (Omega Bio-tek) following the manufacturer's protocol. Purified PCR products, ∼40 ng, were sent to GenScript for sequencing following their guidelines.

RNA Extraction, cDNA Synthesis, and qPCR
Quantitative PCR (qPCR) was conducted to measure if differential regulation of the HvWRKY6 gene occurred upon interaction with Ptt isolates in compatible or incompatible interactions. The isolates Ptt 0-1, SM36-2, and SM36-3 were used to inoculate CI5791 and the barley line Tifang was inoculated with Ptt 0-1 only. Primers were designed across exons 1 and exon 2 (Figure 4): wrky6-qpcr-F2 (5 -GTTCCTGCCGTTACTGTCCTCATC-3 ) and wrky6qpcr-R2 (5 -TCGCCATCAAGAAGGAGGACCTCAC-3 ), that specifically amplify ∼120 bp from cDNA and ∼270 bp from gDNA. At least three biological replications were collected from each mock (water + tween 20) and Ptt inoculated plants. Tissues from the first leaves were collected at time point 0 (non-inoculated control), 5, 30 min, 1, 2, 4, 6, 12, 24, 48, 72, 96, 120, 144, and 168 h post inoculation. Tissue samples were immediately flash frozen in liquid nitrogen and stored at −80 • C for further processing. Total RNA was extracted from the collected tissue using an RNeasy Plant Mini Kit (Qiagen) following the manufacturer's instructions. The total RNA was quantified using the Qubit TM Fluorometer and the Qubit TM RNA BR assay kit (Thermo Fisher Scientifics) per the manufacturer's instructions. To ensure RNA integrity and that the RNA was free of gDNA contamination, 1 µl of total RNA was denatured in 4 volumes of denaturing buffer (Formaldehyde Load Dye, Ambion) at 80 • C for 5 min and visualized on a 1% agarose gel with GelRed R (Biotium). RNA samples with the four-intact ribosomal RNA (rRNA) bands at the expected molecular weights of ∼ 3.4, 1.8, 1.5, and 1.1 kb corresponding to the nuclear 28S and 18S rRNAs and the 23S and 16S plastid rRNAs, respectively, without high molecular weight gDNA contamination were considered as quality RNA and used for cDNA synthesis. The GoScript TM Reverse Transcription System (Promega) was used to synthesize cDNA following the manufacturer's protocol. Briefly, ∼1 µg of total RNA was mixed with oligo(dT) 15 primer (0.5 µg) and incubated at 70 • C for 5 min. The RNA sample was then mixed with 15 µl of reverse transcription reaction mix (GoScript TM Reaction Buffer (5X), MgCl 2 (1.5 mM), PCR Nucleotide Mix (0.5 mM each dNTP), Recombinant RNasin Ribonuclease Inhibitor (20 units), and Reverse Transcriptase, and incubated at 25 • C for 5 min followed by 42 • C for 60 min and inactivated at 70 • C for 15 min. The 20 µl cDNA synthesis reactions were diluted with 80 µl H 2 O (1:5). A 10 µl qPCR reaction was prepared by mixing 4 µl of diluted cDNA, 5 µl of SsoAdvanced Universal SYBR Green Supermix (Bio-Rad), and 0.5 µl of each forward and reverse primer (10 µM). The qPCR was conducted in a CFX96 Real-time system thermocycler (Bio-Rad) with cycling parameters of 95 • C for 30 s followed by 40 cycles of 95 • C for 15 s and 60 • C for 30 s; 65 • C for 30 s; and 60 cycles of temperature increasing from 60 to 95 • C with fluorescence readings acquired at 0.5 • C increments per cycle. Three technical replications were used for each biological rep. The barley HvSnor14 gene was used as the reference (Ferdous et al., 2015) to normalize HvWRKY6 gene expression. The efficiency of qPCR for the HvWRKY6 and Snor14 primers were calculated by generating a standard curve by running qPCR on a 10-fold serial dilution starting from 200 pg of the PCR amplified template of HvWRKY6 and Snor14. Differential expression was calculated by using the CT method on the Bio-Rad CFX Manager 3.1 software. A t-test was performed to check the significance of difference at p < 0.05 using a standard error of mean of 1.

BSMV-VIGS
The barley stripe mosaic virus-virus induced gene silencing (BSMV-VIGS) system was exploited to functionally validate the HvWRKY6 as required for resistance in the barley line CI5791. A unique 65 bp sequence was selected from the HvWRKY6 gene by performing a BLASTn search against the low and high confidence gene list in the IPK barley database 1 to reduce the cross amplification and off target silencing of other WRKY TF homologs in the barley genome. Two primer pairs based on the 5 and 3 termini of this unique sequence were designed with NotI and PacI adaptor sequences attached to the 5 ends of the respective primers. These adaptors were reciprocally utilized in order to develop sense and antisense constructs. The first primer set was designed with a NotI adaptor on the forward primer and PacI adaptor on the reverse primer and the second set with the PacI adaptor on the forward primer and NotI adaptor on the reverse primer.

First Primer Set
These two primer sets were used to clone the unique 65 bp HvWRKY6 fragment into the γRNA strand of the BSMV-VIGS infectious cDNA clone PSL38.1 in both sense and antisense orientations. First, the two primer sets were used to produce the gene specific amplicon from CI5791 cDNA in 20 µl PCR reactions consisting of 2 µl of cDNA template, 0.5 µl of each forward and reverse primers (10 µM), 0.3 µl of dNTPs (500 µM), 0.2 µl of GoTaq R (1.25 units), 4 µl of GoTaq R buffer (10x), and 12.5 µl of H 2 O. The PCR cycle parameters had an initial denaturation at 95 • C for 5 min, followed by 35 cycles of 95 • C for 30 s, 60 • C for 30 s, and 72 • C for 30 s followed by a final extension of 72 • C for 5 min. The amplicon was purified using an E.Z.N.A R Cycle Pure centrifugation column (Omega Bio-tek). The purified PCR products was digested in a 30 µl reaction consisting of 0.5 µl of NotI HF (NEB), 0.5 µl of PacI (NEB), 3 µl of Cut Smart Buffer (NEB), 11 µl of H 2 O, and 15 µl of PCR product. The digestion reaction was allowed to incubate at 37 • C for 2 h followed by inactivation at 65 • C for 20 min. The BSMV vector PSL38.1-MCS for cloning the target amplicon was also digested with 3 units of NotI and PacI double digestion reactions using 5 µg of plasmid in a 30 µl reaction. Digested PCR product (2 µl) was mixed thoroughly in an 8 µl ligation reaction mix comprised of 1 µl of predigested vector (∼80 ng), 1 µl of ligation buffer (10X), 1 µl of T4 DNA ligase, and 5 µl of H 2 O and incubated at 4 • C for 24 h. Chemically competent One Shot R TOP10 E. coli cells (Thermo Fisher Scientific) according to the manufacturer's protocol were then transformed with the ligation mix and inoculated into 250 µl of Luria Broth (LB) liquid media and incubated at 37 • C with 230 rpm shaking for 1 h. A total of 100 µl of each transformation was plated onto LB agar plates with 100 µg/ml of ampicillin and incubated overnight (∼12 h) at 35 • C. Ten random colonies were picked from each transformation and inoculated into 2 ml of LB broth [5 g of NaCl, 5 g of tryptone, 2.5 g of yeast extract, and 500 ml of H 2 O and ampicillin (100 µg/ml)] in 12 ml borosilicate culture tubes and incubated overnight with shaking at 230 rpm at 37 • C. The cell cultures were transferred to a 2ml microcentrifuge tube and centrifuged at 12,000 rcfs for 5 min to pellet the cells and the waste supernatant was discarded. The plasmid DNA was extracted from the pelleted cells using the PureYield TM Plasmid miniprep System (Promega) following the manufacturer's protocol.
The BSMV tripartite viral genomic RNAs (α, β, and γ genomes) with the γ fragment containing the unique 65 bp fragments of the CI5791 HvWRKY6 allele cloned in both the sense and antisense orientations, were synthesized via in vitro transcription using the mMESSAGE mMACHINE TM T7 Transcription Kit (Thermo Fisher Scientific) according to the manufacturer's protocol. Twenty microliter reactions of each of the α genome, β genome, γ-HvWRKY6 sense, and γ-HvWRKY6 antisense genomes were combined with 370 µl of FES buffer (100 ml of GP buffer, 5 g of sodium pyrophosphate decahydrate, 5 g of bentonite, 5 g of celite, and up to 500 ml of H 2 O) as the BSMV-VIGS inoculum. A total of 20 µl of each of the α genome, β genome, and γ genome were combined with 390 µl FES buffer for the BSMV-VIGS control inoculum.
Single seeds of the barley line CI5791 were planted per conetainer and placed in racks. Newly emerged secondary leaves still whorled at the ∼10-11 days old seedling stage were inoculated with either 5 µl of each tripartite RNAs or BSMV-VIGS control virus (both in FES buffer). Approximately 40 individual plants were inoculated for each BMSV-VIGs experimental RNAs and control RNAs. Plants were first misted heavily and then inoculated by gently rubbing the leaves with 5 µl of each BSMV-VIGs construct. After incubation in the mist chamber for 24 h at 100% humidity, inoculated plants were moved back to the growth chamber set at 21 • C with a 12 h photoperiod. Once typical BSMV symptoms, mottling and striping, appeared on the expanded or expanding tertiary leaves, plants were inoculated with Ptt isolate 0-1 as previously described in Friesen et al. (2006). Inoculum preparation, inoculation, and disease reading were performed as described before. CI5791 and Robust were used as a resistance and a susceptible check, respectively. Disease reading was performed 7 and 12 days after Ptt inoculation (Figure 4 and Supplementary Tables 3, 4) using the 1-10 scale developed by Tekauz (1985).

mRNA Extraction, RNAseq Library Preparation, and Sequencing
Three leaf samples of equal size (∼2 cm) from each replicate per treatment, non-inoculates and inoculated with Ptt isolate LDN at 3, 21, and 45 h post pathogen inoculation (hpi), were combined in a single tube and flash frozen in liquid nitrogen then stored at −80 • C until further processing for total RNA extraction. The total RNA was extracted from the frozen leaf samples using the RNeasy mini kit (Qiagen) following the manufacturer's standard protocol. RNA concentrations were measured using the Qubit R Broad Range RNA kit on a Qubit TM 2.0 Fluorometer (Thermo Fisher Scientific), and RNA samples were visualized on 1% agarose gels stained with GelRed R (Biotium) to confirm the integrity of the RNA samples. RNA samples with four sharp ribosomal RNA (rRNA) bands; approximate molecular weights of 3.4, 1.8, 1.5, and 1.1 kb corresponding to nuclear 28S and 18S rRNAs and 23S and 16S plastid rRNAs, respectively, without a high molecular weight genomic DNA contamination band were considered quality RNA. One microgram of total RNA was used for RNA sequencing (RNAseq). The library construction was performed using the TruSeq RNA Library Prep Kit v2 (Illumina) following the manufacturer's standard protocol. The final library was validated and quantified on the Agilent 2100 Bioanalyzer. The cDNA libraries from four different samples were pooled into one single tube and were normalized according to the manufacturer's protocol. Each of the library pools were diluted to a concentration of 1.8 pm and sequenced on the Illumina NextSeq R 500 sequencer on a single flow cell at the USDA Cereal Genotyping Centre, Fargo, ND, United States. The NextSeq R 500/550 High Output Kit v2 (150 cycles) was used for the generation of 150 bp single end reads. The raw sequence reads were demultiplexed and converted into individual FASTQ files using bcl2fastq software v2.17.1.14 (Illumina). The FASTQ reads were quality trimmed in CLC Genomics Workbench v8.0.3 (Qiagen) using default settings.

Expression Analysis
The analysis pipeline for mapping reads to the reference genome, quality check, and for expression analyses was performed by first mapping the high quality trimmed sequencing reads to the barley RefSeq v1.0 (see text footnote 1) in CLC Genomics Workbench v8.0.3 (Qiagen). Gene specific and transcript specific reads were obtained from reference genes as well as from the gene track and mRNA tract information. This enabled reads to align to both intronic and intergenic regions. Reads less than 90% identical for 90% of the read length and that mapped to more than 10 positions were discarded. The total reads mapped for each gene model were normalized to obtain reads per kilobase of exon model per million mapped reads (RPKM) values for each sample. In all the comparisons, the false discovery rate (FDR)-corrected p-values were calculated by the exact test using the EdgeR bioconductor package in CLC Genomics Workbench. Analyses were based on a threshold of 0.05 for FDR-corrected p-value and a fold change of 3. All treatments were compared with 0 h control (no Ptt inoculation).

Mutant Identification and Validation
To identify mutants compromised for Ptt resistance seeds of the highly resistant barley line CI5791 was γ-irradiated and ∼10,000 M 2 seedlings derived from ∼1,400 M 1 individuals were used in a forward genetics screen for their reaction to Ptt isolate LDN. Two mutant individuals designated CI5791-γ3 and CI5791-γ8 were identified and advanced to the M 3 generation. The two mutant lines were confirmed by phenotyping M 3 individuals in replicated trials and shown to express similar phenotypes resembling susceptible SFNB reactions when inoculated with Ptt isolate 0-1 the causal agent of NFNB (Figures 1A,B). The CI5791 WT and the resistant variety Heartland, which was used to develop the mutant mapping populations described below, exhibited highly resistant reactions (pin point necrotic lesions) to Ptt isolate 0-1 with an average infection type (IT) of 1.5 based on the 1-10 NFNB rating scale developed by Tekauz (1985) (Figure 1A and Supplementary Tables 1, 5). Phenotyping CI5791-γ3 and CI5791-γ8 showed average ITs of 6.5 and 6.0, respectively with Ptt isolate 0-1 using the NFNB 1-10 scale (Tekauz, 1985; Figure 1A and Supplementary Table 5). The susceptible barley variety Robust exhibited typical NFNB susceptible infection types ( Figure 1A). Inoculation of CI5791 WT with the Moroccan Ptt isolates SM36-2 and SM36-3 exhibited moderately susceptible ITs averaging 4.5 and 3.5, but the mutants were susceptible showing average disease scores of 6.5 and 6 for CI5791-γ3 and 7.0 and 6.5 for CI5791-γ8, respectively ( Figure 1C and Supplementary Table 6). The variety Hockett, exhibited highly resistant reactions to Ptt isolates SM36-2 and SM36-3, with average ITs of 1 ( Figure 1C and Supplementary Table 6). The susceptible barley variety Hector exhibited highly susceptible ITs with an average disease score of 9.5 and 8.0 to isolates SM36-2 and SM36-3, respectively ( Figure 1C and Supplementary Table 6).
CI5791-γ3 and CI5791-γ8 were crossed with cultivar Heartland and F 1 individuals were allowed to self-pollinate to produce the F 2 populations. The CI5791-γ8 × Heartland F 2 populations contained 116 individuals. The F 2 population was inoculated with Ptt isolate 0-1 and an IT of 4 was used as the cutoff for resistance/susceptibility with the CI5791-γ8 × Heartland F 2 s on the 1-10 NFNB scale (Tekauz, 1985). The F 2 population showed a segregation ratio not significantly different from a 3 resistant: 1 susceptible ratio as would be expected for a single recessive mutant gene ( Supplementary  Tables 1, 5, 7). Since the F 2 phenotyping data determined that a single recessive mutation was responsible for the susceptible phenotype in the CI5791-γ8 mutant, reciprocal crosses between CI5791-γ3 and -γ8 were made to determine if the mutations were allelic. Six CI5791-γ3 × -γ8 and ten CI5791-γ8 × -γ3 F 1 individuals were challenged with Ptt isolate 0-1 and all the F 1 individuals showed susceptible reactions that were similar to each of the mutant parental lines with average scores of 6.45 and 6.25, respectively (Supplementary Table 3).

Mutant Mapping
Thirty-four F 2 individuals from a CI5791-γ3 × Heartland population showing the characteristic SFNB-like susceptible lesions indicating that they were homozygous for the mutated gene, were genotyped using PCR-GBS to identify 123 polymorphic SNP markers spread across the seven barley chromosomes (Figure 2). The QTL mapping utilizing the genotyping of the 34 CI5791-γ3 × Heartland F 2 homozygous susceptible lines, representing 68 recombinant gametes, localized the mutation to chromosome 3H within a high confidence ∼75 cM interval flanked by the SNP markers 11_20742 (POPSEQ position; chr = 3H cM = 15.15) and 11_21493 (POPSEQ position; chr = 3H cM = 90.33) (Figure 2). The most significant marker 11_10444 (POPSEQ position; chr = 3H cM = 74.99) had a LOD score of 71 (Figure 2).

Exome Capture Sequencing and Analysis
Sequencing of CI5791 WT, CI5791-γ3, and CI5791-γ8 gDNA enriched via exome capture on an Illumina NextSeq R flow cell resulted in a total of 111,251,482, 103,796,564, and 120,530,567 reads, respectively. Thus, the parallel sequencing of the three exome-captured genotypes represented a balanced sequencing library. Utilizing the gene models underlying the ∼75 cM interval containing the mutant gene, deletion variant analysis identified a 1 bp cytosine deletion in the predicted coding region of the HORVU3Hr1G033740.2 gene model in CI5791-γ8 (Figure 3). Of the ∼120 million sequence reads from the CI5791-γ8 exome capture, 2,085 reads were aligned to the HORVU3Hr1G033740.2 gene model. The single base deletion in gene model HORVU3Hr1G033740.2 was the only deletion identified within the chromosome 3H region containing the mutant gene (Figure 2). The single base deletion in HORVU3Hr1G033740.2 is predicted to be within the second predicted exon, of the barley ortholog of the Arabidopsis WRKY transcription factor 6 gene, designated HvWRKY6. Coverage analysis showed that HvWRKY6 is completely deleted from the CI5791-γ3 mutant as only three sequence reads from the CI5791-γ3 mutant were mapped to the HORVU3Hr1G033740.2 reference sequence gene model (Figure 3). Considering the balance of the multiplexed sequencing library which yielded 111,251,482, 103,796,564, and 120,530,567 reads for CI5791 WT, CI5791-γ3, and CI5791-γ8, respectively, and the numbers of reads that mapped to HvWRKY6 for CI5791 WT (1,983) and CI5791-γ8 (2,085), the three reads that mapped to CI5791-γ3 is well below the threshold of contamination (Figure 3).

Characterization of the Candidate Gene
The variant analysis of the exome capture data pinpointed a single cytosine base deletion at nucleotide position 545 in relation to the adenosine of the start methionine as base 1 of the barley HvWRKY6 gene model (Figure 3). The mutation in the CI5791-γ8 mutant occurs in the second exon, which resulted in a frame shift and the predicted translation of a non-functional 141 amino acid (aa) truncated protein (Figure 3). Read depth analysis showed that HvWRKY6 was completely deleted from the CI5791-γ3 mutant as only three reads out of ∼104 million mapped to the HORVU3Hr1G033740.2 gene model compared to 1,983/∼111 million for CI5791 WT and 2,085/∼120 million for CI5791-γ8 (Figure 3). The HvWRKY6 gene spans 8,026 bp of gDNA localized to barley chromosome 3H at ∼50.03 cM (Figure 2) based on POPSEQ positions (Mascher et al., 2013a). HvWRKY6 is predicted to transcribe a 1,710 nucleotides mRNA consisting of six exons (Figure 3) predicted to encode a 569 aa functional protein (∼59.67 kDa) containing WRKY transcription factor domains including the highly conserved WRKYGQK DNA binding motif (Figures 3, 4). Homology searches utilizing NCBI BLASTp identified the candidate HORVU3Hr1G033740.2 predicted protein as an ortholog of the Arabidopsis WRKY transcription factor 6, thus, was designated HvWRKY6. The predicted HvWRKY6 protein has 50% aa identity and 59% aa similarity with the Arabidopsis WRKY6 protein (query cover 89% and e-value 4e −130 ) (Figure 5). A reciprocal result was obtained when the AtWRKY6 protein was used as the query in a BLASTp search using the IPK barley blast server, identifying only one matching WRKY protein in the barley genome corresponding to HORVU3Hr1G033740.2. Thus, HvWRKY6 represents the only known AtWRKY6 ortholog in the barley genome. InterProScan SMART domain identified a conserved WRKY domain at 300-360 aa in HvWRKY6 with high confidence prediction 2 (Letunic et al., 2015). Analysis of the full length HvWRKY6 gene sequence from CI5791, Morex, and Tifang were identical (Supplementary Figure 1) suggesting that the gene is conserved for its primary aa sequence across both resistant and susceptible barley genotypes. BLAST analysis of the barley pan-genome (Jayakodi et al., 2020) showed 100% sequence identity for the lines Akashinriki, Barke, Golden Promise, Hockett, HOR 3081, HOR 3365, HOR 7552, HOR 8148, HOR 9043, HOR 10350, HOR 13821, HOR 13942, HOR 22559, Igri, OUN333, RGT Planet, ZDM02064 (Jayakodi et al., 2020;Schreiber et al., 2020), Bowman (IBGSC, 2012), Haruna Nijo (Sato et al., 2016), and Lasa Goumang (Zeng et al., 2020). Whereas, ZDM01467 (Jayakodi et al., 2020) contained two concatenated hits of highly identical sequences suggesting a misassembly of the region. Zangqing320, a Tibeten hulless variety (Dai et al., 2018) and B1K-04-12, a wild barley (H. vulgare subsp. spontaneum, Jayakodi et al., 2020) both show 99.9% sequence identity. Two additional wild hordeum lines H. bulbosum and H. pubiflorum (Mascher et al., 2013b) showed, 96.8 and 95.4% sequence identity, respectively, further exemplifying the high conservation of HvWRKY6 in hordeum. In addition, WRKY6 appears to be more highly conserved amongst monocots (wheat 2 http://smart.embl.de/ and Brachypodium distachyyon) than dicots (Arabidopsis and soybean, Figure 5).
All homozygous susceptible F 2 individuals from both the CI5791-γ3 × Heartland and CI5791-γ8 × Heartland populations showed a mutant HvWRKY6 genotype with CI5791-γ8 mutant specific primers (WRKY6-F1 + WRKY6-Mt-R2; Figure 3) or WT specific primers. This was determined by no observed amplification with either primer pair on the CI5791-γ3 × Heartland homozygous susceptible F 2 individuals, which is consistent with the entire gene deletion detected with the exome capture experiment (Figure 3). With the CI5791-γ8 × Heartland homozygous susceptible F 2 individuals there were amplicons produced with the mutant specific primer pair (WRKY6-F1 + WRKY6-Mt-R2) but no observed amplification with the WT specific primer pair, which is consistent with the 1 bp deletion detected with the exome capture experiment showing that they were homozygous mutant individuals (Figure 3). Fifteen randomly selected resistant F 2 individuals from the CI5791-γ8 × Heartland showed a 1 homozygous: 2 heterozygous genotype segregation (data not shown). This genotyping perfectly linked the genetic mutation with the mutant phenotype in this small F 2 population representing 68 recombinant gametes. Also, the genotypes of all the reciprocal F 1 (CI5791-γ3 × -γ8 or CI5791-γ8 × -γ3) individuals had a CI5791-γ8 mutant like genotype, lacking a WT allele, providing further evidence that the two mutants CI5791-γ3 and CI5791-γ8 are allelic.

BSMV-VIGS
To validate the function of HvWRKY6 in NFNB resistance the barley stripe mosaic virus (BSMV) tripartite genome was utilized for post-transcriptional gene silencing constructs (Figure 4). The disease reactions of the BSMV-WRKY6 inoculated plants targeted for the post-transcriptional gene silencing of HvWRKY6 were significantly more susceptible than the BSMV-pBS virus inoculated controls at both 7 and 12 days post inoculation when inoculated with Ptt isolate 0-1 (Figure 4 and Supplementary  Tables 4, 8, 9). The BSMV-pBs virus control inoculations did not show the shift from resistance toward susceptibility.

qPCR and RNAseq
The qPCR experiment conducted on WT CI5791 inoculated with the Ptt isolate 0-1 showed that HvWRKY6 is upregulated at 4 h post inoculation (hpi) at least 5-fold until 6 hpi, then it gradually decreases and maintains a level of ∼1-fold upregulation until 168 hpi (Figure 6). With the Moroccan Ptt isolate SM36-3, which is moderately virulent on CI5791, HvWRKY6 was upregulated 1 hpi by 1.6-fold and increased to 12.6-fold upregulation at 4 hpi, 4-fold upregulation at 6 hpi, 16-fold at 12 hpi, and maintained at least 5-fold upregulation after 96 through 168 hpi (Figure 6). When susceptible cultivar Tifang was challenged with Ptt isolate 0-1, HvWRKY6 was upregulated 2.8-fold at 30 min post inoculation to 22-fold at 6 hpi, 4.8-fold at 12 hpi, 9-fold at 24 hpi, and maintained at least 8-fold upregulation after 96 through 168 hpi (Figure 6). The susceptible line Tifang was utilized in this study due to the fact that the 3H locus encompassing HvWRKY6 (Figure 2) was identified as a resistance QTL in a CI5791 x Tifang cross (Koladia et al., 2017). The qPCR analyses confirmed that the expression of HvWRKY6 in CI5791 was significantly higher with the moderately virulent Ptt isolate SM36-3 than the avirulent isolate 0-1 from 1 to 4 h, and 96 to 168 h. Similarly, the expression of HvWRKY6 in the susceptible line Tifang challenged with the virulent Ptt isolate 0-1 was much higher and significantly different than SM36-3 and 0-1 on CI5791 at the times between 30 min through 2 and 96 h through 168 h.
RNA sequencing was utilized for comparative global transcriptomics between non-inoculated and LDN-inoculated CI5791 at 3, 21, and 45 h post-inoculation. An upregulation of HvWRKY6 (HORVU3Hr1G033740) in the inoculated samples was observed at all of the aforementioned time points (Supplementary Table 10) with HvWRKY6 upregulated 7. 9-, 4. 6-, and 8.3-fold at 3, 21, and 45 h post-inoculation, respectively, in the CI5791 resistance response. As expected the gene enrichment analysis of differentially expressed genes at all the time points showed a significant enrichment of several classes of genes that are involved in the response to fungal pathogens (Supplementary Table 11).

DISCUSSION
Allelic deletion mutations that compromise NFNB resistance in barley line CI5791 were genetically mapped to an approximate 75 cM region of barley chromosome 3H. An exome capture mapping by sequencing approach identified independent mutations of the barley ortholog of the Arabidopsis WRKY6 TF, designated HvWRKY6, in the delimited region. A 1 bp and whole HvWRKY6 gene deletion in the two mutant lines CI5791-γ3 and CI5791-γ8 are consistent with the reports that γ-irradiation induce 1-10 kb deletions (Morita et al., 2009). The HvWRKY6 gene was validated as playing a role in Ptt resistance via virus-induced gene silencing (VIGS) utilizing the barley stripe mosaic virus. Thus, we hypothesize that the HvWRKY6 FIGURE 6 | Relative expression of HvWRKY6 in resistant CI5791 inoculated with Pyrenophora teres f. teres (Ptt) isolate 0-1 (blue) and moderately virulent isolate SM36-3 (red); and the susceptible cv. Tifang inoculated with Ptt isolate 0-1 (green). The Y -axis represents the fold change relative to the non-inoculated control (time point 0) and the X-axis represents the time point at which the leaf samples were collected.
TF plays an important role in defense signaling that results in NFNB resistance.
The two independent mutants CI5791-γ3 and CI5791-γ8 exhibited susceptible symptoms to Ptt isolates that were not typical NFNB symptoms, but rather resemble SFNB lesions when inoculated with the Ptt isolates LDN and 0-1. The symptoms exhibited by the mutants were dark brown elliptical necrotic lesions surrounded by an extensive expanding yellow chlorotic margin (Figure 1A). The chlorosis expands and eventually coalesces with other lesions suggesting underlying pathogen growth, yet the necrotic regions remain relatively confined and elliptical, resembling an SFNB type of susceptible reaction (Figures 1A,B). Therefore, we initially phenotyped the CI5791-γ3 x Heartland F 2 population using a 1-5 SFNB rating scale as described in Neupane et al. (2015). However, the CI5791-γ8 × Heartland F 2 population, BSMV-VIGS experiment, and CI5791-γ3/CI5791-γ8 reciprocal allelism crosses were phenotyped at the F 1 stage using the NFNB 1-10 rating scale as described by Tekauz (1985) as this mutant compromises NFNB disease resistance yet results in SFNB-like symptoms. Interestingly, the mutant's symptoms when inoculated with the two Moroccan Ptt isolates, SM36-2 and SM36-3, resembled typical net type lesions with enlarged chlorosis that coalesced and has visible longitudinal and vertical striation but not as prominent as seen in a typical NFNB susceptible reaction. Thus, we speculate that the HvWRKY6 transcription factor underlying the 3H QTL may function in regulating genes involved in restricting the growth of the pathogen. In a typical CI5791 resistance reaction the pathogen apparently penetrates the host as indicated by the formation of the pin point lesions, yet the pathogen growth is arrested early in the infection process and the lesion growth is effectively stopped at this early stage in the infection process. The CI5791 major dominant resistance gene responsible for this early resistance mechanism maps to the centromeric region of chromosome 6H (Koladia et al., 2017), which was the initial gene targeted in the mutant screening and was suspected to represent an immunity receptor. However, the first two mutants identified were allelic and mapped to the 3H QTL identified in the CI5791 x Tifang biparental mapping by Koladia et al. (2017) (Figure 2; Koladia et al., 2017) and disrupted the HvWRKY6 transcription factor that plays a role in arresting pathogen colonization and spread after penetration. Additional studies have also reported large QLT intervals over chromosome 3H that spanned from 36.26 to 76.56 cM (Dinglasan et al., 2019). However, a smaller interval (51.27-51.77 cM) was consistently reported and designated as qPttCLS, that maps physically to 398.2-435.5 Mbp on the Morex reference genome (Mascher et al., 2017) and therefore is not suspected to encompass HvWRKY6. Several QTLs have been mapped to the 3H region, one at 28.7-36.6 cM (Steffenson et al., 1996), QTLUHs-3H at 29-31 cM (König et al., 2014), QTL UH -3H at 45-51 cM (König et al., 2013), and Rpt-3H-4 at 57.0-66.6 cM (Yun et al., 2005) that may all be in close proximity, although these cannot be directly compared due to differing populations and absence of common marker sets. The locus, NBP_QRptt3-2, identified using GWAS spans 49.65-52.03 cM using POPSEQ positions and physically maps to 160.7-491.8 Mbp (Wonneberger et al., 2017) and therefore includes HvWRKY6. Additional studies have also reported significant markers in close physical and genetic proximity to HvWRKY6 (Burlakoti et al., 2017;Islamovic et al., 2017;Koladia et al., 2017;Daba et al., 2019;Novakazi et al., 2019;Rozanova et al., 2019;Vatter et al., 2017).
Necrotrophic pathogens often produce several host specific necrotrophic effectors (NE) including low molecular weight metabolites and small secreted proteins that interact with dominant host susceptibility genes (Wolpert et al., 2002;Liu et al., 2011Liu et al., , 2012Liu et al., , 2015Stergiopoulos et al., 2013;Shjerve et al., 2014). These interactions often follow the inverse-genefor-gene model (Friesen et al., 2007) triggering programmed cell death (PCD) to facilitate necrotrophic fungal growth resulting in compatible interactions or a susceptible reaction called necrotrophic-effector triggered susceptibility (NETS) Friesen and Faris, 2010;Liu et al., 2015;Faris and Friesen, 2020). Ptt is a necrotrophic pathogen that has been shown to produce a proteinaceous effector designated as PttNE1 that targets dominant susceptibility gene/s on chromosome 6H in barley (Liu et al., 2015) in an inverse gene-for-gene manner resulting in NETS. However, the CI5791 dominant resistance mechanism appears to follow the gene-for-gene model and possibly represents an R-gene-Avr gene interaction that results in an early dominant resistance response. The HvWRKY6 gene appears to be a highly conserved TF that is most likely required for regulating early defense response genes post recognition that function to arrest pathogen spread after penetration. Thus, we hypothesize that it may be activated early in the response to the pathogen providing early resistance, which translates into preventing further proliferation of the fungus after penetration and thereby limiting the growth of the lesions. Ptt isolates SM36-2 and SM36-3 may express other effectors involved in virulence or have a variable Avr gene that evades early recognition and therefore activation of the early resistance signaling pathway, resulting in a moderately susceptible reaction in CI5791. Yet, this more prolific early pathogen growth and the chlorosis results in a higher level of susceptibility in the HvWRKY6 mutants than in the WT, suggesting that the HvWRKY6 transcription factor plays a role in regulating the early defense response genes that function to sequester pathogen spread after penetration thus inhibiting pathogen establishment in the host.
The qPCR analysis was performed because there appeared to be no polymorphism in the primary aa sequence of the HvWRKY6 protein from a small number of resistant and susceptible genotypes analyzed, and BLAST analysis of 25 barley genotypes sequenced as part of a barley pan-genome project also showed 100% aa sequence identity for all but two of the Hordeum vulgare accessions sequenced suggesting that functional polymorphism between incompatible (resistant) and compatible (susceptible) interactions may be due to differential transcription. The qPCR analysis of HvWRKY6 was performed during infection processing. The qPCR data showed that the differential expression of HvWRKY6 occurs in the barley line CI5791 in response to pathogen challenge as early as 1 hpi with the moderately virulent Ptt isolates SM36-3 and 4 hpi with the avirulent Ptt isolate 0-1. Similarly, in the susceptible line Tifang, the differential expression occurred as early as 30 min after inoculation and reached a maximum of 22-fold upregulation at 6 hpi. The line Tifang was utilized as the susceptible line in these analyses due to its use in the CI5791 x Tifang biparental population used to map the 3H QTL. Overall, the expression level in the susceptible cultivar Tifang with Ptt isolate 0-1 was significantly higher than CI5791 with 0-1 and SM36-3. The BSMV-VIGS experiment showed that the specific silencing of the candidate HvWRKY6 gene results in a susceptible phenotype when inoculated with Ptt isolate 0-1. The time course RNA sequencing between non-inoculated and Ptt isolate LDN inoculated CI5791 seedlings showed upregulation of HvWRKY6 in all the inoculated samples at all the time points (Supplementary Table 10). The gene enrichment analysis of differentially expressed genes at all the time points showed a significant enrichment and upregulation of several classes of genes that are involved in response to fungal pathogens including chitin responsive genes, salicylic acid (SA) and jasmonic acid (JA) response genes, and interestingly the upregulation of positive regulators of leaf senescence genes, which could play a role in suppressing the inverse gene-for-gene induction of PCD by this necrotrophic pathogen to facilitate disease (Supplementary Table 11).
The WRKY TFs are one of the largest groups of plant transcription regulators with protein domain architecture consisting of a highly conserved amino acid sequence (WRKYGQK) at the N-terminus and a zinc-finger-motif (C-C-H-H/C) at the C-terminus (Eulgem et al., 2000;Eulgem, 2006;Bakshi and Oelmüller, 2014). WRKY proteins bind specific W-box elements with the consensus sequence (TTGAC/T) at the promotor regions of targeted genes, resulting in transcriptional activation, and for some genes, repression (Eulgem et al., 2000;Yu et al., 2001;Rushton et al., 2010;Agarwal et al., 2011). WRKY TFs are important in diverse plant physiological activities such as pathogen defense responses and abiotic stress responses such as wounding, nutrient deficiency, salt stress (Chen et al., 2009;Kasajima et al., 2010;Cai et al., 2017;Hichri et al., 2017;Li et al., 2017), and developmental processes including senescence, and root growth Somssich, 2001, 2002;Skibbe et al., 2008).
The AtWRKY6 gene regulates both plant defense responses against Pseudomonas syringae pv. tomato as well as senescence in Arabidopsis (Robatzek and Somssich, 2002). WRKY6 and WRKY3 were also shown to regulate defense response in Nicotiania attenuate against the larvae of the insect herbivore Manduca sexta (Skibbe et al., 2008). In wheat, the TaWRKY70 TF underlying the fusarium head blight QTL, QTL-2DL, governs resistance against Fusarium graminearum by regulating the downstream genes, TaACT, TaDGK, and TaGLI, which are involved in resistance responses (Kage et al., 2017). Other WRKY gene families have also been reported to play vital roles in defense responses in rice against Magnaporthe grisea (rice blast) and Xanthomonas oryzea (bacterial leaf blight) (Liu et al., 2007;Shimono et al., 2007;Wang et al., 2007). Li et al. (2004) reported enhanced resistance to the biotrophic fungal pathogen Erysiphe cichoracearum, whereas an increase in susceptibility to the bacterial necrotroph Erwinia carotovora subsp. carotovora occurred upon the upregulation of WRKY70 in Arabidopsis.
The WRKY TFs also negatively regulate plant defense responses (Robatzek and Somssich, 2002;Eulgem and Somssich, 2007) as the overexpression of WRKY38 and/or WRKY62 were found to compromise immunity to the bacterial pathogen P. syringae and appeared to be negative regulators of plant basal defense responses (Mao et al., 2007;Kim et al., 2008). Grunewald et al. (2008) identified WRKY23 as the negative regulator of plant defense responses against the cyst nematode Heterodera schachtii. The WRKY11 and WRKY17 TFs were also shown to function as negative regulators of basal defense responses in Arabidopsis (Journot-Catalino et al., 2006).
The data generated in this study show that the HvWRKY6 gene functions in NFNB resistance and likely plays a role in the activation of defense genes that are required to restrict lesion growth once the pathogen attempts to penetrate or after it has entered the host. This question will be answered once microscopy analysis is performed on the mutant vs. wild-type CI5791 across the infection process. The HvWRKY6 gene is expressed at higher levels during later time points in compatible interactions showing that the upregulation of the TF does not correlate with resistance. However, earlier upregulation of HvWRKY6 in the incompatible interaction may be the key to resistance, suggesting that the TF may mediate a defense response that is only effective when induced early in the host-pathogen interaction. Also, the later induced expression of HvWRKY6 is not deterministic of resistance as this later upregulation also occurs in the compatible interactions which actually had higher induced levels of HvWRKY6 expression across most of the time points tested (Figure 6). Interestingly, some transcription factors do act as a negative regulator of the plant basal defense response when highly expressed, which has been shown in previous studies (Li et al., 2004;Journot-Catalino et al., 2006;Mao et al., 2007;Kim et al., 2008;Xing et al., 2008). However, the loss of Ptt resistance in the mutants suggested a positive role of HvWRKY6 as the disruption of the gene produces a predicted non-functional protein in one mutant and the complete deletion of the gene in the other showing a predominantly positive role of HvWRKY6 in the NFNB resistance responses.
Interestingly, the HvWRKY6 gene falls directly under the CI5791 dominant resistance QTL identified by Koladia et al. (2017). However, allele analysis of Morex, CI5791, and Tifang alleles did not reveal any primary functional polymorphism (Supplementary Figure 1) and the expression analysis data did not convincingly determine any expression polymorphism that would explain the differential resistances and susceptibilities. Thus, it is likely that the HvWRKY6 TF is involved in a basal resistance and may not represent the gene that underlies the 3H dominant resistance QTL reported by Koladia et al. (2017) as there was little information explaining functional polymorphism in a CI5791 × Tifang population. However, the possibility cannot be ruled out that slight expression differences from a finely tuned level between genotypes could result in the mapping of the resistance/susceptibility QTL on chromosome 3H. Polymorphisms could exist within additional regulatory regions that result in altered expression levels between the resistant and susceptible host genotypes and represent the functional polymorphisms that were segregated in the CI5791 × Tifang population. The fact that all currently analyzed barley lines are predicted to be translationally identical but have diverse phenotypic responses suggest that an upstream signaling component that upregulates HvWRKY6 to a threshold for this resistance that if exceeded or not induced at all results in susceptibility. However, further genetic and functional analyses will be required to answer these questions.
Several studies have provided evidence that WRKY TFs are an integral part of the plant immune system including roles in PAMP-triggered immunity, effector-triggered immunity, and systemic acquired resistance (Li et al., 2004;Eulgem and Somssich, 2007;Rushton et al., 2010). Certain WRKY DNAbinding factors serve as components of signal transduction pathways in plant cells in response to pathogens and regulate the expression of certain plant defense genes (Riechmann et al., 2000). AtWRKY6 regulates both plant defense responses against P. syringae pv. tomato as well as senescence in Arabidopsis, which regulates the SIRK gene (Senescence-Induced Receptor like serine/threonine protein Kinase) that encodes a receptorlike kinase that is exclusively localized to the plant cell nucleus (Robatzek and Somssich, 2002). The qPCR and RNAseq data reported here determined that HvWRKY6 was upregulated in response to pathogen challenge and similar to Arabidopsis may regulate senescence-induced genes as this class of genes was shown to also be upregulated in response to the pathogen in the highly resistant barley line CI5791. Since the TaWRKY70 TF was identified as a strong candidate gene that regulates TaACT, TaDGK, and TaGLI in wheat resistance to Fusarium graminearum (Kage et al., 2017), we hypothesize that HvWRKY6 may regulate other defense related genes that are required to restrict pathogen/lesion growth in line CI5791.
Technological advances in genomic tools and methodology allows for the more accurate characterization of complex traits including quantitative disease resistances. Resistance or lack of susceptibility to the necrotrophic specialist pathogens like P. teres are typically quantitative in nature and the genes underlying the resistance loci have been difficult to localize and isolate. Utilizing forward genetics, mapping by sequencing, exome capture, and next generation sequencing data, we identified the HvWRKY6 gene located under a chromosome 3H resistance QTL that is required for the broad and remarkable NFNB resistance in barley line CI5791. To the best of our knowledge this is the first resistance gene or component required for NFNB resistance to be identified. We hypothesize that the HvWRKY6 transcription factor positively functions to regulate defense response genes, which are required for resistance to Ptt in the barley line CI5791 by limiting the growth of the pathogen in the host after initial entry. The work here exemplifies the powerful new molecular tools that will help generate knowledge and resources to genetically improve crops, through marker assisted and genomic selection strategies in a more intelligent manner. Thus, expediting agricultural productivity in a sustainable manner.

DATA AVAILABILITY STATEMENT
The exome capture data is publicly available and has been deposited on NCBI under the BioProject ID PRJNA686994. The RNAseq data has been deposited under the BioProject ID PRJNA683717. Further inquires can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
RB, TF, RH, and AB designed the experiments. XL and JH generated the mutant population. PT identified the mutants.
PT, JR, SS, GA, PD, SC, and KE performed the wet lab experiments and analysis. JR and RS carried out the bioinformatics analysis. RB, PT, SC, KE, and AB wrote the manuscript. TF, JR, SS, SC, and GA proofread the manuscript. RB and AB obtained funding for research. All authors read and approved the final manuscript.

FUNDING
The research presented in this manuscript was supported by the USDA National Institute of Food and Agriculture Hatch project 1014919, Crop Improvement and Sustainable Production Systems (WSU reference 00011), the United States Department of Agriculture National Institute of Food and Agriculture award # 2018-67014-28491, and the National Science Foundation grant award # 1759030. Any opinions, findings, conclusions, or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the United States Department of Agriculture National Institute of Food and Agriculture or the National Science Foundation.