Genome-wide association study of salicylic acid provides genetic insights for tea plant selective breeding

Abstract

Salicylic acid (SA) is a phenolic phytohormone widely believed to regulate plant growth and stress response. Despite its significance, the genetic basis of SA-mediated resistance to biotic stressors in tea plants is little understood. Our study investigated the genetic diversity, population structure, and linkage disequilibrium (LD) patterns of 299 tea accessions using 79 560 high-quality single nucleotide polymorphisms (SNPs) obtained from genotyping-by-sequencing (GBS) data. Our genome-wide association study identified CSS0033791.1, an essential gene encoding 9-cis-epoxycarotenoid dioxygenase (CsNCED1), which catalyzes a vital step in abscisic acid (ABA) biosynthesis. Exogenous ABA treatment and transgenic overexpression of the CsNCED1 gene lowered SA content in the respective tea plants by inhibiting the expression of the ICS gene. Further analysis revealed that ABA could reduce the expression levels of the SA receptor gene (NPR1) and NPR1 target genes (PR1 and WRKY18), increasing the plant’s susceptibility to biotic stressors. Furthermore, the feeding behavior of Spodoptera litura revealed that the insect bite area on transgenic leaves was substantially more extensive than that in wild type (WT), implying that the CsNCED1 gene had a negative regulatory role in SA-mediated immune response. This study thus provides the foundation for future insect resistance breeding, sustainable tea plant resource usage, and molecular marker-assisted (MAS) tea plant breeding.

Acknowledgements

This study was supported by the program of NSFC [32060700], the National Guidance Foundation for Local Science and Technology Development of China [[2023] 009] and [2022-1-52], Guiyang Science and Technology Plan Project [Construction Technology Contract [2023] 48-21], the Science and Technology Project of Guizhou Province, China (Qiankehe Foundation-ZK [2022]) (Grant No.702076222101).

Author contributions

N.S. performed the experimental design; D.X. conducted the experiments and analyzed experiment data; D.X., W.X., and N.S. wrote and revised the MS; D.X., A.H., H.D., and B.D. collected samples and analyzed data. All authors of this MS have approved the final MS version.

Data availability

The raw sequence data reported in this study are available at the Genome Sequence Archive in BIG Data Center, Beijing Institute of Genomics (BIG) (http://bigd.big.ac.cn/gsa) under accession number CRA001438.

Conflict of interest

The authors claim no conflict of interest.

Supplementary data

Supplementary data is available at Horticulture Research online.

References

[1] Xu Y, Dong Y, Cheng W. et al. Characterization and phyloge-netic analysis of the complete mitochondrial genome sequence of Diospyros oleifera, the first representative from the family Ebenaceae. Heliyon. 2022; 8:09870

[2] Wang J, Zhang T, Ren L. et al. Establishment of a direct somatic embryogenesis regeneration system using immature cotyledon explants in Camellia sinensis cv. Shuchazao. IndCropProd. 2024; 210:118076

[3] Chen L, Yang YJ, Yu FL. Tea germplasm research in China: recent progresses and prospects. J Plant Genet Resour. 2004; 5:389-92

[4] Hazarika L, Bhuyan M, Hazarika B. Insect pests of tea and their management. Annu Rev Entomol. 2009; 54:267-84

[5] Song J, Yan J, Sun B. et al. Abscisic acid regulates Cl(-) efflux via the ABI5-ZAT10-SLAH 3 module in chloride-stressed Malus hupehensis. Hortic Res. 2024; 11:uhae200

[6] Feng S, Li N, Chen H. et al. Large-scale analysis of the ARF and Aux/IAA gene families in 406 horticultural and other plants. Mol Hortic. 2024; 4:13

[7] Hong J, Meng K, Rae TH. et al. Reframing agriculture by light: the role of light-mediated jasmonates/salicylic acid reg-ulation in plant defense, development and beyond. Veg Res. 2024; 4:e027

[8] Zhang Z, Chen C, Jiang C. et al. VvWRKY5 positively regulates wounding-induced anthocyanin accumulation in grape by inter-playing with VvMYBA1 and promoting jasmonic acid biosynthe-sis. Hortic Res. 2024; 11:uhae083

[9] Gah-Hyun L. Regulation of salicylic acid and N-hydroxy-pipecolic acid in systemic acquired resistance. Plant Pathol J. 2023; 39:21-7

[10] Duan D, Yang J, Yi R. et al. MdVQ17 negatively regulates apple resistance to Glomerella leaf spot by promoting MdWRKY17-mediated salicylic acid degradation and pectin lyase activity. Hortic Res. 2024; 11:uhae159

[11] Jiang CJ, Shimono M, Sugano S. et al. Abscisic acid inter-acts antagonistically with salicylic acid signaling pathway in rice-Magnaporthe grisea interaction. Mol Plant-Microbe Interact. 2010; 23:791-8

[12] Zeng W, Melotto M, He SY. Plant stomata: a checkpoint of host immunity and pathogen virulence. Curr Opin Biotechnol. 2010; 21: 599-603

[13] Cao FY, Yoshioka K, Desveaux D. The roles of ABA in plant-pathogen interactions. JPlant Res. 2011; 124:489-99

[14] Du M, Zhai Q, Deng L. et al. Closely related NAC transcription factors of tomato differentially regulate stomatal closure and reopening during pathogen attack. Plant Cell. 2014; 26:3167-84

[15] Ding Y, Dommel M, Mou Z. Abscisic acid promotes proteasome-mediated degradation of the transcription coactivator NPR1 in Arabidopsis thaliana. Plant J. 2016; 86:20-34

[16] Baloch FS, Altaf MT, Liaqat W. et al. Recent advancements in the breeding of sorghum crop: current status and future strategies for marker-assisted breeding. Front Genet. 2023; 14:1150616

[17] Wang Y, Niu S, Deng X. et al. Genome-wide association study, population structure, and genetic diversity of the tea plant in Guizhou plateau. BMC Plant Biol. 2024; 24:79

[18] Huang F, Duan J, Lei Y. et al. Genetic diversity, population struc-ture and core collection analysis of Hunan tea plant germplasm through genotyping-by-sequencing. Beverage Plant Res. 2022; 2: 1-7

[19] Liu Z, Zhao Y, Yang P. et al. Population whole-genome rese-quencing reveals the phylogenetic relationships and population structure of four Hunan typical tea landraces. Beverage Plant Res. 2023; 3:1-8

[20] Wang X, Liu Z, Zhang F. et al. Integrative genomics reveals the polygenic basis of seedlessness in grapevine. Curr Biol. 2024; 34:3763-3777.e5

[21] Gai W, Yang F, Yuan L. et al. Multiple-model GWAS identifies optimal allelic combinations of quantitative trait loci for malic acid in tomato. Hortic Res. 2023; 10:uhad021

[22] Zhang Z, Peng C, Xu W. et al. Genome-wide association study of agronomic traits related to nitrogen use efficiency in Henan wheat. BMC Genomics. 2024; 25:7

[23] Li D, Zhang F, Pinson SRM. et al. Assessment of rice sheath blight resistance including associations with plant architecture, as revealed by genome-wide association studies. Rice (N Y). 2022; 15:31

[24] Fang K, Xia Z, Li H. et al. Genome-wide association analysis identified molecular markers associated with important tea flavor-related metabolites. Hortic Res. 2021; 8:42

[25] Rui L, Su J, Li T. et al. Molecular regulation of immunity in tea plants. Mol Biol Rep. 2022; 50:2883-92

[26] Zhang Y, Zhang J, Yan C. et al. Metabolome and microbiome signatures in the leaves of wild tea plant resources resistant to Pestalotiopsis theae. Front Microbiol. 2022; 13:907962

[27] Evans BA, Rozen DE. Significant variation in transformation frequency in Streptococcus pneumoniae. ISME J. 2013; 7:791-9

[28] Wang H, Tang X, Yang X. et al. Exploiting natural variation in crown root traits via genome-wide association studies in maize. BMC Plant Biol. 2021; 21:346-6

[29] Wu B, Ren W, Zhao L. et al. Genome-wide association study of root system architecture in maize. Genes. 2022; 13:181-1

[30] Laura TC, Zhang Z, Yu J. Status and prospects of genome-wide association studies in plants. Plant Genome. 2021; 14:e20077

[31] Vandana J, Tirthankar B, Kumar SR. et al. Multi-environment GWAS identifies genomic regions underlying grain nutrient traits in foxtail millet (Setaria italica). Plant Cell Rep. 2023; 43:6-6

[32] Guang Y, Pan Y, Pan W. et al. Combined GWAS and eGWAS reveals the genetic basis underlying drought tolerance in emmer wheat (Triticum turgidum L.). New Phytol. 2024; 242: 2115-31

[33] Xu X, Wang Z, Xu S. et al. Identifying loci controlling total starch content of leaf in Nicotiana tabacum through genome-wide association study. Funct Integr Genomics. 2022; 22:537-52

[34] Hifzur R, Prashant V, Hu Y. et al. Mining genomic regions associ-ated with agronomic and biochemical traits in quinoa through GWAS. Sci Rep. 2024; 14:9205-5

[35] Maruki T, Lynch M. Genome-wide estimation of linkage dise-quilibrium from population-level high-throughput sequencing data. Genetics. 2014; 197:1303-13

[36] Zhu X, Dong L, Jiang L. et al. Constructing a linkage-linkage disequilibrium map using dominant-segregating markers. DNA Res. 2016; 23:1-10

[37] Hirayama T, Shinozaki K. Perception and transduction of abscisic acid signals: keys to the function of the versatile plant hormone ABA. Trends Plant Sci. 2007; 12:343-51

[38] Yuan H, Yang J, Ningkun X. et al. OsNCED5, a 9-cis-epoxycarotenoid dioxygenase gene, regulates salt and water stress tolerance and leaf senescence in rice. Plant Sci. 2019; 287:110188

[39] Nidhi S, Barkha R, Saini LK. et al. Voltage-dependent anion channel 3 (VDAC3) mediates P. syringae induced ABA-SA sig-naling crosstalk in Arabidopsis thaliana. Plant Physiol Biochem. 2024; 206:108237

[40] Jiang Y, Yue Y, Wang Z. et al. A novel ABA structural analogues enhanced plant resistance by inducing the plant immunity and inactivating ABA signaling pathway. Adv Agrochem. 2024; 3:64-73

[41] Xiang N, Zhang B, Hu J. et al. Modulation of carotenoid biosyn-thesis in maize (Zea mays L.) seedlings by exogenous abscisic acid and salicylic acid under low temperature. Plant Cell Rep. 2023; 43:1

[42] Du Y, Zhang H, Jia K. et al. Role of abscisic acid-mediated stomatal closure in responses to pathogens in plants. Physiol Plant. 2024; 176:e14135

[43] Liang X, Liao G, Li J. et al. Exogenous ABA promotes resistance to Sitobion avenae (Fabricius) in rice seedlings. Pest Manag Sci. 2024; 80:3389-400

[44] Allasia V, Industri B, Ponchet M. et al. Quantification of salicylic acid (SA) and SA-glucosides in Arabidopsis thaliana. Bio Protoc. 2018; 8:2844

[45] Niu S, Song Q, Hisashi K. et al. Genetic diversity, linkage dise-quilibrium, and population structure analysis of the tea plant (Camellia sinensis) from an origin center, Guizhou plateau, using genome-wide SNPs developed by genotyping-by-sequencing. BMC Plant Biol. 2019; 19:328

[46] Xia E, Tong W, Hou Y. et al. The reference genome of tea plant and resequencing of 81 diverse accessions provide insights into its genome evolution and adaptation. Mol Plant. 2020; 13:1013-26

[47] McKenna A, Hanna M, Banks E. et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20:1297-303

[48] Slifer SH. PLINK: key functions for data analysis. Curr Protoc Hum Genet. 2018; 97:59

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

[50] Mussmann SM, Douglas MR, Chafin TK. et al. ADMIXPIPE: pop-ulation analyses in ADMIXTURE for non-model organisms. BMC Bioinformatics. 2020; 21:337

[51] Kumar S, Stecher G, Li M. et al. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018; 35:1547-9

[52] Zhang C, Dong SS, Xu JY. et al. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics. 2019; 35:1786-8

[53] Huang L, Min Y, Schiessl S. et al. Integrative analysis of GWAS and transcriptome to reveal novel loci regulation flowering time in semi-winter rapeseed. Plant Sci. 2021; 310:110980

[54] Deng J, Jiao Q, Wang Y. et al. Transcriptomic and metabolomic effects of exogenous ABA application on tobacco seedling growth. Plant Growth Regul. 2023; 101:399-414

[55] Li H, Guo L, Yan M. et al. A rapid and efficient transient expression system for gene function and subcellular localiza-tion studies in the tea plant (Camellia sinensis) leaves. Sci Hortic. 2022; 297:110927

[56] Qin L, Li X, Ren S. et al. AcCNGC2, coding a cyclic nucleotide-gated channel protein from onion, enhances Phytophthora nicotianae resistance in transgenic Nicotiana benthamiana. Sci Hortic. 2024; 326:112779

[57] Weigel D, Glazebrook J. Transformation of Agrobacterium using the freeze-thaw method. CSH Protoc. 2006;2006:pdb.prot4666

[58] Qin L, Mo N, Zhang Y. et al. CaRDR1, an RNA-dependent RNA polymerase plays a positive role in pepper resistance against TMV. Front Plant Sci. 2017; 8:1068

[59] Cao Y, Zhao X, Tang M. et al. GbKCS14 synthesizes VLCFAs to respond to low temperatures during fuzz initial stage of Gossypium barbadense. IndCropProd. 2024; 217:118875

[60] Zhang Q, Wang L, Wang Z. et al. The regulation of cell wall ligni-fication and lignin biosynthesis during pigmentation of winter jujube. Hortic Res. 2021; 8:238-8