Evolution of the MLO gene families in octoploid strawberry (Fragaria ×ananassa) and progenitor diploid species identified potential genes for strawberry powdery mildew resistance

Abstract

Powdery mildew (PM) caused by Podosphaera aphanis is a major fungal disease of cultivated strawberry. Mildew Resistance Locus O (MLO) is a gene family described for having conserved seven-transmembrane domains. Induced loss-of-function in specific MLO genes can confer durable and broad resistance against PM pathogens. However, the genomic structure and potential role of MLO genes for PM resistance have not been characterized yet in the octoploid cultivated strawberry. In the present study, MLO gene families were characterized in four diploid progenitor species (Fragaria vesca, F. iinumae, F. viridis, and F. nipponica) and octoploid cultivated (Fragaria ×ananassa) strawberry, and potential sources of MLO-mediated susceptibility were identified. Twenty MLO sequences were identified in F. vesca and 68 identified in F. ×ananassa. Phylogenetic analysis divided diploid and octoploid strawberry MLO genes into eight different clades, in which three FveMLO (MLO10, MLO17, and MLO20) and their twelve orthologs of FaMLO were grouped together with functionally characterized MLO genes conferring PM susceptibility. Copy number variations revealed differences in MLO composition among homoeologous chromosomes, supporting the distinct origin of each subgenome during the evolution of octoploid strawberry. Dissecting genomic sequence and structural variations in candidate FaMLO genes revealed their potential role associated with genetic controls and functionality in strawberry against PM pathogen. Furthermore, the gene expression profiling and RNAi silencing of putative FaMLO genes in response to the pathogen indicate the function in PM resistance. These results are a critical first step in understanding the function of strawberry MLO genes and will facilitate further genetic studies of PM resistance in cultivated strawberry.

References

[1] Devoto, A. et al. Topology, subcellular, localization, and sequence diversity of the Mlo family in plants. J. Biol. Chem. 274, 13 (1999).

[2] Kim, M. C. et al. Mlo, a modulator of plant defense and cell death, is a novel calmodulinbinding protein. Isolation and characterization of a rice Mlo homologue. J. Biol. Chem. 277, 19304-19314 (2002).

[3] Ablazov, A. & Tombuloglu, H. Genome-wide identification of the mildew resistance locus O (MLO) gene family in novel cereal model species Brachypodium distachyon. Eur. J. Plant Pathol. 145, 239-253 (2015).

[4] Liu, Q. & Zhu, H. Molecular evolution of the MLO gene family in Oryza sativa and their functional divergence. Gene 409, 1-10 (2008).

[5] Pessina, S. et al. Characterization of the MLO gene family in Rosaceae and gene expression analysis in Malus domestica. BMC Genomics 15, 618 (2014).

[6] Kusch, S., Pesch, L. & Panstruga, R. Comprehensive phylogenetic analysis sheds light on the diversity and origin of the MLO family of integral membrane proteins. Genome Biol. Evol. 8, 878-895 (2016).

[7] Acevedo-Garcia, J., Kusch, S. & Panstruga, R. Magical mystery tour: MLO proteins in plant immunity and beyond. New Phytol. 204, 273-281 (2014).

[8] Feechan, A., Jermakow, A., Torregrosa, L., Panstruga, R. & Dry, I. Identification of grapevine MLO gene candidates involved in susceptibility to powdery mildew. Funct. Plant Biol. 35, 1255-1266 (2008).

[9] Deshmukh, R., Singh, V. K. & Singh, B. D. Comparative phylogenetic analysis of genome-wide Mlo gene family members from Glycine max and Arabidopsis thaliana. Mol. Genet Genomics 289, 345-359 (2014).

[10] Kusch, S. & Panstruga, R. mlo-Based resistance: an apparently universal “Weapon” to defeat powdery mildew disease. Mol. Plant Microbe Interact. 30, 179-189 (2017).

[11] Pastruga, R. Serpentine plant MLO proteins as entry portals for powdery mildew fungi. Biochem. Soc. Trans. 33, 389-392 (2005).

[12] Jørgensen, I. H. Discovery, characterization and exploitation of Mlo powdery mildew resistance in barley. Euphytica 63, 141-152 (1992).

[13] Consonni, C. et al. Conserved requirement for a plant host cell protein in powdery mildew pathogenesis. Nat. Genet 38, 716-720 (2006).

[14] Zheng, Z. et al. Loss of function in Mlo orthologs reduces susceptibility of pepper and tomato to powdery mildew disease caused by Leveillula taurica. PLoS ONE 8, e70723 (2013).

[15] Qiu, X. et al. Antisense RhMLO1 genetransformation enhances resistance to the powdery mildew pathogen in Rosa multiflora. Plant Mol. Biol. Report. 33, 1659-1665 (2015).

[16] Pessina, S. et al. Knockdown of MLO genes reduces susceptibility to powdery mildew in grapevine. Hortic. Res. 3, 16016 (2016).

[17] Hong, C., Weiping, K., Junfen, L. & Jiping, L. Analysis of powdery mildew resistance in wild melon MLO mutants. Horticultural Plant J. 1, 165-171 (2015).

[18] Nie, J. et al. Loss-of-function mutations in CsMLO1 confer durable powdery mildew resistance in cucumber (Cucumis sativus L.). Front. Plant Sci. 6, 1155- 1155 (2015).

[19] Pavan, S. et al. Pea powdery mildew er1 resistance is associated to loss-of-function mutations at a MLO homologous locus. Theor. Appl. Genet. 123, 1425-1431 (2011).

[20] Appiano, M. et al. Identification of candidate MLO powdery mildew susceptibility genes in cultivated Solanaceae and functional characterization of tobacco NtMLO1. Transgenic Res. 24, 847-858 (2015).

[21] Bai, Y. et al. Naturally occurring broad-spectrum powdery mildew resistance in a Central American tomato accession is caused by loss of mlo function. Mol. Plant Microbe Interact. 21, 30-39 (2008).

[22] 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. 15, 1069-1077 (2002).

[23] Devoto, A. et al. Molecular phylogeny and evolution of the plant-specific seventransmembrane MLO family. J. Mol. Evol. 56, 77-88 (2003).

[24] Acevedo-Garcia, J. et al. mlo-based powdery mildew resistance in hexaploid bread wheat generated by a non-transgenic TILLING approach. Plant Biotechnol. J. 15, 367-378 (2017).

[25] Nekrasov, V. et al. Rapid generation of a transgene-free powdery mildew resistant tomato by genome deletion. Sci. Rep. 7, 482 (2017).

[26] Peries, O. S. Studies on strawberry mildew, caused by Sphaerotheca macularis (Wallr. ex Fries) Jaczewski*. Ann. Appl. Biol. 50, 211-224 (1962).

[27] Xiao, C. L. et al. Comparison of epidemics of botrytis fruit rot and powdery mildew of strawberry in large plastic tunnel and field production systems. Plant Dis. 85, 901-909 (2001).

[28] Paulus, A. Fungal disease of strawberry. HortScience 25, 885-889 (1990).

[29] Peres, N. A. & Mertely. J. C. Powdery Mildew of Strawberries. Electronic Data Information Source UF/IFAS Extension (2018).

[30] Nelson, M. D., Gubler, W. D. & Shaw, D. V. Relative resistance of 47 strawberry cultivars to powdery mildew in california greenhouse and field environments. Plant Dis. 80, 326-328 (1996).

[31] Njuguna, W., Liston, A., Cronn, R., Ashman, T.-L. & Bassil, N. Insights into phylogeny, sex function and age of Fragaria based on whole chloroplast genome sequencing. Mol. Phylogenet. Evol. 66, 17-29 (2013).

[32] Shulaev, V. et al. The genome of woodland strawberry (Fragaria vesca). Nat. Genet 43, 109-116 (2011).

[33] Anciro, A. et al. FaRCg1: a quantitative trait locus conferring resistance to Colletotrichum crown rot caused by Colletotrichum gloeosporioides in octoploid strawberry. Theor. Appl Genet 131, 2167-2177 (2018).

[34] Salinas, N., Verma, S., Peres, N. & Whitaker, V. M. FaRCa1: a major subgenome-specific locus conferring resistance to Colletotrichum acutatum in strawberry. Theor. Appl. Genet. 132, 1109-1120 (2019).

[35] Roach, J. A. et al. FaRXf1: a locus conferring resistance to angular leaf spot caused by Xanthomonas fragariae in octoploid strawberry. Theor. Appl. Genet. 129, 1191-1201 (2016).

[36] Hirakawa, H. et al. Dissection of the octoploid strawberry genome by deep sequencing of the genomes of fragaria species. DNA Res. 21, 169-181 (2013).

[37] Edger, P. P. et al. Reply to: revisiting the origin of octoploid strawberry. Nat. Genet. 52, 5-7 (2020).

[38] Feng, C. et al. Tracing the diploid ancestry of the cultivated octoploid strawberry. Mol. Biol. Evol. 38, 478-485 (2020).

[39] Edger, P. P. et al. Origin and evolution of the octoploid strawberry genome. Nat. Genet 51, 541-547 (2019).

[40] Edger, P. P. et al. Single-molecule sequencing and optical mapping yields an improved genome of woodland strawberry (Fragaria vesca) with chromosome-scale contiguity. Gigascience 7, 1-7 (2018).

[41] Miao, L. X. et al. Genomic identification, phylogeny, and expression analysis of MLO genes involved in susceptibility to powdery mildew in Fragaria vesca. Genet. Mol. Res. 15, gmr.15038400 (2016).

[42] Shitsukawa, N. et al. Genetic and epigenetic alteration among three homeologous genes of a class E MADS box gene in hexaploid wheat. Plant Cell 19, 1723 (2007).

[43] Mutti, J. S., Bhullar, R. K. & Gill, K. S. Evolution of gene expression balance among homeologs of natural polyploids. G3 Genes Genomes Genet. 7, 1225-1237 (2017).

[44] Dobson, L., Reményi, I. & Tusnády, G. E. CCTOP: a consensus constrained topology prediction web server. Nucleic acids Res. 43, W408-W412 (2015).

[45] Horton, P. et al. WoLF PSORT: protein localization predictor. Nucleic Acids Res. 35, W585-W587 (2007).

[46] Panstruga, R. Discovery of novel conserved peptide domains by ortholog comparison within plant multi-protein families. Plant Mol. Biol. 59, 485-500 (2005).

[47] Appiano, M. et al. Monocot and dicot MLO powdery mildew susceptibility factors are functionally conserved in spite of the evolution of class-specific molecular features. BMC Plant Biol. 15, 257 (2015).

[48] Jiao, Y. & Paterson, A. H. Polyploidy-associated genome modifications during land plant evolution. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 369, 20130355 (2014).

[49] Hardigan, M. A. et al. Genome Synteny has been conserved among the octoploid progenitors of cultivated strawberry over millions of years of evolution. Front. Plant Sci. 10, 1789 (2020).

[50] Jiwan, D., Roalson, E. H., Main, D. & Dhingra, A. Antisense expression of peach mildew resistance locus O (PpMlo1) gene confers cross-species resistance to powdery mildew in Fragaria xananassa. Transgenic Res. 22, 1119-1131 (2013).

[51] Tao, Y., Zhao, X., Mace, E., Henry, R. & Jordan, D. Exploring and exploiting pangenomics for crop improvement. Mol. Plant 12, 156-169 (2019).

[52] Chen, Z. et al. Expression analysis of the AtMLO gene family encoding plant-specific seven-transmembrane domain proteins. Plant Mol. Biol. 60, 583-597 (2006).

[53] Leister, D. Tandem and segmental gene duplication and recombination in the evolution of plant disease resistance gene. Trends Genet. 20, 116-122 (2004).

[54] Panchy, N., Lehti-Shiu, M. & Shiu, S. H. Evolution of gene duplication in plants. Plant Physiol. 171, 2294-2316 (2016).

[55] Zheng, Z. et al. Genome-wide study of the tomato SlMLO Gene family and its functional characterization in response to the powdery mildew fungus Oidium neolycopersici. Front Plant Sci. 7, 380 (2016).

[56] Andolfo, G. et al. Evolutionary conservation of MLO gene promoter signatures. BMC Plant Biol. 19, 150 (2019).

[57] Voorrips, R. E. MapChart: software for the graphical presentation of linkage maps and QTLs. J. Hered. 93, 77-78 (2002).

[58] Guo, A.-Y. et al. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics 31, 1296-1297 (2014).

[59] Miller, M. A., Pfeiffer, W. & Schwartz, T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. in 2010 Gateway Computing Environments Workshop (GCE) 1-8 (2010).

[60] Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792-1797 (2004).

[61] Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312-1313 (2014).

[62] Kearse, M. et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647-1649 (2012).

[63] Suyama, M., Torrents, D. & Bork, P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 34, W609-W612 (2006).

[64] Yang, Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586-1591 (2007).

[65] Sarkar, D. Lattice: Multivariate Data Visualization with R. (Springer, 2008).

[66] Brors, B., Gu, L., Schlesner, M., Eils, R. & Gu, Z. circlize implements and enhances circular visualization in R. Bioinformatics 30, 2811-2812 (2014).

[67] Marchler-Bauer, A. et al. CDD: a Conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 39, D225-D229 (2011).

[68] Liu, W. et al. IBS: an illustrator for the presentation and visualization of biological sequences. Bioinformatics 31, 3359-3361 (2015).

[69] Mitchell, C. MultAlin-multiple sequence alignment. Bioinformatics 9, 614- 614 (1993).

[70] Sánchez-Sevilla, J. F. et al. Gene expression atlas of fruit ripening and transcriptome assembly from RNA-seq data in octoploid strawberry (Fragaria ×ananassa). Sci. Rep. 7, 13737 (2017).

[71] Porebski, S., Bailey, L. G. & Baum, B. R. Modification of a CTAB DNA extraction protocol for plants containing high polysaccharide and polyphenol components. Plant Mol. Biol. Report. 15, 8-15 (1997).

[72] Koren, S. et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722-736 (2017).

[73] Calis, O., Cekic, C., Soylu, S. & Tör, M. Identification of new resistance sources from diploid wild strawberry against powdery mildew pathogen. Pak. J. Agric. Sci. 52, 677-683 (2015).

[74] Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using realtime quatitative PCR and the 2^−ΔΔCT method . Methods 25, 402-408 (2001).

[75] Ahmed, F. et al. A web server for designing effective and specific plant siRNAs with genome-wide off-target assessment. Plant Physiol. 184, 65 (2020). pssRNAit.

[76] Karimi, M., Inzé, D. & Depicker, A. GATEWAY™ vectors for Agrobacterium-mediated plant transformation. Trends Plant Sci. 7, 193-195 (2002).

[77] Wei, W. et al. Ectopic expression of FvWRKY42, a WRKY transcription factor from the diploid woodland strawberry (Fragaria vesca), enhances resistance to powdery mildew, improves osmotic stress resistance, and increases abscisic acid sensitivity in Arabidopsis. Plant Sci. 275, 60-74 (2018).

[78] Kennedy, C. et al. Evaluation of strawberry species and cultivars for powdery mildew resistance in open-field and high tunnel production systems. HortScience 48, 5 (2013).

[79] Chaube, H. S. & Singh, U. S. Plant Disease Management: Principles and Practices (CRC Press, Boca Raton, 1991).