Mining of CtOPR2 as an essential gene regulating jasmonic acid mediated insect tolerance in Camellia tachangensis through genome-wide association studies

Abstract

Jasmonic acid (JA), a pivotal lipid-derived phytohormone, serves as a critical regulator in plant growth and defense mechanisms. However, the genetic mechanisms of OPR2 gene in JA-dependent biotic defenses of Camellia tachangensis have rarely been investigated. In this study, we performed a genome-wide association study to analyze 100,720 high-quality single nucleotide polymorphisms (SNPs) among 350 tea accessions from Guizhou province to identify genetic variations associated with JA. Analysis showed C. tachangensis displayed higher levels of JA content, further analysis identified 60 high-quality SNPs and nine candidate genes related to JA. Among them, CtOPR2 encoding 12-oxophytodienoic acid reductase 2 is responsible for catalyzing the conversion of 4,5-didehydro-JA (4,5-ddh-JA) to JA. The expression level of CtOPR2 in three tea accessions with different JA content was consistent with the dynamic changes of JA content. The expression level of synthetic (AOS, AOC, and ACX) and responsive (WRKY18 and MYC2) genes were significantly decreased and increased in asODN- CtOPR2-treated shoot tips and transgenic tobacco lines overexpressing CtOPR2, respectively, which were consistent with the JA content. These results further revealed that CtOPR2 gene played essential roles in promoting JA biosynthesis. A significant reduction in insect bite area was observed on transgenic tobacco leaves compared to wild-type leaves in feeding experiments with Spodoptera litura, highlighting that the positive regulatory function of CtOPR2 gene in JA-mediated immune responses. This study provides a robust theoretical foundation for marker-assisted selection breeding in tea, aimed at developing high-JA germplasm with potentially enhanced pest resistance for cultivation in Guizhou.

References

[1] Zhang Z, Feng X, Wang Y, et al. Advances in research on functional genes of tea plant. Gene. 2019; 711: 143940

[2] Pan S, Nie Q, Tai H, et al. Tea and tea drinking: China’s outstanding contributions to the mankind. Chin Med. 2022; 17: 27

[3] Gao J, Chen D, Lin Z, et al. Research progress on the antidiabetic activities of tea and its bioactive components. Bev Plant Res. 2023; 3: 32

[4] Yang M, Zhou L, Kan Z, et al. Beneficial health effects and possible health concerns of tea consumption: a review. Bev Plant Res. 2025; 5: e035

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

[6] Zhang Y, Quan J, Dong Y, et al. Genome-wide identification, molecular evolution, and expression analysis of the bHLH gene family in Mesembryanthemum crystallinum. Veg Res. 2025; 5: e045

[7] Li J, Jian G, Qian J, et al. Emission pattern and anti-insect function of indole from tea plant (Camellia sinensis) attacked by tea geometrids . Bev Plant Res. 2024; 4: e003

[8] Zainab A, Chen M . Indirect plant defense against insect herbivores: a review. Insect Sci. 2018; 25: 2-23

[9] Wu J, Baldwin IT . New insights into plant responses to the attack from insect herbivores. Annu Rev Genet. 2010; 44: 1-24

[10] Major IT, Guo Q, Zhai J, et al. A Phytochrome B-independent pathway restricts growth at high levels of jasmonate defense. Plant Physiol. 2020; 183: 733-49

[11] Han Z, Zhang H, Li P, et al. Genome-wide characterization and identification of candidate CsRAV genes involved in JA pathway responses to cold stress in tea plant (Camellia sinensis) . Bev Plant Res. 2023; 3: 27

[12] Shi R, Wang Y, Zhou F, et al. Nitrogen fertilizer reduction based on bioorganic fertilizer improves the yield and quality of fresh leaves of alpine tea in summer. Bev Plant Res. 2024; 4: e035

[13] Turner JG, Ellis C, Devoto A . The jasmonate signal pathway. Plant Cell. 2002; 14: S153-64

[14] Chini A, Monte I, Zamarreño AM, et al. An OPR3-independent pathway uses 4,5-didehydrojasmonate for jasmonate synthesis. Nat Chem Biol. 2018; 14: 171-8

[15] Schaller F, Weiler EW . Molecular cloning and characterization of 12-oxophytodienoate reductase, an enzyme of the octadecanoid signaling pathway from Arabidopsis thaliana. Structural and functional relationship to yeast old yellow enzyme . J Biol Chem. 1997; 272: 28066-72

[16] Li W, Liu B, Yu L, et al. Phylogenetic analysis, structural evolution and functional divergence of the 12-oxo-phytodienoate acid reductase gene family in plants. BMC Evol Biol. 2009; 9: 90

[17] Nie W, Chen Y, Tao J, et al. Identification of the 12-oxophytoeienoic acid reductase (OPR) gene family in pepper (Capsicum annuum L.) and functional characterization of CaOPR6 in pepper fruit development and stress response . Genome. 2022; 65: 537-45

[18] Tani T, Sobajima H, Okada K, et al. Identification of the OsOPR7 gene encoding 12-oxophytodienoate reductase involved in the biosynthesis of jasmonic acid in rice. Planta. 2008; 227: 517-26

[19] Pigolev AV, Miroshnichenko DN, Pushin AS, et al. Overexpression of Arabidopsis OPR3 in Hexaploid wheat (Triticum aestivum L.) alters plant development and freezing tolerance . Int J Mol Sci. 2018; 19: 3989

[20] Stintzi A, Browse J . The Arabidopsis male-sterile mutant, opr3, lacks the 12-oxophytodienoic acid reductase required for jasmonate synthesis. Proc Natl Acad Sci USA. 2000; 97: 10625-30

[21] Xin Z, Zhang J, Ge L, et al. A putative 12-oxophytodienoate reductase gene CsOPR3 from Camellia sinensis, is involved in wound and herbivore infestation responses . Gene. 2017; 615: 18-24

[22] Verhoeven A, Finkers-Tomczak A, Prins P, et al. The root-knot nematode effector MiMSP32 targets host 12-oxophytodienoate reductase 2 to regulate plant susceptibility. New Phytol. 2023; 237: 2360-74

[23] Tibbs Cortes L, Zhang Z, Yu J . Status and prospects of genome-wide association studies in plants. Plant Genome. 2021; 14: e20077

[24] Xin W, Chen N, Wang J, et al. Candidate gene analysis of rice grain shape based on genome-wide association study. Theor Appl Genet. 2024; 137: 241

[25] Xu R, Yang Q, Liu Z, et al. Genome-wide association analysis and genomic prediction of salt tolerance trait in soybean germplasm. Front Plant Sci. 2024; 15: 1494551

[26] Zeng Y, Xu X, Jiang J, et al. Genome-wide association analysis and genomic selection for leaf-related traits of maize. PLoS One. 2025; 20: e0323140

[27] Liang S, Wang H, Yamashita H, et al. Genome-wide identification and expression analysis of sucrose phosphate synthase and sucrose-6-phosphate phosphatase family genes in Camellia sinensis . Bev Plant Res. 2024; 4: 0

[28] Deng X, Huang D, Wang Y, et al. Genome-wide association study of salicylic acid provides genetic insights for tea plant selective breeding. Hortic Res. 2025; 12: uhae362

[29] OuYang Y, Rao K, He L, et al. Identification of CsRNF144 as a key gene regulating di-hydroxy catechins accumulation in Camellia tachangensis through genome-wide association studies . Hortic Plant J. 2025; 2468-0141

[30] Niu S, Song Q, Koiwa H, et al. Genetic diversity, linkage disequilibrium, 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

[31] Kottawa-Arachchi JD, Gunasekare MTK, Ranatunga MAB . Biochemical diversity of global tea [ Camellia sinensis (L.) O. Kuntze] germplasm and its exploitation: a review . Genet Resour Crop Evol. 2019; 66: 259-73

[32] Pandey J, Scheuring DC, Koym JW, et al. Genetic diversity and population structure of advanced clones selected over forty years by a potato breeding program in the USA. Sci Rep. 2021; 11: 8344

[33] Alatawi I, Xiong H, Mou B, et al. Genomic insights of leafminer resistance in spinach through GWAS approach and genomic prediction. Hortic Plant J. 2024; 12: 2468-0141

[34] 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

[35] Bayer PE, Gill M, Danilevicz MF, et al. Producing high-quality single nucleotide polymorphism data for genome-wide association studies. Methods Mol Biol (Clifton, NJ). 2022; 2481: 153-9

[36] Borthakur D, Wang B, Meinhardt LW, et al. Multiple distinctive lineages of Camellia sinensis var. assamica from India and China revealed by Single Nucleotide Polymorphism markers . Bev Plant Res. 2025; 5: e026

[37] Huang J, Huang X, Li M . Wild tea species as a genetic resource for future breeding programs. Tree Genet Mol Breed. 2024; 14: 277-85

[38] Kong W, Kong X, Xia Z, et al. Genomic analysis of 1,325 Camellia accessions sheds light on agronomic and metabolic traits for tea plant improvement. Nat Genet. 2025; 57: 997-1007

[39] Wu Q, Song L, Bai D, et al. Genome-wide association studies reveal potential candidate genes associated with amino acid in tea plants. PeerJ. 2025; 13: e18969

[40] Zhang Z, Ersoz E, Lai C, et al. Mixed linear model approach adapted for genome-wide association studies. Nat Genet. 2010; 42: 355-60

[41] Liu X, Huang M, Fan B, et al. Iterative usage of fixed and random effect models for powerful and efficient genome-wide association studies. PLoS Genet. 2016; 12: e1005767

[42] Proietti S, Caarls L, Coolen S, et al. Genome-wide association study reveals novel players in defense hormone crosstalk in Arabidopsis. Plant Cell Environ. 2018; 41: 2342-56

[43] Zhang N, Zhou S, Yang D, et al. Revealing shared and distinct genes responding to JA and SA signaling in Arabidopsis by meta-analysis. Front Plant Sci. 2020; 11: 908

[44] Shui J, Liu X, Chen Q, et al. Differential expression analysis of jasmonic acid pathway genes between resistant and susceptible pepper cultivars against Myzus persicae infestation. Chin J Trop Crops. 2023; 4: 774-84

[45] Chen Y, Jin G, Liu M, et al. Multiomic analyses reveal key sectors of jasmonate-mediated defense responses in rice. Plant Cell. 2024; 36: 3362-77

[46] Wang Y, Ma C, Kong X, et al. Identification of a DUF538 gene conferring resistance to tea green leafhopper (Empoasca onukii) in Camellia sinensis . Hortic Res. 2025; 13: uhaf264

[47] Zhang W, Zhang Y, Qiu H, et al. Genome assembly of wild tea tree DASZ reveals pedigree and selection history of tea varieties. Nat Commun. 2020; 11: 3719

[48] Fu J, Ren R, Jin S, et al. Overexpression of a putative 12-oxophytodienoate reductase gene, EpOPR1, enhances acetyl-shikonin production in Echium plantagineum . In Vitro Cell Dev Biol Plant. 2022; 58: 311-20

[49] Xie Q, Dong W, Wang M, et al. BpWRKY6 regulates insect resistance by affecting jasmonic acid and terpenoid synthesis in Betula platyphylla. Plant Biotechnol J. 2025; 23: 3682-96

[50] Aliakbari M, Tahmasebi S, Sisakht JN . Jasmonic acid improves barley photosynthetic efficiency through a possible regulatory module, MYC2-RcaA, under combined drought and salinity stress. Photosynth Res. 2024; 159: 69-78

[51] Goodin MM, Zaitlin D, Naidu RA, et al. Nicotiana benthamiana: its history and future as a model for plant-pathogen interactions. Mol Plant Microbe Interact. 2008; 21: 1015-26

[52] Uyehara AN, Del Valle-Echevarria AR, Hunter CT, et al. Cytokinin promotes jasmonic acid accumulation in the control of maize leaf growth. Plants (Basel, Switz). 2023; 12: 3014

[53] Shan X, Zhang Y, Peng W, et al. Molecular mechanism for jasmonate-induction of anthocyanin accumulation in Arabidopsis. J Exp Bot. 2009; 60: 3849- 60

[54] Huang H, Liu B, Liu L, et al. Jasmonate action in plant growth and development. J Exp Bot. 2017; 68: 1349-59

[55] Stepanenko AA, Heng HH . Transient and stable vector transfection: pitfalls, off-target effects, artifacts. Mutat Res Rev Mutat Res. 2017; 773: 91-103

[56] Feng S, Liu Z, Chen H, et al. PHGD: an integrative and user-friendly database for plant hormone-related genes. Imeta. 2024; 3: e164

[57] Liu H, Li X, Xiao J, et al. A convenient method for simultaneous quantification of multiple phytohormones and metabolites: application in study of rice-bacterium interaction. Plant Methods. 2012; 8: 2

[58] Xia E, Li F, Tong W, et al. Tea plant information archive: a comprehensive genomics and bioinformatics platform for tea plant. Plant Biotechnol J. 2019; 17: 1938-53

[59] Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. 2013; arXiv preprint arXiv:1303.3997

[60] 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

[61] Chen W, Hou L, Zhang Z, et al. Genetic diversity, population structure, and linkage disequilibrium of a core collection of Ziziphus jujuba assessed with genome-wide SNPs developed by genotyping-by-sequencing and SSR markers . Front Plant Sci. 2017; 8: 575

[62] Hussain W, Baenziger PS, Belamkar V, et al. Genotyping-by-sequencing derived high-density linkage map and its application to QTL mapping of flag leaf traits in bread wheat. Sci Rep. 2017; 7: 16394

[63] Eltaher S, Sallam A, Belamkar V, et al. Genetic diversity and population structure of F3:6 Nebraska winter wheat genotypes using genotyping-by-sequencing. Front Genet. 2018; 9: 76

[64] Danecek P, Auton A, Abecasis G, et al. The variant call format and VCFtools. Bioinf (Oxf). 2011; 27: 2156-8

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

[66] Wang J, Zhang Z, Gong Z, et al. Analysis of the genetic structure and diversity of upland cotton groups in different planting areas based on SNP markers. Gene. 2022; 809: 146042

[67] Khan MMH, Rafii MY, Ramlee SI, et al. DNA fingerprinting, fixation-index (Fst), and admixture mapping of selected Bambara groundnut (Vigna subterranea [L.] Verdc.) accessions using ISSR markers system . Sci Rep. 2021; 11: 14527

[68] Bradbury P, Zhang Z, Kroon DE, et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinf (Oxf). 2007; 23: 2633-5

[69] Wang J, Yu J, Lipka AE, et al. Interpretation of Manhattan plots and other outputs of genome-wide association studies. Methods Mol Biol (Clifton, NJ). 2022; 2481: 63-80

[70] Yu T, Zhang J, Cao J, et al. Identification of multiple genetic loci related to low-temperature tolerance during germination in maize (Zea maize L.) through a genome-wide association study . Curr Issues Mol Biol. 2023; 45: 9634-55

[71] Chen YJ, Niu SZ, Deng XY, et al. Genome-wide association study of leaf-related traits in tea plant in Guizhou based on genotyping-by-sequencing. BMC Plant Biol. 2023; 23: 196

[72] Rao K, OuYang Y, Chen Y, et al. Identification of CsAK as a critical caffeine-related upstream gene in tea accessions through genome-wide association study. J Integr Agric. 2025; 25: 2095-3119

[73] Xie Z, Sundström JF, Jin Y, et al. A selection strategy in plant transformation based on antisense oligodeoxynucleotide inhibition. Plant J. 2014; 77: 954-61

[74] Xiang F, Su Y, Zhou L, et al. Gibberellin promotes theanine synthesis by relieving the inhibition of CsWRKY71 on CsTSI in tea plant (Camellia sinensis) . Hortic Res. 2024; 12: uhae317

[75] Sunilkumar G, Vijayachandra K, Veluthambi K . Preincubation of cut tobacco leaf explants promotes agrobacterium-mediated transformation by increasing Vir gene induction. Plant Sci. 1999; 141: 51-8

[76] Chandra A, Hazra KK, Parihar AK, et al. ‘ImageJ’ based screening method for field pea germplasm against pea leaf miner (Phytomyza horticola Goureau) . Nat Acad Sci Lett. 2025; 48: 7-11