# Towards reducing the immunogenic potential of wheat flour: omega gliadins encoded by the D genome of hexaploid wheat may also harbor epitopes for the serious food allergy WDEIA

Susan B. Altenbach, Han-Chang Chang, Annamaria Simon-Buss, You-Ran Jang, Sandra Denery Papini, Florence Pineau, Yong Q. Gu, Naxin Huo, Sun-Hyung Lim, Chon-Sik Kang, et al.

► To cite this version:

Susan B. Altenbach, Han-Chang Chang, Annamaria Simon-Buss, You-Ran Jang, Sandra Denery Papini, et al.. Towards reducing the immunogenic potential of wheat flour: omega gliadins encoded by the D genome of hexaploid wheat may also harbor epitopes for the serious food allergy WDEIA. BMC Plant Biology, 2018, 18, pp.1-12. 10.1186/s12870-018-1506-z. [https://hal.inrae.fr/hal-02626394](https://hal.inrae.fr/hal-02626394)

HAL Id: hal-02626394

https://hal.inrae.fr/hal-02626394

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH ARTICLE

Open Access

<figure><img src="image_2.png" /></figure>

CrossMark

# Towards reducing the immunogenic potential of wheat flour: omega gliadins encoded by the D genome of hexaploid wheat may also harbor epitopes for the serious food allergy WDEIA

Susan B. Altenbach¹, Han-Chang Chang¹, Annamaria Simon-Buss¹, You-Ran Jang<fnref n="2" />, Sandra Denery-Papini³, Florence Pineau³, Yong Q. Gu¹, Naxin Huo¹, Sun-Hyung Lim<fnref n="2" />, Chon-Sik Kang⁴ and Jong-Yeol Lee<fnref n="2" />*

## Abstract

**Background:** Omega-5 gliadins are a group of highly repetitive gluten proteins in wheat flour encoded on the 1B chromosome of hexaploid wheat. These proteins are the major sensitizing allergens in a severe form of food allergy called wheat-dependent exercise-induced anaphylaxis (WDEIA). The elimination of omega-5 gliadins from wheat flour through biotechnology or breeding approaches could reduce the immunogenic potential and adverse health effects of the flour.

**Results:** A mutant line missing low-molecular weight glutenin subunits encoded at the *Glu-B3* locus was selected previously from a doubled haploid population generated from two Korean wheat cultivars. Analysis of flour from the mutant line by 2-dimensional gel electrophoresis coupled with tandem mass spectrometry revealed that the omega-5 gliadins and several gamma gliadins encoded by the closely linked *Gli-B1* locus were also missing as a result of a deletion of at least 5.8 Mb of chromosome 1B. Two-dimensional immunoblot analysis of flour proteins using sera from WDEIA patients showed reduced IgE reactivity in the mutant relative to the parental lines due to the absence of the major omega-5 gliadins. However, two minor proteins showed strong reactivity to patient sera in both the parental and the mutant lines and also reacted with a monoclonal antibody against omega-5 gliadin. Analysis of the two minor reactive proteins by mass spectrometry revealed that both proteins correspond to omega-5 gliadin genes encoded on chromosome 1D that were thought previously to be pseudogenes.

**Conclusions:** While breeding approaches can be used to reduce the levels of the highly immunogenic omega-5 gliadins in wheat flour, these approaches are complicated by the genetic linkage of different classes of gluten protein genes and the finding that omega-5 gliadins may be encoded on more than one chromosome. The work illustrates the importance of detailed knowledge about the genomic regions harboring the major gluten protein genes in individual wheat cultivars for future efforts aimed at reducing the immunogenic potential of wheat flour.

**Keywords:** Food allergy, Omega-5 gliadins, Proteomics, Reduced immunogenic potential, Wheat flour

<footnote>* Correspondence: jy0820@korea.kr</footnote>

<footnote n="2">²National Institute of Agricultural Sciences, RDA, Jeonju 54874, South Korea
Full list of author information is available at the end of the article</footnote>
However, the analysis revealed that several minor but highly immunogenic omega-5 gliadins remained in this line. The work highlights some of the challenges faced in using breeding approaches to reduce the immuno- genic potential of wheat flour.

## Results

### 2-DE analysis of total flour proteins

Analyses of total flour proteins from the parental lines Keumkang and Olgeuru and the mutant line DH20 by two-dimensional gel electrophoresis (2-DE) revealed that a number of gluten proteins between ~34 and 55 kDa are missing in the mutant (Fig. 1). As documented previously, these include a number of LMW-GS encoded by the *Glu-B3* locus [18]. Spots 7–11 in Keumkang (Fig. 1d), identified by MS/MS as s-type LMW-GS ACA63869 and ACA63865 (Table 1, Additional file 1), and spots 8–10 in

Olgeuru (Fig. 1e), identified as s-type LMW-GS ACA63864 and ACA63868 (Table 1, Additional file 1), are missing in the mutant line DH20 (Fig. 1f). These proteins have the N-terminal sequences SHIPGLERP and differ by only a few amino acids. In addition, DH20 is missing an m-type LMW-GS also encoded by the *Glu-B3* locus. The LMW-GS ACA63871 was identified in spot 17 in Keumkang (Fig. 1d) and the LMW-GS ACP27643 was identified in spot 16 in Olgeuru (Fig. 1e). These two proteins begin with the N-terminal sequences METSHIP and differ by only 5 amino acids. The DH20 mutant contained the i-type LMW-GS AAS10189 (spot 1 in Fig. 1f) contributed by the Keumkang parent (spot 6 in Fig. 1d), confirming the presence of the *Glu-A3c* allele in the mutant line. AAS10189 was also identified in spot 2 in DH20 and was a minor component of Spot 7 in Keumkang (Fig. 1d, Additional file 1). All three lines

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Total flour proteins from parental lines Keumkang (a, d) and Olgeuru (b, e) and mutant line DH20 (c, f). Protein spots within the black boxes in panels (a, b and c) were identified by MS/MS. Proteins circled in red, blue, yellow, purple, green and black in panels (d, e and f) were identified as omega-5 gliadins, s-type LMW-GS, i-type LMW-GS, m-type LMW-GS, gamma gliadins and triticins, respectively. Dashed circles denote proteins identified in at least one of the parental lines but absent in the mutant line. MS/MS data are summarized in Table 1 and Additional file 1</figcaption></figure>
<table><caption>**Table 1** Identification of 2-DE spots from parental and mutant lines by MS/MS</caption><thead><tr><th>Cultivar</th><th>Spot #</th><th>Predominant Protein</th><th>Accession Number</th><th># Exclusive Unique Peptides</th><th># Spectra</th><th>% Coverage</th></tr></thead><tbody><tr><td>Keumkang</td><td>1</td><td>omega-5 gliadin</td><td>CAR82267</td><td>16</td><td>47</td><td>29</td></tr><tr><td>Keumkang</td><td>2</td><td>omega-5 gliadin</td><td>BAE20328</td><td>16</td><td>48</td><td>25</td></tr><tr><td>Keumkang</td><td>3</td><td>omega-5 gliadin</td><td>BAE20328</td><td>28</td><td>59</td><td>38</td></tr><tr><td>Keumkang</td><td>4</td><td>omega-5 gliadin</td><td>BAE20328</td><td>6</td><td>12</td><td>19</td></tr><tr><td>Keumkang</td><td>5</td><td>omega-5 gliadin</td><td>BAE20328</td><td>15</td><td>41</td><td>40</td></tr><tr><td>Keumkang</td><td>6</td><td>LMW-GS (i-type)<sup>a</sup></td><td>AAS10189</td><td>59</td><td>180</td><td>81</td></tr><tr><td>Keumkang</td><td>7</td><td>LMW-GS (s-type)<sup>b</sup></td><td>ACA63869</td><td>67</td><td>191</td><td>73</td></tr><tr><td>Keumkang</td><td>8</td><td>LMW-GS (s-type)<sup>b</sup></td><td>ACA63869</td><td>64</td><td>211</td><td>72</td></tr><tr><td>Keumkang</td><td>9</td><td>LMW-GS (s-type)<sup>b</sup></td><td>ACA63869</td><td>90</td><td>305</td><td>76</td></tr><tr><td>Keumkang</td><td>10</td><td>LMW-GS (s-type)<sup>b</sup></td><td>ACA63869</td><td>130</td><td>443</td><td>88</td></tr><tr><td>Keumkang</td><td>11</td><td>LMW-GS (s-type)<sup>b</sup></td><td>ACA63865</td><td>90</td><td>313</td><td>78</td></tr><tr><td>Keumkang</td><td>12</td><td>triticin</td><td>ACB41345</td><td>30</td><td>74</td><td>29</td></tr><tr><td>Keumkang</td><td>13</td><td>gamma gliadin</td><td>Bu-5<sup>d</sup></td><td>53</td><td>177</td><td>60</td></tr><tr><td>Keumkang</td><td>14</td><td>gamma gliadin</td><td>Bu-5<sup>d</sup></td><td>47</td><td>164</td><td>61</td></tr><tr><td>Keumkang</td><td>15</td><td>gamma gliadin</td><td>AGJ50340</td><td>46</td><td>164</td><td>62</td></tr><tr><td>Keumkang</td><td>16</td><td>LMW-GS (m-type)<sup>c</sup></td><td>ABC84361</td><td>66</td><td>179</td><td>74</td></tr><tr><td>Keumkang</td><td>17</td><td>LMW-GS (m-type)<sup>c</sup></td><td>ACA63871</td><td>70</td><td>357</td><td>84</td></tr><tr><td>Olgeuru</td><td>1</td><td>omega-5</td><td>AIU64835</td><td>7</td><td>26</td><td>48</td></tr><tr><td>Olgeuru</td><td>2</td><td>omega-5</td><td>All26682</td><td>6</td><td>15</td><td>17</td></tr><tr><td>Olgeuru</td><td>3</td><td>omega-5</td><td>All26682</td><td>4</td><td>15</td><td>16</td></tr><tr><td>Olgeuru</td><td>4</td><td>omega-5</td><td>CAR82266</td><td>6</td><td>13</td><td>25</td></tr><tr><td>Olgeuru</td><td>5</td><td>omega-5</td><td>BAE20328</td><td>36</td><td>75</td><td>45</td></tr><tr><td>Olgeuru</td><td>6</td><td>omega-5</td><td>All26682</td><td>9</td><td>25</td><td>34</td></tr><tr><td>Olgeuru</td><td>7</td><td>LMW-GS (i-type)<sup>a</sup></td><td>BAB78764</td><td>60</td><td>175</td><td>77</td></tr><tr><td>Olgeuru</td><td>8</td><td>LMW-GS (s-type)<sup>b</sup></td><td>ACA63864</td><td>40</td><td>110</td><td>54</td></tr><tr><td>Olgeuru</td><td>9</td><td>LMW-GS (s-type)<sup>b</sup></td><td>ACA63868</td><td>44</td><td>124</td><td>62</td></tr><tr><td>Olgeuru</td><td>10</td><td>LMW-GS (s-type)<sup>b</sup></td><td>ACA63864</td><td>92</td><td>310</td><td>73</td></tr><tr><td>Olgeuru</td><td>11</td><td>LMW-GS (i-type)<sup>c</sup></td><td>ACY08811</td><td>47</td><td>111</td><td>69</td></tr><tr><td>Olgeuru</td><td>12</td><td>triticin</td><td>ACB41345</td><td>24</td><td>60</td><td>27</td></tr><tr><td>Olgeuru</td><td>13</td><td>gamma gliadin</td><td>AGO17716</td><td>23</td><td>51</td><td>45</td></tr><tr><td>Olgeuru</td><td>14</td><td>gamma gliadin</td><td>AGO17716</td><td>31</td><td>74</td><td>57</td></tr><tr><td>Olgeuru</td><td>15</td><td>LMW-GS (m-type)<sup>c</sup></td><td>ABC84361</td><td>60</td><td>164</td><td>74</td></tr><tr><td>Olgeuru</td><td>16</td><td>LMW-GS (m-type)<sup>c</sup></td><td>ACP27643</td><td>23</td><td>118</td><td>51</td></tr><tr><td>DH20</td><td>1</td><td>LMW-GS (i-type)<sup>a</sup></td><td>AAS10189</td><td>104</td><td>310</td><td>83</td></tr><tr><td>DH20</td><td>2</td><td>LMW-GS (i-type)<sup>a</sup></td><td>AAS10189</td><td>30</td><td>98</td><td>70</td></tr><tr><td>DH20</td><td>3</td><td>triticin</td><td>ACB41345</td><td>39</td><td>129</td><td>37</td></tr><tr><td>DH20</td><td>4</td><td>gamma gliadin</td><td>AGJ50340</td><td>49</td><td>161</td><td>58</td></tr><tr><td>DH20</td><td>5</td><td>LMW-GS (m-type)<sup>c</sup></td><td>ABC84361</td><td>60</td><td>164</td><td>74</td></tr></tbody></table>

contained the same m-type LMW-GS ABC84361 beginning with the N-terminal sequence METSRV that is contributed by the *Glu-D3a* allele (spot 16 in

Keumkang (Fig. 1d), spot 15 in Olgeuru (Fig. 1e) and spot 5 in DH20 (Fig. 1f)). These spots serve as references among the gels.

<footnote><sup>a</sup>refers to LMW-GS that begin with the N-terminal sequence ISQQQ</footnote>

<footnote><sup>b</sup>refers to LMW-GS that begin with the N-terminal sequence SHIP</footnote>

<footnote><sup>c</sup>refers to LMW-GS that begin with M</footnote>

<footnote><sup>d</sup>best match is to a sequence from cv. Butte 86 reported in [19]</footnote>
Most notably, the omega-5 gliadins were absent from the mutant line. In Keumkang, spots 1–5 were identified as omega-5 gliadins CAR82267 or BAE20328 (Fig. 1d, Table 1, Additional file 1) while spots 1–6 in Olgeuru were identified as CAR82267, BAE20328 or the partial sequences AIU64835 or AII26682 (Fig. 1e, Table 1, Additional file 1). All proteins in these regions of the 2-D gels were missing in DH20 (Fig. 1f).

In addition, several gamma gliadins were missing from DH20, including spots 13 and 14 from Keumkang (Fig. 1d) and 13 and 14 from Olgeuru (Fig. 1e). Both spots in Keumkang were identified as Bu-5, a gamma gliadin from the U.S. cultivar Butte 86 [19], and both spots in Olgeuru were identified as AGO17716. The two proteins differ by only 8 amino acids. Spots 15 in Keumkang and 4 in DH20, identified as gamma gliadin AGJ50340 along with spots 12 in both Keumkang and Olgeuru and spot 3 in DH20, identified as triticin, serve as reference spots among the gels.

## PCR analysis of parental and mutant lines

PCR analysis confirmed that a number of genes were deleted in DH20 (Fig. 2). PCR using primers specific for the *Glu-B3-4* LMW-GS gene (Table 2) yielded a fragment of ~1600 bp from genomic DNA of Keumkang, Olgeuru, and the reference wheat Chinese Spring, but not from the DH20 mutant or the N1BT1A and N1BT1D nullisomic tetrasomic lines of Chinese Spring (Fig. 2a). Likewise, the expected 1200 bp fragment was obtained from Keumkang, Olgeuru and Chinese Spring using primers specific for the omega-5 gliadin gene (Table 2), but not from DH20 or the N1BT1A and N1BT1D nullisomic tetrasomic lines of Chinese Spring (Fig. 2b), confirming the location of the omega-5 gliadin

genes on chromosome 1B and the deletion of these genes in DH20.

To determine the extent of the deletion in DH20, repeat junction primers were selected based on the recently published genomic sequence of a 6.5 Mb region of chromosome 1B from Chinese Spring that contains LMW-GS, omega-5 gliadin and gamma gliadin genes [11] (Table 2). These primers amplified the expected size fragments from the parental lines Keumkang and Olgeuru as well as Chinese Spring, but did not amplify DH20 or the N1BT1A and N1BT1D nullisomic tetrasomic Chinese Spring lines, suggesting that, at a minimum, a 5.8 Mb region of the 1B chromosome had been deleted in DH20.

## Immunogenic potential of parental and mutant lines

Since the omega-5 gliadins are the major sensitizing proteins in the serious food allergy WDEIA, the immunogenic potential of DH20 was evaluated by 2-D immunoblot analysis using sera from patients with verified cases of WDEIA. Patient sera were characterized previously by ELISA using individual gluten protein fractions and by 2D immunoblot analysis [16]. In the two parental lines, proteins identified as omega-5 gliadins showed the greatest amounts of IgE reactivity (Fig. 3a, b, d, e). In addition, there was some reactivity with two protein spots with molecular weights slightly less than the omega-5 gliadins as well as with some proteins in the LMW-GS and alpha gliadin regions of the 2D gels. IgE reactivity was reduced markedly in DH20 (Fig. 3c, f). There was no reactivity in the region of the omega-5 gliadins. However, the two protein spots just below the omega-5 gliadin region, indicated with arrows in Fig. 3, showed reactivity with the sera from both WDEIA patients. These spots correspond to very minor proteins on the Coomassie-stained gels (Fig. 1).

<figure><img src="image_4.png" /><figcaption>**Fig. 2** PCR analysis with primers specific for LMW-GS encoded at the *Glu-B3* locus (a), omega-5 gliadins (b), and repeat junction primers 19S-1.3-2 (c) and 143E-1-600 (d) located at the ends of a 5.8 Mb region of chromosome 1B in Chinese Spring. In each panel, genomic DNA from Keumkang (1), Olgeuru (2), DH20 (3), N1BT1A (4) and N1BT1D (5) nullisomic tetrasomic lines of Chinese Spring or Chinese Spring (6) was amplified. The sizes of molecular weight markers in kb are shown in lane M in each panel. Primer sequences are shown in Table 2</figcaption></figure>
<table><caption>**Table 2** PCR primers and amplification conditions used in this study</caption><thead><tr><th>Target</th><th>Primer set</th><th>Sequence (5' - 3')</th><th>Amplicon size</th><th>PCR conditions</th><th>Reference</th></tr></thead><tbody><tr><td rowspan="2">GluB3-4<br>LMW-GS</td><td>LB4F</td><td>CACCCTATACAAGGTCCAAAAT</td><td rowspan="2">1636</td><td rowspan="2">94 °C 60s, 60 °C 60s, 72 °C 60s</td><td rowspan="2">[27]</td></tr><tr><td>LB4R</td><td>TATTTCCATAATTTAAACTAGTTTGT</td></tr><tr><td rowspan="2">Omega-5 gliadin</td><td>O5GF</td><td>AGTAGGCTGCTAAGCCCTAGA</td><td rowspan="2">1212</td><td rowspan="2">94 °C 30s, 60 °C 45s, 72 °C 60s</td><td rowspan="2">This study</td></tr><tr><td>O5GR</td><td>ATATTGTTGGTATGGGGAAGG</td></tr><tr><td rowspan="2">Chromosome 1B Repeat<br>Junction Primers</td><td>19S-1.3-2F</td><td>CGGGACTCATTGAAGAATCC</td><td rowspan="2">418</td><td rowspan="2">94 °C 30s, 55 °C 45s, 72 °C 60s</td><td rowspan="2">This study</td></tr><tr><td>19S-1.3-2R</td><td>ACTAAAGTAGCGATTCAAATCCTC</td></tr><tr><td rowspan="2">Chromosome 1B Repeat<br>Junction Primers</td><td>143E-1-600F</td><td>CAACGACATATTCTAACCTCCACA</td><td rowspan="2">482</td><td rowspan="2">94 °C 30s, 60 °C 45s, 72 °C 60s</td><td rowspan="2">This study</td></tr><tr><td>143E-1-600R</td><td>ACATACATACATCACTGACCGACTG</td></tr><tr><td rowspan="2">Omega-D4</td><td>OD4-F1</td><td>ACTAGGCAACTAAGCCCTAGA</td><td rowspan="2">1037</td><td rowspan="2">95 °C 30s, 60 °C 30s, 72 °C 45s</td><td rowspan="2">This study</td></tr><tr><td>OD4-R1</td><td>GCTTCTTGCGATTGTTGTTGG</td></tr></tbody></table>

Immunoblots were also reacted with a monoclonal antibody made against the N-terminal sequence of omega-5 gliadin (Fig. 4). As expected, the antibody reacted strongly with the proteins identified as omega-5 gliadins in Keumkang (spots 1–5 in Fig. 1d) and Olgeuru (spots 1–6 in Fig. 1e). However, the antibody also reacted with the two minor protein spots just below the identified omega-5 gliadins in the parental lines as well as in the mutant line (Fig. 4), suggesting that these spots

also correspond to omega-5 gliadins. Minor reactivity with a more basic protein of about 35 kDa was also observed in all of the lines, but was difficult to align with a specific spot on the Coomassie stained gels.

**Arrangement of omega gliadin genes in cv. Chinese Spring**

To investigate the origin of the minor reactive proteins in DH20, the arrangement of omega gliadin genes in the reference wheat Chinese Spring was examined in detail.

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /></figure>

Of the nineteen omega gliadin genes annotated in Chinese Spring [11], four are located on chromosome 1A, eight on 1B and seven on 1D (Fig. 5). Repetitive motifs of seven genes are characteristic of omega-1,2 gliadins while those of 12 genes are characteristic of omega-5 gliadins. A cluster of eight omega-5 gliadin genes are located within a 162 kb region of chromosome 1B that lies between two ancestor genes common to all three genomes as well as *Brachypodium*, rice and sorghum [20]. Six of these are pseudogenes with various mutations that introduce frameshifts and/or premature stop codons into the coding regions. Interspersed among the pseudogenes are omega-B3 and omega-B6 encoding complete omega-5 gliadins. Interestingly, four other omega-5 gliadin genes are located on a 51 kb region between the orthologous ancestor genes on chromosome 1D (Fig. 5). Three of these are obvious pseudogenes while one, omega-D4, contains a single base deletion 50 bp

upstream from the stop codon. If expressed, omega-D4 would encode a protein of 44.4 kDa with an altered C-terminus.

The amino acid sequences of the omega-5 gliadin proteins encoded by the B genome in Chinese Spring (omega-B3 and omega-B6) were compared to that of omega-D4 (Fig. 6). The only other full-length omega-5 gliadin sequence in NCBI, BAE20328, was also included in the analysis. Omega-B3 differs from BAE20328 by two amino acid substitutions and a 28 amino acid deletion while omega-B6 has six amino acid substitutions and a five amino acid insertion. In comparison, the sequence of omega-D4 is quite different from BAE20328. The N-terminal sequences of the B-encoded proteins and BAE20328 are SRL while that of omega-D4 is TRQ. Additionally, the last 16 amino acids at the C-termini of omega-B3, omega-B6 and BAE20328 are identical while nine of the last 16 amino acids are altered in omega-D4

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /><figcaption>**Fig. 6** Sequence comparison of omega-5 gliadins encoded by the B and D genomes in Chinese Spring with the full-length omega-5 gliadin BAE20328. N- and C-terminal sequences are shown in black and blue boxes, respectively. Within each protein sequence, the dominant WDEIA epitopes QQFPQQQ and QQIPQQQ are shown in red and blue, respectively. Lines above the protein sequence highlight peptides unique to omega-D4 that were identified by MS/MS in protein spots from DH20 reacting with sera from WDEIA patients. Red lines highlight peptides identified in the more basic protein while blue lines highlight peptides identified in the more acidic protein. MS/MS data are shown in Additional file 1</figcaption></figure>

and the last six amino acids are missing. Despite the smaller size of omega-D4, the protein contains 6 to 8 more copies of the dominant WDEIA epitopes than BAE20328 or the Chinese Spring B-encoded omega-5 gliadins (Table 3). Omega-D4 contains 27 copies of QQFPQQQ and 6 copies of QQIPQQQ. Notably, more than 90% of these epitopes overlap with other epitopes. In fact, there are two regions of the protein where two epitopes overlap, three regions where four epitopes overlap, one region where five epitopes overlap, and one region where nine epitopes overlap. In comparison, only 52–56% of epitopes in B-encoded omega-5 gliadins overlap and the regions of overlap are shorter. BAE20328 and omega-B3 each have only two regions where two epitopes overlap, two regions

where 3 epitopes overlap and one region where four epitopes overlap while omega-B6 has four regions where two epitopes overlap, one region where three epitopes overlap and one region where four epitopes overlap (Fig. 6, Table 3).

## Identification of minor immunogenic proteins in DH20

The two minor unidentified proteins that reacted with WDEIA sera in DH20 were excised from 2-D gels, digested with chymotrypsin, thermolysin or trypsin and analyzed by MS/MS. For analysis of spectral data, the sequences of gluten protein genes from Chinese Spring were added to the database used for interrogation. Peptides unique to omega-D4 were identified in spots corresponding to both of the reactive proteins (Fig. 6,

<table><caption>**Table 3** Characteristics of omega-5 gliadins</caption><thead><tr><th rowspan="2">Omega-5 Gliadin</th><th rowspan="2">Predicted MW (kDa)</th><th rowspan="2">Predicted pI</th><th colspan="2"># WDEIA epitopes</th><th rowspan="2"># Overlapping epitopes</th></tr><tr><th>QQFPQQQ</th><th>QQIPQQQ</th></tr></thead><tbody><tr><td>BAE20328</td><td>51.0</td><td>6.26</td><td>23</td><td>4</td><td>14</td></tr><tr><td>Omega-B3 (AWK59777)</td><td>47.7</td><td>6.16</td><td>22</td><td>3</td><td>14</td></tr><tr><td>Omega-B6 (AWK59773)</td><td>51.5</td><td>6.01</td><td>23</td><td>4</td><td>15</td></tr><tr><td>Omega-D4</td><td>44.4</td><td>6.34</td><td>27</td><td>6</td><td>30</td></tr></tbody></table>
<figure><img src="image_9.png" /><figcaption>**Fig. 7** PCR analysis with primers specific for the omega-D4 gene from Chinese Spring. Genomic DNA from Keumkang (1), Olgeuru (2), DH20 (3), Chinese Spring (4) and nullisomic tetrasomic lines of Chinese Spring N1AT1B (5), N1AT1D (6), N1BT1A (7), N1BT1D (8), N1DT1A (9) and N1DT1B (10) was amplified. The sizes of molecular weight markers in kb are shown in lane M</figcaption></figure>

Additional file 1), suggesting that these proteins are encoded by the D genome.

PCR analysis of genomic DNA with primers specific for the Chinese Spring omega-D4 gene (Table 2) demonstrated the presence of omega-D4 in both of the parental lines as well as in the mutant DH20 (Fig. 7). Amplification of aneuploid lines of Chinese Spring further confirmed the location of omega-D4 on the 1D chromosome. Direct sequencing of the PCR product from DH20 using OD4-F1, OD4-R1 and three internal primers (OD4-F3, CAATTCCCCGAACAACAATTC; OD4-F4, AACAGCAATTCCCTCAACAGC; OD4-F5, CAACAACACCAGTTACCG) yielded a single sequence that differed from omega-D4 by an A to C transition 968 bp from the initiation codon. Since this sequence encodes a protein with a pI of 6.84, it is likely that it corresponds to the more basic reactive protein in the immunoblots. It is not known whether DH20 expresses another omega-5 gliadin gene corresponding to the more acid protein.

## Discussion

Analysis of flour proteins by 2-DE combined with MS/MS revealed that s-type and m-type LMW-GS encoded by the B genome as well as omega-5 gliadins and several gamma gliadins are absent in the mutant DH20. It is well known that the LMW-GS genes are tightly linked to gamma and omega gliadin genes on the short arms of the group 1 homoeologous chromosomes in hexaploid wheat. A relatively new group of gluten proteins, the delta gliadins, are also encoded in this region [21]. Recent efforts in genome sequencing revealed the structural organization of these groups of genes in the reference wheat Chinese Spring and identified the sequences of 4 LMW-GS, 4 gamma gliadins, 2 delta gliadins and 3 omega gliadins located on a 5.3 Mb region of the 1A

chromosome, 4 LMW-GS, 6 gamma gliadins, 1 delta gliadin and 8 omega gliadins located on a 6.5 Mb region of the 1B chromosome, and 7 LMW-GS, 4 gamma gliadins, 2 delta gliadins and 7 omega gliadins located on a 5.6 Mb region of the 1D chromosome [11]. Amplification of genomic DNA from the parental lines and the mutant line DH20 with primers for omega-5 gliadin and allele-specific LMW-GS genes confirmed that these genes from the 1B chromosome are missing in DH20. The use of repeat junction primers that amplify sequences near the ends of the 6.5 Mb region of the 1B chromosome further indicated that at least a 5.8 Mb region of the 1B chromosome was deleted in the mutant line. This region of the chromosome also contains multiple resistance genes [11]. Despite the sizeable deletion and the absence of multiple prolamin genes, a previous study found that the flour from the mutant line has reasonable breadmaking quality [18].

The loss of the omega-5 gliadin genes in DH20 is of particular interest given the importance of the omega-5 gliadins in food allergies. As expected, DH20 showed reduced levels of reactivity with IgE from sera of WDEIA patients. However, it was curious that several highly immunoreactive proteins remained in this line. These proteins correspond to minor spots in 2-D gels that were found to react with a monoclonal antibody specific for omega-5 gliadins, suggesting that they are omega-5 gliadins encoded outside the deleted region of chromosome 1B.

It has been generally accepted that omega-5 gliadins are encoded by genes at the *Gli-B1* locus on the short arm of the 1B chromosome. Indeed, primers specific for the known full-length omega-5 gliadin genes amplified DNA from Chinese Spring but not from the N1BT1A or N1BT1D nullisomic tetrasomic lines. However, over the years there also have been a number of reports of omega gliadins encoded by minor loci [22, 23]. Recent genome sequence analyses in the reference wheat Chinese Spring has made it possible to survey the entire collection of gliadin and LMW-GS genes in a hexaploid wheat cultivar [11, 24]. This work uncovered four omega-5 gliadin genes on the D genome in addition to the eight on the B genome. All of the D genome genes were classified as pseudogenes because of frameshift mutations in the coding regions and low FPKM values in transcriptome data. However, one gene, omega-D4, contained only a single mutation near the 3' end of the coding region that, if expressed, would yield a truncated protein with an altered C-terminus. It is interesting that there are a number of sequences for omega-5 gliadin genes in NCBI that encode proteins with altered C-termini. Most were from a study in which 66 omega-5 gliadins were amplified from cDNA prepared from developing grain from a variety of diploid, tetraploid and hexaploid wheats [9]. All were missing portions of the central repetitive regions
and thus would encode proteins of only 4 to 38 kDa. Seventeen of the 66 sequences encoded proteins with altered C-termini similar to that of omega-D4, including 15 from the diploid species *Ae. tauschii*, *Ae. speltoides*, *Ae. searsii* and *T. monococcum*. It is also interesting that 23 of 29 sequences amplified from hexaploid wheat and 14 of 17 sequences amplified from tetraploid wheat encoded proteins with the N-terminal sequence SRL while all 20 of the sequences amplified from diploid wheat encoded proteins with the N-terminal sequence TRQ that is also found in omega-D4. Additionally, a gene with only 4 bp mismatch to omega-D4 was identified in *Ae. tauschii*, the D progenitor of hexaploid wheat [20].

Omega-5 gliadins are unusual in that each protein contains many copies of IgE-reactive epitopes. Omega-D4 contains a greater number of the major WDEIA epitopes than the other full-length omega-5 gliadins deposited in the NCBI database. Additionally, there is considerable overlap among these epitopes in many regions of the protein. The fact that these proteins are very minor species in both the parental lines and DH20, yet the proteins are nearly as IgE reactive with patient sera as the more abundant omega-5 gliadins suggests that both the number and the arrangement of the epitopes in the protein are likely important in determining its immunoreactivity.

In previous studies, RNA interference was used to reduce the levels of omega-5 gliadins in wheat flour, resulting in transgenic plants with reduced IgE reactivity to sera from WDEIA patients without adverse effects on flour end-use functional properties [15-17]. This homology-based approach used a 153 bp target region corresponding to the 3' end of the coding region and effectively eliminated all omega-5 gliadins from the flour with relatively few other changes on the proteome. However, the resulting plants were transgenic and not likely to be incorporated into breeding programs because of consumer acceptance issues. The present work demonstrates that breeding approaches also can be used to reduce the levels of omega-5 gliadins and hence the immunogenic potential of the flour. However, such approaches are complicated by the genetic linkage of omega gliadin, gamma gliadin and LMW-GS genes and the finding that minor but highly immunoreactive omega-5 gliadins may be encoded on chromosomes other than 1B.

Currently, the exact sequences of omega-5 gliadins in the parental lines Keumkang and Olgeuru and the mutant line DH20 are not known. It is also not known whether other hexaploid wheat cultivars contain active omega-5 gliadin genes on either their A or D chromosomes. However, it is clear that future efforts to reduce the levels of omega-5 gliadins and other immunogenic proteins in wheat flour, either by breeding approaches that rely on induced mutations or biotechnology approaches that utilize homology-based methods such as

gene silencing or genome editing, would benefit from a detailed understanding of the complete complement of gluten protein genes in individual cultivars. Clearly, it is important to know the numbers of genes and pseudogenes in these complex families as well as their precise sequences and genomic locations. To this end, the reference genome sequence from Chinese Spring now makes it possible to design gene capture methods to select genomic regions containing gluten protein genes from individual cultivars for long-read sequencing. Thus, it should be possible to obtain the complete sequences of highly repetitive omega-5 gliadin genes that have not been easily cloned in the past and survey the entire collection of these genes in different cultivars of hexaploid wheat. Such studies will also provide insight into the evolution of this interesting family of gluten protein genes.

## Conclusions

Breeding or biotechnology approaches to reduce the immunogenic potential of wheat flour require that multiple gluten proteins be eliminated without impacting the functional properties of the flour. However, the lack of complete information on the sequences of the gluten proteins and the locations of the genes that encode them complicates these efforts. In this study, a mutant wheat line missing the major omega-5 gliadins encoded on chromosome 1B was shown to have reduced IgE reactivity with sera from WDEIA patients. However, several minor but highly immunogenic proteins remained in the flour. Further analysis revealed that these proteins were omega-5 gliadins encoded by genes located on chromosome 1D that had not been described previously and that the proteins contained a greater number of WDEIA epitopes and more overlapping epitopes than those encoded on chromosome 1B. The work illustrates the importance of detailed knowledge about the complete sets of gluten protein genes and proteins in individual wheat cultivars for future efforts aimed at reducing the immunogenic potential of wheat flour.

## Methods

### Plant material

*Triticum aestivum* L. cv. Keumkang, a leading hard wheat cultivar in Korea that is used for both breadmaking and noodle production, and cv. Olgeuru, a soft wheat used for noodle production, were the parental lines used to produce a doubled haploid population from which the mutant line DH20 was selected [18]. Keumkang was shown previously to contain the *Glu-A3c*, *Glu-B3h* and *Glu-D3a* alleles while Olgeuru contained the *Glu-A3d*, *Glu-B3d* and the *Glu-D3a* alleles [25]. The reference cultivar Chinese Spring, containing *Glu-A3a*, *Glu-B3a* and *Glu-D3a* alleles, and nullisomic tetrasomic lines of Chinese Spring, kindly provided by National
BioResource Project-Wheat (Kyoto, Japan), served as controls for PCR analyses. Plant material was grown in the field as described previously [18].

## Protein extraction, 2-DE analysis and mass spectrometry

Total flour proteins from Keumkang, Olgeuru and DH20 were extracted with SDS buffer (2% SDS, 10% glycerol, 50 mM DTT, 40 mM Tris-Cl, pH 6.8), quantified and analyzed by 2-DE as described previously [19]. Individual protein spots were excised from Coomassie-stained gels and transferred to 96-well plates for digestion with thermolysin, chymotrypsin or trypsin using a DigestPro (Intavis, Koeln, Germany). Plates containing the peptides were placed in the autosampler of a nanoflow HPLC interfaced to an Orbitrap Elite mass spectrometer (Thermo Scientific, San Jose, CA). Details of data collection are described in [26]. MS/MS spectra were searched against a database containing *Triticeae* protein sequences downloaded from NCBI on November 1, 2016, cultivar-specific sequences from *Triticum aestivum* cv. Butte 86 reported in [19] and a database of common contaminants downloaded on January 30, 2015. Two search engines, Mascot and X!Tandem, were used for the analyses and the results were compiled in Scaffold Version 4.7.5 (Proteome Software, Inc., Portland, OR). Data extracted from Scaffold are summarized in Additional file 1.

An updated database was used for the identification of minor IgE reactive proteins in DH20. This database included *Triticeae* protein sequences downloaded from NCBI on February 7, 2018, cultivar-specific gluten protein sequences from cv. Butte 86 [19], sequences of gliadins and LMW-GS from cv. Chinese Spring [11, 24] and a database of common contaminants downloaded on January 30, 2015.

## PCR analysis

Genomic DNA was isolated from young leaf tissue of Keumkang and Olgeuru parental lines, DH20 mutant line, Chinese Spring, and the nullisomic tetrasomic lines of Chinese Spring using the DNeasy Plant mini kit (Qiagen, Hilden, Germany) according to protocols provided by the manufacturer. PCR was performed using Extaq polymerase (Takara Bio Inc., Kyoto, Japan). PCR primers and amplification conditions are shown in Table 2. Primers specific for LMW-GS *GluB3-4* gene found in Keumkang, Olgeuru and Chinese Spring were described previously [27]. Forward and reverse primers for omega-5 gliadin were based on the sequence of AB181300 and would be expected to amplify a 1212 bp fragment. The design of repeat junction primers is described in [28]. These primers were based on the assembled sequence of a 6.5 Mb region of chromosome 1B from cv. Chinese Spring reported in [11], NCBI Accession Number MG560141. Primer pair 19S-1.3-2 amplifies a 418 bp region that is 1.74 Mb upstream of the first ancestral gene shown in Fig. 2 in [11]. The forward primer corresponds to the sequence between position

700,904 and 700,923 and the reverse primer corresponds to the sequence between 701,298 and 710,321 in NCBI Accession Number MG560141. Primer pair 143E-1-600 amplifies a 482 bp region that is just before the last two ancestral genes shown in Fig 2 in [11]. The forward primer corresponds to the sequence between 6,512,265 and 6,512,288 and the reverse primer corresponds to the sequence between 6,512,722 and 6,512,746 in NCBI Accession Number MG560141.

PCR primers were also designed to amplify the omega-D4 gene from cv. Chinese Spring reported in [11]. The forward primer corresponds to the sequence from 58 to 78 bp downstream from the initiation codon and the reverse primer corresponds to the sequence from 1074 to 1094 downstream from the initiation codon. The forward and reverse primers have three and seven base mismatches with AB181300, the gene sequence encoding omega-5 gliadin BAE20328.

## Immunoblot analysis

For immunoblot analysis of the parental and mutant lines, 7.5 µg total SDS-extracted protein were separated by 2-DE and transferred onto nitrocellulose membranes as described in [29]. Immunoblot analysis was performed as described in [30] using sera from WDEIA patients that was obtained with the informed consent of patients and approval from the Ethics Committee of Ile de France III and the French Agency for the Safety of Health Products (AFSSAPS) (authorization number 2008-A01565-50) and characterized previously [16]. IgE-binding proteins were visualized using peroxidase-labeled rabbit anti-human IgE (9160-05 SouthernBiotech) diluted 1:500,000 and chemiluminescent substrate (WesternBright Quantum, Advansta K-120420D20) according to the manufacturer's instructions. Alternately, membranes were reacted with a monoclonal antibody against the peptide SRLLSPRGKELG (mono ONT 18A5) that was diluted 1/80 and visualized using peroxidase-labeled goat anti mouse IgG (H, L)-HRP conjugate (170-6516 Biorad) diluted 1:50,000.

## Additional file

**Additional file 1:** Summary of MS/MS data from protein spots shown in Fig. 1. Data was extracted from Scaffold. (XLSX 3032 kb)

### Abbreviations

2-DE: Two-dimensional gel electrophoresis; MS/MS: Tandem mass spectrometry; WDEIA: Wheat dependent exercise-induced anaphylaxis

### Acknowledgements

The authors thank Leslie Harden and Melissa Erickson-Beltran (USDA-ARS) for assistance with mass spectrometry work and Olivier Tranquet (INRA BIA) for providing the monoclonal antibody ONT18A5.

### Funding

This work was supported by grants from the National Institute of Agricultural Science (RDA PJ012458) and the Next-Generation BioGreen 21 Program (RDA
PJ013149 and RDA PJ013159), South Korea and by USDA Agricultural Research Service Research Project 5325-43000-028-00D. Mention of a specific product by the United States Department of Agriculture does not constitute an endorsement and does not imply a recommendation over other suitable products. USDA is an equal opportunity provider and employer.

**Availability of data and materials**

The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.

**Authors' contributions**

SA and J-YL designed the study and managed the research project. H-CC, AS-B, Y-RJ and FP performed the experiments. SA wrote the manuscript. SA, J-YL, SD-P, YG, NH, S-HL and C-SK participated in data analysis and/or review of the manuscript. All authors read and approved the manuscript.

**Ethics approval and consent to participate**

Not applicable.

**Consent for publication**

Not applicable.

**Competing interests**

The authors declare that they have no competing interests.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Author details**

¹USDA-ARS Western Regional Research Center, 800 Buchanan Street, Albany, CA 94710, USA. ²National Institute of Agricultural Sciences, RDA, Jeonju 54874, South Korea. ³UR1268 Biopolymers, Interactions, Assemblies, Institut National de la Recherche Agronomique, Rue de la Géraudière, F-44316 Nantes, France. ⁴National Institute of Crop Science, RDA, Jeonju 55365, South Korea.

Received: 24 July 2018 Accepted: 26 October 2018

Published online: 21 November 2018

**References**

1. Battais F, Courcoux P, Popineau Y, Kanny G, Moneret-Vautrin DA, Denery-Papini S. Food allergy to wheat: differences in immunoglobulin E-binding proteins as a function of age or symptoms. J Cereal Sci. 2005;42:109-17.

2. Morita E, Matsuo H, Mihara S, Morimoto K, Savage AWJ, Tatham AS. Fast ω-gliadin is a major allergen in wheat-dependent exercise-induced anaphylaxis. J Derm Sci. 2003;33:99-104.

3. Tye-Din JA, Stewart JA, Dromey JA, Beissbarth T, van Heel DA, Tatham A, Henderson K, Mannering SI, Gianfrani C, Jewell DP, Hill AVS, McCluskey J, Rossjohn J, Anderson RP. Comprehensive, quantitative mapping of T cell epitopes in gluten in celiac disease. Sci Transl Med. 2010;4(41):41a51.

4. Altenbach SB, Kothari KM. Omega gliadin genes expressed in *Triticum aestivum* cv. Butte 86: effects of post-anthesis fertilizer on transcript accumulation during grain development. J Cereal Sci. 2007;46:169-77.

5. Altenbach SB, Tanaka CK, Hurkman WJ, Whitehand LC, Vensel WH, Dupont FM. Differential effects of a post-anthesis fertilizer regimen on the wheat flour proteome determined by quantitative 2-DE. Proteome Sci. 2011;9:46.

6. Hurkman WJ, Tanaka CK, Vensel WH, Thilmony R, Altenbach SB. Comparative proteomic analysis of the effect of temperature and fertilizer on gliadin and glutenin accumulation in the developing endosperm and flour from *Triticum aestivum* L. cv. Butte 86. Proteome Sci. 2013;11:8.

7. Wan Y, Gritsch CS, Hawkesford MJ, Shewry PR. Effects of nitrogen nutrition on the synthesis and deposition of the ω-gliadins of wheat. Ann Botany. 2014;113:607-15.

8. Dupont FM, Vensel WH, Chan R, Kasarda DD. Characterization of the 1B-tpe ω-gliadins from *Triticum aestivum* cultivar Butte 86. Cereal Chem. 2000;77:607-14.

9. Du X, Tang H, Li M, Ma X, Yin H, Wang H, Zhang X, Qiao X, Li A, Kong L. Molecular characterization of the IgE-binding epitopes in the fast ω-gliadins of *Triticaceae* in relation to wheat-dependent, exercise-induced anaphylaxis. Gene. 2016;591:27-33.

10. Anderson OD, Gu YQ, Kong X, Lazo GR, Wu J. The wheat ω-gliadin genes: structure and EST analysis. Funct Integr Genomics. 2009;9:397-410.

11. Huo N, Zhang S, Zhu T, Dong L, Mohr T, Hu T, Liu Z, Dvorak J, Luo M-C, Wang D, Lee J-Y, Altenbach S, Gu YQ. Gene duplication and evolution dynamics in the homeologous regions harboring multiple prolamin and resistance gene families in hexaploid wheat. Front Plant Sci. 2018;9:673.

12. Matsuo H, Morita E, Tatham AS, Morimoto K, Horikawa T, Osuna H, Ikezawa Z, Kaneko S, Kohno K, Dekio S. Identification of the IgE-binding epitope in ω-5 gliadin, a major allergen in wheat-dependent exercise-induced anaphylaxis. J Biol Chem. 2004;279:12135-40.

13. Denery-Papini S, Lauriere M, Branlard G, Morisset M, Pequet C, Choudat D, Merlino M, Pineau F, Popineau Y, Boulenc E, Bouchez-Mahiout I, Bodinier M, Moneret-Vautrin D-A. Influence of the allelic variants encoded at the *Gli-B1* locus, responsible for a major allergen of wheat, on IgE reactivity for patients suffering from food allergy to wheat. J Agric Food Chem. 2007;55:799-805.

14. Waga J, Skoczowski A. Development and characteristics of ω-gliadin-free wheat genotypes. Euphytica. 2014;195:105-16.

15. Altenbach SB, Allen PV. Transformation of the US bread wheat 'Butte 86' and silencing of omega-5 gliadin genes. GM Crops. 2011;2(1):67-74.

16. Altenbach SB, Tanaka CK, Pineau F, Lupi R, Drouet M, Beaudoin E, Morisset M, Denery-Papini S. Assessment of the allergenic potential of transgenic wheat (*Triticum aestivum*) with reduced levels of ω5-gliadins, the major sensitizing allergen in wheat-dependent exercise-induced anaphylaxis. J Agric Food Chem. 2015;63:9323-32.

17. Altenbach SB, Tanaka CK, Seabourn BW. Silencing of omega-5 gliadins in transgenic wheat eliminates a major source of environmental variability and improves dough mixing properties of flour. BMC Plant Biol. 2014;14:1.

18. Lee J-Y, Kang C-S, Beom H-R, Jang Y-R, Altenbach SB, Lim S-H, Kim Y-M, Park C-S. Characterization of a wheat mutant missing low-molecular-weight glutenin subunits encoded by the B-genome. J Cereal Sci. 2017;73:158-64.

19. DuPont FM, Vensel WH, Tanaka CK, Hurkman WJ, Altenbach SB. Deciphering the complexities of the wheat flour proteome using quantitative two-dimensional electrophoresis, three proteases and tandem mass spectrometry. Proteome Sci. 2011;9:10.

20. Dong L, Huo N, Wang Y, Deal K, Wang D, Hu T, Dvorak J, Anderson OD, Luo M-C, Gu YQ. Rapid evolutionary dynamics in a 2.8-Mb chromosomal region containing multiple prolamin and resistance gene families in *Aegilops tauschii*. Plant J. 2016;87:495-506.

21. Anderson OD, Dong L, Huo N, Gu YQ. A new class of wheat gliadin genes and proteins. PLoS One. 2012;7(12):e52139.

22. Pogna NE, Metakovsky EV, Redaelli R, Raineri F, Dachkevitch T. Recombination mapping of *Gli 5*, a new gliadin-coding locus on chromosomes 1A and 1B in common wheat. Theor Appl Genet. 1993;87:113-21.

23. Ruiz M, Carrillo JM. Linkage relationships between prolamin genes on chromosomes 1A and 1B of durum wheat. Theor Appl Genet. 1993;87:353-60.

24. Huo N, Zhu T, Altenbach S, Dong L, Wang Y, Mohr T, Liu Z, Dvorak J, Luo M-C, Gu YQ. Dynamic evolution of α-gliadin prolamin gene family in homeologous genomes of hexaploid wheat. Sci Rep. 2018;8:5181.

25. Lee J-Y, Jang Y-R, Beom H-R, Altenbach SB, Lim S-H, Lee CK. Allelic analysis of low molecular weight glutenin subunits using 2-DGE in Korean wheat cultivars. Breed Sci. 2017;67:398-407.

26. Vensel WH, Tanaka CK, Altenbach SB. Protein composition of wheat gluten polymer fractions determined by quantitative two-dimensional gel electrophoresis and tandem mass spectrometry. Proteome Sci. 2014;12:8.

27. Wang LH, Zhao XL, He ZH, Ma W, Appels R, Pena RJ, Xia XC. Characterization of low-molecular-weight glutenin subunit *Glu-B3* genes and development of STS markers in common wheat (*Triticum aestivum* L.). Theor Appl Genet. 2009;118:525-39.

28. You FM, Wanjugi H, Huo N, Lazo GR, Luo M-C, Anderson OD, Dvorak J, Gu YQ. RJPrimers: unique transposable element insertion junction discovery and PCR primer design for marker Development Nucleic Acids Res 38, Web Server issue W313-W320. https://doi.org/10.1093/nar/gkq425. 2010.

29. Huebener S, Tanaka CK, Uhde M, Zone JJ, Vensel WH, Kasarda DD, Beams L, Briani C, Green PHR, Altenbach SB, Alaedini A. Specific nongluten proteins of wheat are novel target antigens in celiac disease humoral response. J Proteome Res. 2015;14:503-11.

30. Lupi R, Denery-Papini S, Rogniaux H, Lafandra D, Rizzi C, De Carli M, Moneret-Vautrin DA, Masci S, Larré C. How much does transgenesis affect wheat allergenicity? Assessment in two GM lines over-expressing endogenous genes. J Proteomics. 2013;80:281-91.