Open Access Research article

Assessment of allelic diversity in intron-containing Mal d 1 genes and their association to apple allergenicity

Zhongshan Gao123, Eric W van de Weg23*, Catarina I Matos2, Paul Arens3, Suzanne THP Bolhaar4, Andre C Knulst4, Yinghui Li35, Karin Hoffmann-Sommergruber6 and Luud JWJ Gilissen23

Author Affiliations

1 Department of Horticulture/Allergy Research Center, Zhejiang University, Hangzhou 310029, PR China

2 Allergy Consortium Wageningen, Wageningen University and Research Centre, P.O. Box 16, 6700AA, Wageningen, the Netherlands

3 Plant Research International, Wageningen University and Research Centre, P.O. Box 16, 6700AA, Wageningen, the Netherlands

4 Department of Dermatology/Allergology, University Medical Center Utrecht, P.O. Box 85500, 3508GA Utrecht, the Netherlands

5 The National Key Facility for Crop Gene Resources and Genetic Improvement (NFCRI)/Key Lab of Germplasm & Biotechnology (MOA), Institute of Crop Science, China Academy of Agricultural Science, Beijing, 100081, PR China

6 Department of Pathophysiology, Medical University of Vienna, AKH-EBO-3Q, Währinger Gürtel 18-20, A-1090 Vienna, Austria

For all author emails, please log on.

BMC Plant Biology 2008, 8:116  doi:10.1186/1471-2229-8-116

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2229/8/116


Received:28 November 2007
Accepted:13 November 2008
Published:13 November 2008

© 2008 Gao et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Mal d 1 is a major apple allergen causing food allergic symptoms of the oral allergy syndrome (OAS) in birch-pollen sensitised patients. The Mal d 1 gene family is known to have at least 7 intron-containing and 11 intronless members that have been mapped in clusters on three linkage groups. In this study, the allelic diversity of the seven intron-containing Mal d 1 genes was assessed among a set of apple cultivars by sequencing or indirectly through pedigree genotyping. Protein variant constitutions were subsequently compared with Skin Prick Test (SPT) responses to study the association of deduced protein variants with allergenicity in a set of 14 cultivars.

Results

From the seven intron-containing Mal d 1 genes investigated, Mal d 1.01 and Mal d 1.02 were highly conserved, as nine out of ten cultivars coded for the same protein variant, while only one cultivar coded for a second variant. Mal d 1.04, Mal d 1.05 and Mal d 1.06 A, B and C were more variable, coding for three to six different protein variants. Comparison of Mal d 1 allelic composition between the high-allergenic cultivar Golden Delicious and the low-allergenic cultivars Santana and Priscilla, which are linked in pedigree, showed an association between the protein variants coded by the Mal d 1.04 and -1.06A genes (both located on linkage group 16) with allergenicity. This association was confirmed in 10 other cultivars. In addition, Mal d 1.06A allele dosage effects associated with the degree of allergenicity based on prick to prick testing. Conversely, no associations were observed for the protein variants coded by the Mal d 1.01 (on linkage group 13), -1.02, -1.06B, -1.06C genes (all on linkage group 16), nor by the Mal d 1.05 gene (on linkage group 6).

Conclusion

Protein variant compositions of Mal d 1.04 and -1.06A and, in case of Mal d 1.06A, allele doses are associated with the differences in allergenicity among fourteen apple cultivars. This information indicates the involvement of qualitative as well as quantitative factors in allergenicity and warrants further research in the relative importance of quantitative and qualitative aspects of Mal d 1 gene expression on allergenicity. Results from this study have implications for medical diagnostics, immunotherapy, clinical research and breeding schemes for new hypo-allergenic cultivars.

Background

Many birch pollen sensitised patients (50–70%) in Central and Northern Europe suffer from oral allergy symptoms after eating fresh apples [1]. The prevalence of apple allergic individuals mounts up to ~3% in Central and Northern Europe. This type of apple allergy is caused by cross reactivity of IgE antibodies against the major and sensitizing birch pollen allergen Bet v 1 with Mal d 1, the major allergen of apple. Bet v 1 and Mal d 1 are both pathogenesis-related (PR) proteins. They belong to the PR-10 family and share a high degree of homology [2-6].

From patients' experience it is known for a long time that the severity of allergic reactions to apple was not only related to the specific sensitivity of the individual, but also largely depended on the apple cultivar. This cultivar dependent allergenicity has also been described in literature. For instance, Mal d 1 from the cultivar Golden Delicious was found highly reactive to specific IgE antibodies from allergic patients' sera, whereas Mal d 1 from the cultivar Gloster generally showed much less reactivity [7,8]. In addition, skin prick testing (SPT) with 21 different apple cultivars and confirmations for specific cultivars in double-blind placebo-controlled food challenges (DBPCFC) and oral challenges of whole apples, revealed a wide range of allergenic reactivity from very high to very low [9,10]. As a result from these studies, the new cultivar Santana was identified as hypo- allergenic for 75% of the patients with a mild apple allergy [10], which is usually assumed to be Mal d 1 based. In The Netherlands (where birch pollen-related apple allergy is by far the most common form of apple allergy), this cultivar has recently been marketed as 'suited for individuals with mild apple allergy' in order to meet the general desire of apple allergic persons to be able to add this common fruit to their daily diet.

The differences in allergenicity among cultivars raised a crucial question on the origin of this cultivar-specific degree of allergenicity. Allergenicity may depend on the total amount of Mal d 1 proteins, as suggested by Son et al. [11] from their observed ten-fold difference in Mal d 1 amount between the high-allergenic cultivar Golden Delicious and the low-allergenic cultivar Gloster. However, there is little evidence supporting this hypothesis because only very few cultivars have been studied and, for these, a linear response between Mal d 1 protein content and allergenicity estimates is lacking. On the other hand, qualitative characteristics of the Mal d 1 proteins could be involved too, as can be argued from the differences in binding capacity of birch pollen-specific IgE to two protein variants of Mal d 1 [11-13]. To elaborate this latter issue, research on the genetic variation of Mal d 1 and its expression pattern in the different cultivars is required and should be compared to allergenicity data. It is known that Mal d 1 is coded by a large gene family of 18 members mapped on three linkage groups of the apple genome [14,15]. Not all of these members are likely to be involved in allergenicity since only a limited number of different Mal d 1 proteins and mRNAs have been traced back in apple fruit so far [16-18].

Research towards the relative importance of the quality and quantity of Mal d 1 proteins on the allergenicity of apple cultivars is relevant for designing apple breeding programs for low-allergenic apple cultivars of high quality and healthiness. In this paper, we focused on the genetic diversity of Mal d 1 genes. The Mal d 1 gene family can be subdivided into two major categories: genes with and genes without an intron. Preliminary genetic analyses revealed that the genetic diversity was by far larger in the intron-containing genes. Furthermore, the intron-containing genes cover all three linkage groups that Mal d 1 loci [15]. Therefore, this category of genes has been chosen to start looking for putative qualitative effects of Mal d 1 proteins in cultivar specific allergenicity.

Allelic diversity of the seven intron-containing Mal d 1 genes was assessed among a set of cultivars chosen for their importance in breeding programs and apple production. In order to find putative associations with allergenicity, the presence of alleles coding for different protein variants was subsequently compared with the degree of allergenicity for a subset of cultivars for which allergenicity data from SPT or DBPCFC tests were available.

Results

Diversity of Mal d 1 genes and deduced proteins

The observed DNA polymorphisms in the 10 studied cultivars resulted in a total of 46 different Mal d 1 sequences over seven genes (Table 1, 2). These sequences were denoted according to the occurrence of 1) DNA polymorphisms in the coding region of the gene leading to different protein variants; 2) polymorphism in the coding region that did not affect the protein sequence (silent mutations), and 3) polymorphism in the intron (Table 1, 2). Although the latter two differences are of minor importance with respect to allergenicity, they provided additional landmarks for the development of sequence specific molecular markers.

Table 1. Genetic variation in the intron containing Mal d 1 genes of linkage groups 6 (Mal d 1.05) and 13 (Mal d 1.01) among ten apple cultivars.

Table 2. Genetic variation in the intron containing Mal d 1 genes of linkage group 16 among ten apple cultivars.

Mal d 1.01 and -1.02 showed to be highly conserved at the protein level. The related genes coded for only two variants each that both differed in just a single amino acid (pos. 135 V/A for Mal d 1.01, pos. 56 N/K for Mal d 1.02) and for which the second variant was found only once. The other genes were more variable, coding for three (Mal d 1.05) up to six variants (Mal d 1.06C). Mal d 1.04 showed to be special in that two (out of the three) sequences contained a stop codon in the coding region and were therefore regarded as pseudo alleles (ps1 and ps2). Interestingly, the pseudo alleles occurred frequently as for seven out of the ten cultivars at least one of the alleles was a pseudo allele whereas cultivars Priscilla and Fuji only contained pseudo alleles of Mal d 1.04 (Table 2).

Allergenicity scores of 14 apple cultivars by skin prick test (SPT)

Relative SPT responses of 14 apple cultivars are given in Table 3. Fiesta, Delblush, Pinova and Golden Delicious were ranked in the high allergenic group (83–100%). Priscilla and Santana showed low SPT responses, with wheal areas 30–35% of that of Golden Delicious. Nine cultivars were intermediate (48–72%) allergenic. Santana was also identified as low-allergenic in comparison to Golden Delicious in DBPCFC tests [9] and oral provocation tests [10].

Table 3. SPT responses of apple cultivars relative to Golden Delicious (in %) for 4 experiments.

General associations from the sequenced cultivars

From the cultivars used to sequence the intron containing Mal d 1 genes and to perform SPT on allergenicity, Golden Delicious was ranked as the highest allergenic cultivar whereas Priscilla was ranked as the lowest (Table 3). This difference was not related to protein variant composition of Mal d 1.01, 1.05, 106B and 1.06C because these cultivars have identical protein variants (Table 2). In contrast, both Mal d 1.04 and Mal d 1.06A have different protein variants in Golden Delicious and Priscilla. Mal d 1.04 may contribute to the observed difference in allergenicity between the two cultivars since protein variant 04 is present in Golden Delicious while Priscilla is heterozygous for the two pseudo-alleles ps1 and ps2. Similarly, Mal d 1.06A protein variant -02 was found in Priscilla whereas variants -01 and -03 were found in Golden Delicious. The intermediate cultivars, Jonathan, Prima, and Fuji, have the putatively high allergenic alleles -01 or -03 in single dosage, whereas high allergenic cultivars either have both alleles -01 and -03, or have -01 in homozygous state. This suggests that allergenicity might be determined by both the protein variant and the gene dosage of highly allergenic alleles.

Association analysis by pedigrees: from Golden Delicious to Santana

The identity and origin of genomic alleles and thus protein variants in additional cultivars (not sequenced for Mal d 1) could be traced by developed allele specific SNAP and SSR markers and the use of pedigree information [19] (see Methods). For instance, the deduced flow of protein variants over the pedigree of cultivar Santana is presented in Figure 1. Santana and Priscilla are low allergenic whereas Golden Delicious is high allergenic [9,10]. For Mal d 1.01, 1.02 and 1.05 the same protein variants were found for Golden Delicious and Santana. In contrast, Golden Delicious and Santana differ in their protein variant composition of Mal d.1.04, -1.06A and -1.06B (Fig. 1), indicating a possible involvement of these proteins in the observed difference in allergenicity between these cultivars. Santana, like Priscilla, has only pseudo alleles for Mal d 1.04 that do not result in protein production, while Golden Delicious has one functional allele. The Mal d 1.06A02 variant of Santana is not present in Golden Delicious.

thumbnailFigure 1. Flow of putative protein variants over the pedigree of Santana for seven Mal d 1 genes. Abbreviated cultivar names are according to Table 1. GD is high allergenic whereas PS and ST are low allergenic [9,10].

Associations for Mal d 1.04 and -1.06A

Similarly to Santana, pedigree information and sequence specific markers were used to assess the putative protein variant composition for other additional cultivars for which we have SPT response data. This allowed a further validation of the association of the Mal d 1.04 and 1.06A with allergenicity in a total of 14 cultivars (Table 4). Cultivars Delblush and Pinova evoked an SPT response that was similar to that of Golden Delicious and also had the same variant compositions as Golden Delicious for Mal d 1.04 and Mal d 1.06A, namely variant 04 for Mal d 1.04 (coded for by a single allele) and the two variants 01 and 03 for Mal d 1.06. The remaining five cultivars evoking an intermediate allergenic response also had a similar variant composition as the previously identified intermediate cultivars, namely variants 02 and 03 for Mal d 1.06A, and variant 04 of Mal d 1.04 (coded for by a single allele dosage), except for Prima, that has two alleles coding for Mal d 1.0404 and Fuji that contains both pseudo alleles. The other intron containing genes did not show any association (data not shown).

Table 4. SPT responses of 14 apple cultivars and the putative protein constitutions of these cultivars for two Mal d 1 iso-allergens.

Discussion

Birch pollen induced oral allergy for apple has been the subject in a considerable number of studies. One of the prominent results has been the presence of cultivar-specific differences in allergenicity. Unfortunately, evidence regarding the causes of cultivar-specific allergenicity is still lacking. One of the knowledge gaps concerned the number and identity of Mal d 1 genes and the amount of variation within these genes. Recently, Gao et al. [15] have shown that Mal d 1 genes are members of a large gene family by identifying 18 different loci that are located in three clusters. Based on sequence identity, these 18 genes could be subdivided into intron containing and intronless genes. In order to create a basis for a better understanding of the genetics of Mal d 1 genes and their impact on allergenicity, we have studied the allelic diversity of the intron containing genes in 10 cultivars that are often used in breeding. Development of sequence specific markers and pedigree information enabled the assessment of putative Mal d 1 constitutions of other cultivars. Using this information, we assessed the different Mal d 1 isoforms that cultivars are able to produce and found associations between their putative protein constitutions and SPT-responses.

Allelic diversity and validity of database sequences

Cloning and sequencing of the seven intron-containing Mal d 1 genes in 10 cultivars revealed 46 different alleles that coded for 25 different Mal d 1 isoforms. The variation per gene varied with regard to the number of alleles and deduced proteins. Mal d 1.01 and Mal d 1.02 were diverse at the gDNA level but conserved at the protein level. For both genes only two protein variants were found, of which the second variant, differing in one amino acid only, was found in just a single cultivar. For Mal d 1.04, Mal d 1.05 and the three Mal d 1.06 genes, gDNAs often coded for different proteins and these genes were therefore more variable at the protein level.

Because the examined cultivars are important in the breeding of many modern apple varieties, the set of alleles found in this study likely represents a considerable part of the total variation present in intron containing Mal d 1 genes of common apple varieties.

Although other Mal d 1 sequences are known from public databases, we suspect that many of these sequences may be artefacts derived through strand switching and PCR mutations. The problem associated with PCR amplification of a group of closely related sequences, such as the Mal d 1 gene family, is that besides PCR induced single base pair mutations, in vitro strand switching or re-annealing of incompletely amplified fragments can lead to artefacts as was exemplified by Schenk et al. [20] for birch Bet v1 sequences. For instance, for Mal d 1.01 one of the most studied Mal d 1 genes, over 13 DNA sequences from previous studies are known from public databases (Table 5) indicating the presence of 9 putative protein isoforms. We know Mal d 1.01 is a single locus gene with maximum two alleles present in a cultivar [15], but four sequences from Golden Delicious can be found in databases. From these, only sequence accession AF124830 was identical to one of our two sequences (Mal d 1.0105.01b). The other sequences may be due to artefacts. Firstly, accession AF126402 had one SNP at position 11 (G→A) compared to AF124830, which is due to the cloning primer used. Similarly, sequences from a number of other cultivars showed this 11A mutation too. In our study, the cloning primers used were positioned in the 5'-untranslated region thus avoiding this problem. Secondly, Accession Z48969[21] that was classified as Mal d 1c [10] is quite similar to Mal d 1.01, while the SNPs found are unique to three other Mal d 1 genes. This may be due to an artefact from strand switching during PCR which could explain why we, like other researchers could not retrieve this sequence from Golden Delicious [[10,17]; this study]. Similarly, the first 84 nucleotides of AF124829 and AF124832 match to Mal d 1.02 (typical SNPs are shown at the bottom of Table 5), while the remaining sequence matches to Mal d 1.01 making these accessions also putative erroneous sequences due to strand switching. Therefore, the actual number of Mal d 1.01 variants present among apple cultivars is probably not as high as nine, because at least five of them are likely PCR artefacts. Two protein variants (1.0105 and 1.0109) of Mal d 1.01 could be confirmed in this study. Variants Mal d 1.0104, 1.0107 and 1.0108 are likely artefacts of 1.0105 whereas 1.0101 and 1.0102 are likely artefacts of 1.0109.

Table 5. SNPs among Mal d 1.01 sequences as found in databases.

The occurrence of PCR recombination and mutations in sequences from gene families warrants scrutinised assessment of sequences. The use of two independent PCR-cloning steps for each cultivar may effectively filter out most of these erroneous sequences before database donation since the probability of isolating identical artefacts in independent PCRs is low [20]. Sequence specific markers may be used to validate newly found isoforms that have passed this first sifting. Many of the sequences found in this study were either confirmed by identical sequences retrieved from the other cultivars used or by identical sequences previously donated in the databases as well as through the use of sequence specific markers [[15], this study].

Cultivar specific allergenicity and its relation to quantitative and qualitative differences in Mal d 1

Levels of total Mal d 1 protein can not fully explain allergenicity of cultivars, indicating that other factors are involved as well. For instance, Son et al. [10] found a fourfold higher level of Mal d 1 protein in Granny Smith compared to Gloster. Although this suggests a correlation between total amount and allergenicity, this relation was contradicted by their results on Golden Delicious and Granny Smith, which cultivars were similar allergenic despite of a threefold difference in total Mal d 1 content. Furthermore, it was demonstrated that different Mal d 1 isoforms [10], as well as mutants of specific isoforms [13], have different binding affinities to IgE indicating that the relative Mal d 1 protein composition as well as the presence of specific alleles may have a significant effect on allergenicity.

This research is the first to show an association between the genetic Mal d 1 constitution of apple cultivars and birch pollen related allergenicity. If the allergenicity is a matter of the variants of a single gene, then Mal d 1.06A will be the most outstanding candidate. In case of the involvement of multiple genes, also Mal d 1.04 may play a role. A putative involvement of the intronless genes can not yet be excluded, as their allelic variation has not yet been surveyed. Such a survey will be a challenge considering the high sequence similarity among alleles of different loci, which can be over 98% [15]. With regard to Mal d 1.04, two of the sequenced cultivars showed variant compositions for Mal d 1.04 that did not exactly fit with the relative SPT data found. Prima has the functional allele of Mal d 1.04 in duplex but has an intermediate allergenicity, whereas Fuji is the only intermediate-allergenic cultivar that had two pseudo-alleles for Mal d 1.04 (Table 2). This might indicate that Mal d 1.04 does not show dosage allele effects and cannot explain allergenicity without considering other factors. Alternatively, these two exceptions infer variant 03 of Mal d 1.06A to be more allergenic than variant 01, in which case the stronger effect of 03 in Fuji is counteracted by the absence of a functional allele for Mal d 1.04. Similarly, the milder effect of variant 01 of Prima is then compensated by the double dosage of a functional Mal d 1.04 allele. A stronger effect of variant 03 fits with the tendency of higher responses of the three heterozygous 01/03 cultivars compared to the two 01 homozygous cultivars.

The above described associations could be found due to the presence of allelic variation among the examined apple cultivars and due to performing a complexity reduction of the human variation by only analysing patients with mild SPT responses, thus reducing the effect of variation among humans for sensitivity to different allergens. Studies with larger patient sizes may probably benefit from further grouping to also account for genetically determined human variation in sensitivity to different (iso) allergen variants. Such grouping has probably to be based on allergy responses as no knowledge exists on the involved human genes neither on their allelic composition.

The finding that allergenicity depends on the presence and amount of some specific Mal d 1 isoforms is highly relevant for diagnostics tests and immunotherapy, and justifies additional research on a larger number of apple cultivars as well as atopic individuals. Since the first Mal d 2 and Mal d 4 genes have also been recently mapped [22] and the mapping of additional genes of these allergens is in progress, it will become possible to also investigate the effects of allelic composition of these Mal d allergens on the allergenicity of cultivars by association studies.

Location of amino acid polymorphism in a 3D structure model

For Mal d 1.06A, high-allergenic cultivars have two putative genotypes, homozygous variant 01 or heterozygous variant 01 together with variant 03, whereas low-allergenic cultivars are homozygous for variant 02. The intermediate-allergenic cultivars contained the low allergenic variant 02 in combination with one of the high allergenic variants 01 or 03.

The three Mal d 1.06A variants differ at two amino acids: 13 V/I and 135V/A. Considering the three dimensional structure model of Mal d 1 [23], the first polymorphism is located in the first loop between the β1-strand and the α1-helix, the second is located in the α3-helix structure motif. The amino acid changes are all between hydrophobic amino acids but they have different side chains that may have an effect on the 3D structure of the protein and thus on epitope conformation.

Expression of Mal d 1 genes in fruit

For specific Mal d 1 genes to be involved in allergenicity, expression in apple fruit is a prerequisite. Until now, mRNA expression for five genes was observed in mature fruit through both rtPCR [16,24] and EST sequences (unpublished), representing two genes (Mal d 1.01 and -03E) on LG 13 and three genes (Mal d 1.02, -06A and -06B) on LG 16., whereas mRNA based EST sequences of Mal d 1.04 were only found in mature leaves. Because the number of rtPCR studies on and EST sequences derived from fruit is still limited and usually only assessed at a single time moment, definitive mRNA based conclusions regarding Mal d 1.04 expression in fruit can not yet be made. Recently obtained proof for its expression in the epidermis of fruit [24] is not conclusive as the primers used can probably also amplify mRNA of Mal d 1.05 and Mal d 1.06B. From ESTs identified in cDNA libraries of ripe Gala fruit, tissue specific expression patterns were found. Mal d 1.01 and Mal d 1.02 ESTs were found in both skin and cortex, while Mal d 1.06A, 1.06B and 1.03E ESTs were only found in the skin of apple fruit. At this moment, there is evidence for the presence of two Mal d 1 proteins in apple fruit [17,18], the majority of Mal d 1 protein is Mal d 1.02 (Mal d 1b) and a minor part is Mal d 1.06A [15]. Interestingly, both genes are located on linkage group 16 where also Mal d 1.04 is located. These mRNA and protein data thus allow Mal d 1.06A to be involved in differences in allergenicity among cultivars. The current lack of support for the presence of Mal d 1.04 in fruit might indicate that the observed association is coincidental, but may as well be due to lack of extensive expression studies.

Genotyping for Mal d 1 haplotypes

The Mal d 1 genes in LG 16 are tightly linked to each other [15]. This tight linkage can simplify the genotyping of additional cultivars, at least if their pedigree and the linkage phase of their parental alleles are known. In these cases, genotyping can be performed by a single representative, multi-allelic marker such as the Mal d 1.06A SSR marker. As linkage phases of the Mal d 1 genes of LG 16 are known for all 10 cultivars of our reference set but Discovery (Table 2), this simple and efficient approach was performed in this study for certain cultivars (Figure. 1).

Conclusion

We have shown that differences in allergenicity among apple cultivars are associated with the allelic composition of two specific genes, Mal d 1.04 and Mal d 1.06A, which are both located on linkage 16 of Malus domestica. Furthermore, allele dosage effects are found relevant for Mal d 1.06A. Our findings indicate the need to reconsider the relevance of merely assessing total amounts of Mal d 1 protein in allergy research and diagnostic tests and warrant further research on the association of specific Mal d 1 isoforms and allergenicity among a larger group of cultivars and allergy sufferers.

Methods

Cultivars for cloning and sequencing

Eight cultivars were used for cloning and sequencing of Mal d 1 isoallergen genes: Golden Delicious (GD), Priscilla (PS), Ingrid Marie (IM), Cox (CO), Jonathan (JO), Red Delicious (RD), Fuji (FJ) and Discovery (DS). Including the two parental cultivars Prima (PM) and Fiesta (FS) used in our previous mapping study [15], a total of ten cultivars were under investigation. They were chosen for three reasons: (1) GD, JO, CO and RD are founders in many breeding programmes world wide, and RD, GD and FJ are main cultivars for apple production in the world; (2) GD, IM and PS are members in the pedigree from GD (hyper-allergenic) to Santana (ST) = [PS × Elstar (= GD × IM)] a hypo-allergenic cultivar; (3) JO, DS, PM and FS are the parents of three mapping progenies.

Genomic cloning and sequencing of intron-containing Mal d 1 genes

For each cultivar, all seven intron-containing Mal d 1 genes were cloned and sequenced: Mal d 1.05 of linkage group (LG) 6; Mal d 1.01 of LG 13; and Mal d 1.02, -1.04, -1.06A, -1.06B, -1.06C of LG 16. Six primer pairs were used for cloning [15], of which those for Mal d 1.04 and Mal d 1.05 were newly designed (Table 6). The PCR amplification, cloning and sequencing procedures were described previously [15,22]. For all 10 cultivars, eight to ten clones for each gene were sequenced in both directions. Next, sequences were aligned and putative Single Nucleotide Polymorphisms (SNPs) were identified using the Seqman program (DNAstar, Madison, WI). The coding sequences were deduced and translated into amino acid sequences for alignments and assessment of protein variant with the GeneDoc program http://www.psc.edu/biomed/genedoc webcite. New protein variants or gDNA alleles were named according to Gao et al. [15,22], following a modification of the allergen nomenclature guidelines [25].

Table 6. Cloning primers and PCR conditions.

Cultivars for association studies

Allergenicity data were available for 6 out of the 10 cultivars for which we assessed allelic diversity [9]. Besides these six cultivars, eight additional cultivars were included in the association study. For these eight cultivars allergenicity data were available [9] and their allelic constitutions of the intron-containing Mal d 1 genes could be assessed by their pedigree relationships to the set of 10 sequenced cultivars.

For this, sequence-specific SNAP markers [23,27] were developed (Table 7) and applied to these cultivars. Primer design and test procedures were described previously [15,27]. Primers for the Mal d 1.06A SSR marker [15] were redesigned according to new sequence data obtained. The PCR reaction mixture for this SSR marker consisted of 2 μl 10× buffer, 1.2 μl MgCl2 (25 mM), 0.4 μl dNTPs (10 mM), 1μl of each primer (2 μM), 0.06 μl Taq polymerase (5 U/μl) and 1μl DNA (10 ng/μl) in a total volume of 20 μl. After an initial denaturation at 94°C for 2.5 min, the amplification was carried out for 34 cycles at 94°C for 30 s, 60°C for 30 s and 72°C for 1 min, and a final extension at 72°C for 5 min. PCR products were analysed on an ABI 377 (Applied Biosystems, Foster City, Calif.).

Table 7. Description of nine sequence-specific SNAP markers and one SSR marker used for the genotyping of cultivars (in bold: SNP, underlined are deliberate introduced SNP to increase specificity).

Using the SNAP makers, the pedigree structures allowed us to follow the flow of the Mal d 1 alleles over generations by applying the Identity by Descent principle in the genotyping of cultivars [19]. In total 14 cultivars were thus available for association of protein variant composition with SPT responses.

Allergenicity data

In this study, Skin Prick Test (SPT) responses were used to evaluate allergenicity. The SPT procedure and the history of the patients have been described previously [9]. In short, patients were recruited from the outpatient clinic of the department of Dermatology/Allergology of the UMCU. They all had birch pollinosis manifesting with rhinoconjunctivitis during the birch pollen season (April and May), as well as a positive SPT to fresh apple of at least half the diameter of the positive histamine control. All patients had a typical history of apple allergy, with oral allergy syndrome (OAS) symptoms like itching and mild swelling of the mouth, throat and sometimes rhinoconjuctivitis after eating an apple. SPT were performed on the flexor surface of the forearm using the prick-to-prick-technique according to Dreborg [28,29].

Histamine dihydrochloride (10 mg/ml) was used as a positive control, and the glycerol diluents of the SPT-extracts were used as negative control (ALK-ABELLO, Nieuwegein, The Netherlands). The wheal reaction (a small, itching elevation of the skin, as from the bite of an insect) was marked and transferred with transparent adhesive tape to a record sheet. The skin wheal areas were measured by computer scanning [30]. SPT responses for each cultivar were standardized by dividing the original wheal area of the prick by that obtained from the reference cultivar for high allergenicity Golden Delicious and multiplied by 100. Data have been derived from four experiments, three of which had been published previously [9]. For each experiment we used only a fraction of the data, this is only the data of patients with mild symptoms as preliminary experiments indicated that Mal d 1 is the major allergen to these patients, while other Mal d proteins seem to be major allergens to patients displaying more severe symptoms (Van de Weg, unpublished). Consequently, only 25%–50% of the patients of the previous experiments [9] and only 50% of the patients (4 out of 8) of the forth, new experiment were included. Finally, 11 different patients were involved. In order to combine data from these different experiments, responses of cultivars were expressed as a percentage of the response against the reference cultivar Golden Delicious. The final ranking results were obtained by averaging the responses from four different experiments. This study was reviewed and approved by the Ethics Committee of the University Medical Center Utrecht under document number 01–050. All patients provided written informed consent before enrolment in the study.

Authors' contributions

LJWJG and KHS initiated this study. LJWJG, KHS, EVW and ZSG designed the experiment. ZSG, CIM and PA performed the gene cloning and sequences analysis. STHPB and ACK conducted the skin prick experiments. YHL tested and analysed the SSR markers. EVW performed the association analyses. ZSG, EVW, PA and LJWJG drafted the paper.

Acknowledgements

This study has been carried out with financial support from the Ministry of Agriculture, Nature and Food Quality of The Netherlands, from the China 111 Project (B06014) and from the Commission of the European Communities for the projects SAFE (QLK1-CT-2000-01394) and ISAFRUIT. ISAFRUIT is funded under the Thematic Priority 5-Food Quality and Safety of the 6th Framework Programme of RTD (Contract no. FP6-FOOD-CT-2006-016279). The views and opinions expressed in this publication are purely those of the writers and may not in any circumstances be regarded as stating an official position of the European Commission and in no way anticipates its future policy in this area.

References

  1. Eriksson NE, Formgren H, Svenonius E: Food hyper-sensitivity in patients with pollen allergy.

    Allergy 1982, 37:437-443. PubMed Abstract | Publisher Full Text OpenURL

  2. Ebner C, Birkner T, Valenta R, Rumpold H, Breitenbach M, Scheiner O, Kraft D, Sakamoto T, Hayashi Y, Yamada M, Torii S, Urisu A, Vieths S, Schoning B, Petersen A, Janek K, Aulepp H, Rudeschko O, Fahlbusch B, Henzgen M, Schlenvoigt G, Herrmann D, Jager L: Common epitopes of birch pollen and apples: Studies by Western and Northern blot.

    J Allergy Clin Immunol 1991, 88:588-594. PubMed Abstract | Publisher Full Text OpenURL

  3. Vieths S, Janek K, Aulepp H, Petersen A: Isolation and characterization of the 18-kDa major apple allergen and comparison with the major birch pollen allergen (Bet v I).

    Allergy 1995, 50:421-430. PubMed Abstract | Publisher Full Text OpenURL

  4. Schoning B, Ziegler WH, Vieths S, Baltes W: Apple allergy: the cDNA sequence of the major allergen of apple, determined by performing PCR with a primer based on the N-terminal amino acid sequence, is highly homologous to the sequence of the major birch pollen allergen.

    J Sci Food Agri 1996, 71:475-482. Publisher Full Text OpenURL

  5. Vieths S: Allergenic cross-reactivity, food allergy and pollen.

    Environ Toxicol Chem 1997, 4:61-70. OpenURL

  6. Holm J, Baerentzen G, Gajhede M, Ipsen H, Larsen JN, Lowenstein H, Wissenbach M, Spangfort MD: Molecular basis of allergic cross-reactivity between group 1 major allergens from birch and apple.

    J Chromatogr B 2001, 756:307-313. Publisher Full Text OpenURL

  7. Vieths S, Jankiewicz A, Schoning B, Aulepp H: Apple allergy: the IgE-binding potency of apple strains is related to the occurrence of the 18-kDa allergen.

    Allergy 1994, 49:262-271. PubMed Abstract | Publisher Full Text OpenURL

  8. Hsieh LS, Moos M, Lin Y: Characterization of apple 18 and 31 kd allergens by microsequencing and evaluation of their content during storage and ripening.

    J Allergy Clin Immunol 1995, 96:960-970. PubMed Abstract | Publisher Full Text OpenURL

  9. Bolhaar STHP, Weg WE, Van Ree R, Gonzalez-Mancebo E, Zuidmeer L, Bruijnzeel-Koomen CAFM, Fernandez-Rivas M, Jansen J, Hoffmann-Sommergruber K, Knulst AC, Gilissen LJWJ: In vivo assessment with prick-to-prick testing and double-blind, placebo-controlled food challenge of allergenicity of apple cultivars.

    J Allergy Clin Immunol 2005, 116(5):1080-1086. PubMed Abstract | Publisher Full Text OpenURL

  10. Kootstra HS, Vlieg-Boerstra BJ, Dubois AEJ: Assessment of the reduced allergenic properties of the Santana apple.

    Ann Allergy Asthma Immunol 2007, 99:522-525. PubMed Abstract OpenURL

  11. Son DY, Scheurer S, Hoffmann A, Haustein D, Vieths S: Pollen-related food allergy: Cloning and immunological analysis of isoforms and mutants of Mal d 1, the major apple allergen, and Bet v 1, the major birch pollen allergen.

    Eur J Nutri 1999, 38:201-215. Publisher Full Text OpenURL

  12. Ferreira F, Hirtenlehner K, Jilek A, Godnik CJ, Breiteneder H, Grimm R, Hoffmann Sommergruber K, Scheiner O, Kraft D, Breitenbach M, Rheinberger Hans J, Ebner C: Dissection of immunoglobulin E and T lymphocyte reactivity of isoforms of the major birch pollen allergen bet v 1: Potential use of hypoallergenic isoforms for immunotherapy.

    J Exp Med 1996, 183:599-609. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Ma Y, Gadermaier G, Bohle B, Bolhaar S, Knulst A, Markovic-Housley Z, Breiteneder H, Briza P, Hoffmann-Sommergruber K, Ferreira F: Mutational analysis of amino acid positions crucial for IgE-binding epitopes of the major apple (Malus domestica) allergen, Mal d 1.

    Int Arci Allergy 2006, 139:53-62. Publisher Full Text OpenURL

  14. Atkinson RG, Perry J, Matsui T, Ross GS, MacRae EA: A stress-, pathogenesis-, and allergen-related cDNA in apple fruit is also ripening-related.

    NZ J Crop Hort Sci 1996, 24:103-107. OpenURL

  15. Gao ZS, Weg WE, Schaart JG, Schouten HJ, Tran DH, Kodde L, Meer IM, Geest AHM, Kodde J, Breiteneder H, Hoffmann-Sommergruber H, Bosch D, Gilissen LJWJ: Genomic cloning and linkage mapping of the Mal d 1 (PR-10) gene family in apple (Malus domestica).

    Theor Appl Genet 2005, 111:171-183. PubMed Abstract | Publisher Full Text OpenURL

  16. Puehringer H, Zinoecker I, Marzban G, Katinger H, Laimer M: MdAP, a novel protein in apple, is associated with the major allergen Mal d 1.

    Gene 2003, 321:173-183. PubMed Abstract | Publisher Full Text OpenURL

  17. Beuning LL, Bowen JH, Persson HA, Barraclough D, Bulley S, MacRae EA: Characterisation of Mal d 1-related genes in Malus.

    Plant Mol Biol 2004, 55:369-388. PubMed Abstract | Publisher Full Text OpenURL

  18. Helsper JPFG, Gilissen LJWJ, van Ree R, America Antoine HP, Cordewener Jan HG, Bosch D: Quadrupole time-of-flight mass spectrometry: A method to study the actual expression of allergen isoforms identified by PCR cloning.

    J Allergy Clin Immunol 2002, 110:131-138. PubMed Abstract | Publisher Full Text OpenURL

  19. Weg WE, Voorrips RE, Finkers R, Kodde LP, Jansen J, Bink MCAM: Pedigree genotyping: A new pedigree-based approach of QTL identification and allele mining.

    Acta Hortic 2004, 663:45-50. OpenURL

  20. Schenk MF, LJWJ Gilissen, GD Esselink, MJM Smulders: Seven different genes encode a diverse mixture of isoforms of Bet v 1, the major birch pollen allergen.

    BMC Genomics 2006, 7:168. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  21. Vanek-Krebitz M, Hoffmann-Sommergruber K, Laimer da Camara Machado M, Susani M, Ebner C, Kraft D, Scheiner O, Breiteneder H: Cloning and sequencing of Mal d 1, the major allergen from apple (Malus domestica), and its immunological relationship to Bet v 1, the major birch pollen allergen.

    Biochem Biophys Res Commun 1995, 214:538-551. PubMed Abstract | Publisher Full Text OpenURL

  22. Gao ZS, Weg WE, Schaart JG, Van Arkel G, Breiteneder H, Hoffmann-Sommergruber K, Gilissen LJWJ: Genomic characterization and mapping of the apple allergen genes Mal d 2 (thaumatin-like protein) and Mal d 4 (profilin).

    Theor Appl Genet 2005, 111:1087-1097. PubMed Abstract | Publisher Full Text OpenURL

  23. Neudecker P, Schweimer K, Nerkamp J, Scheurer S, Vieths S, Sticht H, Rösch P: Allergic cross-reactivity made visible: solution structure of the major cherry allergen Pru av 1.

    J Biol Chem 2001, 276:22756-22763. PubMed Abstract | Publisher Full Text OpenURL

  24. Botton A, Lezzer P, Dorigoni A, Barcaccia G, Ruperti B, Ramina A: Genetic and environmental factors affecting allergen-related gene expression in apple fruit (Malus domestica L. Borkh).

    J Agric Food Chem 2008, 56:6707-6716. PubMed Abstract | Publisher Full Text OpenURL

  25. King TP, Hoffman D, Lowenstein H, Marsh David G, Platts Mills Thomas AE, Thomas W: Allergen nomenclature.

    J Allergy Clini Immunol 1995, 96:5-14. Publisher Full Text OpenURL

  26. Drenkard E, Richter BG, Rozen S, Stutius LM, Angell NA, Mindrinos M, Cho RJ, Oefner PJ, Davis RW, Ausubel FM: A simple procedure for the analysis of single nucleotide polymorphisms facilitates map-based cloning in Arabidopsis.

    Plant Physiol 2000, 124:1483-1492. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Gao ZS, Weg WE, Schaart JG, Meer IM, Kodde L, Laimer M, Breiteneder H, Hoffmann-Sommergruber K, Gilissen LJWJ: Linkage map positions and allelic diversity of two Mal d 3 (non-specific lipid transfer protein) genes in the cultivated apple (Malus domestica).

    Theor Appl Genet 2005, 110:479-491. PubMed Abstract | Publisher Full Text OpenURL

  28. Dreborg S, Foucard T: Allergy to apple, carrot and potato in children with birch pollen allergy.

    Allergy 1983, 38:167-172. PubMed Abstract | Publisher Full Text OpenURL

  29. Dreborg S: EAACI position paper.

    Allergy 1993, 48(suppl 14):49-82. Publisher Full Text OpenURL

  30. Poulsen LK, Liisberg C, Bindslev-Jensen C, Malling HJ: Precise determination of skin-prick tests: validation of a scanning device and software for a personal computer.

    Clin Exp Allergy 1993, 23:61-68. PubMed Abstract | Publisher Full Text OpenURL

  31. Brownstein MJ, Carpten JD, Smith JR: Modulation of nontemplated nucleotide addition by tag DNA polymerase: primer modifications that facilitate genotyping.

    Biotechniques 1996, 20:1004-1010. PubMed Abstract OpenURL