Abstract: 【Objective】This study aims at analysing both expression patterns and regulatory pathways of tea plants in response to glyphosate stressing, through which could revealed the effect of glyphosate herbicides on tea plants at transcriptional level and identify key genes of tea plants.【Method】C. sinensis cv ‘Jin-guanyin’ was applied as material plant. A recommended concentration of glyphosate was irrigated to test plants.Samples were collected at different time intervals (0, 0.25, 1, 3 and 7 d). The samples were sequenced by transcriptome, the content of shikimic acid was also quantified.The WGCNA method was used to jointly analyze transcriptome and shikimic acid content data, identify co-expressed gene modules related to glyphosate response, and screen out key regulatory genes.【Results】the content of shikimic acid in tea leaves reduced gradually during first 3 days. However, it suddenly reached a peak on the 7th day( 6.99 times compared to no glyphosate treated sample).A total of 12568 differential expression genes (DEGs) were also identified,whichmainly enriched in phenylpropane, flavonoid biosynthesis and plant hormone signal transduction pathways. In addition, glyphosate treatment induced 24, 52, 31 and 69 genes respectively which related to shikimic acid metabolism, phenylpropane, flavonoid biosynthesis and hormone signal transduction pathways. A total of 19 modules were screened out by WGCNA method. Correlate analysis of transcriptome and shikimic acid content indicated two key modules which including 2024 and 2305 genes, respectively. The top 50 genes with the highest connectivity in the key modules were selected for co-expression analysis, and 6 key response genes were obtained, including 2 resistance genes (SHMTandRPM), 1 drug resistance gene (PDR), 1 ion transport gene (At), 1 membrane transport gene (GPT) and 1 transcription factor gene (ERF). 【Conclusion】Glyphosate could affect downstream genes transcription of phenylpropane, flavonoid biosynthesis and hormone signal transduction pathways by interfering shikimic acid metabolism of tea plants.In addition, the study also identified two co-expression modules closely related to glyphosate response, and found that multiple potential candidate genes and transcription factors can resist glyphosate stress, such asSHMT, RPM, At, PDR, ERFandGPT.

Key words: Camellia sinensis, glyphosate, Shikimic acid, transcriptome, WGCNA