- Research note
- Open Access
Characterization of genomic DNA sequence of the candidate gene for FB_Mfu10 associated with fire blight resistance in Malus species
BMC Research Notes volume 14, Article number: 291 (2021)
The proposed candidate gene underlying the Malus fusca fire blight resistance locus on chromosome 10 was previously predicted to possess 880 amino acids and 8 exons. Eight base pair (8 bp) insertion/deletion in the first exon potentially distinguished resistant genotypes from susceptible ones. This study aimed at analyzing the candidate gene sequence in another set of original resistant and susceptible progeny, characterizing the sequence in a transgenic line transformed with the candidate gene under its own native promoter, as well as deciphering the potential genomic differences between this candidate gene and its homolog in the ‘Golden Delicious’ doubled haploid genome (GDDH13).
Sequences of amplicons of part of the candidate gene amplified in two progenies that showed resistant and susceptible fire blight phenotypes, confirmed the 8 bp insertion that distinguishes susceptible and resistant progenies. The transgenic line was positive for the candidate gene sequence, confirming a successful transfer into the background of apple cultivar ‘Pinova’, and possessed the same genomic sequence as the progeny with a resistant phenotype. Sequence analysis showed that the homolog gene on GDDH13 possesses a significant 18 bp deletion in exon 1 leading to a difference of 15 amino acid from the protein sequence of the candidate gene.
Several Malus species, especially cultivars of the domesticated apple (Malus domestica Borkh.) are susceptible to fire blight disease—caused by the bacterium Erwinia amylovora [1,2,3]. Nevertheless, strong resistance is found in individual accessions of different wild Malus species —for example in M. ×robusta 5 [5, 6], M. floribunda 821, ‘Evereste’ , M. ×arnoldiana MAL0004  and M. fusca MAL0045 . Only in some of these genotypes have candidate genes been detected, which underlie mapped fire blight resistance loci on different linkage groups [10,11,12,13].
The candidate gene proposed to underlie the fire blight resistance locus of Malus fusca MAL0045 possesses B-lectin, PAN-apple and serine/threonine kinase domains, 8 exons and 880 amino acids . This candidate gene was predicted on the sequences of bacterial artificial chromosomes (BACs)—46H22/24N24, which spanned the resistance locus on linkage group 10 (LG10) . The donor of this gene, M. fusca accession MAL0045, is highly resistant to E. amylovora, with the resistance locus detected on LG10 following artificial inoculation of MAL0045 × ‘Idared’ progeny with different wild types and mutant strains of the bacterium . One of the apple reference genomes—the ‘Golden Delicious’ doubled haploid (GDDH13) sequence , was integral in characterizing the fire blight region on LG10 by aligning sequences of closely linked markers as well as BAC clones in order to characterize the region relative to this apple reference genome . Emeriewen et al.  amplified the open reading frame (ORF) of the gene in clones 46H22/24N24 spanning the region for resistance and in clone 94B13, spanning the homolog region, with sequencing of the amplicons detecting an 8 bp insertion present in clone 94B13.
To characterize the candidate gene sequence for FB_Mfu10, PCR was applied on M. fusca (MAL0045), ‘Idared’, and clones 46H22 and 94B13 to prove amplification of part of the ORF (3539 bp) as well as the entire ORF and its border sequences upstream and downstream (6366 bp) as described in . PCR was also applied on two progeny individuals of MAL0045 × ‘Idared’, 05210-147, which displays a resistant phenotype to E. amylovora strain Ea222, and 05210-89, which displays a susceptible phenotype to the same strain. Briefly, PCR was performed using the DreamTaq Hot Start PCR Master Mix Kit (Thermo Scientific, Berlin, Germany) in a 12 μl volume with the following conditions: 95 °C for 5 min, followed by 34 cycles of 95 °C for 1 min, 68 °C for 90 s and 72 °C for 6 min, and an extension of 72 °C for 10 min. Table 1 shows the primers used in this study. 20–40 ng of DNA were used for PCR . Purification of PCR products was performed using the MSB Spin PCRapace kit (Invitek GmbH, Berlin, Germany) according to the manufacturer’s protocol. Eurofins MWG Operon (Ebersberg, Germany) conducted sequencing of the amplification products. The resulting sequences were analysed using the National Centre for Biotechnology (NCBI) BLAST® program or by aligning two sequences with each other. Gene prediction analyses with sequences was done using FGENESH with algorithms for dicot Arabidopsis. Predicted proteins were analyzed as described in .
To obtain a transgenic line, fire blight susceptible apple cultivar ‘Pinova’ was transformed with a construct carrying the candidate gene sequence. This construct was obtained by cloning the ORF and its border sequences into a D034p9-Dao-FLPi_STPK (FB_Mfu10) vector (Additional file 1: File S1). Cloning was performed by DNA cloning service (Hamburg, Germany). Agrobacterium-mediated transformation and regeneration were carried out according to the established protocols of Julius Kühn Institute, Dresden-Pillnitz, Germany, as described by Flachowsky et al. .
DNA was isolated from desired samples using Qiagen DNeasy Plant Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. The transgenic line was tested with nptII primer pair and with candidate gene-specific primers (Table 1). PCR analyses using nptII primer pair and FB_Mfu10-Rseq-ORF primer pair were performed in a 20-µl volume consisting of 1 × Dream Taq buffer (ThermoFisher Scientific, Darmstadt, Germany), 0.2 mM dNTPs, 1 µM each of forward and reverse primers, 0.5 U Dream Taq DNA polymerase (ThermoFisher Scientific, Darmstadt, Germany) and 2-µl of DNA (20 ng). Running PCR conditions for nptII test were 94 °C for 4 min, followed by 25 cycles of 94 °C for 30 s, 56 °C for 1 min and 72 °C for 1 min and an extension of 72 °C for 7 min. Running PCR conditions for FB_Mfu10-Rseq-ORF were 94 °C for 4 min, followed by 32 cycles of 94 °C for 30 s, 67 °C for 1 min and 72 °C for 1 min and an extension of 72 °C for 5 min. To amplify the ORF in the transgenic line with MVH_FB_Mfu10 primer pair, CloneAmp™ Hifi PCR Premix kit (Takara, Saint-Germain-en-Laye, France) was used according to the manufacturer’s instruction with a running PCR program of 98 °C for 30 s, followed by 34 cycles of 98 °C for 10 s, 65 °C for 15 s and 72 °C for 30 s and an extension of 72 °C for 10 min.
Total RNA was extracted from transgenic line and M. fusca in vitro seedling using InviTrap® Spin Plant RNA Mini Kit according to the manufacturer’s protocol (Stratec, Berlin, Germany). DNA contamination removal and cDNA synthesis were performed using the Invitrogen DNA-free™ Kit (Thermo Scientific, Berlin, Germany) and the RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific, Berlin, Germany), respectively, and according to the manufacturers’ protocols. The conditions for amplifying other gene-specific primers are as described earlier .
Only in M. fusca and clone 46H22 was the ORF together with its border regions (6366 bp) amplified. A fragment of the ORF (~ 3539 bp) was amplified in clone 94B13, whereas the susceptible parent ‘Idared’ failed to amplify any gene specific fragments tested (Fig. 1a). Similarly, the 6366 bp was amplified in the resistant progeny, 05210-147, but not in the susceptible progeny, 05210-89 (Fig. 1b). However, similar ORF fragments were amplified in both progenies (Fig. 1c). Sequencing of amplicons amplified by 46H22 and 94B13 confirmed the 8 bp insertion, which is possessed by the susceptible clone. In parallel, amplicon sequencing confirmed that the susceptible progeny, 05210-89 also possesses this 8 bp insertion but not the resistant progeny, 05210-147 (Fig. 1d).
The obtained sequences from amplicons of the aforementioned samples were aligned against each other and against the predicted sequence of the candidate gene. In addition, the predicted sequence of the candidate gene, 3555 bp in total, was aligned against the GDDH13 genome using BLASTN  via the Rosaceae website available at https://www.rosaceae.org/. The output from this alignment showed a best hit of 97% identity (3456/3554 bp) and gaps of 1% (23/3554). The best-hit alignment of these sequences (Additional file 2: File S2) were in the region between 30,708,564 and 30,712,095 on chromosome 10 of the GDDH13 genome . In comparison to the candidate gene sequence, there is a significant 18 bp deletion in the GDDH13 homolog sequence (Additional file 2: File S2). Gene prediction analyses with the homolog sequence from GDDH13 predicted an ORF with the same domains as the candidate gene but with 15 amino acids difference (i.e. 865 aa—Additional file 3: File S3).
A transgenic line, M2117, was successfully generated following agrobacterium-mediated transformation of apple cultivar ‘Pinova’ with the construct carrying the candidate gene with its border sequences. The correct nptII fragment (780 bp) in the construct was amplified in M2117 as well as the candidate gene-specific primer, which amplifies a section of the border sequence upstream and part of the ORF, and the 3539 bp of almost the entire ORF (Fig. 2a). Sequencing of amplicons amplified by M2117 showed the same characteristics of 05210-147 unique to resistant progenies (Fig. 2b). RT-PCR also showed M2117 was positive for the candidate gene (Fig. 2c).
It is quite clear from these results and as initially reported by Emeriewen et al.  that 8 bp differentiates fire blight susceptible individuals from resistant individuals within the MAL0045 × ‘Idared’ progeny. However, the failure of the susceptible progeny to amplify the ORF plus border sequences strongly suggests that the border (promoter) regions of this gene might play a key role in fire blight resistance. It was reported that this 8 bp causes a difference of 28 amino acids in the protein sequence, which is potentially sufficient to affect the corresponding phenotype. Nonetheless, we cannot rule out the role of the promoter. The apple cultivar ‘Pinova’ used for genetic engineering in this study is susceptible to fire blight . The transfer of FB_Mfu10 candidate gene sequence into ‘Pinova’ genetic background led to the generation of M2117 transgenic line, who’s genomic DNA amplified gene-specific fragments, which ‘Pinova’ itself failed to amplify. It is logical that sequences of amplicons amplified by M2117 are the same as the 05210-147, which possesses a resistant phenotype to E. amylovora strain Ea222. However, it would be hasty to suggest or jump to the conclusion that M2117 will possess a resistant phenotype. For this, a phenotypic evaluation of this transgenic line is necessary.
Furthermore, the apple reference genome, GDDH13  showed strong identity with the candidate gene sequence, as shown in this study. However, noteworthy is the 18 bp deletion and several other SNPs within the homolog gene sequence. It is in fact interesting that this indel also situates in the first exon of the homolog gene. The 18 bp deletion in GDDH13, which causes 15 amino acid difference in the protein level, is however potentially sufficient to affect the corresponding phenotypes, or render the homolog gene inactive or non-functional. One base pair switch in the Pi-d2 gene was sufficient to distinguish resistance and susceptibility to rice blast  and a single nucleotide mutation in E. amylovora was sufficient to breakdown the fire blight resistance of M. ×robusta 5 . Although, there is sequence identity between the gene and the homolog on GDDH13, the apple cultivar ‘Golden Delicious’ is highly susceptible to fire blight  whereas MAL0045—the source of FB_Mfu10, is highly resistant to different strains of the pathogen .
Although the aforementioned results confirm the results of Emeriewen et al.  and reinforces this receptor kinase gene as a strong fire blight candidate gene, the generation of a single transgenic line so far and the failure to report its phenotypic evaluation is a limitation of the study. This is because the line is still too small for phenotypic evaluation. Moreover, this study primarily aimed at characterizing the genomic sequences of this line and MAL0045-derived F1 progeny, to confirm the initial study . A larger scale study to generate more transgenic lines is underway. Further studies will aim at generating more transgenic lines overexpressing this candidate gene in ‘Pinova’ as well as in ‘Gala’—another susceptible cultivar, and more importantly, evaluate whether it is functional. These experiments are in progress and we intend to publish the results regardless of the outcome. These results are highly anticipated as to date the only functionally proven fire blight resistance gene in Malus is FB_MR5 , isolated from M. ×robusta 5 fire blight resistance locus region on LG3 [5, 11], and which is already broken down by highly virulent and mutant strains of E. amylovora .
Availability of data and materials
Data generated from this study are published within this article. Sequences obtained in this study were aligned against the GDDH13 genome . Further materials can be provided on request from the corresponding author, Ofere Francis Emeriewen.
Golden delicious doubled haploid genome
Bacterial artificial chromosome
Open reading frame
Polymerase chain reaction
Reverse transcription PCR
Norelli JL, Jones AL, Aldwinckle HS. Fire blight management in the twenty-first century: using new technologies that enhance host resistance in apple. Plant Dis. 2003;87(7):756–65.
Malnoy M, Martens S, Norelli JL, Barny MA, Sundin GW, Smits TH, Duffy B. Fire blight: applied genomic insights of the pathogen and host. Annu Rev Phytopathol. 2012;50:475–94.
Peil A, Emeriewen OF, Khan A, Kostick S, Malnoy M. Status of fire blight resistance breeding in Malus. J Plant Pathol. 2020. https://doi.org/10.1007/s42161-020-00581-8.
Emeriewen OF, Wöhner T, Flachowsky H, Peil A. Malus hosts–Erwinia amylovora interactions: strain pathogenicity and resistance mechanisms. Front Plant Sci. 2019;10:551.
Peil A, Garcia-Libreros T, Richter K, Trognitz FC, Trognitz B, Hanke MV, Flachowsky H. Strong evidence for a fire blight resistance gene of Malus robusta located on linkage group 3. Plant Breed. 2007;126(5):470–5.
Peil A, Hübert C, Wensing A, Horner M, Emeriewen OF, Richter K, Wöhner T, Chagné D, Orellana-Torrejon C, Saeed M, Troggio M, et al. Mapping of fire blight resistance in Malus ×robusta 5 flowers following artificial inoculation. BMC Plant Biol. 2019;19(1):1–2.
Durel CE, Denance C, Brisset MN. Two distinct major QTL for resistance to fire blight co-localize on linkage group 12 in apple genotypes ‘Evereste’ and Malus floribunda clone 821. Genome. 2009;52(2):139–47.
Emeriewen OF, Peil A, Richter K, Zini E, Hanke MV, Malnoy M. Fire blight resistance of Malus ×arnoldiana is controlled by a quantitative trait locus located at the distal end of linkage group 12. Eur J Plant Pathol. 2017;148(4):1011–8.
Emeriewen O, Malnoy M, Richter K, Kilian A, Hanke MV, Peil A. Evidence of a major QTL for fire blight resistance in the apple wild species Malus fusca. In: XIII international workshop on fire blight 1056. 2013. p. 289–93.
Parravicini G, Gessler C, Denance C, Lasserre-Zuber PA, Vergne E, Brisset MN, Patocchi A, Durel CE, Broggini GA. Identification of serine/threonine kinase and nucleotide-binding site–leucine-rich repeat (NBS-LRR) genes in the fire blight resistance quantitative trait locus of apple cultivar ‘Evereste.’ Mol Plant Pathol. 2011;12(5):493–505.
Fahrentrapp J, Broggini GA, Kellerhals M, Peil A, Richter K, Zini E, Gessler C. A candidate gene for fire blight resistance in Malus ×robusta 5 is coding for a CC–NBS–LRR. Tree Genet Genomes. 2013;9(1):237–51.
Emeriewen OF, Richter K, Flachowsky H, Malnoy M, Peil A. Genetic analysis and fine mapping of the fire blight resistance locus of Malus ×arnoldiana on linkage group 12 reveal first candidate genes. Front Plant Sci. 2021. https://doi.org/10.3389/fpls.2021.667133.
Emeriewen OF, Richter K, Piazza S, Micheletti D, Broggini GA, Berner T, Keilwagen J, Hanke MV, Malnoy M, Peil A. Towards map-based cloning of FB_Mfu10: identification of a receptor-like kinase candidate gene underlying the Malus fusca fire blight resistance locus on linkage group 10. Mol Breed. 2018;38(8):1–4.
Emeriewen O, Richter K, Kilian A, Zini E, Hanke MV, Malnoy M, Peil A. Identification of a major quantitative trait locus for resistance to fire blight in the wild apple species Malus fusca. Mol Breed. 2014;34(2):407–19.
Emeriewen OF, Richter K, Berner T, Keilwagen J, Schnable PS, Malnoy M, Peil A. Construction of a dense genetic map of the Malus fusca fire blight resistant accession MAL0045 using tunable genotyping-by-sequencing SNPs and microsatellites. Sci Rep. 2020;10(1):1.
Daccord N, Celton JM, Linsmith G, Becker C, Choisne N, Schijlen E, Van de Geest H, Bianco L, Micheletti D, Velasco R, Di Pierro EA. High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development. Nat Genet. 2017;49(7):1099–106.
Flachowsky H, Peil A, Sopanen T, Elo A, Hanke V. Overexpression of BpMADS4 from silver birch (Betula pendula Roth) induces early-flowering in apple (Malus × domestica Borkh). Plant Breed. 2007;126(2):137–45.
Zhang Z, Schwartz S, Wagner L, Miller W. A greedy algorithm for aligning DNA sequences. J Comput Biol. 2000;7(1–2):203–14.
Vogt I, Wöhner T, Richter K, Flachowsky H, Sundin GW, Wensing A, Savory EA, Geider K, Day B, Hanke MV, Peil A. Gene-for-gene relationship in the host–pathogen system Malus × robusta 5–E rwinia amylovora. New Phytol. 2013;197(4):1262–75.
Chen X, Shang J, Chen D, Lei C, Zou Y, Zhai W, Liu G, Xu J, Ling Z, Cao G, Ma B. AB-lectin receptor kinase gene conferring rice blast resistance. Plant J. 2006;46(5):794–804.
Broggini GA, Wöhner T, Fahrentrapp J, Kost TD, Flachowsky H, Peil A, Hanke MV, Richter K, Patocchi A, Gessler C. Engineering fire blight resistance into the apple cultivar ‘Gala’using the FB _ MR 5 CC-NBS-LRR resistance gene of Malus × robusta 5. Plant Biotechnol J. 2014;12(6):728–33.
We acknowledge DFG for providing funding. We also thank Ute Hille, Astrid Paufler, and Ines Hiller for excellent technical assistance.
Open Access funding enabled and organized by Projekt DEAL. This research is funded by DFG (Deutsche Forschungsgemeinschaft) research grant received by Ofere Francis Emeriewen—Project Number AOBJ 661177.
Ethics approval and consent to participate
Consent for publication
We confirm that all authors read and approved this manuscript for publication.
The authors’ declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: File S1.
Binary plasmid vector harbouring the candidate gene under its own promoter sequence used for transformation. STPK = serine/threonine protein kinase – candidate gene sequence.
Additional file 2: File S2.
Alignment of FB_Mfu10 candidate gene sequence against the GDDH13 genome. The significant 18 bp indel is highlighted in red. Two nucleotides ‘GA’ after position 3553 shown in this alignment completes FB_Mfu10 sequence (i.e. 3555 bp).
Additional file 3: File S3.
Predicted 865 amino acid from GDDH13 homolog. See Emeriewen et al.  for the predicted proteins of this candidate gene.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Emeriewen, O.F., Flachowsky, H. & Peil, A. Characterization of genomic DNA sequence of the candidate gene for FB_Mfu10 associated with fire blight resistance in Malus species. BMC Res Notes 14, 291 (2021). https://doi.org/10.1186/s13104-021-05709-2
- Malus fusca fire blight resistance locus