Skip to main content

Celiac disease T-cell epitopes from gamma-gliadins: immunoreactivity depends on the genome of origin, transcript frequency, and flanking protein variation



Celiac disease (CD) is caused by an uncontrolled immune response to gluten, a heterogeneous mixture of wheat storage proteins. The CD-toxicity of these proteins and their derived peptides is depending on the presence of specific T-cell epitopes (9-mer peptides; CD epitopes) that mediate the stimulation of HLA-DQ2/8 restricted T-cells. Next to the thoroughly characterized major T-cell epitopes derived from the α-gliadin fraction of gluten, γ-gliadin peptides are also known to stimulate T-cells of celiac disease patients. To pinpoint CD-toxic γ-gliadins in hexaploid bread wheat, we examined the variation of T-cell epitopes involved in CD in γ-gliadin transcripts of developing bread wheat grains.


A detailed analysis of the genetic variation present in γ-gliadin transcripts of bread wheat (T. aestivum, allo-hexaploid, carrying the A, B and D genome), together with genomic γ-gliadin sequences from ancestrally related diploid wheat species, enabled the assignment of sequence variants to one of the three genomic γ-gliadin loci, Gli-A1, Gli-B1 or Gli-D1. Almost half of the γ-gliadin transcripts of bread wheat (49%) was assigned to locus Gli-D1. Transcripts from each locus differed in CD epitope content and composition. The Gli-D1 transcripts contained the highest frequency of canonical CD epitope cores (on average 10.1 per transcript) followed by the Gli-A1 transcripts (8.6) and the Gli-B1 transcripts (5.4). The natural variants of the major CD epitope from γ-gliadins, DQ2-γ-I, showed variation in their capacity to induce in vitro proliferation of a DQ2-γ-I specific and HLA-DQ2 restricted T-cell clone.


Evaluating the CD epitopes derived from γ-gliadins in their natural context of flanking protein variation, genome specificity and transcript frequency is a significant step towards accurate quantification of the CD toxicity of bread wheat. This approach can be used to predict relative levels of CD toxicity of individual wheat cultivars directly from their transcripts (cDNAs).


Gluten proteins present in wheat can cause celiac disease (CD) [1]. Wheat gluten consists of a complex mixture of α/β-, γ- and ω-gliadins, and high and low molecular weight (HMW; LMW) glutenins, which are all encoded by medium to large multigene families. In contrast to other dietary proteins, gluten proteins are minimally digested by gastrointestinal proteases, resulting in relatively large peptides that accumulate in the small intestine [2]. In CD patients, during passage through the mucosal tract these gluten fragments can bind to the HLA-DQ2/8 receptors present on antigen presenting cells (APCs) and trigger T-cell responses. This T-cell response to HLA-DQ2/8 restricted gluten peptides occurs only in celiac patients and not in healthy controls. The HLA-DQ2/8 receptor favors peptides that contain one or more amino acids with a negative charge [310]. These are not present in gluten peptides but can be introduced due to the activity of human tissue transglutaminase (TG2) [1114] an enzyme that converts glutamine (Q) residues into negatively charged glutamic acid (E) residues. Gliadin peptides, with their high proline and glutamine content, are perfect substrates for the transglutaminase reaction of TG2, which is critical for the creation of active T-cell epitopes involved in CD (CD epitopes). Since target sites for TG2 may be determined by sequences flanking the core motif, the amino acid composition of the surrounding of the epitope cores can influence the immunogenic capacity of a peptide regarding CD [8, 13, 14].

Several HLA-DQ2 and -DQ8 restricted T-cell epitopes have been identified in wheat gluten and in homologous proteins from barley (hordeins) and rye (secalins) ([7, 8, 1526]. A 33-mer peptide, resistant to the action of proteases, derived from α-gliadins is regarded as one of the most CD-immunodominant gluten peptides [2, 27]. However, there is clear evidence that γ-gliadins also contain T-cell stimulatory peptides [7, 8, 19, 22, 24, 2729]. In several specific cohorts a high frequency of CD patients was observed that mainly reacted to γ-gliadin peptides. Vader et al. [22] and Camarca et al. [24] found that half of the CD patients (respectively 10 of 20 children and 7 of 14 adults from Southern Europe) did not respond to peptides derived from α-gliadins, but instead reacted to peptides derived from γ-gliadins and LMW glutenins. A major CD epitope derived from γ-gliadins is the DQ2-γ-I epitope (PQQSFPQQQ[17, 19]), recognized by one third of the patients [22, 24]. Among pasta and bread wheat varieties differences have been observed in the capacity to stimulate in vitro DQ2-γ-I specific T-cell clones [30, 31]. These facts justify an elaborated investigation of the γ-gliadins with special regard to the natural variation in CD-epitope content.

The γ-gliadins account for approximately 12% of the flour proteins of the hexaploid wheat cultivar Butte 86 [32]. They are encoded by clusters of tightly linked genes present at the Gli-1 loci (Gli-A1, Gli-B1 or Gli-D1) that are located on the short arms of the respective homoeologous group 1 chromosomes (1AS, 1BS and 1DS) of hexaploid bread wheat (T. aestivum L.) [33] and are tightly linked to the Glu-3 (LMW glutenins) and Gli-3 (ω-gliadins) loci [[34], and references therein]. The number of different γ-gliadin genes in the genome of bread wheat was estimated at 15–40 [35, 36] and these can be clustered into four up to eleven groups [3641]. This diversity is qualitatively reflected at the protein level [42], but not all genes are expressed evenly. Epitope content of the genes varies [39], just as it varied among α-gliadins [31, 4345] according to the genome of origin of the loci.

In the present study CD epitopes from γ-gliadins of bread wheat are considered in their natural context of deduced protein variation, TG2 deamidation sites, transcript frequency and genome specificity. For this purpose, we sequenced genomic γ-gliadin sequences from ancestrally related diploid Triticum/Aegilops species and analyzed these together with expressed γ-gliadin sequences from the NCBI database of hexaploid bread wheat (Triticum aestivum subspec. aestivum, with the A, B and D genome) to accomplish the assignment of γ-gliadin transcripts contigs to one of the three genomic γ-gliadin loci, Gli-A1, Gli-B1 or Gli-D1, as was previously done for the α-gliadins [43, 46]. Furthermore, natural occurring variants of a major γ-gliadin CD epitope, DQ2-γ-I, were analysed using T-cell assays.


Plant material

The taxonomy of Triticum and Aegilops is complex, and different taxonomic systems are being used. In this paper we have followed the classification according to Van Slageren [47]. Bread wheat (T. aestivum subsp. aestivum; 2n = 6x = 42; ABD genome) is an allohexaploid that was formed through hybridization and successive chromosome doubling of three diploid Triticum/Aegilops species. The diploid ancestors of the D genome and the A genome of bread wheat are respectively Ae. tauschii (D) and T. urartu (Au). Ae. speltoides var. speltoides (S) species or closely related species have been suggested as the ancestor of the B genome, but the exact diploid progenitor remains uncertain [4850].

We included the following accessions in this study: Ae. speltoides var. speltoides (CGN 10682, according to the passport data originating from Israel, and CGN10684, originating from Turkey) is related to the B genome, Ae. tauschii (к-1368, Uzbekistan) as a D genome representative, and T. monococcum subspec. monococcum (CGN10542, Ab genome) as a representative of the A genome. The accessions were obtained from the Centre of Genetic Resources (CGN), The Netherlands, except K-1368, which is from N.I. Vavilov All-Russian Scientific Research Institute of Plant Industry (VIR), St. Petersburg, Russia.

At 21 days after anthesis, wheat kernels in the late milk ripening stage (250 mg) were grinded in liquid nitrogen and subsequently DNA was extracted with chloroform/isoamylalcohol as described by Van Herpen et al. [43].

Amplification, cloning and sequencing of genomic γ-gliadin sequences

Primers to amplify γ-gliadin genes from genomic DNA using PCR were designed on the conserved sequences at the 5’ and 3’ end of the coding region of γ-gliadin gene sequences obtained from the EMBL database (forward primer, γ1F: 5’-atgaagaccttactcatcc-3’, and reverse γ11R: 5’-ggacaWagacRttgcacatg-3’). Primers were degenerated to cover as many different γ-gliadin sequences as possible. Amplification was performed in a 25 μl reaction volume, containing 0.2 μM reverse and 0.2 μM forward primer, dNTP mix (0.25 mM each), 1 x Pfu buffer (Stratagene) and a 1:4 (v/v) mixture of Pfu DNA polymerase (Stratagene) (2.5 U/μl) and Goldstar DNA polymerase (Eurogentec) (5 U/μl). Twenty ng of genomic DNA was used as a template for PCR. The PCR programs consisted of 5 min at 94°C followed by 24 cycles 94 °C for 1 min, 53°C for 1 min and 72°C for 2 min with a final extension at 72°C for 10 min. The PCR products were ligated into the pCRII-TOPO vector (Invitrogen) and subsequently transferred into E. coli-XL1-blue cells (Stratagene). Recombinants were identified using blue-white color selection. Positive colonies were picked and grown overnight at 37°C in freeze media (36 mM K2HPO4, 13.2 mM KH2PO4, 1.7 mM trisodium citrate, 0.4 mM MgSO4, 6.8 mM (NH4)2SO4, 4.4% v/v glycerol, 100 μg/ml ampiciline, 10 g/l tryptone, 5 g/l yeast extract and 5 g/l NaCl). The cloned insert was amplified directly from the culture in a PCR reaction using the M13 forward (5’-cgc cag ggt ttt ccc agt cac gac-3’) and reverse primer (5’-agc gga taa caa ttt cac aca gga -3’) in 20 μl reaction volume containing 2 μl of culture. The reaction mixture consisted of the same components as well as concentrations, and utilized the same PCR program as described before (annealing at 55°C). The amplified product was used in sequencing reactions using γ1F or γ6R(new) primers. Additional primers were designed on two internal, conserved regions of the γ-gliadin gene to generate additional sequences: one internal forward primer, γFi2: 5’-ccc(ac)tgcaagaat(at)t(ct)c-3’, and one internal reverse primer, γRi2: 5’-g(ag)a(at)attcttgca(gt)ggg-3’. The sequence data were manually checked using SeqMan (DNASTAR). In total, 69 unique γ-gliadin sequences have been deposited in the EMBL database (Table 1).

Table 1 Genomic γ-gliadin sequences

Neighbor-joining analysis of genomic γ-gliadin clones

The 69 genomic γ-gliadin sequences derived from the diploid wheat species T. monococcum subspec. monococcum (Ab), Ae. speltoides var. speltoides (S) and Ae. tauschii (D), were included in alignment and distance analysis. All analyses were conducted in MEGA4 [51]. Complete sequences were trimmed to include a region spanning sequence domain I to V of the γ-gliadin protein. Due to the typical sequence repeats in the repetitive domains II and IV, the γ-gliadin sequences differed in length, ranging from 722 to 775 base pairs for Ae. tauschii, from 637 to 820 base pairs for T. monococcum subspec. monococcum, and from 725 to 968 base pairs for Ae. speltoides var. speltoides. Next, the trimmed sequences were analyzed using the Neighbor-Joining method [52]. The Neighbor-Joining method returned an optimal tree with the sum of branch length = 4.85106104. The sequence distances were computed using the Maximum Composite Likelihood method [53] in units of the number of base substitutions per site. All positions containing alignment gaps and missing data were eliminated only in pair wise sequence comparisons (Pairwise deletion option). There were a total of 968 positions in the final dataset. Codon positions included were 1st + 2nd + 3rd + non-coding. The number of base substitutions per site within species and between species was using respectively the Within-group function and Between-group function in MEGA4.

Analysis and genomic assignment of γ-gliadin transcripts from bread wheat

A total of 939 γ-gliadin transcripts originating from T. aestivum subspec. aestivum were extracted from the NCBI database (on 13 November 2007). The γ-gliadin transcripts were assembled in 201 contigs of sequences with 98% to 99% sequence identity (SeqMan, Lasergene 8, DNAStar). The consensus nucleotide sequences of those contigs that were composed of at least four transcripts were further analyzed, resulting in 717 transcripts arranged in 47 contigs (ranging from 306 to 1038 nucleotides in length), coding for 26 different γ-gliadin protein isoforms, 10 of which are full length γ-gliadins (Additional file 1).

The γ-gliadins are encoded by clusters of tightly linked genes present at the Gli-1 loci (Gli-A1, Gli-B1 or Gli-D1) located on the short arms of the homoeologous group 1 chromosomes of bread wheat [33]. To determine whether actively expressed γ-gliadin sequences from bread wheat can be assigned to a genomic locus, based on sequence relationship to one of the diploid wheat species representing the ancestors of the bread wheat genomes, as was possible for α-gliadins [43], the sequences of the 26 expressed γ-gliadin isoforms detected in bread wheat were included in a Neighbor-Joining analysis together with different genomic γ-gliadins cloned from the diploid wheat species used in these studies as representatives for the ancestors of cultivated wheat. As the 26 different γ-gliadin isoforms were deduced from contigs of transcripts which are mostly partial sequences of expressed genes, not all of them covered the complete γ-gliadin gene. Therefore, a region of 333 to 384 base pairs in length, covered by all expressed isoforms, starting with domain III (NPC motif) and ending with domain V (MCN motif) was created and compared for sequence distances in a Neighbor-Joining analysis together with the same nucleotide region of genomic sequences derived from diploid Triticum and Aegilops genomes. The nucleotide sequence distances were computed using the Maximum Composite Likelihood method [53] and are in the units of the number of base substitutions per site. All positions containing gaps and missing data were eliminated from the dataset (Complete deletion option). There were in total 193 positions in the final dataset. The tree optimal (sum of branch length = 0.739407630) is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree.

CD epitope screening

The deduced amino acid sequences of 26 γ-gliadin contigs, encompassing the CD epitope containing regions (domain II and V), representing a total of 717 transcripts, were explored for epitopes and surrounding sequence regions using a text explorer program (PatternResearch, in-house developed tool for extracting sequences with specific sequence motifs from large datasets; I Padioleau, EMJ Salentijn and MJM Smulders, unpublished) after which the output file was analysed in Excel for further calculations and sorting. In addition, the regions of domain II and domain IV that harbor CD epitopes were aligned (amino acid sequences, MEGA4) to detect sequence variation and to reconstruct the genomic organization of γ-gliadins expressed from the respective loci, Gli-A1, Gli-B1 and Gli-D1. Sites for tryptic and peptic enzymes were identified with PeptideCutter (at

T-cell proliferation assay

A DQ2-γ-I epitope specific T clone was generated from small intestinal biopsies of a celiac disease patient (DQ2.5) as described before [15, 16, 22]. The patient signed an informed consent form which was approved by the hospital ethics committee. Proliferation assays were performed in triplicate in 150 μl Iscove’s Modified Dulbecco’s Medium (Bio Whittaker, Verviers, Belgium) with 10% pooled normal human serum in 96 well flat-bottom plates using 2x104 gluten specific T-cells stimulated with 105 irradiated HLA-DQ2-matched allogeneic peripheral blood mononuclear cells (3000 rad) in the presence or absence of antigen (1–10 μg/ml, TG2 treated 17 mer pepides). After 2 days 3 H-thymidine (0.5 μCi/well) was added to the cultures, and 18–20 hours thereafter the cells were harvested. 3 H-thymidine incorporation in the T-cell DNA was determined with a liquid scintillation counter (1205 Betaplate Liquid Scintillation Counter, LKB Instruments, Gaithersburg, MD).


γ-gliadins from the Ab, S and D genome

To analyze the genetic diversity of the γ-gliadin gene family a total of 69 different γ-gliadin sequences were cloned from diploid Triticum and Aegilops species containing the Ab, S or D genome (Table 1). These diploid genomes were studied because they are related to the ancestors of respectively the A, B and D genome of hexaploid bread wheat (T. aestivum). γ-Gliadins contain an N-terminal signal peptide to direct them into the lumen of the endoplasmic reticulum. The mature protein consists of a unique N-terminal domain (Figure 1, domain I) followed by an alternation of two variable domains (Figure 1, domain II and IV) and two more conserved sequence domains (Figure 1, domain III and V). In the conserved domains eight conserved cysteine residues are present that are involved in the formation of intra-molecular disulphide bounds, six in domain III and two in domain V [54]. The six conserved cysteine residues of domain III were present in all sequences whereas no conclusions could be drawn on the cysteine’s in domain V, because this domain was only partly sequenced. Most γ-gliadins from Ae. speltoides harboured a ninth conserved cysteine residue on position 14 of domain II [40]. Overall, 6.9% of the clones contained an internal stop codon. The highest percentage of internal stop codons occurred in sequences from T. monococcum, in which a quarter (6 out of 16) of all genomic γ-gliadin sequences had a stop codon in domain II of the γ-gliadin sequence due to C to T transitions, generating an in frame truncation after the sequence motif QPQ(Q/L)QFPQPQQPstop.

Figure 1

Schematic structure of a γ-gliadin protein. Structure of a γ-gliadin protein according to Anderson et al.[36] with the location of the CD epitopes. The protein consists of a short N-terminal signal peptide (S) followed by a unique N-terminal domain (I), two repetitive domains, R1 and R2 (II and IV) and two non-repetitive domains, NR1 and NR2 (III and V). In the first non-repetitive domain and in the second non-repetitive domain respectively six and two conserved cystein residues are present (indicated with vertical lines) that form four interchain disulfide bonds (indicated with arrows). The minimal T-cell epitopes are shown and their approximate position is indicated.

The genetic variation among the genomic clones was studied in more detail using the Neighbor-Joining method, which resulted in four significantly different groups of genomic γ-gliadin sequences. Two groups contained only sequences derived from Ae. speltoides, another group consisted mainly of sequences derived from Ae. tauschii together with a few sequences (JQ269782-JQ269785) from Ae. speltoides, and the fourth group contained the T. monococcum sequences. The number of base substitutions per site among genomic γ-gliadin sequences depended on the species analyzed. For example, all sequences from Ae. tauschii showed an average number of base substitutions of 0.216 per site (21.6 base alterations per 100 base pairs) whereas the level of genetic diversity of γ-gliadin sequences derived from Ae. speltoides and T. monococcum was higher with respectively 0.423 and 0.789 base substitutions per site.

γ-gliadin transcripts assigned to the A, B and D genome of bread wheat

To analyze the genetic diversity of actively expressed γ-gliadins in bread wheat (T. aestivum), γ-gliadin transcripts (939; cDNAs and expressed sequence tags) were extracted from the NCBI database. Of the γ-gliadin transcripts thus obtained 717 were found to code for 26 different γ-gliadin protein isoforms, each represented by contigs with at least four transcripts at a homology level of at least 98%. The number of transcripts for each isoform ranged from 4 to 220, indicating large differences in relative expression rates and redundancy as a part of the differences can be related to different alleles coding for the same amino acid sequence (protein isoform).

Subsequently, the 26 γ-gliadin isoforms of modern bread wheat were assigned to a genomic locus (Gli-A1, Gli-B1 or Gli-D1). To realize this, the consensus nucleotide sequences were studied for sequence distances by Neighbor-Joining analysis (domain III, IV and V; 333 to 384 base pairs) together with the ancestrally related genomic γ-gliadin sequences from diploid Triticum and Aegilops species obtained in these studies, which are carrying the Ab, S or D genome (related to respectively the A, B and D genome from T. aestivum). In the resulting Neighbor-Joining tree 12 topology groups of γ-gliadins could be distinguished (Figure 2). The very center of the tree has low bootstrap values, due to the fact that several branches are joined here at almost the same point, but in varying topology. The different γ-gliadin transcript contigs of bread wheat showed up in seven topology groups (Figure 2, group 3, 5, 6, 7, 8, 9 and 10, red dots). Taking into account the number of transcripts represented by a consensus sequence, almost half (49%, 354 ESTs out of 717 ESTs) of the γ-gliadin transcriptome of bread wheat clustered with D-genomic sequences (Figure 2, group 9 and 10) of Ae. tauschii and they were therefore assigned to the γ-gliadin locus Gli-D1. The other transcripts were assigned evenly to the other homoeologous loci: 25% to locus Gli-A1 on basis of their clustering with Ab-genome sequences derived from T. monococcum (Figure 2, group 6 and 7), and 25% to locus Gli-B1 either by grouping with genomic sequences of Ae. speltoides (S genome, 8%, group 3 and 8) or by BlastN homology to Ae. searsii (Ss genome, 17%, group 5).

Figure 2

Neighbor-Joining relationships of γ-gliadin sequences. Neighbor-Joining relationship of γ-gliadin nucleotide sequences including transcript-contigs of bread wheat (T. aestivum, ABD genome, red circles) and genomic γ-gliadin sequences derived from diploid wheat species; T. monococcum (Ab, squares), Ae. speltoides (S, black circles), Ae. tauschii (D, triangles). In all cases, a region of the γ-gliadin gene, shared by all transcripts sequences and genomic sequences, encompassing domain III, IV and V was analyzed. Numbers on branch points indicate the percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (500 replicates).

CD epitopes of T.aestivum γ-gliadins in their natural context

In γ-gliadins several CD-immunogenic peptides have been identified (Table 2) and an analysis on the basis of sequence identity showed that all groups of γ-gliadin transcripts of bread wheat contained CD epitope cores but with a different frequency (Table 3). With a frequency of 10.1 CD epitope-cores per transcript, the Gli-D1 transcripts were coding for the highest number of γ-gliadin CD epitopes followed by the Gli-A1 transcripts (8.6) and the Gli-B1 transcripts (5.4). As the Gli-D1 transcripts were also the highest expressed ones, most CD epitopes found in the γ-gliadin transcriptome originated from Gli-D1 transcripts (in total 3588 out of 6006 epitope cores or 59%) followed by Gli-A1 (1533 epitope cores or 25%) and Gli-B1 (1005 epitope cores or 16%) (Table 3).

Table 2 CD-epitope cores derived from γ-gliadins
Table 3 The number of CD-epitope cores in γ-gliadin transcripts of T. aestivum

By analyzing the CD epitopes in their natural context, distinct differences in genomic organization and epitope composition were observed among bread wheat γ-gliadin transcripts derived from the three γ-gliadin loci, Gli-A1, Gli-D1 and Gli-B1 (Additional file 2, Figure 3). In many cases, repeat motifs contained overlapping CD epitopes. For instance, transcripts expressed from Gli-D1 contained a repeat motif of 17 amino acids (Q QPQ QPFPQ QPQ QPFPQ; underlined Q residues are target sites for the enzyme tissue transglutaminase) that included three overlapping CD epitopes (DQ2-γ-VIIb, DQ2-Glia-γ2a, DQ2-γ-VI). Because predicted tryptic cleavage sites were absent in these repeats they are most likely resistant to enzymatic degradation after ingestion. Indeed, a sequence of 26 amino acids present in a part of the Gli-D1 transcripts (Figure 3), 9% of the overall γ-gliadin transcriptome and harbouring four CD-epitope cores (FLQPQ QPFPQ QPQ QPYPQ QPQ QPFPQ) has been reported to be resistant to proteolysis [56]. Such CD epitope-rich repeat motifs were also present in transcripts expressed from Gli-A1 and Gli-B1, despite the shorter repetitive domain of Gli-B1 transcripts (Additional file 2). The genomic structure of the C-terminal domain V is less complex and harbours a single CD epitope, DQ2-γ-II (glia-γ30) (IQPQ QPAQL) (Table 4). Only two variants representing 9% of the transcripts displayed an amino acid change within the 9-mer core (L9 to Y9 and P6 to L6).

Figure 3

CD-epitopes of γ-gliadin transcripts of T. aestivum in their natural context. The deduced aminoacid sequences of a T.aestivum γ-gliadin transcripts (N = 61) assigned to locus Gli-D, Neighbor Joining topology group 10, spanning a part of the first repetitive domain of the γ-gliadin sequence. CD T-cell epitopes are depicted: γ-I (PQQSFPQQQ), γ-III (QQPQQPYPQ), γ-IV (SQPQQQFPQ), γ-VI (QQPFPQQPQ), γ-VIIb (QQPQQPFPQ), 26-mer (FLQPQQPFPQQPQQPYPQQPQQPFPQ). Glutamine residues that are a primary targets for the enzyme tissue transglutaminase are underlined (Q) in QxP target sites whereas moderate target sites are depicted in italics (Q) [13, 14]. Cleavage sites: in grey with white letters, chymotrypsin-high specificity; in grey with black letters, chymotrypsin-low specificity; cleavage occurs at the right side (C-terminal direction) of the marked amino acid.

Table 4 The DQ2- γ-II epitope in transcripts of T.aestivum

Immunogenicity of DQ2-γ-I epitope variants

The DQ2-γ-I epitope (PQ QSFPQQQ) found in the variable domain II of γ-gliadins is known as a major CD epitope derived from the γ-gliadin fraction of gluten [17, 22, 24]. In the present study thirteen natural variants of the DQ2-γ-I epitope (including the 9-mer core motif and four flanking amino acid positions on both sides of the epitope core) were found among genomic sequences of T. monococcum subspec. monococcum (Ab genome), Ae. speltoides var. speltoides (S-genome), Ae. tauschii (D-genome), and T. aestivum transcripts (ABD genome) (Table 5). According to the rules for enzymatic TG2 deamidation [13, 14], position Q2 of the DQ2-γ-I epitope core (PQ QSFPQQQ) is a moderate target for TG2 whereas the genetic variation at the C-terminal flanking positions +1 to +3 can influence the deamidation pattern of Q7 to Q9. Based on the natural variation in this flanking region several different deamidation profiles of the DQ2-γ-I epitope core are evident (Table 5 and 6). These variants were deamidated by TG2 and tested for their capacity to stimulate a DQ2-γ-I epitope specific T-cell clone (Table 6). The results show that variants in which Q2 and Q9 are deamidated stimulate the T-cell clone while variants in which Q8 instead of Q9 is deamidated or the F5 is replaced by a S do not display such activity. Thus, natural variants of the DQ2-γ-I epitope exist that influence the T-cell stimulatory capacity.

Table 5 Natural variants of the DQ2-γ-I epitope
Table 6 T-cell stimulating capacity of DQ2-γ-I is influenced by deamidation pattern and flanking amino acids


Wheat gluten is the source of a large repertoire of T-cell epitopes involved in the pathogenesis of CD. It is clear that most patients react to epitopes derived from α-gliadins. However, there is increasing evidence that more T-cell epitopes from other gluten genes are involved in CD [24, 57]. It was shown that next to α-gliadins also T-cell epitopes from γ-gliadins play a significant role in the development of CD [2, 1719, 22, 24]. Tye-Din et al. [57] found that gluten specific polyclonal T-cells in the peripheral blood of CD patients were specific for the same set of gluten peptides after feeding the patients with the same cereal but after the patients ingested respectively wheat, barley and rye different T-cells, recognizing different sets of gluten peptides, were found. Also, it was hypothesized that a T-cell response to gluten can be initiated against a relatively large number of peptides, and when the response evolves it focuses on the most immuno-dominant sequences [22]. So, there are some indications that the composition of different cereals and wheat varieties in the diet [24, 57] as well as the disease status [22] may play a role in the defining the pattern of T cell reactivity towards different gluten peptides in a celiac patient.

Using polyclonal, gliadin reactive T-cell lines, Camarca et al. [24] showed that T-cell reactivity towards the γ-gliadin derived peptides is more heterogeneous compared to the reaction to α-gliadin derived peptides; reaction to the latter was focused towards a few immunodominant peptides whereas a larger repertoire of different γ-gliadin derived peptides was found and patients were recognizing varying sets of peptides. They suggested that this less focussed behaviour of γ-gliadin peptides in celiac disease may reflect the high genetic diversity of γ-gliadins.

To obtain insight in CD epitopes derived from γ-gliadins we performed a detailed study of the presence and natural variation of CD epitopes derived from the γ-gliadin gluten fraction of bread wheat (T. aestivum L.) which is an allo-hexaploid (2n = 6x = 42) carrying three homoeologous genomes A, B and D. This study included two datasets: (1) 69 novel genomic γ-gliadin clones from the diploid wheat species T. monococcum, Ae. speltoides and Ae. tauschii. These diploid wheat species are carrying respectively genome Ab, S and D which are ancestrally related to the three respective homoeologous genomes of bread wheat, A, B and D; and (2) γ-gliadin transcripts from bread wheat (T. aestivum) derived from the public database (Genbank, NCBI). By similarity analysis the transcripts of bread wheat could be assigned to locus Gli-A1, Gli-B1, or Gli-D1 (on the respective homoeologous chromosomes 1A, 1B and 1D). To reduce the complexity of the multigene γ-gliadin family, transcripts of bread wheat were analyzed (instead of genomic sequences). As such, the γ-gliadin transcriptome is a minimal estimate for the γ-gliadin variants in gluten of bread wheat. The γ-gliadin transcripts of bread wheat (N = 717) encode for 26 different γ-gliadin protein isoforms of which ten represented full length γ-gliadins (contigs >300 bp, ≥ four transcripts, 98% to 99% sequence identity). This figure is in keeping with the 15–40 estimated γ-gliadin genes in hexaploid wheat estimated before [35, 36, 41]. Half of the γ-gliadin transcripts from bread wheat grouped with Ae. tauschii sequences (D genome, 6 contigs, 354 transcripts) and were assigned to the γ-gliadin locus on chromosome 1D (Gli-D1), indicating that the majority of the γ-gliadins are expressed from the D genome. The remainder of the transcripts were equally assigned to Gli-A1 (9 contigs, 178 transcripts) and Gli-B1 (11 contigs, 185 transcripts). Similarly, for the α-gliadins it was found that the highest expression occurred from the D genome (Gli-D2) [58].

With regard to CD epitope content, our analysis shows that γ-gliadins from bread wheat contain on average between 4 and 10 potential CD epitopes in the first variable domain (Table 3, range 4.6 to10.9 epitope cores per transcript). The 26-mer γ-gliadin peptide that harbours four distinct CD epitopes [56] is only present in D-genome γ-gliadins. It is present in 9% of the γ-gliadin transcripts of bread wheat. Thus, the γ-gliadins from all three genomes encode a large number of different immunogenic peptides which exceeds the number of identified immunogenic peptides in the α-gliadins. Perhaps this is the basis of the observation that immuno-dominant responses to α-gliadin derived CD epitopes are generally found in patients, while responses to the γ-gliadins appear less consistent as was previously suggested by Camarca et al. [24]. If patients make frequently T-cell responses to a limited number α-gliadin peptides but can respond to a large number of different γ-gliadin derived peptides, T-cells for some specific α-gliadin epitopes may prevail numerically but collectively γ-gliadin derived T-cell epitopes are probably similarly effective in causing CD.

The DQ2-γ-I epitope (PQ QSFPQQQ) is regarded as a major CD epitope as it is frequently recognized by CD patients [22, 24]. We observed a high level of natural genetic variation in the C-terminal flanking sequence of this CD epitope. According to the rules for deamidation by TG2 [13, 14] the variation in the C-terminal flanking sequences determines the deamidation pattern of this epitope core. Among diploid and hexaploid wheat six variants of the 9-mer DQ2-γ-I epitope core were evident with respect to variation in amino acid composition and deamidation pattern. Testing of these DQ2-γ-I variants for their capacity to trigger proliferation of DQ2-γ-I specific T-cells confirmed that TG2 deamidation of glutamine (Q) at position 9 is essential for T-cell stimulation, most likely because it is providing a negative charge which is essential for HLA-DQ2 binding. Variant peptides that displayed a serine (S) at position 5 instead of an phenylalaninine (F) were found to have lost T-cell stimulatory capacity. Moreover, the presence of a positively charged arginine (R) residue at position +1 diminished T-cell proliferation (Additional file 3), most likely because it influenced the deamidation of Q9 as previously observed by Dørum et al. [28]. Similarly, a tryptophan (W) at position +2 inhibited the T-cell response (Additional file 3). Thus, both substitutions within and outside the epitope core can influence the T-cell stimulatory properties of the DQ2-γ-I epitope.

Taken together, the γ-gliadins from all three wheat genomes appear to be a significant source of CD epitopes. Similar to the α-gliadins, the highest number of potential immunogenic γ-gliadin peptides are encoded by the D genome of bread wheat which is thus the most critical genome with regard to CD toxicity. However, gluten derived from tetraploid wheat varieties lacking the D genome is not tolerated by CD patients as well, indicating that the mere elimination of the D genome is not sufficient for the generation of safe wheat. Regarding the A-genome of wheat, there are indications for the existence of T. monococcum spp. with a level or type of gliadin that is unable to induce IFN-γ production and histologic damage as was observed for instance in duodenal biopsy specimens from patients with celiac disease [59]. So, the high level of genetic variation among wheat lines and the presence of genetic variation influencing the immunogenicity of the major CD epitopes may offer possibilities to generate wheat varieties with a reduced CD-immunogenicity, tailored to major CD epitopes (i.e. DQ2-α-I, DQ2-α-II; [46] and/or DQ2-γ-I, this study). Such varieties would help to reduce the presence of immunogenic CD epitopes in wheat flour and, while not safe for consumption by patients, might help to prevent the onset of CD in people that carry genetic risk factors [44, 45]. However, introgression from the diploid to the hexaploid level is a time-consuming process, and consequently it will take many years before the product of such a synthetic hexaploidisation has been bred to sufficient agronomic quality. Alternatively, one could screen a large number of existing varieties for differences in their toxicity for CD patients as the analysis of CD epitope regions in transcript sequences does provide an accurate and quick method for screening [46, 58]. Consequently, we are currently analyzing CD epitope regions in the gluten transcriptome of tetraploid wheat cultivars in a medium-throughput way by employing next-generation sequencing technology.


The results presented here for γ-gliadins form a significant contribution to the approach whereby the CD toxicity of wheat cultivars is quantified directly from protein-coding mRNA sequences.


  1. 1.

    Dicke WK, Weijers NA, van der Kramer JH: Coeliac disease 2: The presence in wheat of a deleterious effect in cases of coeliac disease. Acta Paediatr. 1953, 42: 34-

    Article  CAS  PubMed  Google Scholar 

  2. 2.

    Shan L, Molberg Ø, Parrot I, Hausch F, Filiz F, Gray GM, Sollid LM, Khosla C: structural basis for gluten intolerance in celiac sprue. Science. 2002, 297: 2275-2279. 10.1126/science.1074129.

    Article  CAS  PubMed  Google Scholar 

  3. 3.

    Van de Wal Y, Kooy YMC, Drijfhout JW, Amons R, Koning F: Peptide binding characteristics of the coeliac disease-associated DQ (α1*0501, β1*0201) molecule. Immunogenetics. 1996, 44: 246-253. 10.1007/BF02602553.

    Article  CAS  PubMed  Google Scholar 

  4. 4.

    Vartdal F, Johansen BH, Friede T, Thorpe CJ, Stevanovic S, Eriksen JE, Sletten K, Thorsby E, Rammensee HG, Sollid LM: The peptide binding motif of the disease associated HLA-DQ (α1*0501, β1*0201) molecule. Eur J Immunol. 1996, 26: 2764-2772. 10.1002/eji.1830261132.

    Article  CAS  PubMed  Google Scholar 

  5. 5.

    Godkin A, Friede T, Davenport M, Stevanovic S, Willis A, Jewell D, Hill A, Rammensee HG: Use of eluted peptide sequence data to identify the binding characteristics of peptides to the insulin-dependent diabetes susceptibility allele HLA-DQ8 (DQ 3.2). Int Immunol. 1997, 9: 905-911. 10.1093/intimm/9.6.905.

    Article  CAS  PubMed  Google Scholar 

  6. 6.

    Kim C-Y, Quarsten H, Bergseng E, Khosla C, Sollid LM: Structural basis for HLA-DQ2-mediated presentation of gluten epitopes in celiac disease. Proc Natl Acad Sci U S A. 2004, 101: 4175-4179. 10.1073/pnas.0306885101.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  7. 7.

    Qiao S-W, Bergseng E, Molberg Ø, Jung G, Fleckenstein B, Sollid LM: Refining the rules of gliadin T-cell epitope binding to the disease-associated DQ2 molecule in celiac disease: importance of proline spacing and glutamine deamidation. J Immunol. 2005, 175: 254-261.

    Article  CAS  PubMed  Google Scholar 

  8. 8.

    Tollefsen S, Arentz-Hansen H, Fleckenstein B, Molberg Ø, Ráki , Kwok WW, Jung G, Lundin KEA, Sollid LM: HLA-DQ2 and -DQ8 signatures of gluten T-cell epitopes in celiac disease. J Clin Invest. 2006, 116: 2226-2236. 10.1172/JCI27620.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  9. 9.

    Henderson KN, Tye-Din JA, Reid HH, Chen Z, Borg NA, Beissbarth T, Tatham A, Mannering SI, Purcell AW, Dudek NL, Van Heel DA, McCluskey J, Rossjohn J, Anderson RP: A structural and immunological basis for the role of human leukocyte antigen DQ8 in celiac disease. Immunity. 2007, 27: 23-34. 10.1016/j.immuni.2007.05.015.

    Article  CAS  PubMed  Google Scholar 

  10. 10.

    Stepniak D, Wiesner M, de Ru AH, Moustakas AK, Drijfhout JW, Papadopoulos GK, van Veelen PA, Koning F: Large-scale characterization of natural ligands explains the unique gluten-binding properties of HLA-DQ2. J Immunol. 2008, 180: 3268-3278.

    Article  CAS  PubMed  Google Scholar 

  11. 11.

    Molberg Ø, Mcadam SN, Körner R, Quarsten H, Kristiansen C, Madsen L, Fugger L, Scott H, Norén O, Roepstorff P, Lundin KEA, Sjöström H, Sollid LM: Tissue transglutaminase selectively modifies gliadin peptides that are recognized by gut-derived T-cells in celiac disease. Nat Med. 1998, 4: 713-717. 10.1038/nm0698-713.

    Article  CAS  PubMed  Google Scholar 

  12. 12.

    Van de Wal Y, Kooy YMC, van Veelen P, Pena AS, Mearin ML, Papadopoulos GK, Koning F: Selective deamidation by tissue transglutaminase strongly enhances gliadin-specific T-cell reactivity. J Immunol. 1998, 161: 1585-1588.

    CAS  PubMed  Google Scholar 

  13. 13.

    Vader LW, De Ru A, Van der Wal Y, Kooy YMC, Benckhuijsen W, Mearin ML, Drijfhout JW, Van Veelen PA, Koning F: Specificity of tissue transglutaminase explains cereal toxicity in celiac disease. J Exp Med. 2002, 195: 643-649. 10.1084/jem.20012028.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  14. 14.

    Fleckenstein B, Molberg Ø, Qiao SW, Schmid DG, von der MF Mülbe, Elgstøen K, Jung G, Sollid LM: Gliadin T-cell epitope selection by tissue transglutaminase in celiac disease. J Biol Chem. 2002, 277: 34109-34116. 10.1074/jbc.M204521200.

    Article  CAS  PubMed  Google Scholar 

  15. 15.

    Van de Wal Y, Kooy YMC, Van Veelen P, Pẽna AS, Mearin ML, Molberg Ø, Lundin KE, Sollid LM, Mutis T, Benckhuijsen WB, Drijfhout JW, Koning F: Small intestinal T-cells of celiac disease patients recognize a natural pepsin fragment of gliadin. Proc Natl Acad Sci USA. 1998, 95: 10050-10054. 10.1073/pnas.95.17.10050.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  16. 16.

    de Wal Van, Kooy YMC, Van Veelen P, Vader W, August SA, Drijfhout JW, Pẽna SA, Koning F: Glutenin is involved in the gluten-driven mucosal T-cell response. Eur J Immunol. 1999, 29: 3133-3139. 10.1002/(SICI)1521-4141(199910)29:10<3133::AID-IMMU3133>3.0.CO;2-G.

    Article  Google Scholar 

  17. 17.

    Sjöström H, Lundin KE, Molberg O, Körner R, McAdam SN, Anthonsen D, Quarsten H, Norén O, Roepstorff P, Thorsby E, Sollid LM: Identification of a gliadin T-Cell epitope in coeliac disease: general importance of gliadin deamidation for intestinal T-cell recognition. Scand J Immunol. 1998, 48: 111-115. 10.1046/j.1365-3083.1998.00397.x.

    Article  PubMed  Google Scholar 

  18. 18.

    Arentz-Hansen H, Körner R, Molberg Ø, Quarsten H, Vader W, Kooy YMC, Lundin KEA, Koning F, Roepstorff P, Sollid LM, McAdam SN: The intestinal T-cell response to α-gliadin in adult celiac disease is focused on a single deamidated glutamine targeted by tissue transglutaminase. J Exp Med. 2000, 191: 603-612. 10.1084/jem.191.4.603.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  19. 19.

    Arentz-Hansen H, McAdam SN, Molberg Ø, Fleckenstein B, Lundin KEA, Jørgensen TJD, Jung G, Roepstorff P, Sollid LM: Celiac lesion T-cells recognize epitopes that cluster in regions of gliadins rich in proline residues. Gastroenterology. 2002, 123: 803-809. 10.1053/gast.2002.35381.

    Article  PubMed  Google Scholar 

  20. 20.

    Arentz-Hansen H, Fleckenstein B, Molberg Ø, Scott H, Koning F, Jung G, Roepstorff P, Lundin KEA, Sollid LM: The molecular basis for oat intolerance in patients with celiac disease. PLoS Med. 2004, 1: 84-92.

    Article  CAS  Google Scholar 

  21. 21.

    Stepniak D, Vader WL, Kooy Y, van Veelen PA, Moustakas A, Papandreou NA, Eliopoulos E, Drijfhout JW, Papadopoulos GK, Koning F: T-cell recognition of HLA-DQ2-bound gluten peptides can be influenced by an N-terminal proline at p-1. Immunogenetics. 2005, 57: 8-15. 10.1007/s00251-005-0780-8.

    Article  CAS  PubMed  Google Scholar 

  22. 22.

    Vader LW, Kooy Y, Van Veelen PA, De Ru A, Harris D, Benckhuijsen W, Pen S, Mearin ML, Drijfhout JW, Koning F: The gluten response in children with celiac disease is directed toward multiple gliadin and glutenin peptides. Gasteroenterology. 2002, 122: 1729-1737. 10.1053/gast.2002.33606.

    Article  CAS  Google Scholar 

  23. 23.

    Vader LW, Stepniak , Bunnik EM, Kooy YMC, De Haan W, Drijfhout JW, Van Veelen PA, Koning F: Characterization of cereal toxicity for celiac disease patients based on protein homology in grains. Gasteroenterology. 2003, 125: 1105-1113. 10.1016/S0016-5085(03)01204-6.

    Article  CAS  Google Scholar 

  24. 24.

    Camarca A, Anderson RP, Mamone G, Fierro O, Facchiano A, Costantini S, Zanzi D, Sidney J, Auricchio S, Sette A, Troncone R, Gianfrani C: Intestinal T-cell responses to gluten peptides are largely heterogeneous: implications for a peptide-based therapy in celiac disease. J Immunol. 2009, 182: 4158-4166. 10.4049/jimmunol.0803181.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  25. 25.

    Anderson RP, van Heel DA, Tye-Din JA, Jewell DP, Hill AVS: In vivo antigen challenge in celiac disease identifies a single transglutaminase-modified peptide as the dominant A-gliadin T-cell epitope. Nat Med. 2000, 6: 337-342. 10.1038/73200.

    Article  CAS  PubMed  Google Scholar 

  26. 26.

    Anderson RP, van Heel DA, Tye-Din JA, Jewell DP, Hill AVS: Antagonists and non-toxic variants of the dominant wheat gliadin T-cell epitope in coeliac disease. Gut. 2006, 55: 485-491. 10.1136/gut.2005.064550.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  27. 27.

    Qiao S-W, Bergseng E, Molberg Ø, Xia J, Fleckenstein B, Khosla C, Sollid LM: Antigen presentation to celiac lesion-derived T-cells of a 33-mer gliadin peptide naturally formed by gastrointestinal digestion. J Immunol. 2004, 173: 1757-1762.

    Article  CAS  PubMed  Google Scholar 

  28. 28.

    Dørum S, Qiao SW, Sollid LM, Fleckenstein B: A quantitative analysis of transglutaminase 2-mediated deamidation of gluten peptides: Implications for the T-cell response in celiac disease. J Proteome Res. 2009, 8: 1748-1755. 10.1021/pr800960n.

    Article  PubMed  Google Scholar 

  29. 29.

    Dørum S, Arntzen MØ, Qiao S-W, Holm A, Koehler CJ, Thiede B, Sollid LM, Fleckenstein B: The preferred substrates for transglutaminase 2 in a complex wheat gluten digest are peptide fragments harboring celiac disease T-cell epitopes. PLoS One. 2010, 5: e14056-10.1371/journal.pone.0014056.

    PubMed Central  Article  PubMed  Google Scholar 

  30. 30.

    Spaenij-Dekking L, Kooy-Winkelaar Y, Van Veelen P, Drijfhout JW, Jonker H, van Soest L, Smulders MJM, Bosch D, gilissen LJWJ, Koning F, Bosch D: Natural variation in toxicity of wheat: potential for selection of nontoxic varieties for celiac disease patients. Gasteroenterology. 2005, 129: 797-806. 10.1053/j.gastro.2005.06.017.

    Article  CAS  Google Scholar 

  31. 31.

    Molberg O, Uhlen AK, Jensen T, Solheim Flaete N, Fleckenstein B, Arentz-Hanzen H, Raki M, Lundin KEA, Sollid LM: Mapping of gluten T-cell epitopes in the bread wheat ancestors: implications for celiac disease. Gasteroenterology. 2005, 128: 393-401. 10.1053/j.gastro.2004.11.003.

    Article  CAS  Google Scholar 

  32. 32.

    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-10.1186/1477-5956-9-10.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  33. 33.

    Payne PI: Genetics of wheat storage proteins and the effect of allelic variation on bread making quality. Annu Rev Plant Physiol. 1987, 38: 141-153. 10.1146/annurev.pp.38.060187.001041.

    Article  CAS  Google Scholar 

  34. 34.

    Gao S, Gu Y-Q, Wu J, Coleman-Derr D, Huo N, Crossman C, Jia J, Zuo Q, Ren Z, Anderson OD, Kong X: Rapid evolution and complex structural organization in genomic regions harboring multiple prolamin genes in the polyploid wheat genome. Plant Mol Biol. 2007, 65: 189-203. 10.1007/s11103-007-9208-1.

    Article  CAS  PubMed  Google Scholar 

  35. 35.

    Sabelli P, Shewry PR: Characterization and organisation of gene families at the Gli-1 loci of bread and durum wheats by restriction fragment analysis. Theor Appl Genet. 1991, 83: 209-216.

    Article  CAS  PubMed  Google Scholar 

  36. 36.

    Anderson OD, Hsia CC, Torres V: The wheat γ-gliadin genes: characterization of ten new sequences and further understanding of γ-gliadin gene family structure. Theor Appl Genet. 2001, 103: 323-330. 10.1007/s00122-001-0551-3.

    Article  CAS  Google Scholar 

  37. 37.

    Pistón F, Dorado G, Martín A, Barro F: Cloning of nine γ-gliadin mRNAs (cDNAs) from wheat and the molecular characterization of comparative transcript levels of γ-gliadin subclasses. J Cereal Sci. 2006, 43: 120-128. 10.1016/j.jcs.2005.07.002.

    Article  Google Scholar 

  38. 38.

    Zhang W, Gianibelli MC, Ma W, Rampling L, Gale KR: Identification of SNPs and development of allele-specific PCR markers for gamma-gliadin alleles in Triticum aestivum. Theor Appl Genet. 2003, 107 (1): 130-138.

    CAS  PubMed  Google Scholar 

  39. 39.

    Qi P-F, Wei Y-M, Ouellet T, Chen Q, Tan X, Zheng Y-L: The gamma-gliadin multigene family in common wheat (Triticum aestivum) and its closely related species. BMC Genomics. 2009, 10: 168-10.1186/1471-2164-10-168.

    PubMed Central  Article  PubMed  Google Scholar 

  40. 40.

    Huang Z, Long H, Wei YM, Qi PF, Tan ZH, Zheng YL: Characterization and classification of γ-gliadin multigene sequences from Aegilops section Sitopsis. Cereal Res Commun. 2010, 38 (1): 1-14. 10.1556/CRC.38.2010.1.1.

    Article  Google Scholar 

  41. 41.

    Altenbach SB, Vensel WH, DuPont FM: Analysis of expressed sequence tags from a single wheat cultivar facilitates interpretation of tandem mass spectrometry data and discrimination of gamma gliadin proteins that may play different functional roles in flour. BMC Plant Biology. 2010, 10: 7-10.1186/1471-2229-10-7.

    PubMed Central  Article  PubMed  Google Scholar 

  42. 42.

    Metakovsky EV: Gliadin allele identification in common wheat II. Catalogue of gliadin alleles in common wheat. J Genet Breed. 1991, 45: 325-344.

    Google Scholar 

  43. 43.

    Van Herpen TWJM, Goryunova SV, Schoot J, Mitreva M, Salentijn EMJ, Vorst O, Schenk MF, Van Veelen PA, Koning F, Van Soest LJM, Vosman B, Bosch D, Hamer RJ, Gilissen LJWJ, Smulders MJM: Alpha-gliadin genes from the A, B, and D genomes of wheat contain different sets of celiac disease epitopes. BMC Genomics. 2006, 7: 1-10.1186/1471-2164-7-1.

    PubMed Central  Article  PubMed  Google Scholar 

  44. 44.

    Van den Broeck HC, de Jong HC, Salentijn EMJ, Dekking L, Bosch D, Hamer RJ, Gilissen LJWJ, van der Meer IM, Smulders MJM: Presence of celiac disease epitopes in modern and old hexaploid wheat varieties: wheat breeding may have contributed to increased prevalence of celiac disease. Theor Appl Genet. 2010, 121: 1527-1539. 10.1007/s00122-010-1408-4.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  45. 45.

    Van den Broeck H, Chen H, Lacaze X, Dusautoir J-C, Gilissen L, Smulders M, Van der Meer I: In search of tetraploid wheat accessions reduced in celiac disease-related gluten epitopes. Mol Biosyst. 2010, 6: 2206-2213. 10.1039/c0mb00046a.

    Article  CAS  PubMed  Google Scholar 

  46. 46.

    Mitea C, Salentijn EMJ, van Veelen P, Goryunova SV, van der Meer IM, van den Broeck HC, Mujico JR, Monserrat V, Gilissen LJWJ, Drijfhout JW, Dekking L, Koning F, Smulders MJM: A universal approach to eliminate antigenic properties of alpha-gliadin peptides in celiac disease. PLoS One. 2010, 5 (12): e15637-10.1371/journal.pone.0015637.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  47. 47.

    Van Slageren MW: Wild wheats: a monograph of Aegilops L. and Amblyopyrum (Jaub. & Spach) Eig (Poaceae). Wageningen Agriculture University Papers, ICARDA, Wageningen. 1994, 94 (7): 1-513.

    Google Scholar 

  48. 48.

    Huang S, Sirikhachornkit A, Su X, Faris J, Gill B, Haselkorn R, Gornicki P: Genes encoding plastid acetyl-CoA carboxylase and 3-phosphoglycerate kinase of the Triticum/Aegilops complex and the evolutionary history of polyploid wheat. Proc Natl Acad Sci U S A. 2002, 99 (12): 8133-8138. 10.1073/pnas.072223799.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  49. 49.

    Feldman M, Lupton FGH, Miller TE: Wheats. Evolution of Crop Plants. Edited by: Smart J, Simmonds NW. 1995, Longman Group, London, 184-192.

    Google Scholar 

  50. 50.

    Salse J, Chagué V, Bolot S, Magdelenat G, Huneau C, Pont C, Belcram H, Couloux A, Gardais S, Evrard A, Segurens B, Charles M, Ravel C, Samain S, Charmet G, Boudet N, Chalhoub B: New insights into the origin of the B genome of hexaploid wheat: evolutionary relationships at the SPA genomic region with the S genome of the diploid relative Aegilops speltoides. BMC Genomics. 2008, 9: 555-567. 10.1186/1471-2164-9-555.

    PubMed Central  Article  PubMed  Google Scholar 

  51. 51.

    Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.

    Article  CAS  PubMed  Google Scholar 

  52. 52.

    Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.

    CAS  PubMed  Google Scholar 

  53. 53.

    Tamura K, Nei M, Kumar S: Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci U S A. 2004, 101: 11030-11035. 10.1073/pnas.0404206101.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  54. 54.

    Cassidy BG, Dvorak J, Anderson OD: The wheat low-molecular-weight glutenin genes: characterization of six new genes and progress in understanding gene family structure. Theor Appl Genet. 1998, 96: 743-750. 10.1007/s001220050797.

    Article  CAS  Google Scholar 

  55. 55.

    Sollid LM, Qiao SW, Anderson RP, Gianfrani C, Koning F: Nomenclature and listing of celiac disease relevant gluten T-cell epitopes restricted by HLA-DQ molecules. Immunogenetics. 64 (6): 455-460.

  56. 56.

    Shan L, Qiao S-W, Arentz-Hansen H, Molberg Ø, Gray GM, Sollid LM, Khosla C: Identification and analysis of multivalent proteolytically resistant peptides from gluten: implications for celiac sprue. J Proteome Res. 2005, 4 (5): 1732-1741. 10.1021/pr050173t.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  57. 57.

    Tye-Din JA, Steward 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, 2: 1-14.

    Article  Google Scholar 

  58. 58.

    Salentijn EMJ, Goryunova SV, Bas N, Van der Meer IM, Van den Broeck HC, Bastien T, Gilissen LJWJ, Smulders MJM: Tetraploid and hexaploid wheat varieties reveal large differences in expression of alpha-gliadins from homoeologous Gli-2 loci. BMC Genomics. 2009, 10: 48-10.1186/1471-2164-10-48.

    PubMed Central  Article  PubMed  Google Scholar 

  59. 59.

    Pizzuti D, Buda A, D'Odorico A, D'Incà R, Chiarelli S, Curioni A, Martines D: Lack of intestinalmucosal toxicity of Triticum monococcum in celiac disease patients. Scand J Gastroenterol. 2006, 41: 1305-1311. 10.1080/00365520600699983.

    Article  CAS  PubMed  Google Scholar 

Download references


This research was funded by the Celiac Disease Consortium, an Innovative Cluster approved by the Netherlands Genomics Initiative and partially funded by the Dutch Government (BSIK03009), and the DLO program 'Plant and Animal for Human Health' (KB-05-001-019-PRI).

Author information



Corresponding author

Correspondence to Elma MJ Salentijn.

Additional information

Competing interest

The authors declare that they have no competing interest.

Authors’ contributions

EMJS and MJMS designed the studies. EMJS carried out sequence analysis and wrote the manuscript. DCM carried out T-cell tests. SVG performed cloning and sequence analysis. ISMP and carried out sequence analysis and bioinformatics. FK and MJMS helped to coordinate the study and draft the manuscript. LJWJG and IMM helped to draft the manuscript. All authors read and approved the final manuscript. ISMP is currently employed at Université de Genève as a bio-informatician. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Ten full-length γ-gliadin sequences. Ten full length γ-gliadin nucleotide sequences obtained after the assemblage of 717 γ-gliadins transcripts at 98% homology. In brackets, the number of sequences in a contig. (TXT 9 KB)

CD-epitopes of γ-gliadin transcripts of

Additional file 2: T. aestivum in their natural context. Alignments of the deduced aminoacid sequences of T.aestivum γ-gliadin transcript contigs (717 transcripts from the Genbank NCBI) spanning a part of the first repetitive domain and a part of the γ-gliadin sequence. CD T-cell epitopes are depicted in bold: γ-I (PQQSFPQQQ), γ-III (QQPQQPYPQ), γ-IV (SQPQQQFPQ), γ-VI (QQPFPQQPQ), γ-VIIa (PQPQQQFPQ), γ-VIIb (QQPQQPFPQ), Glia-γ2a (FPQQPQQPF), 26-mer FLQPQQPFPQQPQQPYPQQPQQPFPQ. Gli-A: T.aestivum γ-gliadin transcripts expressed from locus Gli-A1, Neighbor Joining topology group 6 (N = 140 transcripts) and 7 (N = 38 transcripts). Gli-B: T.aestivum γ-gliadin transcripts expressed from locus Gli-B1, Neighbor Joining topology group 3 (N = 41 transcripts), 5 (N = 125 transcripts) and 8 (N = 19 transcripts). Gli-D: T.aestivum γ-gliadin transcripts expressed from locus Gli-D1, Neighbor Joining topology group 9 (N = 293 transcripts) and 10 (N = 61 transcripts). Alignment gaps are indicated with dashes (−). Shorter sequences, not connected to domain I are marked with #. Glutamine residues that are a primary targets for the enzyme tissue transglutaminase are underlined (Q) in QxP target sites whereas moderate target sites are depicted in italics (Q) [13, 14]. Variants of DQ2-γ-I are indicated. Cleavage sites: In black, trypsin; in grey with white letters, chymotrypsin-high specificity; in grey with black letters, chymotrypsin-low specificity; cleavage occurs at the right side (C-terminal direction) of the marked amino acid. (PDF 164 KB)

T-cell stimulating capacity of natural variants of celiac disease epitope, DQ2-γ-I.

Additional file 3: In vitro T-cell stimulating capacity of natural variants of celiac disease epitope, DQ2-γ-I. Stimulation of a T-cell clone specific for CD epitope DQ2-γ-I, with natural occurring variants of DQ2-γ-I (9-mer epitope core PQQSFPQQQ and residues in the positions −1 to −4 and +1 to +4). Glutamine residues that are a primary targets for the enzyme tissue transglutaminase are underlined (Q) in QxP target sites whereas moderate target sites are depicted in italics (Q) [13, 14]. SI = Stimulation Index = cpm of the stimulated culture/cpm unstimulated culture (with APC only); ++ = SI ≥ 50, + = 20 ≤ SI <50, ± =10 ≤ SI <20, - = SI <10. The stimulation is mediated by HLA-DQ2 carrying antigen presenting cells (APC). Peptide no. corresponds to the peptide numbering in Table 5. (DOCX 14 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Salentijn, E.M., Mitea, D.C., Goryunova, S.V. et al. Celiac disease T-cell epitopes from gamma-gliadins: immunoreactivity depends on the genome of origin, transcript frequency, and flanking protein variation. BMC Genomics 13, 277 (2012).

Download citation


  • Wheat
  • Gluten
  • γ-gliadins
  • Celiac disease
  • T-cell epitopes