Integrative GWAS and transcriptomic analyses reveal markers and candidate genes associated with resistance to Botrytis cinerea fruit rot in blueberry

Abstract

Botrytis cinerea is a broad host range fungal pathogen causing gray mold disease and crop losses worldwide. In blueberries, symptoms include blossom blight in the field and postharvest fruit rot, affecting the entire supply chain. With control options constrained by regulatory restrictions and fungicide resistance, the dissection of the genetic and molecular basis of blueberry response to B. cinerea can accelerate breeding for resistance. In this study, we phenotyped 354 blueberry selections using a high-throughput Botrytis infection fruit assay. The same population was genotyped by targeted sequencing for genome-wide association study (GWAS). In addition, we performed RNA-seq time-course (0–96 hours post-inoculation) for resistant and susceptible genotypes. Our results showed a continuum of tolerance levels and moderate narrow-sense heritability estimates for the disease-related traits (0.46–0.61). GWAS identified small-effect loci, consistent with quantitative resistance observed in other host plant species. Intersecting differentially expressed genes with GWAS intervals revealed eight candidate genes. Transcriptomic analyses showed that, at early stages, the resistant genotype upregulated components of basal innate immunity, including wax and cutin biosynthesis, responses to wounding and fungal-derived molecules, the MAPK cascade, and ethylene and jasmonate signaling. In contrast, susceptible genotypes displayed delayed activation of these defense pathways and altered cell wall-related processes. The moderate correlation between disease traits and wax bloom further supported a role for wax in disease response. Together, our findings provide molecular markers and candidate genes for Botrytis fruit rot resistance in blueberry with significant applications in breeding programs and opportunities to future validation studies.

References

[1] Williamson B, Tudzynski B, Tudzynski P, et al. Botrytis cinerea: the cause of grey mould disease . Mol Plant Pathol. 2007; 8: 561-80

[2] Elad Y, Pertot I, Cotes Prado AM, et al. Plant hosts of Botrytis spp. In: Fillinger S, Elad Y, eds. Botrytis-The Fungus, the Pathogen and Its Management in Agricultural Systems. Springer: Cham, 2016, 413-86

[3] Weiberg A, Wang M, Lin F-M, et al. Fungal small RNAs suppress plant immunity by hijacking host RNA interference pathways. Science. 2013; 342: 118-23

[4] Barrau C, Santos B, Romero F . Susceptibility of southern highbush and rabbiteye blueberry cultivars to postharvest diseases in Huelva, Spain. Acta Hortic. 2006; 715: 525-30

[5] Smith BJ . Botrytis blossom blight of southern blueberries: cultivar susceptibility and effect of chemical treatments. Plant Dis. 1998; 82: 924-7

[6] Vasquez P, Baldomá JA, Wright ER, et al. First report of blueberry Botrytis blight in Buenos Aires, Entre Ríos, and Córdoba, Argentina. Plant Dis. 2007; 91: 639

[7] Bristow PR, Milholland RD, Schilder AC. Botrytis blight and fruit rot. In: Polashock JJ, Caruso FL, Averill AL, et al., eds. Compendium of Blueberry, Cranberry, and Lingonberry Diseases and Pests.2nd ed. APS Press: St. Paul, MN, 2017, 17-9

[8] Harmon PF . Botrytis Blossom Blight of Southern Highbush Blueberry . Gainesville, FL: University of Florida IFAS Extension; 2018

[9] Rivera SA, Zoffoli JP, Latorre BA . Infection risk and critical period for the postharvest control of gray mold (Botrytis cinerea) on blueberry in Chile . Plant Dis. 2013; 97: 1069-74

[10] Bestfleisch M, Luderer-Pflimpfl M, Höfer M, et al. Evaluation of strawberry (Fragaria L.) genetic resources for resistance to Botrytis cinerea . Plant Pathol. 2015; 64: 396-405

[11] Droby S, Lichter A . Post-harvest Botrytis infection: etiology, development and management. In: Elad Y, Williamson B, Tudzynski P, Delen N, eds. Botrytis: Biology, Pathology and Control.Springer: Dordrecht, 2007, 349-67

[12] Gabler FM, Fassel R, Mercier J, et al. Influence of temperature, inoculation interval, and dosage on biofumigation with Muscodor albus to control postharvest gray mold on grapes. Plant Dis. 2006; 90: 1019-25

[13] Kwon JH, Cheon MG, Choi O, et al. First report of Botrytis cinerea as a postharvest pathogen of blueberry in Korea. Mycobiology. 2011; 39: 52-3

[14] Romanazzi G, Lichter A, Gabler FM, et al. Recent advances on the use of natural and safe alternatives to conventional methods to control postharvest gray mold of table grapes. Postharvest Biol Technol. 2012; 63: 141-7

[15] Sautua FJ, Baron C, Pérez-Hernández O, et al. First report of resistance to carbendazim and procymidone in Botrytis cinerea from strawberry, blueberry and tomato in Argentina. Crop Prot. 2019; 125: 104879-20

[16] Cantu D, Vicente AR, Greve LC, et al. The intersection between cell wall disassembly, ripening, and fruit susceptibility to Botrytis cinerea . Proc Natl Acad Sci USA. 2008; 105: 859-64

[17] Tani H, Koshino H, Sakuno E, et al. Botcinins A, B, C, and D, metabolites produced by Botrytis cinerea, and their antifungal activity against Magnaporthe grisea, a pathogen of rice blast disease . J Nat Prod. 2005; 68: 1768-72

[18] Bessire M, Chassot C, Jacquat AC, et al. A permeable cuticle in Arabidopsis leads to a strong resistance to Botrytis cinerea . EMBO J. 2007; 26: 2158-68

[19] Elad Y, Evensen K . Physiological aspects of resistance to Botrytis cinerea . Phytopathology. 1995; 85: 637-43

[20] Sun J, Zhang X, Zheng J, et al. Importance of cell wall permeability and cell wall degrading enzymes during infection of Botrytis cinerea in hazelnut. Forests. 2023; 14: 565

[21] Bi K, Liang Y, Mengiste T, et al. Killing softly: a roadmap of Botrytis cinerea pathogenicity. Trends Plant Sci. 2023; 28: 211-22

[22] Thomma BPHJ, Eggermont K, Penninckx IAMA, et al. Separate jasmonate-dependent and salicylate-dependent defense-response pathways in Arabidopsis are essential for resistance to distinct microbial pathogens. Proc Natl Acad Sci USA. 1998; 95: 15107-11

[23] Thomma BPHJ, Eggermont K, Tierens KF-J, et al. Requirement of functional ethylene-insensitive 2 gene for efficient resistance of Arabidopsis to infection by Botrytis cinerea . Plant Physiol. 1999; 121: 1093-101

[24] Boudsocq M, Sheen J . CDPKs in immune and stress signaling. Trends Plant Sci. 2013; 18: 30-40

[25] Prusky D, Lichter A . Activation of quiescent infections by postharvest pathogens during transition from the biotrophic to the necrotrophic stage. FEMS Microbiol Lett. 2007; 268: 1-8

[26] Su K, Zhao W, Lin H, et al. Candidate gene discovery of Botrytis cinerea resistance in grapevine based on QTL mapping and RNA-seq. Front Plant Sci. 2023; 14: 1127206

[27] Goyal N, Bhatia G, Garewal N, et al. Genome-wide identification, characterization and expression analysis of the effector-triggered immunity (ETI) pathway-mediated defense-responsive genes in grapes against powdery and downy mildew infections. Res Sq. 2020; v2: 1-34

[28] Denby KJ, Kumar P, Kliebenstein DJ . Identification of Botrytis cinerea susceptibility loci in Arabidopsis thaliana . Plant J. 2004; 38: 473-86

[29] Fordyce RF, Soltis NE, Caseys C, et al. Digital imaging combined with genome-wide association mapping links loci to plant-pathogen interaction traits. Plant Physiol. 2018; 178: 1406-22

[30] Corwin JA, Copeland D, Feusier J, et al. The quantitative basis of the Arabidopsis innate immune system to endemic pathogens depends on pathogen genetics. PLoS Genet. 2016; 12: e1005789

[31] Rowe HC, Kliebenstein DJ . Complex genetics control natural variation in Arabidopsis thaliana resistance to Botrytis cinerea . Genetics. 2008; 180: 2237-50

[32] Finkers R, Bai Y, Van Den Berg P, et al. Quantitative resistance to Botrytis cinerea from Solanum neorickii . Euphytica. 2008; 159: 83-92

[33] Davis J, Yu D, Evans W, et al. Mapping of loci from Solanum lycopersicoides conferring resistance or susceptibility to Botrytis cinerea in tomato . Theor Appl Genet. 2009; 119: 305-14

[34] Petrasch S, Mesquida-Pesci SD, Pincot DDA, et al. Genomic prediction of strawberry resistance to postharvest fruit decay caused by the fungal pathogen Botrytis cinerea . G3 (Bethesda). 2022; 12: jkab378

[35] Zhang W, Kwon ST, Chen F, et al. Isolate dependency of Brassica rapa resistance QTLs to Botrytis cinerea . Front Plant Sci. 2016; 7: 1-13

[36] Anuradha C, Gaur PM, Pande S, et al. Mapping QTL for resistance to Botrytis grey mould in chickpea. Euphytica. 2011; 182: 1-9

[37] Fu Y, van Silfhout A, Shahin A, et al. Genetic mapping and QTL analysis of Botrytis resistance in Gerbera hybrida. Mol Breed. 2017; 37: 13

[38] Maggiolini FAM, Prencipe A, Bergamini C, et al. A comparative transcriptomic study reveals temporal and genotype-specific defense responses to Botrytis cinerea in grapevine. J Fungi (Basel). 2025; 11: 124

[39] Lee K, Lee JG, Min K, et al. Transcriptome analysis of the fruit of two strawberry cultivars “Sunnyberry” and “Kingsberry” that show different susceptibility to Botrytis cinerea after harvest. Int J Mol Sci. 2021; 22: 1518

[40] Sniezko RA, Nelson CD . Resistance breeding against tree pathogens. In: Asiegbu FO, Kovalchuk A, eds. Forest Microbiology, Vol. 2: Forest Tree Health. Cambridge, MA, USA: Academic Press, 159-75

[41] Poland JA, Balint-Kurti PJ, Wisser RJ, et al. Shades of gray: the world of quantitative disease resistance. Trends Plant Sci. 2009; 14: 21-9

[42] St. Clair DA . Quantitative disease resistance and quantitative resistance loci in breeding. Annu Rev Phytopathol. 2010; 48: 247-68

[43] Gajewska KA, Lescesen H, Ramialison M, et al. Nuclear transporter Importin-13 plays a key role in the oxidative stress transcriptional response. Nat Commun. 2021; 12: 5904

[44] Audenaert K, De Meyer GB, Höfte MM . Abscisic acid determines basal susceptibility of tomato to Botrytis cinerea and suppresses salicylic acid-dependent signaling mechanisms. Plant Physiol. 2002; 128: 491-501

[45] Asselbergh B, Curvers K, França SC, et al. Resistance to Botrytis cinerea in sitiens, an abscisic acid-deficient tomato mutant, involves timely production of hydrogen peroxide and cell wall modifications in the epidermis. Plant Physiol. 2007; 144: 1863-77

[46] Liu S, Kracher B, Ziegler J, et al. Negative regulation of ABA signaling by WRKY33 is critical for Arabidopsis immunity towards Botrytis cinerea 2100 . Elife. 2015; 4: e07295

[47] Song X, Zhang B, Zhou Y . Golgi-localized UDP-glucose transporter is required for cell wall integrity in rice. Plant Signal Behav. 2011; 6: 1097-100

[48] Nohales MÁ, Flores R, Daròs JA . Viroid RNA redirects host DNA ligase 1 to act as an RNA ligase. Proc Natl Acad Sci USA. 2012; 109: 13805-10

[49] Goggin FL, Shah J, Gillaspy G . Editorial: lipid metabolism and membrane structure in plant biotic interactions. Front Plant Sci. 2022; 13: 1096268

[50] Lu J, Xu Y, Wang J, et al. The role of triacylglycerol in plant stress response. Plants (Basel). 2020; 9: 472

[51] Seth T, Asija S, Umar S, et al. The intricate role of lipids in orchestrating plant defense responses. Plant Sci. 2024; 338: 111904

[52] Liu R, Zhang L, Xiao S, et al. Ursolic acid, the main component of blueberry cuticular wax, inhibits Botrytis cinerea growth by damaging cell membrane integrity. Food Chem. 2023; 415: 135753

[53] Gabler FM, Smilanick JL, Mansour M, et al. Correlations of morphological, anatomical, and chemical features of grape berries with resistance to Botrytis cinerea . Phytopathology. 2003; 93: 1263-73

[54] Jiang B, Liu R, Fang X, et al. Botrytis cinerea infection affects wax composition, content and gene expression in blueberry fruit. Postharvest Biol Technol. 2022; 192: 112020

[55] Mittler R, Zandalinas SI, Fichman Y, et al. Reactive oxygen species signalling in plant stress responses. Nat Rev Mol Cell Biol. 2022; 23: 663-79

[56] Martin RE, Postiglione AE, Muday GK . Reactive oxygen species function as signaling molecules in controlling plant development and hormonal responses. Curr Opin Plant Biol. 2022; 69: 102293

[57] Ali M, Kaderbek T, Khan MA, et al. Biosynthesis and multifaceted roles of reactive species in plant defense mechanisms during environmental cues. Plant Stress. 2025; 18: 101102

[58] Govrin EM, Levine A . The hypersensitive response facilitates plant infection by the necrotrophic pathogen Botrytis cinerea . Curr Biol. 2000; 10: 751-7

[59] Blanco-Ulate B, Vincenti E, Powell ALT, et al. Tomato transcriptome and mutant analyses suggest a role for plant stress hormones in the interaction between fruit and Botrytis cinerea . Front Plant Sci. 2013; 4: 142

[60] Jin Y, Zhang Y, Lin L, et al. Cucumber PGIP2 is involved in resistance to gray mold disease. Gene. 2024; 923: 148588

[61] Yuan S, Gou X, Hu J, et al. Mutation of tomato xyloglucan transglucosylase/hydrolase5 increases fruit firmness and contributes to prolonged shelf life. J Plant Physiol. 2024; 303: 154350

[62] Meuwissen THE, Hayes BJ, Goddard ME . Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001; 157: 1819-29

[63] Ferrão LFV, Amadeu RR, Benevenuto J, et al. Genomic selection in an outcrossing autotetraploid fruit crop: lessons from blueberry breeding. Front Plant Sci. 2021; 12: 676326

[64] Benito EP, Ten Have A, van’t Klooster JW, et al. Fungal and plant gene expression during synchronized infection of tomato leaves by Botrytis cinerea . Eur J Plant Pathol. 1998; 104: 207-20

[65] Krasnow C, Ziv C . Non-chemical approaches to control postharvest gray mold disease in bell peppers. Agronomy. 2022; 12: 216

[66] Salotti I, Ji T, Rossi V . Temperature requirements of Colletotrichum spp. belonging to different clades. Front Plant Sci. 2022; 13: 953760

[67] Casorzo G, Ferrão LF, Adunola P, et al. Understanding the genetic basis of blueberry postharvest traits to define better breeding strategies. G3. 2024; 14: jkae163

[68] Colle M, Leisner CP, Wai CM, et al. Haplotype-phased genome and evolution of phytonutrient pathways of tetraploid blueberry. Gigascience. 2019; 8: giz012

[69] Lee W-P, Stromberg MP, Ward A, et al. MOSAIK: a hash-based algorithm for accurate next-generation sequencing short-read mapping. PLoS One. 2014; 9: e90581

[70] Garrison E, Marth G . Haplotype-based variant detection from short-read sequencing. arXiv:1207.3907, 2012.

[71] Danecek P, Auton A, Abecasis G, et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8

[72] Gerard D, Ferrão LFV, Garcia AAF, et al. Genotyping polyploids from messy sequencing data. Genetics. 2018; 210: 789-807

[73] Rosyara UR, De Jong WS, Douches DS, et al. Software for genome-wide association studies in autopolyploids and its application to potato. Plant Genome. 2016; 9: 1-10

[74] Ghimire L, Adunola P, Harmon PF, et al. Genome-wide association study reveals candidate loci for resistance to anthracnose in blueberry. G3 (Bethesda). 2025; 15: jkaf187

[75] Amadeu RR, Cellon C, Olmstead JW, et al. AGHmatrix: R package to construct relationship matrices for autotetraploid and diploid species: a blueberry example. Plant Genome. 2016; 9: 1-10

[76] Pérez P, de los Campos G . BGLR: a statistical package for whole genome regression and prediction. Genetics. 2014; 198: 483-95

[77] Andrews S. FastQC: A Quality Control Tool for High-throughput Sequence Data . Cambridge, UK: Babraham Bioinformatics; 2010

[78] Kim D, Paggi JM, Park C, et al. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019; 37: 907-15

[79] Anders S, Pyl PT, Huber W . HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015; 31: 166-9

[80] Love MI, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15: 550

[81] Ferrão LFV, Benevenuto J, Oliveira IB, et al. Insights into the genetic basis of blueberry fruit-related traits using diploid and polyploid models in a GWAS context. Front Ecol Evol. 2018; 6: 107

[82] Törönen P, Holm L . PANNZER-a practical tool for protein function prediction. Protein Sci. 2022; 31: 118-28

[83] Huerta-Cepas J, Szklarczyk D, Heller D, et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 2019; 47: D309-14

[84] Yu G, Wang L-G, Han Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16: 284-7

[85] Supek F, Bošnjak M, Skunca N, et al. Revigo summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011; 6: e21800