Rare allele of HvLox-1 associated with lipoxygenase activity in barley (Hordeum vulgare L.)

Key message Identification and allele-specific marker development of a functional SNP of HvLox - 1 which associated with barley lipoxygenase activity. Abstract Improving the stability of the flavor of beer is one of the main objectives in breeding barley for malting, and lipoxygenase-1 (LOX-1) is a key enzyme controlling this trait. In this study, a modified LOX activity assay was used for null LOX-1 mutant screening. Four barley landraces with no detected level of LOX-1 activity were screened from 1,083 barley germplasm accessions from China. The genomic sequence diversity of the HvLox-1 gene of the four null LOX-1 Chinese landraces was compared with that of a further 76 accessions. A total of 104 nucleotide polymorphisms were found, which contained 83 single-nucleotide polymorphisms (SNPs), 7 multiple-nucleotide polymorphisms, and 14 insertions and deletions. Most notably, we found a rare C/G mutation (SNP-61) in the second intron which led to null LOX-1 activity through an altered splicing acceptor site. In addition, an allele-specific polymerase chain reaction marker was developed for the genotyping of SNP-61, which could be used in breeding programs for barley to be used for malting. The objective was to improve beer quality. Electronic supplementary material The online version of this article (doi:10.1007/s00122-014-2362-3) contains supplementary material, which is available to authorized users.


Introduction
The metabolism products of lipid degradation seriously affect the quality of cereals and their processed products during storage (Garbus et al. 2009;Wang et al. 2008;Zhang et al. 2007). In the lipid degradation pathway, the polyunsaturated fatty acid substrate is converted to 9-hydroperoxide (9-HPOD) and 13-hydroperoxide (13-HPOD) by the dioxygenase effect of lipoxygenase (Kuroda et al. 2003;Kuroda et al. 2005). Subsequently, 9-HPOD is converted to trans-2-nonenal (T2N) and trihydroxy octadecenoic acids (THODs) by HvHPL2 and its isoenzymes. At the same time, 13-HPOD is metabolized to hexanal (Feussner and Wasternack 2002). T2N was identified as being responsible for the cardboard-like flavor in beer (Hambraeus and Nyberg 2005), and it is generated and accumulated during beer storage Vanderhaegen et al. 2006). In addition, THOD also has an adverse effect on foam stability and the flavor of beer (Kaneda et al. 2001; 1 3 Kobayashi et al. 2002). Therefore, there has been extensive screening of artificial or natural variations for null-or reduced-LOX barley mutants (Hirota et al. 2005;Oozeki et al. 2007; Patent US7420105 B2).
As a key enzyme in the conversion of linoleic acid to THODs in the lipid degradation pathway, LOX is encoded by a gene family in plants. For example, three Lox genes were identified in barley-Lox-1 (LoxA), Lox-2 (LoxC), and Lox-3 (LoxB) (Rouster et al. 1997(Rouster et al. , 1998van Mechelen et al. 1999;van Mechelen et al. 1995). In beer production, LOX-1 catalyzes the formation of 9-HPOD while LOX-2 catalyzes the formation of 13-HPOD. LOX-3 is considered irrelevant to beer quality given its low expression level in barley kernels and unclear product specificity (van Mechelen et al. 1999). Collectively, LOX-1 is considered the main contributor to malt LOX activity (Doderer et al. 1992;Kuroda et al. 2003;Yang and Schwarz 1995). To date, several artificial and natural mutants of null LOX-1 activity have been screened out and used for the breeding of flavor-stable malting barley. Forced aging tests on the null LOX beer produced also demonstrated its improved flavor stability through reduced T2N content (Hirota et al. 2005(Hirota et al. , 2006. Genetic improvement, extending the stability of beer flavor and its shelf life has been one of the main targets in breeding barley for malting. We sought to identify the related gene and develop the gene's molecular markers for the quality characteristics useful to barley breeders. Thus, in the present study, 1,083 barley landrace accessions were screened for null LOX-1 activity mutants. Additionally, 77 barley landraces, including undetectable LOX activity mutants and randomly selected genotypes, were sequenced and analyzed for the diverse alleles/haplotypes and key single-nucleotide polymorphisms (SNPs) affecting LOX-1 activity. The mechanism of LOX-1 activity loss in the mutant is also discussed.

Plant materials
In the process of screening for null LOX-1 activity barley, 1,083 barley landraces from the National Genebank of China were purified by a single plant harvest and individually milled with a Perten™ laboratory mill (PLM3100/C, 0.8 mm standard sieve). Crude homogenates were prepared by mixing and vigorously vortexing 0.1 g of barley powder in 1 mL of ice cold water in a 1.5 mL Eppendorf tube, and then centrifuging at 3,000 rpm for 2 min at 4 °C. Then, 50 μL of supernatant was aspirated for the lipoxygenase enzyme activity assay.

Lipoxygenase assay
Chemical preparation and the activity analysis of the linoleic hydroperoxide reaction of LOX were performed using a modified DMAB-MBTH assay as described in Anthon and Barrett (2001). For the modified rapid assay, a fresh working solution A was prepared by mixing 10 mL of DMAB solution (20 mM DMAB, 100 mM phosphate buffer, and pH 6.0), 0.4 mL of linoleic acid substrate (25 mM linoleic acid stock), and 9.6 mL of water. Solution B was also freshly prepared by mixing 0.4 mL of MBTH (10 mM), 0.4 mL of hemoglobin (5 mg/mL), and 19.2 mL of water. All assays of linoleic hydroperoxide were carried out at room temperature (20-25 °C). The sample was incubated initially with 0.5 mL of solution A for 20 min and 0.5 mL of solution B was added subsequently. After an additional 10 min, the reaction was terminated by adding 0.5 mL of 1 % (w/v) sodium lauryl sulfate (SLS). Blank tests were performed by adding SLS before the reaction. For samples with normal LOX-1 activity, stable dark blue products were quickly produced. The presence or absence of dark blue products was used as a rapid visual determination of LOX-1 activity.

Nucleotide extraction
Barley plants were grown in an illuminated incubator under the following conditions: 22 °C for 16 h of light and 18 °C for 8 h of darkness. Two-week-old seedlings were harvested individually and frozen in liquid nitrogen and then stored at −80 °C for DNA and RNA extraction. DNA and RNA were isolated and purified with the DNeasy™ and RNeasy™ plant mini kit (Qiagen) according to the manufacturer's instructions.
A total of 77 barley germplasms were sequenced for allele mining and haplotype analysis. These included 73 LOX-1 enzyme activity positive lines from the major barley-growing regions of China and four null LOX-1 landraces (ZDM00279, ZDM00626, ZDM00677, and ZDM05300). Moreover, all null LOX-1 lines were confirmed with a second sequencing. Details of accession names, their haplotypes, row type, and origin are listed in Table S2.
The amplification of the HvLox-1 gene fragment was performed using Pfu-DNA polymerase (TIANGEN, China). The PCR conditions were 35 cycles at 95 °C for 30 s, 58 °C for 30 s, and 72 °C for 2 min, followed by a 7-min extension at 72 °C. The PCR product was purified and sequenced directly using an ABI3730 sequencer. The haplotypes corresponding to full-length cDNA of HvLox-1 from eight lines (ZDM00470, ZDM00797, ZDM03953, ZDM03595, ZDM00279, ZDM00626, ZDM00677, and ZDM05300) were amplified and sequenced. Nucleotide sequences were assembled with PHRAP software (Bastide and McCombie 2007) and manually edited using BioEdit 7.0 (Hall 1999). Multiple-sequence alignments were carried out with the ClustalW v2.1 program (Larkin et al. 2007). The allelic haplotypes and sliding-window plot (window length 100 bp, step size 25 bp) of nucleotide diversity (π) of the HvLox-1 gene were defined and generated by DnaSP 5.10 (Librado and Rozas 2009). A neighbor-joining (NJ) tree of seven haplotypes of HvLox-1 was constructed using MEGA 6.0 with a 1,000 times bootstrap (Tamura et al. 2013). All nucleotide variations, amino acids, and predicted regulatory motif changes are shown in Table S3.
Allele-specific PCR for functional SNP genotyping An allele-specific PCR fragment length polymorphism analysis (AS-PCR) assay was used for SNP validation. Two pairs of primers were developed for the identification and genotyping of SNP-61 as described in Wangkumhang et al. (2007). A pair of primers, RNF/Lox1.FL.1R, was designed for 'C' allele genotyping with a 320 bp specific PCR amplicon. The primer pairs Lox1.FL.2F and RNR were able to discriminate the 'G' allele from null LOX-1 barley lines for which a 720-bp specific PCR product was generated. A 10 bp (CACTAGTGAT), randomly selected flanking sequence from a common cloning vector was added to the 5′-end of the RNR primer to improve PCR stability. All these primer pairs were multiplexed in a single-tube PCR assay to assess the allelic status at SNP-61 (Table S1; Fig. 4a). The amplification of the allele-specific gene fragment was performed under PCR conditions of 35 cycles at 95 °C for 30 s, 65 °C for 30 s, and 72 °C for 30 s; followed by a 5-min extension at 72 °C. The PCR products were separated using a 1 % agarose gel.

Characterization of LOX-1 activity in Chinese barley landraces
A total of 1,083 Chinese barley (Hordeum vulgare L.) landraces were analyzed using a rapid visual determination of LOX-1 activity assay. Four barley landraces (ZDM00279, ZDM00626, ZDM00677, and ZDM05300) showed no significant LOX-1 activity in silenced seeds (Fig. 1), which was further confirmed by a malted barley analysis (data not shown). The four barley landraces without LOX-1 activity originated from Shan-Dong and He-Nan Provinces of China (Table S2). The nucleotide variation analysis of HvLox-1 was carried out using the contig sequences of three non-Chinese barley cultivars (cvs Morex, Bowman, and Barke) and 77 sequenced Chinese barley landraces (including four null LOX-1 activity lines). A total of 104 varied DNA polymorphisms-including 83 SNPs, 7 multiple-nucleotide polymorphisms (MNPs) and 14 insertions and deletions (InDels)-were observed in the 5,884-bp genomic region of HvLox-1, with an average frequency of ~17 DNA polymorphisms for every kilo base-pair (π = 0.01209). The nucleotide diversity in the non-coding region was much higher than in the exons. There were 23 nucleotide polymorphisms in the exons (2,589 bp) and 43 in the introns (1,464 bp) of HvLox-1 (Table S3). Six out of seven multinucleotide polymorphisms (MNPs) and all of the 14 InDels were located in the non-coding region. Also, nine non-synonymous SNPs were located in three exons (exons 1, 3, and 7) ( Fig. 2a; Table 1). In addition, the second intron of the gene contained 30 nucleotide variations in the 634 bp region, with about half of these being concentrated in a 100 bp region (1,500-1,600 bp) ( Fig. 2a; Table 2). Slidingwindow analysis of the polymorphism pattern across the HvLox-1 gene also revealed that the nucleotide diversity in the coding region was far less than in the exons and promoter region. Notably, the peak of nucleotide polymorphism (pi) was located at the second intron of the HvLox-1 gene, which demonstrated that intron 2 is the high-density area of DNA polymorphisms (Fig. 3a).
A total of 34 nucleotide polymorphisms, containing 29 SNPs, three MNPs and two InDels, were found in the  Table 1). The nucleotide variations of intron 2 are shown below (corresponding to Table 2). b SNP-61 in the second intron causing alternative splicing of HvLox-1. SNP-61 affecting an intron-exon splice junction and resulting in a premature termination (5-8) and frame-shift (6-8) in the coding sequence. The positions of the SNP-61 natural mutant alleles in the genomic DNA and processed cDNA sequence alignment and two types of alternative splicing events are shown below the gene structure. 1-8: as for Fig. 1. c Alternative splicing induced cDNA length polymorphism of HvLox-1. 1-8: as for Fig. 1. In contrast to the normal splicing form of intron 2 in 1-4 (H2-H5), 5 (H6) is the only AltA-1 event shown, and 6-8 (H7) are shown as the two types of alternative splicing events, AltA-1 and AltA-2 (color figure online)     (Table S3).

No (F)
Haplotype analysis and identification of critical variation associated with null LOX-1 activity The 80 barley accessions analyzed contained seven haplotypes according to the detected nucleotide polymorphisms. The American and European cvs Morex, Bowman, and Barke were defined as haplotype H1. Six haplotypes (H2-H7) were constructed from the remaining 77 Chinese barley landrace lines. A total of 51 landraces were identified as haplotype H2; these landraces were widely distributed in the major barley-growing regions of China. Five landraces were assigned to haplotype H3, eight to H4, and nine to H5. Of the four barley accessions with null LOX-1 activity, ZDM00279 was assigned to H6 and ZDM00626, ZDM00677, and ZDM05300 to H7. These latter shared the same SNP (SNP-61, C/G) and originated from He-Nan and Shan-Dong Provinces (Table S2). Haplotypes H2, H3, and H6 partly contained H1 SNP alleles and were quite distinct from haplotypes H4, H5, and H7, which showed completely different SNP alleles to H1 (Table S3) and were grouped into one phylogenetic sub-tree (Fig. 3b). Thus, SNP-61 was considered as the candidate variation that may affect the function of HvLOX-1 protein. Two clades were detected through a phylogenetic analysis of HvLox-1 (Fig. 3b), which corresponded with the promoter analysis result and with the polyphyletic domestication division of cultivated barley (Dai et al. 2012).

Discussion
Barley HvLOX-1 involvement in lipid degradation and the formation of beer-deteriorating substances have been studied previously (Kuroda et al. 2003). Furthermore, null LOX beer has become a highly desirable product, since it has significant stability of both flavor and foam, which allows the beer to remain fresh for longer. So far, several natural Fig. 3 Sliding-window analysis and neighbor-joining tree of the seven haplotypes of HvLox-1. a Sliding-window representation of the HvLox-1 gene polymorphism sites, which correspond to the exonintron structure of the seven haplotypes, was executed using DnaSP 5.10, the window length was 100 bp, and the step size was 25 bp. b A neighbor-joining (NJ) tree of seven haplotypes of the HvLox-1. The NJ tree of the seven haplotypes was developed using the neighborjoining algorithm of MEGA 6.0 and artificial null LOX-1 barley lines have been identified. Six LOX-1 inactive barley landrace lines which shared the unique polymorphism in the HvLox-1 structural gene region were screened from 1,153 lines (Hirota et al. 2005). A LOX-deficient mutant line was identified from a sodium azide (NaN 3 )-treated malting cultivar. Sequence analysis revealed that a C/T substitution introduced an internal stop codon into the third exon of HvLox-1 (Oozeki et al. 2007). In addition, two null LOX-1 (A618 and D112) and one low LOX-1 (line G) mutant lines were also screened from generation M 3 of a malting barley cultivar. It was believed that a splicing error and/or translational stop caused by nucleotide substitution resulted in a null LOX-1 phenotype (Patent US7420105 B2). Both artificial and natural null LOX-1 mutants are used for breeding flavor-stable malting barley for the brewing industries. Therefore, the extensive screening of natural null LOX-1 barley mutants is indispensable to high-quality barley breeding and to understanding the genetic variation of HvLox-1. In the present study, four null LOX-1 landrace lines were screened from 1,083 Chinese barley germplasm accessions. Interestingly, these four lines originated from Shan-Dong and He-Nan Provinces in the lower reaches of the Yellow River.
Sequence diversity analysis showed that there were more nucleotide polymorphisms in the intron regions (43) than in the promoter (34) and exons (23) ones of HvLox-1. Based on the 34 nucleotide polymorphisms discovered in the promoter region of HvLox-1, Chinese barley landraces and American and European barley cultivars are clearly divided into two types. SNP changes between promoters indicated that HvLox-1 evolved from two distinct ancestral gene pools and independent domestication processes as previously reported (Dai et al. 2012). Notably, the second intron contained the most abundant nucleotide variations, and SNP-61 was located in this region (Table 2). A nonsynonymous amino acid substitution in the third exon of the HvLox-1 gene in Line G has contributed to its low LOX-1 enzyme activity phenotype (Patent US7420105 B2). Thus, the unconventional amino acid changes we discovered are likely to be functional properties of the enzyme. Since the modified rapid lipoxygenase assay is quite suitable for the qualitative detection of LOX-1 activity, the potential impact of amino acid substitutions of different haplotypes and LOX-1 activity remains unclear. A further study of a quantitative lipoxygenase assay is required to understand the relationship between non-synonymous mutation and the LOX activity.
Simple, highly reproducible, and co-dominant molecular markers derived from polymorphic sites within functional genes are essential for marker-assisted selection (Geng et al. Fig. 4 Allele-specific PCR marker developing and genotyping. a Schematic illustration of the development of allele-specific PCR primers for SNP-61. A pair of primers (RNF/Lox1.FL.1R) was designed for 'C' allele genotyping with a 320-bp specific PCR amplicon. The primer pairs Lox1.FL.2F and RNR were able to discriminate the 'G' allele from the null-LOX1 barley lines for which the primer pairs were generated as a 720-bp specific PCR product. All these primer pairs were multiplexed in a single-tube PCR assay to assess the allelic status at SNP-61. b Gel analysis can clearly distinguish the allelic status of SNP-61. Lanes 1-8 are the same as Fig. 1. Two arrows indicate the 720-bp (G allele) and 320-bp (C allele) specific PCR products separately; M DL2000 ladder 2012; Hirota et al. 2005). Cleaved amplified polymorphic sequence markers and allele-specific PCR fragment length polymorphism analysis (AS-PCR) are commonly used for locus-specific SNP identification (Su et al. 2011;Yang et al. 2012). In the present study, AS-PCR primers were developed for the identification of null LOX-1 barley lines. This simple method would be effective in direct marker-based selection of the null LOX-1 gene in breeding programs to create advanced null LOX-1 malting barley varieties. Finally, this work demonstrates that quality-related mutations and other valuable rare allele identifications in beer barley germplasm from different sources have considerable potential.
As a common phenomenon and an important post-transcriptional regulatory mechanism in plants, alternative splicing results in a gain or loss of protein function. The different forms of translated protein are often manifested as changes in cellular localization, stability, and activity (Barbazuk et al. 2008;Ner-Gaon et al. 2004;Wang and Brendel 2006). One of the earliest examples of alternative splicing demonstrated in plants was for Rubisco (Werneke et al. 1989). Genome-wide analyses of alternative splicing in plants were performed and reviewed in subsequent studies (Barbazuk et al. 2008;International Barley Genome Sequencing Consortium 2012;Ner-Gaon et al. 2004). Extensive alternative splicing is one notable feature of post-transcriptional processing in the barley genome. About 73 % of the introncontaining high-confidence barley genes show evidence of alternative splicing (International Barley Genome Sequencing Consortium 2012). In addition, it was reported that strict 5′-and 3′-splice site sequences of plant introns were essential to splicing factors (Reddy 2007). Thus, further investigation of the mechanism of SNP-61 inducing re-organization changes in the second intron splice site is required. It was reported that a single-nucleotide mutation altered the splicing donor site (AltD) of the fifth intron, which caused LOX deficiency in six barley lines with no significant LOX activity (Hirota et al. 2005). In the present study, a base substitution (SNP-61, C/G) in the second intron of HvLox-1 was found to be associated with the null LOX-1 phenotype of the four barley landrace lines. However, further analysis of cDNA sequences identified that SNP-61 (C/G) resulted in an alternative splicing acceptor site (AltA) in the second intron of HvLox-1, which caused the LOX-1 deficiency by frame-shift and premature termination in translation (Fig. 2b, c). Since null LOX-1 haplotypes H6 and H7 were not grouped into one phylogenetic sub-tree (Fig. 3b), it was probable that a recombination event happened before and split a single mutation origin into two haplotypes. Although this finding differs from those of previous reports, given the materials used in the studies, the events of post-transcriptional processing observed in the different studies suggest that alternative splicing is an important and extensive phenomenon in barley.