- Research article
- Open Access
Characterization of the MLO gene family in Rosaceae and gene expression analysis in Malus domestica
BMC Genomics volume 15, Article number: 618 (2014)
Powdery mildew (PM) is a major fungal disease of thousands of plant species, including many cultivated Rosaceae. PM pathogenesis is associated with up-regulation of MLO genes during early stages of infection, causing down-regulation of plant defense pathways. Specific members of the MLO gene family act as PM-susceptibility genes, as their loss-of-function mutations grant durable and broad-spectrum resistance.
We carried out a genome-wide characterization of the MLO gene family in apple, peach and strawberry, and we isolated apricot MLO homologs through a PCR-approach. Evolutionary relationships between MLO homologs were studied and syntenic blocks constructed. Homologs that are candidates for being PM susceptibility genes were inferred by phylogenetic relationships with functionally characterized MLO genes and, in apple, by monitoring their expression following inoculation with the PM causal pathogen Podosphaera leucotricha.
Genomic tools available for Rosaceae were exploited in order to characterize the MLO gene family. Candidate MLO susceptibility genes were identified. In follow-up studies it can be investigated whether silencing or a loss-of-function mutations in one or more of these candidate genes leads to PM resistance.
Powdery mildew (PM) is a major fungal disease for thousands of plant species , including cultivated Rosaceae such as apple (Malus domestica), apricot (Prunus armeniaca), peach (Prunus persica), and strawberry (Fragaria x ananassa). Powdery mildew occurs in all major growing regions of Rosaceous crops, leading to severe losses . The major PM causal agents are Podosphaera leucotricha in apple , Sphaerotheca pannosa var. persicae in peach , Podosphaera tridactyla in apricot  and Podosphaera aphanis (syn. Sphaerotheca macularis f. sp. fragariae) in strawberry . Powdery mildew shows similar symptoms in the four species: white spots appear on young green tissues, particularly leaves in the first days after opening, whereas mature leaves show some resistance. Infected leaves crinkle, curl, and prematurely drop. Blossoms and fruits are not the primary targets of PM fungi, but infections of these tissues are possible [2, 3, 5]. In peach, apricot and apple, PM spores overwinter in buds and then in spring, with the reprise of vegetative growth, the spores start a new infection [2, 3].
Cultivars resistant to PM are fundamental in order to reduce the use of pesticides in agricultural practices. The usual strategy in breeding focuses on dominant plant resistance genes (R-genes), however these genes often originate from wild-relatives of the cultivated species, and thus interspecific crossability barriers could prevent their introgression . Moreover, in case of a successful cross, several undesirable traits are incorporated with the R-gene, making extensive backcrossing necessary, which is time-consuming in woody species. Finally, the durability of R-genes is generally limited due to the appearance of virulent strains of the pathogen, which can overcome resistance in a few years . Two examples are Venturia inaequalis race 6, which is able to overcome Rvi6 resistance to scab in apple , and P. leucotricha strains able to breakdown Pl-1 and Pl-2, two major PM R-genes of apple .
An alternative to the use of R-genes is based on plant susceptibility genes (S-genes), defined as genes whose loss-of-function results in recessively inherited resistance . Barley mlo PM resistance, first characterized in 1942, is a remarkable example of immunity due to the absence of an S-gene, as it derives from a loss-of-function mutation in the MLO (M ildew L ocus O) gene, encoding for a protein with seven transmembrane domains [11, 12]. mlo resistance has long been considered as a unique form of immunity, characterized by durability, broad-spectrum effectiveness and recessive inheritance . However, the characterization of the sources of resistance in other plant species, such as Arabidopsis, pea [15, 16] and tomato , has revealed that resistance resulting from loss-of-function mutations in MLO functional orthologs is more common than previously thought. Therefore, it has been suggested that the inactivation of MLO susceptibility genes could represent a valid strategy to introduce PM resistance across a broad range of cultivated species .
Histological characterization of mlo resistance revealed that it is based on a pre-penetration defense system, associated with the formation of cell-wall appositions [14, 18] and at least partially dependent on the actin cytoskeleton . It has been suggested that functional MLO proteins negatively regulate vesicle-associated and actin-dependent defense pathways at the site of attempted PM penetration . MLO proteins are therefore targeted by PM fungi as a strategy to induce pathogenesis. The early stages of PM infection are associated with an increase in transcript abundance of MLO susceptibility genes, showing a peak at six hours after inoculation. This has been shown to occur in tomato , barley , pepper  and grape [23, 24].
MLO susceptibility genes are members of a gene family which shows tissue specific expression patterns and are involved in a variety of physiological processes, besides the response to PM fungi: one of the 15 MLO genes of Arabidopsis, AtMLO7, is involved in pollen tube reception by the embryo sac and its mutation results in reduced fertility . Two other Arabidopsis genes, named AtMLO4 and AtMLO11, are involved in the control of root architecture, as mutants with null alleles of these two genes display asymmetrical root growth and exaggerated curvature .
Previous phylogenetic analysis of the MLO protein family identified six clades . Currently, all MLO proteins functionally related to PM susceptibility in dicot species appear in a single clade, namely Clade V [14, 17, 23, 24]. Similarly, Clade IV harbours all characterized PM susceptibility proteins from monocots [20, 27].
MLO genes have been intensively studied in many monocots and dicots, but very little has been performed in Rosaceae. In this investigation, we characterized the MLO gene family in a number of Rosaceous species, with respect to their structural, genomic and evolutionary features. Moreover, we monitored the transcript abundances of apple MLO homologs following P. leucotricha inoculation in three apple cultivars.
In silico and in vitro characterization of Rosaceae MLO homologs
A database search for Rosaceae MLO homologs produced 21 significant matches in peach, 23 in strawberry and 28 in apple. Of these, six (five from M. domestica and one from F. vesca) showed a very limited alignment region with other MLO genes, whereas eight (two from M. domestica, two from P. persica and four from F. vesca) were characterized by markedly different length with respect to MLO homologs reported in the genomes of Arabidopsis and grapevine [23, 28], i.e. less than 350 amino acids (aa) or more than 700 aa. Details on genomic localization amino acid number, putative transmembrane domains and predicted exon/intron structure of the remaining homologs, together with information about the MLO homologs nomenclature chosen in this study is provided in Tables 1, 2 and 3.
Peach and apricot are evolutionary very close to each other, and show a high degree of homology in DNA sequence. Phylogenetic analysis (see next paragraph) indicated peach homologs PpMLO1, PpMLO3 and PpMLO4 as candidates for being required for PM susceptibility. Therefore, we used the sequences of these genes to design primers to identify full-length apricot MLO genes. This approach resulted in the amplification and the successive characterization of three MLO sequences, which were by analogy named PaMLO1, PaMLO3, and PaMLO4 (deposited in the NCBI database with the accession numbers KF177395, KF177396, and KF177397, respectively).
Phylogenetic relations and inference of orthology
We performed a phylogenetic study on the newly identified Rosaceae MLO proteins. The dataset was completed with four homologs recently characterized in Rosa hybrida (rose)  (RhMLO1, RhMLO2, RhMLO3 and RhMLO4), the complete Arabidopsis thaliana AtMLO protein family , a series of MLO homologs which have been functionally associated with PM susceptibility, namely tomato (Solanum lycopersicum) SlMLO1 , pea (Pisum sativum) PsMLO1 [15, 16], pepper (Capsicum annuum) CaMLO2 , lotus (Lotus japonicus) LjMLO1 , barrel clover (Medicago truncatula) MtMLO1 , barley (Hordeum vulgare) HvMLO , rice (Oryza sativa) OsMLO2 , wheat (Triticum aestivum) TaMLO_B1 and TaMLO_A1b , and grapevine (Vitis vinifera) VvMLO14, the only dicot MLO homolog known to belong to clade IV . Clustering analysis using the UPGMA algorithm resulted in a total of eight distinct clades and no divergent lineage (Figure 1). Clade numbers from I to VI were assigned based on the position of Arabidopsis AtMLO homologs and barley HvMLO, according to the previous study of Feechan et al.. The two additional clades (named VII and VIII) were found to include Rosaceae MLO homologs only, both having one homolog from P. persica, one from F. vesca and one from M. domestica. Further clustering analysis with a Neighbour-Joining algorithm resulted in merging clade VII and VIII (not shown).
Four apple MLO homologs (MdMLO5, MdMLO7, MdMLO11 and MdMLO19) and three MLO homologs from peach (PpMLO1, PpMLO3 and PpMLO4), apricot (PaMLO1, PaMLO3 and PaMLO4) and woodland strawberrry (FvMLO1, FvMLO4 and FvMLO12) were found to cluster together in the phylogenetic clade V, containing all the dicot MLO proteins experimentally shown to be required for PM susceptibility (e.g. [16, 23]. One homolog from strawberry (FvMLO17) and one from peach (PpMLO12) were found to group, together with grapevine VvMLO14, in clade IV, which contains all monocot MLO proteins acting as PM susceptibility factors (Figure 1).
We used the GBrowse-Syn tool to detect syntenic blocks encompassing P. persica, F. vesca and M. domestica MLO genes. As syntenic blocks derive from the evolution of the same chromosomal region after speciation, orthology between MLO genes could be inferred. In total, twelve orthologous relationships were predicted between P. persica and F. vesca, nine between P. persica and M. domestica and eight between F. vesca and M. domestica (Table 4, Figure 2 and Additional file 1).
Transcription of putative apple MLO genes in response to Podosphaera leucotrichainoculation
To identify MLO genes that respond to the PM fungus P. leucotricha, we measured the transcript abundance of 19 out of 21 apple MLO genes in leaves 4, 6, 8 and 24 hours after artificial inoculation with the pathogen, and compared these data with the ones of non-inoculated leaves. Three cultivars, Golden Delicious, Braeburn and Gala, were analysed in order to investigate whether up-regulation was comparable among them and results could therefore be generalized for all apple cultivars. Three genes, namely MdMLO11, MdMLO18 and MdMLO19, were found to be significantly up-regulated after inoculation with the pathogen (Figure 3 and Additional file 2). Up-regulation of these genes was about 2-fold compared to non-inoculated plants, with peaks of 4-fold up-regulation at very early time points (‘Braeburn’- MdMLO11 - 6 hpi; ‘Gala’- MdMLO18 - 4 hpi; ‘Golden Delicious’- MdMLO19 - 6hpi). MdMLO11 and MdMLO18 were up-regulated in all cultivars, while MdMLO19 was only up-regulated in ‘Braeburn’ and ‘Golden Delicious’.
Two of the genes, MdMLO11 and MdMLO19 belong to Clade V, while MdMLO18 belongs to the newly identified Clade VII (Figure 1).
Genomic organization and phylogenetic relations between Rosaceae MLOhomologs
We report here the identification, through an in silico approach, of 19 MLO homologs in the genome of peach and 18 in the genome of strawberry. This is consistent with the results of previous genome-wide studies carried out on dicotyledonous species, indicating the presence of 15 MLO homologs in Arabidopsis, 17 in grapevine and 16 in tomato [9, 13]; Appiano et al., unpublished results; . Conversely, the number of MLO homologs detected in apple (21) is lower than expected, considering that a relatively recent genome-wide duplication event has occurred in the Pyreae tribe .
Most PpMLO, FvMLO and MdMLO homologs appeared to be widely distributed within the respective genomes (Tables 1, 2 and 3), indicating segmental duplication as the prevailing evolutionary mechanism for the Rosaceae MLO gene family. However, we also found cases of clusters of adjacent homologs (PpMLO3, PpMLO8 and PpMLO18, PpMLO12 and PpMLO16, PpMLO1 and PpMLO14, FvMLO3 and FvMLO4, FvMLO6 and FvMLO7, MdMLO2, MdMLO3 and MdMLO8), which are likely the result of tandem duplication events.
Inference of phylogenetic relationships between MLO proteins revealed the presence of apple, strawberry, peach and apricot homologs in the clade V, containing all dicot MLO homologs shown so far to be involved in PM susceptibility, thus making them candidates to act as susceptibility factors. Although the simple clustering in clade V is not enough to recognize a gene as a susceptibility factor, it does provide the first evidence for functionality and allows for the reduction in the number of candidates for further functional analysis. Clade IV, that contains functional MLO susceptibility homologs from monocots, was found to include one homolog from F. vesca (FvMLO17) and one from P. persica (PpMLO12). In accordance with this finding, a MLO homolog from the dicot species V. vinifera also clusters in clade IV [23, 24]. Figure 1). Interestingly, phylogenetic analyses carried out in this study also revealed the presence of one or two additional clades, depending on the type of phylogenetic reconstruction (UPGMA or Neighbour-Joining), which were not reported to occur in earlier investigations. Moreover, they appear to be characteristic of Rosaceae, since they contain only homologs from this family. Clearly, the exclusivity for Rosaceae of these clade(s) needs to be confirmed by further studies containing a larger dataset of MLO proteins. Additional studies could be also addressed to the functional characterization of Rosaceae MLO homologs grouped in clade VII. Indeed, this appears to be basal to both clade IV and clade V (Figure 1), and thus might have contained ancestral proteins which later on evolved into PM susceptibility factors.
Synteny between apple, peach and woodland strawberry MLOgenes
We took advantage of recent developments in Rosaceae genomics in order to detect synteny between P. persica, F. vesca and M. domestica chromosomal regions containing MLO homologs. This permitted the inference of ortholgous relationships between MLO genes in these species. Notably, all predicted MLO orthologs from different Rosaceae species, fell in the same phylogenetic clade (Tables 1, 2 and 3; Figure 1 and Additional file 1). This is to be expected, since orthologs generally share the same function and thus are characterized by a high level of sequence conservation. It is noteworthy that the chromosomal localization of predicted MLO orthologs between P. persica, M. domestica and F. vesca is in accordance with the results of the synteny study performed after the release of the three genomes [33, 34]. In particular, genes situated on peach scaffold 2, 7 and 8 were predicted to have orthologs on strawberry chromosome 7, 1 and 2, respectively, whereas genes on peach scaffold 4 were predicted to have orthologs on strawberry chromosomes 2 or 3 (Table 4). FvMLO3 was predicted to be orthologous to two peach MLO genes, PpMLO8 and PpMLO18, which are localised in close proximity to each other on peach scaffold 6 and grouped together in clade VI. In this case, we hypothesize a relation of co-orthology due to the occurrence of a recent tandem duplication event in the peach genome. Similarly, PpMLO5 and FvMLO16 were predicted to be orthologs of two apple MLO genes, MdMLO10 and MdMLO17, located on chromosomes 3 and 11. This is consistent with indications of duplications of large segments of these two chromosomes during the evolution of the apple genome .
Transcription of apple putative MLO genes in response to P. leucotrichainoculation
In barley, pea and tomato, only one of the clade V MLO homologs seems to be involved in powdery mildew susceptibility, whereas in A. thaliana three MLO genes in Clade V are required to be inactivated in order to achieve a fully resistant phenotype [16, 27]. This implies that, within Clade V MLO genes, a further selection might be required to identify PM susceptibility genes. Accumulating evidence indicates that MLO susceptibility genes are up-regulated upon challenge with powdery mildew fungi . Therefore, we analysed the expression level of apple MLO genes identified in this study in response to the interaction with P. leucotricha. Three pathogen-dependent gene up-regulations were detected. Two up-regulated MLO homologs, MdMLO11 and MdMLO19, encode for proteins falling in clade V, thus making them likely candidates to act as PM susceptibility genes in apple. MdMLO11 and MdMLO19 are located on chromosomes 4 and 12 respectively, and are therefore both generated from a duplication event in the 9-chromosome ancestor of apple . A third pathogen-dependent up-regulated gene, MdMLO18, was found, which encodes a protein grouping in the newly identified Clade VII (Figure 1). The presence of a PM upregulated gene outside clade V is consistent with transcriptome analyses recently performed in tomato (Appiano et al., unpublished results). Apple clade V also contains two genes, MdMLO5 and MdMLO7, which show no significant changes in expression following inoculation. Accordingly, the lack of up-regulation of some clade V MLO genes has been observed in grapevine and tomato [23, 24]; Appiano et al., unpublished results). The possible role of these genes as susceptibility factors has not yet been highlighted.
PpMLO3, PaMLO3 and FvMLO4 are likely to represent true orthologs of MdMLO19 (Table 4). Since orthologs often maintain the same function during evolution, we conjecture that the expression of these genes might also be responsive to PM fungi attacking corresponding species. Moreover, FvMLO15 and PpMLO9 are likely orthologs of MdMLO18, so they should also be considered as putative transcriptionally responsive genes to PM fungi attack. Further studies aimed at the functional characterization of these genes (e.g. through the application of reverse genetic approaches of targeted mutagenesis or gene silencing), in apple but also in peach and strawberry, might lead to the identification of resistant phenotypes, which could be used for the development of PM resistant cultivars. Particularly, studies on MdMLO18 could lead to the characterization of a possible role for clade VII in the interaction with PM fungi.
Our work led to the identification of 19 MLO homologs in peach, 17 in strawberry and 21 in apple. Three, three and four homologs, respectively, belong to clade V and therefore are candidates for being S-genes. Due to the high similarity between peach and apricot, we were able to amplify and characterize three Clade V apricot MLO genes.
The phylogenetic analysis revealed two new Rosaceae specific clades for the MLO family, although this needs to be confirmed by the use of a larger MLO proteins dataset.
Through inoculation of apple with P. leucotrica, we identified three up-regulated genes, i.e. MdMLO11, MdMLO18 and MdMLO19. MdMLO11 and MdMLO19, that belong to Clade V, are positioned in duplicated regions and have high sequence identity, therefore they are likely to be recent paralogs. MdMLO18 belongs to the newly identified Clade VII.
In silicoidentification and comparison of MLO predicted proteins in peach, woodland strawberry and apple
Predicted peptides from the peach genome (v. 1.0) and the strawberry genome (v.1.0) gene model databases, available at the website of the Genomic Database for Rosaceae (http://www.rosaceae.org) , were queried for the presence of MLO homologs protein sequences. First, a BLAST search, using the tomato SlMLO1 amino acid sequence as query was carried out. A further search was performed with the HMMER programme, which uses a method for homolog searches based on the profile hidden Markov probabilistic model . The sequences obtained with the previously mentioned BLAST search, were used together with other known MLO sequences from dicot and monocot species, namely: four RhMLOs from Rosa hybrida, 15 AtMLOs from Arabidopsis thaliana, SlMLO1 from Solanum lycopersicum, CaMLO2 from Capsicum annuum, PsMLO1 from Pisum sativum, MtMLO1 from Medicago truncatula, LjMLO1 from Lotus japonicus, VvMLO14 from V. Vinifera, HvMLO from Hordeum vulgare, TaMLO1_A1b and TaMLO_B1 from Triticum aestivum and OsMLO2 from Oryza sativa. MLO protein sequences from apple (Malus domestica cv. Golden Delicious) were identified by searching for the MLO domain profile (IPR004326) in the apple genome available at FEM-IASMA computational biology web resources (http://genomics.research.iasma.it). The resulting list was integrated with a BLAST search, carried out with the amino acid sequences previously listed for the HMMER search in peach and strawberry.
Chromosomal localization and predicted introns/exons structure of each MLO gene of apple, peach and strawberry was deduced based on the available genomic information at the GDR database. The presence and number of membrane spanning helices was predicted using the online software InterPro (http://www.ebi.ac.uk/interpro). Alignments for conserved amino-acids analysis were performed with the CLC Sequence Viewer v. 6.9 software (http://clcbio.com).
A total of 90 MLO protein sequences, including three apricot MLO sequences isolated in vitro (see next paragraph), were used for Clustal alignment (http://www.ebi.ac.uk/Tools/msa/clustalw2/). UPGMA-based and Neighbour-Joining-based phylogenetic trees were obtained with the CLC sequence viewer software. The UPGMA clustering algorithm was further used as input for the Dendroscope software, suitable for the visualization of large phylogenetic trees .
Relationships of orthology between MLO candidate genes from peach, woodland strawberry and apple were inferred by running the GBrowse-Syn tool available at GDR (http://www.rosaceae.org/gb/gbrowse_syn/peach_apple_strawberry) [35, 38]. This displays syntenic regions among the three available genomes of Rosaceae, as detected by the Mercator programme [35, 39]. For 50 Kb chromosomal stretches flanking each P. persica PpMLO homolog, syntenic regions from F. vesca and M. domestica were searched. Orthology was called upon the identification of F. vesca or M. domestica MLO homologs within syntenic blocks.
In vitro isolation of apricot MLOhomologs
RNA from apricot leaves (cultivar Orange Red) was extracted by using the SV Total RNA Isolation System Kit (Promega), and corresponding cDNA was synthesized by using the QuantiTect Reverse Transcription Kit (Qiagen) with oligo(dT) primers. Sequences of the peach MLO homologs PpMLO1, PpMLO3 and PpMLO4, are phylogenetically close to MLO homologs functionally associated to PM susceptibility, and were therefore used to design the primer pairs 5 ′ -ATGGCAGCCGCAACCTCAGGAAGA-3 ′ /5 ′ -TTATATACTTTGCCTATTGTCAAAC-3 ′ , 5 ′ -ATGGCAGGGGGAAAAGAAGGACG-3 ′ /5 ′ -TCAACTCCTTTCTGATTTCTCAA-3 ′ and 5 ′ -ATGGCCGAACTAAGTAAAGA-3 ′ /5 ′ TCAACTTCTTGATTTTCCTTTGC-3 ′ , respectively. These were employed to amplify full-length cDNA sequences of apricot putative orthologs, by using the AccuPrime Taq polymerase (Invitrogen). Amplicons were purified by using the NucleoSpin Extract II kit (Macherey-Nagel) and ligated (molar ratio 1:1) into the pGEM-T Easy vector (Promega). Recombinant plasmids were cloned in E. coli DH10Î² chemically competent cells and recovered by using the Qiaprep spin miniprep kit (Qiagen). Sequencing reactions were performed twice, by using universal T7 and SP6 primers (Eurofins MWG Operon).
Glasshouse test with apple cultivars
A total of 192 apple plants from three cultivars (Braeburn, Golden Delicious and Gala) were used to measure transcript abundance of MLO genes. Budwoods from these cultivars were grafted on M9 rootstocks in January 2012. The grafts were kept at −1°C for 2 months, and potted at the beginning of March in greenhouse. The plants grew for 6 weeks in the greenhouse at 20°C during the day, 17°C during the night, relative humidity of 70% and natural day/night cycle.
P. leucothrica was collected from apple trees in an unsprayed test orchard and used to infect greenhouse grown apple seedlings from ‘Gala Galaxy’ seeds. Four weeks after inoculation, conidia were used for the inoculation experiment, or transferred to new seedlings, to keep them viable. We inoculated by touching the plants with heavily infected apple seedlings. Control plants were not inoculated and kept separated in the same greenhouse of the inoculated plants. Inoculated and control plants were grown in the greenhouse at the growing conditions previously mentioned. The leaf samples were collected 4, 6, 8 and 24 hours post-inoculation (hpi).
Eight experimental repeats were performed and each sample contained three or four young leaves collected from each single plant. Every plant was used for sampling only once, to avoid any possible effect of wounding on the expression of MLO genes. The smallest statistical unit was a plant. The leaves were flash-frozen and ground in liquid nitrogen, and stored at −80°C until RNA extraction.
qPCR analysis of transcript levels
RNA extraction was carried out with the MagMAX-96 Total RNA isolation kit (Applied Biosystem) that includes DNAse treatment. The kit yielded between 50 and 200 ng/ul, of good quality RNA per sample.
Primers for gene expression analysis were designed with NCBI Primer Designing Tool (http://www.ncbi.nlm.nih.gov/tools/primer-blast/). Four serial dilutions of cDNA (1/5 - 1/25 – 1/125 – 1/625) were used to calculate the efficiency of each primer pair with iCycler software (Biorad). In case of efficiency lower than 1.80 or greater than 2.20, the primer pair was discarded and a new one tested, with the exception of MdMLO9, for which was not possible to design a primer pair with better efficiency. It was only possible to analyse 19 MLO genes because for MdMLO12 and MdMLO16 was not possible to design specific and efficient primer pairs, despite numerous attempts. Presence of a specific final dissociation curve was determined after each qPCR amplification reaction with progressive increment of temperature from 65°C to 95°C (0.5°C each step, 5 sec) and the size of the product was confirmed by agarose gel electrophoresis.
Quantitative Real Time-PCR (qPCR) was performed with SYBR greenER mix (Invitrogen) in a 15-Î¼L reaction volume, using a Bio-Rad iCycler iQ detection system, run by the Bio-Rad iCycler iQ multicolor 3.1 software. The software applies comparative quantification with an adaptive baseline. Samples were run in two technical replicates with the following thermal cycling parameters: 95°C 3 min – 95°C 15 sec, 60°C 1 min (repeated 40 times) – 95°C 10 sec.
Reference genes Î²-actin (NCBI accession number DT002474; Plaza accession number MD00G171330 - http://bioinformatics.psb.ugent.be/plaza/), ubiquitin (Plaza accession number MD05G001920) and elongation factor 1 (Plaza accession number MD09G014760) were used as reference genes (Table 5). All these three genes were used in previous works [40–42]. For additional control, we assessed the stability of our genes with the software geNorm (medgen.ugent.be/~jvdesomp/genorm/). An M-value lower than 1.5 is generally considered as stable enough [43–45] and all three reference genes in all three cultivars considered are within this threshold. We saw differences in stability between cultivars: ‘Golden Delicious’ was the most stable cultivar (actin: 0.824 – ubiquitin: 0,852 – elongation factor 1: 0,926), whereas ‘Braeburn’ was the less stable (actin: 1.246 – ubiquitin: 1,293 – elongation factor 1: 1,369) and ‘Gala’ showed intermediate stability (actin: 1.039 – ubiquitin: 1,152 – elongation factor 1: 1,078).
Each of the biological replicates was analysed in duplicate and the average of these two replicates was used for further analysis. In case of excessive difference between the two replicates (one Ct or more), the run was repeated. Considering the high number of samples and genes of interest, we opted for this approach in order to reduce the number of total runs. Data analysis was performed according to Hellemans et al., using the statistical package SPSS (IBM). This analysis method takes into account the efficiency value of each primer pair. For some genes it was necessary to apply a natural log transformation to the data, in order to obtain normal distribution of residues. To investigate the differences between control and inoculated samples, we used T-test (p ≤ 0.05).
Glawe DA: The powdery mildews: a review of the world most familiar (yet poorly known) plant pathogens. Annu Rev Phytopathol. 2008, 46: 27-51.
Turechek WW, Carroll JE, Rosenberger DA: Powdery Mildew of Apple. 2004, New York State Integrated Pest Management Program: Cornell Universit, (http://www.nysipm.cornell.edu/factsheets/treefruit/diseases/pm/apple_pm.pdf)
Foulongne M, Pascal T, Pfeiffer F, Kervella J: QTLs for powdery mildew resistance in peach × Prunus davidiana crosses: consistency across generations and environments. Mol Breed. 2003, 12: 33-50.
Boesewinkel HJ: Differences between the conidial states of Podosphaera tridactyla and Sphaerotheca pannosa. Annales de Phytopathologie. 1979, 11: 525-527.
Xiao CL, Chandler CK, Price JF, Duval JR, Mertely JC, Legard DE: Comparison of epidemics of botrytis fruit rot and powdery mildew of strawberry in large plastic tunnel and field production systems. Plant Dis. 2001, 85: 901-909.
Fu XL, Lu Y, Liu XD, Li JQ: Crossability barriers in the interspecific hybridization betweenOryza sativaandO. Meyeriana. J Integr Plant Biol. 2009, 51: 21-28.
Parlevliet JE: What is durable resistance, a general outline. Durability of Disease Resistance. Edited by: Jacobs TH, Parlevliet JE. 1993, Dordrecht: Kluwer, 23-29.
Parisi L, Lespinasse Y, Guillaumes SJ, Kruger J: A new race of Venturia inequalis virulent to apples with resistance due to the Vf gene. Phytopathology. 1993, 85: 533-537.
Krieghoff O: Entwicklung einer In-vitro-selektionsmethode auf Resistenz von Malus-Genotypen gegenüber Podosphaera leucotricha (Ell. Et Ev.) Salm. und In-vitron-Differenzierung von Virulenzunterschieden des Erregers. 1995, Berlin: PhD dissertation, Humboldt-Universität
Pavan S, Jacobsen E, Visser RGF, Bai Y: Loss of susceptibility as a novel breeding strategy for durable and broad-spectrum resistance. Mol Breed. 2010, 25: 1-12.
Büschges R, Hollricher K, Panstruga R, Simons G, Wolter M, Frijters A, van Daelen R, van der Lee T, Diergaarde P, Groenendijk J, Töpsch S, Vos P, Salamini F, Schulze-Lefert P: The barley Mlo gene: a novel control element of plant pathogen resistance. Cell. 1997, 88 (5): 695-705.
Devoto A, Piffanelli P, Nilsson I, Wallin E, Panstruga R, Von Heijne G, Schulze-Lefert P: Topology, subcellular localization, and sequence diversity of the Mlo family in plants. J Biol Chem. 1999, 274: 34993-35004.
JÃ¸rgensen JH: Discovery, characterization and exploitation of Mlo powdery mildew resistance in barley. Euphytica. 1992, 63: 141-152.
Consonni C, Humphry ME, Hartmann HA, Livaja M, Durner J, Westphal L, Vogel J, Lipka V, Kemmerling B, Schulze-Lefert P, Somerville SC, Panstruga R: Conserved requirement for a plant host cell protein in powdery mildew pathogenesis. Nat Genet. 2006, 38 (6): 716-720.
Humphry M, Reinstädler A, Ivanov S, Bisseling T, Panstruga R: Durable broad-spectrum powdery mildew resistance in pea er1 plants is conferred by natural loss-of-function mutations in PsMLO1. Mol Plant Pathol. 2011, 12 (9): 866-878.
Pavan S, Schiavulli A, Appiano M, Marcotrigiano AR, Cillo F, Visser RGF, Bai Y, Lotti C, Ricciardi L: Pea powdery mildew er1 resistance is associated to loss-of-function mutations at a MLO homologous locus. Theor Appl Genet. 2011, 123: 1425-1431.
Bai Y, Pavan S, Zheng Z, Zappel N, Reinstädler A, Lotti C, De Giovanni C, Ricciardi L, Lindhout P, Visser RGF, Theres K, Panstruga R: Naturally occurring broad-spectrum powdery mildew resistance in a Central American tomato accession is caused by loss of MLO function. Mol Plant Microbe Interact. 2008, 21: 30-39.
Aist JR, Bushnell WR: Invasion of plants by powdery mildew fungi, and cellular mechanisms of resistance. The Fungal Spore and Disease Initiation in Plants and Animals. Edited by: Cole GT, Hoch HC. 1991, New YorK: plenum press, 321-345.
Miklis M, Consonni C, Bhat RA, Lipka V, Schulze-Lefert P, Panstruga R: Barley MLO modulates actin-dependent and actin-independent antifungal defense pathways at the cell periphery. Plant Physiol. 2007, 144: 1132-1143.
Panstruga R: Serpentine plant MLO proteins as entry portals for powdery mildew fungi. Biochem Soc Transact. 2005, 33 (Pt 2): 389-392.
Piffanelli P, Zhou FS, Casais C, Orme J, Jarosch B, Schaffrath U, Collins NC, Panstruga R, Schulze-Lefert P: The barley MLO modulator of defense and cell death is responsive to biotic and abiotic stress stimuli. Plant Physiol. 2002, 129: 1076-1085.
Zheng Z, Nonomura T, BÃ³ka K, Matsuda Y, Visser RGF, Toyoda H, Kiss L, Bai Y: Detection and quantification of Leveillula taurica growth in pepper leaves. Phytopathology. 2013, 103: 6-
Feechan A, Jermakow AM, Torregrosa L, Panstruga R, Dry IB: Identification of grapevine MLO gene candidates involved in susceptibility to powdery mildew. Funct Plant Biol. 2008, 35: 1255-1266.
Winterhagen P, Howard SF, Qiu W, Kovács LG: Transcriptional up-regulation of grapevine MLO genesin response to powdery mildew infection. Am J Enol Vit. 2008, 59: 2-
Kessler SA, Shimosato-Asano H, Keinath NF, Wuest SE, Ingram G, Panstruga R, Grossniklaus U: Conserved molecular components for pollen tube reception and fungal invasion. Science. 2010, 330: 968-
Chen Z, Noir S, Kwaaitaal M, Hartmann A, Wu MJ, Mudgil Y, Sukumar P, Muday G, Panstruga R, Jones AM: Two seven-transmembrane domain MILDEW RESISTANCE LOCUS O proteins cofunction in Arabidopsis root thigmomorphogenesis. Plant Cell. 2009, 21: 1972-1991.
Reinstädler A, Müller J, Czembor JH, Piffanelli P, Panstruga R: Novel induced mlo mutant alleles in combination with site-directed mutagenesis reveal functionally important domains in the heptahelical barley Mlo protein. BMC Plant Biol. 2010, 10: 31-
Devoto A, Hartmann HA, Piffanelli P, Elliott C, Simmons C, Taramino G, Goh CS, Cohen FE, Emerson BC, Schulze-Lefert P, Panstruga R: Molecular phylogeny and evolution of the plant-specific seven-transmembrane MLO family. J Mol Evol. 2003, 56 (1): 77-88.
Elliott C, Muller J, Miklis M, Bhat RA, Schulze-Lefert P, Panstruga R: Conserved extracellular cysteine residues and cytoplasmic loop-loop interplay are required for functionality of the heptahelical MLO protein. Biochem l J. 2005, 385 (Pt 1): 243-254.
Kaufmann H, Qiu X, Wehmeyer J, Debener T: Isolation, molecular characterization, and mapping of four rose MLO orthologs. Front Plant Sci. 2012, 3: 244-
Elliott C, Zhou F, Spielmeyer W, Panstruga R, Schulze-Lefert P: Functional conservation of wheat and rice Mlo orthologs in defense modulation to the powdery mildew fungus. Mol Plant Microbe Interact. 2002, 15 (10): 1069-1077.
Velasco R, Zharkikh A, Affourtit J, Dhingra A, Cestaro A, Kalyanaraman A, Fontana P, Bhatnagar SK, Troggio M, Pruss D, Salvi S, Pindo M, Baldi P, Castelletti S, Cavaiuolo M, Coppola G, Costa F, Cova V, Dal Ri A, Goremykin V, Komjanc M, Longhi S, Magnago P, Malacarne G, Malnoy M, Micheletti D, Moretto M, Perazzolli M, Si-Ammour A, Vezzulli S, et al: The genome of the domesticated apple (Malus × domestica Borkh). Nat Genet. 2010, 42: 833-839.
Shulaev V, Sargent DJ, Crowhurst RN, Mockler TC, Folkerts O, Delcher AL, Jaiswal P, Mockaitis K, Liston A, Mane SP, Burns P, Davis TM, Slovin JP, Bassil N, Hellens RP, Evans C, Harkins T, Kodira C, Desany B, Crasta OR, Jensen RV, Allan AC, Michael TP, Setubal JC, Celton J, Rees DJG, Williams KP, Holt SH, Ruiz-Rojas JJ, Chatterjee M, et al: The genome of woodland strawberry (Fragaria vesca). Nat Genet. 2011, 43: 109-116.
Jung S, Cestaro A, Troggio M, Main D, Zheng P, Cho I, Folta K, Sosinski B, Abbot B, Celton JM, ArÃ¹s P, Shulaev V, Verde I, Morgante M, Rokhsar D, Velasco R, Sargent DJ: Whole genome comparisons ofFragaria, Prunus and Malus reveal different modes of evolution between Rosaceous subfamilies. BMC Genomics. 2012, 13: 129-
Jung S, Ficklin SP, Lee T, Cheng C, Blenda A, Zheng P, Yu J, Bombarely A, Cho I, Ru S, Evans KM, Peace C, Abbott A, Mueller L, Olmstead M, Main D: The genome database for rosaceae (GDR): year 10 update. Nucleic Acid Res. 2014, 42: D1237-D1244.
Finn RD, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011, 39: W29-W37.
Huson D, Richter D, Rausch C, Dezulian T, Franz M, Rupp R: Dendroscope: an interactive viewer for large phylogenetic trees. BMC Bioinformatics. 2007, 8: 460-
McKay SJ, Vergara IA, Stajich JE: Using the generic synteny browser (GBrowse_syn). Curr Protoc Bioinformatics. 2010, Chapter 9 (Unit 9): 12-
Dewey CN: Aligning multiple whole genomes with Mercator and MAVID. Methods Mol Biol. 2007, 395: 221-236.
Kürkcüoglu S, Degenhardt J, Lensing J, Al-Masri AN, Gaul AE: Identification of differentially expressed genes in Malus domestica after application of the non-pathogenic bacterium Pseudomonas fluorescensBk3 to the phyllosphere. J Exp Bot. 2007, 58 (3): 733-741.
Giorno F, Guerriero G, Baric S, Mariani C: Heat shock transcriptional factors in Malus domestica: identification, classification and expression analysis. BMC Genomics. 2012, 13: 639-
Dal Cin V, Danesin M, Boschetti A, Dorigoni A, Ramina A: Ethylene biosynthesis and perception in apple fruitlet abscission (Malus domestica L. Borck). J Exp Bot. 2005, 56 (421): 2995-3005.
Ling D, Salvaterra PM: Robust RT-qPCR data normalization: validation and selection of internal reference genes during post-experimental data analysis. PLoS One. 2011, 6: 3-
Van Hiel MB, Van Wielendaele P, Temmerman L, Van Soest S, Vuerinckx K, Huybrechts R, Broeck JV, Simonet G: Identification and validation of housekeeping genes in brains of the desert locust Schistocerca gregaria under different developmental conditions. BMC Mol Biol. 2009, 10: 56-
Strube C, Buschbaum S, Wolken S, Schnieder T: Evaluation of reference genes for quantitative real-time PCR to investigate protein disulfide isomerase transcription pattern in the bovine lungworm Dictyocaulus viviparus. Gene. 2008, 425: 36-43.
Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J: qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007, 8: R19-
The authors would like to thank Remmelt Groenwold (WUR) for helping with apple inoculation with P. leucotricha and Lorenza Dalla Costa (FEM) for all the valuable advices about qPCR data analysis.
This work has been founded by Fondazione Edmund Mach and GMPF PhD programme.
The authors declare that they have no competing interests.
SPe carried out the inoculation and gene expression analysis, identified MLO homologs in apple and wrote the major part of the manuscript. SPa contributed to design the experiment, identified MLO homologs in peach, apricot and strawberry, performed the analysis of synteny and contributed to write the manuscript. DC contributed to design the experiment, performed the phylogenetic analysis and revised the manuscript. AG contributed to design the experiment and revised the manuscript. YB contributed to the design of the experiments. RGFV contributed to the design of the experiments and revised the manuscript. MM contributed to the design of the experiments, contributed to the identification of MLO homologs in apple and revised the manuscript. HS contributed to the design of the experiments and revised the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Results of search for F. vesca and M. domestica regions syntenic to 50 kb P. persica chromosomal stretches containing the PpMLO homologs identified in this study. Shaded polygons indicate aligned regions between genomes. Grid lines are drawn to indicate insertions/deletions between the genomes of F. vesca and M. domestica with respect to the P. persica reference sequence. P. persica, F. vesca and M. domestica MLO homologs, named according to the nomenclature of the Genomic Database of Rosaceae, are boxed. (PDF 817 KB)
Additional file 2: P. leucotricha. Transcription abundances of 19 MLO-like genes following powdery mildew (PM) inoculation in ‘Golden Delicious’ (1a), ‘Gala’ (1b) and ‘Braeburn (1c) leaf samples. The graphs show expression values of inoculated samples relative to control samples, averaged from four to eight biological replicate, normalized, that are in turn the average of two experimental replicates. The Ct values have been normalized with three reference genes: actin, ubiquitin and elongation factor 1. Statistical significance was determined with a t-test for each individual pair of inoculated and control samples at each time point (4, 6, 8 and 24 hpi). The error bars show standard errors of the means. Significant differences between inoculated samples and control samples are indicated with a *(P < 0.05). (PDF 1 MB)
About this article
Cite this article
Pessina, S., Pavan, S., Catalano, D. et al. Characterization of the MLO gene family in Rosaceae and gene expression analysis in Malus domestica. BMC Genomics 15, 618 (2014). https://doi.org/10.1186/1471-2164-15-618
- Powdery mildew
- Malus domestica