Integrative genomic resequencing and transcriptome atlas identify InDel markers for identifying tea plants resistant to anthracnose

Abstract

The development of insertion‒deletion (InDel) markers provides a valuable tool for genetic and genomic studies in tea (Camellia sinensis) plants. In this study, 4,370,229 InDel markers were identified through genome resequencing of Lingtou Dancong (BY) and Shuixian 7 (SX7). Transcriptome based on RNA sequencing revealed differential expression of 1,125 genes involved in plant-pathogen interaction pathways. Based on genome resequencing of the two tea varieties with a shared genetic background and RNA-seq dataset for post-anthracnose infection, molecular markers for anthracnose resistance in tea plants were developed by identifying specific InDels within the disease resistance-associated interval. The InDel marker linked to anthracnose resistance was located at 123,508,930 bp on chromosome Chr10. Tea plants with DNA bands of 590-bp or heterozygous bands amplified by the primer were resistant to the disease, while those with 431-bp bands were susceptible. When tested across 206 different tea varieties, the Indel marker CsRc demonstrated an 88.59% genotype–phenotype match, confirming its applicability for detecting anthracnose resistance/susceptibility in tea plants. Additionally, the CsRc marker was found to be located within the serine/threonine protein kinase gene, suggesting that ATP plays a critical role in the metabolic processes of tea plants. Taken together, our findings unravel important Indel markers for efficiently distinguishing valuable tea genotpes with desirable resistance to anthracnose, which holds great potential for facilitating disease resistance breeding in the future.

Keywords

Camellia sinensis / InDel / Genome resequencing / Transcriptome / Metabolome / Colletotrichum

References

[1] An Y, Mi X, Zhao S, Guo R, Xia X, Liu Set al.. Revealing distinctions in genetic diversity and adaptive evolution between two varieties of Camellia sinensis by whole-genome resequencing. Front Plant Sci, 2020, 11 603819

[2] Anarjan MB, Bae I, Lee S. Marker-assisted evaluation of two powdery mildew resistance candidate genes in Korean cucumber inbred lines. Agronomy, 2021, 11 2191

[3] Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JMet al.. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium Nat Genet, 2000, 25: 25-29

[4] Barone A. Molecular marker-assisted selection for potato breeding. Am J Potato Res, 2004, 81: 111-117

[5] Brody H. Tea. Nature, 2019, 566 S1

[6] Chen D, Cao Y, Li H, Kim D, Ahsan N, Thelen Jet al.. Extracellular ATP elicits DORN1-mediated RBOHD phosphorylation to regulate stomatal aperture. Nat Commun, 2017, 8: 13

[7] Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Yet al.. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant, 2020, 13: 1194-1202

[8] Chivasa S, Ndimba BK, Simon WJ, Lindsey K, Slabas AR. Extracellular ATP functions as an endogenous external metabolite regulating plant cell viability. Plant Cell, 2005, 17: 3019-3034

[9] Cho SH, Tóth K, Kim D, Vo PH, Lin CH, Handakumbura PPet al.. Activation of the plant mevalonate pathway by extracellular ATP. Nat Commun, 2022, 13 450

[10] Cho SH, Nguyen CT, Pham AQ, Stacey G. Computational prediction and in vitro analysis of the potential ligand binding site within the extracellular ATP receptor, P2K2. Plant Signal Behav, 2023, 18: 6

[11] Choi J, Tanaka K, Cao Y, Qi Y, Qiu J, Liang Yet al.. Identification of a plant receptor for extracellular ATP. Science, 2014, 343: 290-294

[12] Choi YR, Lee JY, Hwang S, Kim HU. PCR-based InDel marker associated with powdery mildew-resistant MR-1. Agronomy, 2020

[13] Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang Let al.. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly, 2012, 6: 80-92

[14] Clarke C, Richter BS, Rathinasabapathi B. Genetic and morphological characterization of United States tea (Camellia sinensis): insights into crop history, breeding strategies, and regional adaptability. Front Plant Sci, 2023, 14 1149682

[15] Dark A, Demidchik V, Richards SL, Shabala S, Davies JM. Release of extracellular purines from plant roots and effect on ion fluxes. Plant Signal Behav, 2011, 6: 1855-1857

[16] Drew L. The growth of tea. Nature, 2019, 566: S2-4

[17] Eisenstein M. Tea's value as a cancer therapy is steeped in uncertainty. Nature, 2019, 566: S6-7

[18] Foolad MR, Panthee DR. Marker-assisted selection in tomato breeding. Crit Rev Plant Sci, 2012, 31: 93-123

[19] Foria S, Magris G, Copetti D, Coleman C, Morgante M, Di Gaspero G. Indel markers for monitoring the introgression of downy mildew resistance from wild relatives into grape varieties. Mol Breed, 2018, 38 124

[20] Gilbert N. The science of tea's mood-altering magic. Nature, 2019, 566: S8-9

[21] Hanks SK, Hunter T. Protein kinases 6. The eukaryotic protein kinase superfamily: kinase (catalytic) domain structure and classification. FASEB J. 1995;9:576–96. https://doi.org/10.1096/fasebj.9.8.7768349.

[22] Hassan MZ, Rahim MA, Natarajan S, Robin AHK, Kim HT, Park JIet al.. Gummy stem blight resistance in melon: inheritance pattern and development of molecular markers. Int J Mol Sci, 2018, 19 2914

[23] Hassan O, Kim S, Kim K, Chang T. First report of leaf anthracnose caused by Colletotrichum camelliae on tea plants (Camellia sinensis) in South Korea. Plant Dis, 2023, 107: 2881

[24] Hayashi K, Yoshida H, Ashikawa I. Development of PCR-based allele-specific and indel marker sets for nine rice blast resistance genes. Theor Appl Genet, 2006, 113: 251-260

[25] Islam MR, Hossain MR, Jesse DMI, Jung HJ, Kim HT, Park JIet al.. Development of molecular marker linked with bacterial fruit blotch resistance in melon (Cucumis melo L.). Genes, 2020, 11 220

[26] Jander G, Norris SR, Rounsley SD, Bush DF, Levin IM, Last RL. Arabidopsis map-based cloning in the post-genome era. Plant Physiol, 2002, 129: 440-450

[27] Jin JQ, Liu YF, Ma CL, Ma JQ, Hao WJ, Xu YXet al.. A novel f3′5′h allele with 14 bp deletion is associated with high catechin index trait of wild tea plants and has potential use in enhancing tea quality. J Agric Food Chem, 2018, 66: 10470-10478

[28] Jo J, Venkatesh J, Han K, Lee HY, Choi GJ, Lee HJet al.. Molecular mapping of PMR1, a novel locus conferring resistance to powdery mildew in pepper (Capsicum annuum). Front Plant Sci, 2017, 8 2090

[29] Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh Met al.. KEGG for linking genomes to life and the environment. Nucleic Acids Res, 2008, 36: D480-D484

[30] Khan RS, Ntui VO, Chin DP, Nakamura I, Mii M. Production of marker-free disease-resistant potato using isopentenyl transferase gene as a positive selection marker. Plant Cell Rep, 2011, 30: 587-597

[31] Kim S, Kim CW, Choi MS, Kim S. Development of a simple PCR marker tagging the Allium roylei fragment harboring resistance to downy mildew (Peronospora destructor) in onion (Allium cepa L.). Euphytica, 2016, 208: 561-569

[32] Kim B, Hwang IS, Lee HJ, Oh CS. Combination of newly developed SNP and InDel markers for genotyping the Cf-9 locus conferring disease resistance to leaf mold disease in the tomato. Mol Breed, 2017, 37 59

[33] Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol, 2019, 37: 907-915

[34] Kim D, Chen DQ, Ahsan N, Jorge GL, Thelen JJ, Stacey G. The raf-like MAPKKK integrin-linked kinase 5 regulates purinergic receptor-mediated innate immunity in Arabidopsis. Plant Cell, 2023, 35: 1572-1592

[35] Kumar S, Tripathi D, Okubara PA, Tanaka K. Purinoceptor P2K1/DORN1 enhances plant resistance against a soilborne fungal pathogen, Rhizoctonia solani. Front Plant Sci, 2020, 11 572920

[36] Kumari S, Nazir F, Maheshwari C, Kaur H, Gupta R, Siddique KHMet al.. Plant hormones and secondary metabolites under environmental stresses: enlightening defense molecules. Plant Physiol Biochem, 2024, 206 108238

[37] Li X, Li MH, Deng WW, Ahammed GJ, Wei JP, Yan Pet al.. Exogenous melatonin improves tea quality under moderate high temperatures by increasing epigallocatechin-3-gallate and theanine biosynthesis in Camellia sinensis L.. J Plant Physiol, 2020, 253 153273

[38] Liao Y, Smyth GK, Shi W. Featurecounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics, 2014, 30: 923-930

[39] Liu B, Wang Y, Zhai W, Deng J, Wang H, Cui Yet al.. Development of indel markers for Brassica rapa based on whole-genome re-sequencing. Theor Appl Genet, 2013, 126: 231-239

[40] Liu L, Dang PM, Chen CY. Development and utilization of indel markers to identify peanut (Arachis hypogaea) disease resistance. Front Plant Sci, 2015, 6 988

[41] Liu G, Xu JH, Zhang M, Li PF, Yao XF, Hou Qet al.. Exploiting illumina sequencing for the development of InDel markers in watermelon (Citrullus lanatus). J Hortic Sci Biotechol, 2016, 91: 220-226

[42] Liu G, Zhao T, You X, Jiang J, Li J, Xu X. Molecular mapping of the Cf-10 gene by combining SNP/InDel-index and linkage analysis in tomato (Solanum lycopersicum). BMC Plant Biol, 2019, 19 15

[43] Liu S, An Y, Tong W, Qin X, Samarina L, Guo Ret al.. Characterization of genome-wide genetic variations between two varieties of tea plant (Camellia sinensis) and development of InDel markers for genetic research. BMC Genomics, 2019, 20 935

[44] Liu Y, Hao X, Lu Q, Zhang W, Zhang H, Wang Let al.. Genome-wide identification and expression analysis of flowering-related genes reveal putative floral induction and differentiation mechanisms in tea plant (Camellia sinensis). Genomics, 2020, 112(3): 2318-2326

[45] Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2^-ΔΔCt method. Methods, 2001, 25: 402-408

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

[47] Lv HH, Yang LM, Kang JG, Wang QB, Wang XW, Fang ZYet al.. Development of InDel markers linked to fusarium wilt resistance in cabbage. Mol Breed, 2013, 32: 961-967

[48] Ma S, An C, Lawson AW, Cao Y, Sun Y, Tan EYJet al.. Oligomerization-mediated autoinhibition and cofactor binding of a plant NLR. Nature, 2024, 632: 869-876

[49] Mills RE, Luttig CT, Larkins CE, Beauchamp A, Tsui C, Pittard WSet al.. An initial map of insertion and deletion (INDEL) variation in the human genome. Genome Res, 2006, 16: 1182-1190

[50] Mishra A, Wickneswari R, Bhuiyan MAR, Jena KK, Shamsudin NA. Broad spectrum blast resistance alleles in newly developed Malaysian rice (Oryza sativa L.) genotypes. Euphytica. 2021;217:8. https://doi.org/10.1007/s10681-020-02738-z.

[51] Mphangwe NIK, Vorster J, Steyn JM, Nyirenda HE, Taylor NJ, Apostolides Z. Screening of tea (Camellia sinensis) for trait-associated molecular markers. Appl Biochem Biotechnol, 2013, 171: 437-449

[52] Nvsvrot T, Xia W, Xiao Z, Zhan C, Liu M, Yang Xet al.. Combining QTL mapping with genome resequencing identifies an Indel in an R gene that is associated with variation in leaf rust disease resistance in poplar. Phytopathology, 2020, 110: 900-906

[53] Orrock JM, Rathinasabapathi B, Richter BS. Anthracnose in US tea: pathogen characterization and susceptibility among six tea accessions. Plant Dis, 2020, 104: 1055-1059

[54] Pham AQ, Cho SH, Nguyen CT, Stacey G. Arabidopsis lectin receptor kinase P2K2 is a second plant receptor for extracellular ATP and contributes to innate immunity. Plant Physiol, 2020, 183: 1364-1375

[55] Qiu T, Cao FR. Leveraging the advantages of Lingtou single-variety tea industry to promote the rural revitalization of Lingtou Village. Guangdong Tea Industry. 2022;4:32–5(Chinese).

[56] Rao IS, Srikanth B, Kishore VH, Suresh PB, Chaitanya U, Vemireddy LRet al.. Indel polymorphism in sugar translocation and transport genes associated with grain filling of rice (Oryza sativa L.). Mol Breed, 2011, 28: 683-691

[57] Shan Y, Zhang D, Luo Z, Li T, Qu H, Duan Xet al.. Advances in chilling injury of postharvest fruit and vegetable: Extracellular ATP aspects. Compr Rev Food Sci Food Saf, 2022, 21: 4251-4273

[58] Tanaka K, Gilroy S, Jones AM, Stacey G. Extracellular ATP signaling in plants. Trends Cell Biol, 2010, 20: 601-608

[59] Thakur O, Randhawa GS. Identification and characterization of SSR, SNP and InDel molecular markers from RNA-seq data of guar (Cyamopsis tetragonoloba, L. Taub.) roots. BMC Genomics, 2018, 191 951

[60] Tian Z, Yan C, Qian Q, Yan S, Xie H, Wang Fet al.. Development of gene-tagged molecular markers for starch synthesis-related genes in rice. Chin Sci Bull, 2010, 55: 3768-3777

[61] Uzunova M, Ecke W, Weissleder K, Röbbelen G. Mapping the genome of rapeseed (Brassica napus L.). I. Construction of an RFLP linkage map and localization of QTLs for seed glucosinolate content. Theor Appl Genet, 1995, 90: 194-204

[62] Varet H, Brillet-Guéguen L, Coppée JY, Dillies MA. SARTools: a DESeq2- and EdgeR-based R pipeline for comprehensive differential analysis of RNA-seq data. PLoS ONE, 2016, 11 e0157022

[63] Wang L, Wang Y, Cao H, Hao X, Zeng J, Yang Yet al.. Transcriptome analysis of an anthracnose-resistant tea plant cultivar reveals genes associated with resistance to Colletotrichum camelliae. PLoS ONE, 2016, 11 e0148535

[64] Wang Y, Liu H, Li H, Teng RM, Zhuang J. Genome-based identification and analysis of the genes involved in auxin biosynthesis and signal transduction during tea plant leaf development. Sci Hortic, 2020, 261 109030

[65] Wang M, Li J, Liu X, Liu C, Qian J, Yang Jet al.. Characterization of key odorants in lingtou dancong oolong tea and their differences induced by environmental conditions from different altitudes. Metabolites, 2022, 12 1063

[66] Wei C, Yang H, Wang S, Zhao J, Liu C, Gao Let al.. Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality. Proc Natl Acad Sci U S A, 2018, 115: E4151-E4158

[67] Wu PF, Gong HE, Yao XH, Chang WX, Wang KL. Evaluating the resistance of Camellia oleifera clones to Anthracnose. Forestry Science Research. 2023;31:158–63. https://doi.org/10.13275/j.cnki.lykxyj.2018.04.022.

[68] Zhao Y, Wang K, Wang WL, Yin TT, Dong WQ, Xu CJ. A high-throughput SNP discovery strategy for RNA-seq data. BMC Genomics, 2019, 20 160

[69] Zheng Q, Guo L, Huang J, Hao X, Li X, Li Net al.. Comparative transcriptomics provides novel insights into the mechanisms of selenium accumulation and transportation in tea cultivars (Camellia sinensis (L.) O. Kuntze). Front Plant Sci, 2023, 14 1268537

[70] Zhou G, Zhang Q, Tan C, Zhang XQ, Li C. Development of genome-wide indel markers and their integration with SSR, DArT and SNP markers in single barley map. BMC Genomics, 2015, 16 804

Funding

Open Fund of Collaborative Innovation Center of Chinese Oolong Tea Industry(2024W04)

Study on the effect of slope interplanting herbs on tea garden ecology and tea quality improvement(2024I0030)

RIGHTS & PERMISSIONS

The Author(s)