Phenotypically, progeny of luteostrians can be classified into three categories: (1) green-yellow striped, (2) completely yellow, and (3) albino (Figure 1A). The green-yellow striped seedlings were able to complete the life cycle while the yellow and albino plants do not survive beyond the three-leaf vegetative growth stage. Considering the variegated pattern to be analogous to the albostrians mutant of barley (Li et al., 2019b), we named the causal gene underlying the luteostrians striped phenotype as HvLST, representing Hordeum vulgare LUTEOSTRIANS. Selfing of heterozygous plants leads to about a quarter of aborted grains, indicating that homozygosity for the gene (lst/lst) is lethal at post-zygote or early embryonic stage (Table 1). Inheritance of the yellow striped phenotype was gamete-dependent (Figure 1B) since variegation is observed if the lst allele is inherited through the female gamete (i.e., lst/LST). This interpretation is supported by three lines of evidence: (1) when wild type (LST/LST) plants are fertilized with pollen from striped (lst/LST) or green heterozygous plants (LST/lst), in either case, only half of the F1 plants offspring showed phenotypic segregation in F2 generation; (2) a segregation ratio of 2:1 (green:striped) was observed in segregant offspring in F2 generation; and (3) heterozygous F2 plants (LST/lst or lst/LST) segregated with green (LST/LST or LST/lst) and striped (lst/LST) progenies with a ratio of 2:1 in subsequent F3 generation (Figure 1B). Thus, in addition to its essential function during embryogenesis, we postulate that HvLST plays an essential role in plastid differentiation/programming during gametic stage in the egg cell.

Defective chloroplasts do not contain 70S ribosomes in the albostrians mutant (Hess et al., 1993; Li et al., 2019b). In an attempt to check the function of the translation machinery in plastids of the luteostrians mutant, we initially examined accumulation of the rRNAs in wild type and luteostrians mutant. In analogy to the albostrians mutant, the 16S and 23S rRNA species are not observed in defective plastids of the luteostrians mutant (Figure 2), indicating the lack of 70S ribosomes and consequently missing chloroplast translation in plastids of the yellow leaf sectors. Next, we studied chloroplast/plastid ultrastructure by transmission electron microscopy. These analyses revealed normal chloroplast development in wild type and in green leaf sections of the luteostrians mutant; in both cases chloroplasts contained well-developed stroma and grana thylakoids (Figures 3A,B). In contrast, plastids in yellow sectors of the variegated luteostrians leaves contained no grana and only rudimentary stroma lamellae (Figure 3C). Notably, chloroplast 70S ribosomes were not detectable in yellow plastids of the luteostrians mutant (Figure 3). Altogether, based on the absence of chloroplast rRNA species and the lack of 70S ribosomes it can be postulated that chloroplast translation is abolished in the luteostrians mutant.

Two F2 mapping populations, designated as ‘BL’ and ‘ML’, were constructed for the purpose of genetic mapping of the HvLST gene. The genotypic status of the luteostrians locus of F2 plants was determined by phenotypic segregation analysis of their respective F3 progenies. As homozygosity (lst/lst) leads to grain abortion, the BL population contained 95 wild type (LST/LST) and 172 heterozygotes (LST/lst and lst/LST), consistent with the segregation of a single recessive gene (Table 2). Genetic mapping in the ML population was affected by segregation distortion with genotype ratios deviating from the expected Mendelian ratios (Table 2). Genotyping-by-sequencing was performed for 267 and 269 F2 genotypes in BL and ML populations, respectively. Sequencing data was mapped to the reference genome of barley (Monat et al., 2019) for SNP calling. In total, 3,745 and 5,507 SNP markers were obtained genome-wide at a minimum sequencing coverage of six-fold for BL and ML populations, respectively. By applying a permissive threshold of 5% missing data for both molecular marker and F2 genotype, mapping could be performed in 124 F2/3369 SNPs and 146 F2/4854 SNPs for the BL and ML populations, respectively (Supplementary Table 1). Genetic maps for seven linkage groups, representing the seven barley chromosomes, were obtained for BL (LOD ≥ 6) and ML (LOD = 10) populations, comprising mapped markers for 66.4% (2238/3369) and 71.5% (3469/4584) of the originally defined SNPs, respectively (Supplementary Table 1). The accuracy of the linkage maps was consistency-checked by aligning genetic marker positions with their respective physical order in the reference genome of barley (Monat et al., 2019). The HvLST gene was assigned to the region of the genetic centromere region of chromosome 2H with a physical distance between the flanking markers of 461.7 and 499.9 Mbp in BL and ML populations, respectively (Supplementary Figure 1).

Next, PCR-based KASP (Kompetitive Allele-Specific PCR) markers were designed and employed to saturate the identified HvLST intervals in both mapping populations. Initially, genotyping all the 267 (BL population) and 269 (ML population) F2 individuals with two population-specific KASP markers (chr2H_178489464 and chr2H_461744325 for BL population; chr2H_200083049 and chr2H_463274491 for ML population) allowed to confirm and further narrow down the HvLST target region. Notably, the numerical value within each marker designation indicates the physical coordinate of the mapped SNP position on the reference genome of barley (Monat et al., 2019; Figure 4). Subsequently, saturation mapping of the HvLST interval was performed in a total of 18 recombinants (13 in BL/5 in ML) with additional KASP markers that were suitable for both mapping populations (Supplementary Table 2). Finally, we delimited the HvLST target region to a 0.74 centiMorgan (cM) interval between flanking markers chr2H_431057673 and chr2H_458001177, spanning a distance of ∼26.94 Mbp (Figure 5 and Supplementary Table 2). A cluster of markers co-segregated with the HvLST locus, suggesting that maximum genetic resolution was achieved at the given size of both mapping populations (Figure 4).

Based on the annotated reference genome sequence of barley (Monat et al., 2019) the genetic interval of HvLST delimited by the closest flanking molecular markers is annotated with 284 genes (Figures 5A,C), a number too large for realistically spotting directly a candidate gene for HvLST. One option, higher resolution genetic mapping of the gene HvLST by screening for recombinants in a much larger mapping population was rejected due to the unfavorable genetic placement of the gene in a region of very low recombination. Instead, we evaluated the possibility to survey directly for mutations in any of the predicted genes in the HvLST interval by high-throughput re-sequencing.

First, we simulated the feasibility of the approach based on existing data for the variegation mutant albostrians (Li et al., 2019b). Previously, low-resolution mapping with 91 F2 genotypes delimited the gene HvCMF7 to a 6.05 cM genetic interval between flanking markers Zip_2613 and 1_0169. Anchoring the flanking markers to the barley reference genome (Monat et al., 2019) revealed a 17.86 Mbp physical interval comprising 323 genes (Figures 5A,B); comparable to the situation for HvLST. Next, we used whole genome re-sequencing data for the mutant M4205 (original albostrians mutant) and barley cv. Haisa (genetic background used for induction of the albostrians mutant). In an initial screen we filtered M4205-specific SNPs in genic regions revealing nine candidate genes. The second step of filtering for functional SNPs (e.g., leading to non-synonymous exchange of amino acid, change of splice junction, or premature stop) retrieved four candidate genes (Figures 5A,B, Table 3 and Supplementary Table 3). Based on functional annotation of these four genes, a single gene, HORVU.MOREX.r2.7HG0603920.1, homolog to the Arabidopsis gene CIA2, was suggested as the most promising albostrians candidate, supported by photosynthesis-related pale green phenotype of the Arabidopsis cia2 mutant (Sun et al., 2001). HORVU.MOREX.r2.7HG0603920.1 in fact represents the genuine albostrians gene, as was verified by independent mutant analysis (Li et al., 2019b). In conclusion, based on a delimited genetic interval, survey sequencing of wild type and mutant genetic background may provide a direct path for candidate gene identification for photosynthesis related phenotype mutants in barley.

We adopted this approach to identify a HvLST candidate gene. A wild type (line luteostrians-P2_2, LST/LST) and a mutant genotype (line luteostrians-P3_6, lst/LST) were selected from the originally segregating luteostrians mutant, to ensure both analyzed genotypes sharing the same genetic background of line MC20. In contrast to the previously described simulation for the gene albostrians, the initial screening exercise for sequence polymorphisms between mutant and wild type had to consider the fact that homozygous lst/lst is embryo lethal, hence functional polymorphisms were expected to be present at heterozygous state. Consequently, 54 out of the 284 genes, identified for the luteostrians mapping interval, carried heterozygous SNPs specific to the luteostrians mutant. For a shortlist of eleven genes, the observed SNP was predicted to induce a putative functional change (Figures 5A,C, Table 3 and Supplementary Table 3). BLASTp search (Mount, 2007) against the Arabidopsis proteome revealed presence of orthologs for 8 of the 11 candidate genes. None of the remaining three genes likely represented a genuine HvLST candidate as they either encoded for a putative retrotransposon protein in barley or showed similarity to an organelle gene lacking essential function in chloroplast biogenesis in Arabidopsis (Table 3 and Supplementary Table 3). We then inspected functional annotation information and, if available, phenotype information for mutants of the narrowed shortlist of eight Arabidopsis genes. One of these genes AT5G50920 encodes an ATP-dependent Clp protease ATP-binding subunit ClpC1. Mutants of this gene exhibited a pale-yellow phenotype with reduced photosynthetic performance (Sjogren et al., 2004). Notably, inactivation of the ClpP genes (e.g., ClpP4, ClpP5, ClpP6) in Arabidopsis resulted in embryo lethality and antisense repression lines exhibited a variegated ‘yellow-heart’ phenotype (Clarke et al., 2005). In conclusion, HORVU.MOREX.r2.2HG0135340.1 (HvClpC1), homolog (putative ortholog) of Arabidopsis ClpC1, represented the most likely candidate for the gene luteostrians. The predicted heterozygous SNP in HvClpC1 was confirmed by Sanger sequencing (Supplementary Table 4). The luteostrians mutant carries a heterozygous SNP at position 2,078 (i.e., G2078A; coordinate refers to the coding sequence), consequently, changing glycine into aspartic acid at position 693 (i.e., G693D) (Figure 6B).

The HvLST candidate HvClpC1 has two close homologs in barley, HORVU.MOREX.r2.5HG0373350.1 (homolog 1) and HORVU.MOREX.r2.4HG0286520.1 (homolog 2), which share 85.57 and 74.26% amino acid identity with HvClpC1, respectively. All three homologs share the same gene structure with nine exons. In line with the chloroplast localization of ClpC1 in Arabidopsis (Sjogren et al., 2014), ChloroP predicts presence of a chloroplast transit peptide for all three homologs in barley (Figure 6A), potentially suggesting a chloroplast localization.

We looked up the expression profile of HvClpC1 and its two close homologs in BARLEX database (Colmsee et al., 2015). HvClpC1 shows ubiquitous expression in all examined tissues except young developing inflorescences (INF1) (Figure 6C). The expression pattern of homolog 1 resembles that of HvClpC1 but in most tissues at a distinctly lower level. In contrast, expression of homolog 2 is barely detectable across all the samples.

Low-resolution mapping delimited HvLST to a large interval spanning a physical distance of 26.9 Mbp comprising a total of 284 genes. By comparative whole genome re-sequencing of wild type and mutant lines, a homolog of the Arabidopsis gene AT5G50920 was identified as candidate LUTEOSTRIANS gene. AT5G50920 belongs to the ATP-dependent protease family and encodes the ATP-binding subunit ClpC1. The clpC1 null mutant exhibited a homogenous pale-yellow phenotype (Sjogren et al., 2004). The distinct phenotypes of ClpC1 knockout lines in Arabidopsis and luteostrians barley are reminiscent of the recently reported orthologous pair of genes of the variegated albostrians barley mutant and its Arabidopsis homolog pale green mutant cia2 (Sun et al., 2001; Li et al., 2019b).

The CLP protease system is a central component of the chloroplast protease network (Olinares et al., 2011). It plays essential roles in coordinating plastid proteome dynamics during developmental transitions such as embryogenesis and leaf development (Nishimura and van Wijk, 2015). Genetic studies demonstrated that members of the CLP protease family contribute to chloroplast biogenesis (Constan et al., 2004; Sjogren et al., 2004) and embryogenesis (Kovacheva et al., 2007; Kim et al., 2009). ClpC1 plays important roles on chlorophyll biosynthesis as it controls turnover of glutamyl-tRNA reductase (GluTR) and chlorophyllide a oxygenase (CAO); GluTR and CAO catalyze an early step of the tetrapyrrole biosynthesis pathway and the conversion of chlorophyllide a to b, respectively (Nakagawara et al., 2007; Apitz et al., 2016; Rodriguez-Concepcion et al., 2019). Although single mutants were viable, siliques of the double mutant clpC1clpC2 (previously described as hsp93-V/hsp93-III-2) contained aborted seeds (because of a block in the zygote-embryo transition) and failed ovules (because of a moderate defect in female gametophytes) (Kovacheva et al., 2007). This suggested that both homologs ClpC1 and ClpC2 confer important functions to cell viability during gamete stage and early zygotic stage during embryogenesis. Moreover, the proteome phenotype of clpC1clpS1 suggested a ClpS1 and ClpC1 interaction effect on plastid gene expression components and nucleoid interactors, including RNA processing and editing, as well as 70S ribosome biogenesis (Nishimura et al., 2013). So far, biochemical and genetic evidences gathered for Arabidopsis single and double mutants related to clpC1 were in line with the observed yellow-striped, chloroplast ribosome deficient and embryonic lethal phenotype (at homozygous state) of the luteostrians mutant. Mutants of related genes showed also variegation in Arabidopsis, e.g., Arabidopsis mutant yellow variegated 2 (var2) exhibited a patched phenotype (Chen et al., 2000). VAR2 encodes an ATP-dependent metalloprotease FtsH2. Together with the CLP and the LON family proteases, FtsH2 belongs to the ATP-dependent protease AAA + (ATPase associated with various cellular activities) superfamily (van Wijk, 2015). FtsH2 showed functional redundancy with FtsH8 and a threshold model (i.e., level of FtsH protein complexes formed in the thylakoid membrane) was proposed for the underlying mechanism of var2 leaf variegation (Aluru et al., 2006). In analogy to the barley mutant albostrians (Li et al., 2019b), green and yellow sectors of luteostrians have the same nuclear genotype. Furthermore, no differences could be determined for the chloroplast ultrastructure and plastid rRNA content between wild-type leaves and green sectors of striped luteostrians leaves. Therefore, variegation of luteostrians leaves may be caused by a threshold-dependent mechanism.

Though HvClpC1 represents the most promising candidate gene for HvLST, further experimental evidence is required to confirm this hypothesis. As previously demonstrated, reverse genetic approaches like TILLING or site-directed mutagenesis using Cas9 are viable options in barley to further check the identity of the candidate gene with LUTEOSTRIANS (Gottwald et al., 2009; Lawrenson et al., 2015; Li et al., 2019b).

In this study, we demonstrated that a combination of low-resolution genetic mapping and comparative whole-genome skim sequencing analysis of mutant and wild type can serve as an effective strategy for gene cloning in barley, especially for genes that share a highly conserved function in plants, including those involved in photosynthesis and chloroplast differentiation and development. The probability of having available functional characterization data for such genes in model plants is high and there is increased probability for seeing related or conserved phenotypic effects in both model and non-model plants. By adopting this strategy, we pinpointed the albostrians gene HvCMF7 among 323 candidates within a 17.86 Mbp interval delimited with 91 F2 genotypes, circumventing the time-consuming and laborious fine-mapping process as performed previously (Li et al., 2019b). Implementation of the pipeline in case of the gene luteostrians suggested HvClpC1 to be the likely causative gene underlying the yellow-stripe variegation phenotype. Notably, the cloning strategy is independent of the effect of genomic location. Compared to the barley genome-wide average ratio of physical to genetic distance of 4.4 Mb/cM (Kunzel et al., 2000; Mascher et al., 2017), HvAST and HvLST are located in genomic regions with either increased (2.95 Mb/cM) or suppressed (36.35 Mb/cM) rates of recombination, respectively. The success of candidate gene identification in both cases can be attributed to the achieved high-quality and completeness of the reference genome of barley (Mascher et al., 2017; Monat et al., 2019). Whole genome re-sequencing of wild type and mutant parental lines revealed four and eleven candidates for HvAST and HvLST, respectively. To further narrow down the shortlist to a single candidate gene, it was critical that functional analyses for the putative Arabidopsis homologs were available. In both cases, the best candidate gene showed a chlorophyll-deficient phenotype for mutants of the respective Arabidopsis homolog.

Genetic mapping information is available for 881 barley mutants of the so-called Bowman near-isogenic lines (NILs) (Druka et al., 2011); among them, 142 mutants with pigmentation phenotype. Although the original genetic background of the Bowman NILs is represented by over 300 different barley genotypes, approximately one-half of the mutants (426/881) were identified in one of eight barley cultivars (i.e., Bonus, Foma, Betzes, Akashinriki, Morex, Steptoe, Volla, and Birgitta) (Druka et al., 2011). Since leveraging of the high-quality barley reference genome sequence (Mascher et al., 2017; Monat et al., 2019), whole-genome resequencing of the parental lines would help to filter background mutations and facilitate candidate gene identification and isolation. Therefore, it is conceivable that gene identification in barley, given this pilot attempt based on chlorophyll/photosynthesis related mutants, reached to a point to become even faster and more systematic now.

The luteostrians mutant was derived from line MC20 [Institute of Genetics ‘Ewald A. Favret’ (IGEAF), INTA CICVyA, mutant accession number] by chemical mutagenesis using sodium azide. MC20 was derived from the spring barley cv. Maltería Heda by gamma ray mutagenesis. Two mapping populations were constructed. The six-rowed spring barley (Hordeum vulgare) cv. ‘Morex’ and two-rowed spring barley cv. ‘Barke’ were used as maternal parent in crossings with the variegated mutant line either luteostrians-P1 (lst/LST) or luteostrians-P3 (lst/LST), respectively. The derived F1 progeny were self-pollinated to generate F2 mapping populations designated as ‘ML’ (Morex x luteostrians) and ‘BL’ (Barke x luteostrians). All the F2 and F3 individuals were grown under the greenhouse condition with a day/night temperature and photoperiod cycle of 20°C/15°C and 16-h light/8-h darkness, respectively. Supplemental light (at 300 μmol photons m^–2 s^–1) was used to extend the natural light with incandescent lamps (SON-T Agro 400, MASSIVEGROW).

For ribosomal RNA analysis, wild type (luteostrians-P2; genotype LST/LST) and mutant (luteostrians-P1; genotype lst/LST) lines were grown under the greenhouse light condition. For dark treatment, seeds were germinated within a carton box wrapped with aluminum foil under the greenhouse condition mentioned above.

The trait ‘leaf color’ was scored for the F2 mapping populations at first, second and third leaf seedling stages. The phenotype of the seedlings was classified into two categories: green and variegated (i.e., green-yellow striped). Green phenotype defined all the three leaves of each seedling were purely green; and variegated phenotype defined green-white striped pattern can be observed on all or any of the three leaves. According to inheritance pattern of the causal gene luteostrians (Figure 1B), variegated plants are heterozygous for the mutant allele (lst/LST). Green plants, however, can be either wild type (LST/LST) or heterozygous (LST/lst). Therefore, the luteostrians genotype of all the F2 plants was determined by phenotyping 32 seedlings of each corresponding F3 family. A wild type (LST/LST) F2 would produce 100% green progeny, progeny of heterozygous (LST/lst) F2 green plant would follow a Mendelian segregation (green:variegated = 2:1).

Total RNA was extracted using TRIzol reagent (Invitrogen, Braunschweig, Germany) following manufacturer’s instructions. Initially, concentration of the RNA was determined by help of a Qubit^® 2.0 Fluorometer (Life Technologies, Darmstadt, Germany) according to manufacturer’s instructions. RNA samples were further diluted within a quantitative range of 1 - 10 ng/μL. RNA quality and quantity were then measured using Agilent High Sensitivity RNA ScreenTape following the manufacturer’s manual (Agilent, Santa Clara, United States).

Primary leaves were collected from 7-day-old seedlings of wild type and variegated mutant (luteostrians-P1). Preparation of 1-2 mm² leaf cuttings for ultrastructure analysis including aldehyde/osmium tetroxide fixation, dehydration, resin infiltration and ultramicrotomy were performed as previously described (Li et al., 2019a).

Genomic DNA was extracted from primary leaves of 7-day-old seedlings using a GTC-NaCL method in 96-well plate format. Frozen leaf samples were homogeneously crushed by help of a Mixer Mill MM400 (Retsch GmbH, Haan, Germany) for 30 s at 30 Hz. 600 μL of preheated (65°C for 1 h) GTC extraction buffer (1 M guanidine thiocyanate, 2 M NaCl, 30 mM NaCOOH pH = 6.0, 0.2% Tween 20) was added to the frozen powder and mixed thoroughly for 30 s at 30 Hz under Mixer Mill MM400. Centrifugation at 2,500 x g for 10 min at 4°C after incubation the samples at 65°C for 1 h. 480 μL of supernatant was transferred into a 96-well EconoSpin plate (Epoch Life Science, Texas, United States). EconoSpin plate then vacuumed on the vacuum manifold (MACHEREY-NAGEL, Dueren, Germany) until no droplet drop down, followed by washing theDNA samples twice with 880 μL washing buffer (50 mM NaCl, 10 mM Tris-HCl pH = 8.0, 1 mM EDTA, 70% ethanol). The EconoSpin plate was planced onto a 96-well Microtiter^TM microplate (Thermo Fisher Scientific, Braunschweig, Germany) and centrifugation at 2,500 x g for 3 min to remove residual wash solution. DNA was dissolved with 100 μL preheated (65°C for 1 h) TElight buffer (0.1 mM EDTA, 10 mM Tris-HCl pH = 8.0). The isolated DNA were used for downstream analysis or put at −20°C for long term storage.

For the genetic mapping of HvLST, GBS were prepared from genomic DNA extracted as described previously, digested with PstI and MspI (New England Biolabs, Frankfurt am Main, Germany) following published procedure (Wendler et al., 2015). DNA samples were pooled in an equimolar manner per lane and sequenced on Illumina HiSeq2500 for 107 cycles, single read, using a custom sequencing primer. The GBS reads were aligned to the reference genome of barley as described by Milner et al. (2019). Reads were trimmed using cutadapt (Martin, 2011), mapped with BWA-MEM version 0.7.12a (Li and Durbin, 2009) to barley reference genome (Monat et al., 2019) and the resulting bam files were sorted using Novosort (Novocraft Technologies Sdn Bhd, Selangor, Malaysia). Variants were called with samtools version 1.7 and bcftools version 1.6 (Li, 2011) and filtered following the protocol of Milner et al. (2019) for a minimum depth of sequencing of six to accept a genotype call, a maximum fraction of heterozygous call of 60% and a maximum fraction of 20% of missing data. In the case of ‘BL’ population, SNP calls were converted to reflect the polymorphisms between the two parents, using the calls for barley cv. Barke.

Genetic linkage groups were constructed by help of the JoinMap^® 4.1 software (Van Ooijen, 2006) following the instruction manual. Homozygous wild type and heterozygous allele calls were defined as A and H, respectively; missing data was indicated by a dash. A permissive threshold of 5% missing data for both molecular marker and F2 genotype was applied to the datasets. Regression mapping algorithm and Kosambi’s mapping function were chosen for building the linkage maps. Markers were grouped into seven groups based on Logarithm of Odds (LOD ≥ 6) groupings. The seven linkage groups were in corresponding to the barley chromosomes according to the locus coordinates determined during read mapping to the reference genome of barley (Monat et al., 2019). Visualization of maps derived from JoinMap^® 4.1 was achieved by MapChart software (Voorrips, 2002).

Sequence 50 bp upstream and downstream of the SNP was extracted from the barley reference genome (Monat et al., 2019). Allele-specific forward primers and one common reverse primer were designed by help of the free assay design service of 3CR Bioscience³ (Supplementary Table 4). KASP primer mix was prepared in a volume of 100 μL containing 12 μM of each allele-specific forward primer and 30 μM common reverse primer. KASP-PCR reactions were performed in a total volume of 5 μL containing 2.5 μL 2X PACE-IR^TM Genotyping Master Mix (001-0010, 3CR Bioscience, Essex, United Kingdom), 0.07 μL KASP primer mix, and 40 ng of template DNA. PCR program was used with a HydroCycler (LGC, Teddington, United Kingdom): initial denaturation at 94°C for 15 min followed by 10 cycles at 94°C for 20 s, at 65 to 57°C (−0.8°C/cycle) for 1 min, and proceeded for 30 cycles at 94°C for 20 s, at 57°C for 1 min, and followed at 30°C for 30 s. In case of the genotyping clusters not well separated, an additional PCR was performed with 6 cycles at 94°C for 20 s, at 57°C for 1 min. Pre-read and post-read of the fluorescence signals were performed on ABI 7900HT instrument (Applied Biosystems, Thermo Fisher Scientific).

DNA was isolated according to the protocol of Doyle and Doyle (1990) from leaf samples collected from one-week-old, greenhouse-grown seedlings of barley cv. Haisa, albostrians mutant M4205, and lines luteostrians-P2_2 and luteostrians-P3_6. For lines luteostrians-P2_2 and luteostrians-P3_6, DNA was fragmented (200 - 500 bp) with ultrasounds, then 150 bp paired-end libraries were prepared according to Illumina standard protocol and sequenced with Illumina Hiseq X Ten. For Haisa and albostrians mutant M4205, 350 bp paired-end libraries were prepared and sequenced according to the protocol as described previously (International Barley Genome Sequencing Consortium, 2012). Reads mapping and variants calling to the barley reference genome were performed as above described for GBS data. Notably, re-sequencing data of barley cv. Barke, maternal parent of the BL mapping population (for luteostrians) and BM mapping population (for albostrians) (Li et al., 2019b) was also included for SNP calling for identification of mutant-specific SNPs. Functional effect of the mutant-specific SNPs was annotated by help of SNPeff version 4.3 (Cingolani et al., 2012).