Comparative Transcriptome Analysis of Three Artemisia Species in Qinghai Tibet Plateau

doi:10.7525/j.issn.1673-5102.2022.02.005

Abstract

Abstract:

Artemisia dalai-lamae， Artemisia frigida and Artemisia hedinii， the common herbs in arid and semi-arid areas of Northwest China， have the effects of wind prevention and sand fixation， preventing grassland degradation， sterilization and anti-inflammatory， and treating a variety of diseases. The transcriptome of three Artemisia species were sequenced by Illumina method and then 54 268 322， 46 434 864 and 43 971 646 Clean reads were obtained， respectively. SNP and alternative splicing analyses were performed on the sequencing results. Differentially expressed genes were screened from the comparison groups of 6 107（ADL vs AF）， 4 822（ADL vs AH） and 3 755（AF vs AH） respectively. Go enrichment annotated the differentially expressed genes into three categories： biological process， cell component and molecular function. KEGG enrichment annotated the differentially expressed genes into 198， 198 and 197 pathways respectively. The top ten high expressed genes in each group were analyzed. We screened out 25 differentially expressed genes related to the active components of alkaloids，terpenoids and flavonoid. This study provid a scientific basis for species identification， stress resistance research and resource utilization of Artemisia.

Key words: Artemisia, comparative transcriptome analysis, SNP, alternative splicing

CLC Number:

Q949.783.5

Jingya Yu, Mingze Xia, Hao Xu, Faqi Zhang. Comparative Transcriptome Analysis of Three Artemisia Species in Qinghai Tibet Plateau[J]. Bulletin of Botanical Research, 2022, 42(2): 200-210.

Figures/Tables 8

Table 1

Fig.1

Table 2

Fig.2

Fig.3

Fig.4

Table 3

Functional analysis of high expression genes

对比组 Group	基因编号 Gene ID	log2（FC）	GO	KEGG	KOG
ADL vs AF	CTI12_AA551550	23.73	—	—	—
	CTI12_AA109090	23.71	氨酰-tRNA链接酶活性 Aminoacyl-tRNA ligase activity	—	—
	CTI12_AA619470	22.84	防御反应Defense response	—	—
	CTI12_AA626460	22.56	—	—	—
	CTI12_AA388450	22.07	防御反应Defense response	—	—
	CTI12_AA097740	21.79	防御反应Defense response	—	—
	CTI12_AA345130	21.75	—	—	—
	CTI12_AA130990	21.50	DNA结合 DNA binding	—	—
	CTI12_AA254130	21.31	—	—	—
	CTI12_AA047940	21.25	脂质运输Lipid transport	—	—
ADL vs AH	CTI12_AA106760	25.68	铁离子结合Iron ion binding	—	C-4甾醇甲基氧化酶 C-4 sterol methyl oxidase
	CTI12_AA274180	24.85	铁离子结合Iron ion binding	—	C-4甾醇甲基氧化酶 C-4 sterol methyl oxidase
	CTI12_AA551550	23.73	—	—	—
	CTI12_AA619470	22.84	防御反应Defense response	—	—
	CTI12_AA626460	22.56	—	—	—
	CTI12_AA485110	22.39	—	—	—
	CTI12_AA284240	22.35	—	—	—
	CTI12_AA465820	22.12	—	大亚基核糖体蛋白LP1 Large subunit ribosomal protein LP1	60s酸性核糖体蛋白P1 60s acidic ribosomal protein P1
	CTI12_AA001800	22.00	—	—	—
	CTI12_AA097740	21.79	防御反应Defense response	—	—
AF vs AH	CTI12_AA106760	25.47	铁离子结合Iron ion binding	—	C-4甾醇甲基氧化酶 C-4 sterol methyl oxidase
	CTI12_AA274180	23.77	铁离子结合Iron ion binding	—	C-4甾醇甲基氧化酶 C-4 sterol methyl oxidase
	CTI12_AA396830	23.18	DNA结合 DNA binding	—	—
	CTI12_AA161720	23.11	—	—	分子伴侣（小热休克蛋白Hsp26/Hsp42） Molecular chaperone （small heat-shock protein Hsp26/Hsp42）
	CTI12_AA622290	22.68	转运酶活性 Transferase activity	谷胱甘肽转移酶 Glutathione S-transferase	谷胱甘肽转移酶 Glutathione S-transferase
	CTI12_AA056630	22.64	氧化还原酶活性 Oxidoreductase activity	—	铁/抗坏血酸盐氧化还原酶 Iron/ascorbate family oxidoreductases
	CTI12_AA169950	22.55	—	—	—
	CTI12_AA351280	22.51	—	HSP20家族蛋白 HSP20 family protein	分子伴侣（小热休克蛋白Hsp26/Hsp42） Molecular chaperone （small heat-shock protein Hsp26/Hsp42）
	CTI12_AA366420	22.37	核糖体的结构成分 Structural constituent of ribosome	大亚基核糖体蛋白L10 Large subunit ribosomal protein L10	60s亚基核糖体蛋白L10 60s ribosomal protein L10
	CTI12_AA559880	22.27	DNA结合 DNA binding	—	—

Table 3

Table 4

Gene classification of active ingredients

基因编号 Gene ID	log2（FC）			基因名称 Gene name	基因描述 Desription	KO编号 KO_ID
基因编号 Gene ID	ADL vs AF	ADL vs AH	AF vs AH	基因名称 Gene name	基因描述 Desription	KO编号 KO_ID
生物碱Alkaloi
CTI12_AA405270	1.14	-1.43	-1.17	ALDH	醛脱氢酶（NAD+） Aldehyde dehydrogenase（NAD+）	K00128
CTI12_AA407480	-2.76	1.15	3.91	E2.1.1.104	咖啡酰辅酶a甲基转移酶 Caffeoyl-CoA O-methyltransferase	K00588
CTI12_AA512880	-2.47	1.35	3.81
CTI12_AA477420	3.95	2.57	-1.03
CTI12_AA528100	5.42	2.78	-2.64
CTI12_AA466140	-1.67	1.02	2.69	OGDH，sucA	2-氧戊二酸脱氢酶E2成分2-oxoglutarate dehydrogenase E2 component	K00164
CTI12_AA436910	-2.90	-1.54	1.36	GOT2	天冬氨酸转氨酶Aspartate aminotransferase	K14455
CTI12_AA496610	-2.98	-1.45	1.54	TAT	tyrosine aminotransferase（酪氨酸转氨酶）	K00815
萜类Terpenes
CTI12_AA205030	3.20	1.84	-1.36	DXS	1-脱氧-D-木酮糖-5-磷酸合酶 1-deoxy-D-xylulose-5-phosphate synthase	K01662
CTI12_AA270670	-1.68	2.02	3.70	VTE3，APG1	MPBQ/MSBQ甲基转移酶 MPBQ/MSBQ methyltransferase	K12502
CTI12_AA332590	-1.65	-3.32	-1.67	wrbA	NAD（P）H脱氢酶（醌） NAD（P）H dehydrogenase（quinone）	K03809
CTI12_AA374370	1.39	2.99	1.61	E2.1.1.95	生育酚O-甲基转移酶Tocopherol O-methyltransferase	K05928
CTI12_AA473110	2.54	3.78	1.23	GERD	大根香叶烯D合酶 Germacrene D synthase	K15803
CTI12_AA555350	-1.37	16.55	17.92	GERD	大根香叶烯D合酶 Germacrene D synthase	K15803
CTI12_AA495640	3.03	1.80	-1.22	DHDDS，RER2，SRT1	反式，多顺式—聚戊烯基二磷酸合酶Ditrans，polycis-polyprenyl diphosphate synthase	K11778
CTI12_AA496610	-2.98	-1.45	1.54	TAT	酪氨酸转氨酶Tyrosine aminotransferase	K00815
CTI12_AA506410	-2.98	-4.28	-1.30	E1.1.1.208	（+）-新薄荷醇脱氢酶（+）-neomenthol dehydrogenase	K15095
CTI12_AA510120	-1.29	2.78	4.07	CYP82G1	三甲基三癸四烯/二甲基壬三烯合酶Trimethyltridecatetraene/dimethylnonatriene synthase	K17961
类黄酮Flavonoid
CTI12_AA407480	-2.76	1.15	3.91	E2.1.1.104	咖啡酰辅酶a甲基转移酶 Caffeoyl-CoA O-methyltransferase	K00588
CTI12_AA512880	-2.47	1.35	3.81
CTI12_AA528100	5.42	2.78	-2.64
CTI12_AA445950	1.08	2.62	1.54	CYP98A，C3′H	5-O-（4-酰基）-D-奎尼酸3′-单氧酶 5-O-（4-coumaroyl）-D-quinate 3′-monooxygenase	K09754
CTI12_AA426020	5.80	3.30	-2.51	FLS	黄酮醇合成酶Flavonol synthase	K05278
CTI12_AA430920	3.16	4.70	1.53	CHS	查耳酮合酶 Chalcone synthase	K00660
CTI12_AA627070	3.08	4.15	1.07	CHS	查耳酮合酶 Chalcone synthase	K00660

Table 4

References 33

1	Funk V A.Systematics，evolution，and biogeography of compositae［M］.Vienna：IAPT，2009.
2	林有润.论世界蒿属植物区系［J］.植物研究，1995，15（1）：1-37.
	Lin Y R.On the floristics of Artemisia L.in the world［J］.Bulletin of Botanical Research，1995，15（1）：1-37.
3	Lin Y R，Christopher J H，Michael G G.Flora of China vol. 20-21［M］.Beijing：Science Press，2011，2：845-1753.
4	Bora K S，Sharma A.The genus Artemisia：a comprehensive review［J］.Pharmaceutical Biology，2011，49（1）：101-109.
5	贾丽娜，吴斌，丁国栋，等.宁夏盐池县3种常见沙生植物群落特征及阻沙能力［J］.东北林业大学学报，2009，37（11）：56-58，63.
	Jia L N，Wu B，Ding G D，et al.Community characters and sand trapping capability of three common species of psammophyte in Yanchi County，Ningxia［J］.Journal of Northeast Forestry University，2009，37（11）：56-58，63.
6	王静，杨持，王铁娟.冷蒿（Atemisia frigida）种群在放牧干扰下构件的变化［J］.生态学报，2006，26（3）：960-965.
	Wang J，Yang C，Wang T J.Changes of model of Atemisia frigida populations under the disturbance of grazing［J］.Acta Ecologica Sinica，2006，26（3）：960-965.
7	张汝民，王玉芝，侯平，等.几种牧草幼苗对冷蒿茎叶水浸提液化感作用的生理响应［J］.生态学报，2010，30（8）：2197-2204.
	Zhang R M，Wang Y Z，Hou P，et al.Physiological responses to allelopathy of aquatic stem and leaf extract of Artemisia frigida in seedling of several pasture plants［J］.Acta Ecologica Sinica，2010，30（8）：2197-2204.
8	杜品.青藏高原甘南藏药植物志［M］.兰州：甘肃科学技术出版社，2006.
	Du P.Flora of Tibetan medicine in Gannan of Qinghai Tibet Plateau［M］.Lanzhou：Gansu science and Technology Press，2006.
9	崔箭，唐丽，蓝蓉，等.蒙药阿给炭治疗支气管扩张咯血临床观察［J］.中央民族大学学报：自然科学版，2006，15（2）：149-151，155.
	Cui J，Tang L，Lan R，et al.Clinical analysis of 63 cases on curative effect of Arteisia frigida Wild in treatment of hemoptysis caused by bronchiectasis［J］.Journal of the CUN：Natural Sciences Edition，2006，15（2）：149-151，155.
10	Rashid M U，Alamzeb M，Ali S，et al.The chemistry and pharmacology of alkaloids and allied nitrogen compounds from Artemisia species：A review［J］.Phytother Research，2019，33（10）：2661-2684.
11	Shen Q，Zhang L D，Liao Z H，et al.The genome of Artemisia annua provides insight into the evolution of Asteraceae family and artemisinin biosynthesis［J］.Molecular Plant，2018，11（6）：776-788.
12	Bolger A M，Lohse M，Usadel B.Trimmomatic：a flexible trimmer for Illumina sequence data［J］.Bioinformatics，2014，30（15）：2114-2120.
13	Kim D，Langmead B，Salzberg S L.HISAT：a fast spliced aligner with low memory requirements［J］.Nature Methods，2015，12（4）：357-360.
14	Pertea A M，Pertea G M，Antonescu C M，et al.StringTie enables improved reconstruction of a transcriptome from RNA-seq reads［J］.Nature Biotechnology，2015，33（3）：290-295.
15	Li H.A statistical framework for SNP calling，mutation discovery，association mapping and population genetical parameter estimation from sequencing data［J］.Bioinformatics，2011，27（21）：2987-2993.
16	Cingolani P，Platts A，Wang L L，et 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［J］.Fly，2012，6（2）：80-92.
17	Florea L，Song L，Salzberg S L.Thousands of exon skipping events differentiate among splicing patterns in sixteen human tissues［J］.F1000res，2013，2：188.
18	Love M I，Huber W，Anders S.Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2［J］.Genome Biology，2014，15（12）：550.
19	Yu G C，Wang L G，Han Y Y，et al.clusterProfiler：an R package for comparing biological themes among gene clusters［J］.OMICS：A Journal of Integrative Biology，2012，16（5）：284-287.
20	夏铭泽，张雨，余静雅，等.多裂骆驼蓬叶片转录组分析［J］.广西植物，2021，41（04）：503-513.
	Xia M Z，Zhang Y，Yu J Y，et al.Transcriptome analysis for leaves of Peganum mulstecutm［J］.Guihaia，2021，41（04）：503-513.
21	Roberts W R，Roalson E H.Phylogenomic analyses reveal extensive gene flow within the magic flowers（Achimenes）［J］.American Journal of Botany，2018，105（4）：726-740.
22	冯雅岚，熊瑛，张均，等.可变剪切在植物发育和非生物胁迫响应中的作用［J］.核农学报，2020，34（1）：62-70.
	Feng Y L，Xiong Y，Zhang J，et al.Role of alternative splicing in plant development and abiotic stress responses［J］.Journal of Nuclear Agricultural Sciences，2020，34（1）：62-70.
23	Ding Y Q，Wang Y，Qiu C，et al.Alternative splicing in tea plants was extensively triggered by drought，heat and their combined stresses［J］.PeerJ，2020，8：e8258.
24	Singh A，Jindal S，Longchar B，et al.Overexpression of Artemisia annua sterol C-4 methyl oxidase gene，AaSMO1，enhances total sterols and improves tolerance to dehydration stress in tobacco［J］.Plant Cell， Tissue and Organ Culture（PCTOC），2015，121（1）：167-181.
25	Waters E R，Vierling E.Plant small heat shock proteins-evolutionary and functional diversity［J］.New Phytologist，2020，227（1）：24-37.
26	Alhaithloul H A S.Impact of combined heat and drought stress on the potential growth responses of the desert grass Artemisia sieberi alba：relation to biochemical and molecular adaptation［J］.Plants，2019，8（10）：416.
27	Horn P J，Benning C.The plant lipidome in human and environmental health［J］.Science，2016，353（6305）：1228-1232.
28	王垿，孙昕，李鹏飞，等.双对栅藻FACHB-78甘油三酯积累的盐胁迫条件优化［J］.中国环境科学，2019，39（12）：5248-5253.
	Wang X，Sun X，Li P F，et al.Optimization of salt stress condition for accumulation of triglycerides in Scenedesmus bijuga FACHB-78［J］.China Environmental Science，2019，39（12）：5248-5253.
29	王昕.氮磷胁迫对微藻糖脂累积的影响［D］.上海：上海交通大学，2016.
	Wang X.Effects of different nitrate and hydrophosphate concentration on glaycolipid accumulation in microalgae［D］.Shanghai：Shanghai Jiao Tong University，2016.
30	Ling C C，Zheng L J，Yu X R，et al.Cloning and functional analysis of three aphid alarm pheromone genes from German chamomile（Matricaria chamomilla L.）［J］.Plant Science，2020，294：110463.
31	Li X Q，Guo R R，Li J，et al.Genome-wide identification and analysis of the aldehyde dehydrogenase （ALDH） gene superfamily in apple（Malus × domestica Borkh.）［J］.Plant Physiology and Biochemistry，2013，71：268-282.
32	Pellicer J，Saslis-Lagoudakis C H，Carrió E，et al.A phylogenetic road map to antimalarial Artemisia species［J］.Journal of Ethnopharmacology，2018，225：1-9.
33	Ma D M，Wang Z L，Wang L J，et al.A genome-wide scenario of terpene pathways in self-pollinated Artemisia annua［J］.Molecular Plant，2015，8（11）：1580-1598.

[1]	Yibo Yin, Jixiang Li, Yingjie Guo, Ziting Lu, Ying Xiao, Hualing Liu, Yaguang Zhan, Fansuo Zeng. Variation of Lignin Content and Association Analysis of FmPAL Nucleotide Polymorphism in Progenies of Interspecific Hybrids of Fraxinus [J]. Bulletin of Botanical Research, 2022, 42(2): 191-199.
[2]	Yi-Liang LI, Fen-Cheng ZHAO, Yang LIU, Sui-Ying ZHONG, Chang-Ming LIN, Zhi-Qiang TAN, Fu-Ming LI, Fang-Yan LIAO, Hui-Shan WU, Wen-Bin GUO, Zhe WANG. Relationship Between Parental Genetic Distance and Heterosis of Growth Traitsin Pinus elliottii × P.caribaea var. Hondurensis Based on SNP marker [J]. Bulletin of Botanical Research, 2021, 41(5): 738-743.
[3]	Xu ZHANG, Xiu-Jiao ZHANG, Yong-Peng MA, Zheng-Hong LI, You-Ming WAN, Hong MA. Genetic Diversity Assessment of Rhododendron sinofalconeri with Genotyping by Sequencing(GBS) [J]. Bulletin of Botanical Research, 2021, 41(3): 429-435.
[4]	ZHANG Shan-Shan, KANG Hong-Mei, YANG Wen-Zhong. Population Genetic Analysis of Nyssa yunnanensis by Reduced-representation Sequencing Technique [J]. Bulletin of Botanical Research, 2019, 39(6): 899-907.
[5]	YANG Yun, CHEN Meng-Jiao, DU Qing-Xin, ZHU Jing-Le, DU Hong-Yan, YANG Shao-Bin. SNP Sites Developed by Whole Genome Resequencing Analysis in Eucommia ulmoides ‘Hongye’ [J]. Bulletin of Botanical Research, 2019, 39(6): 947-954.
[6]	DUAN Yi-Zhong, WANG Jian-Wu, DU Zhong-Yu, KANG Fu-Ren. SNP Sites Developed by Specific Length Amplification Fragment Sequencing(SLAF-seq) and Genetic Analysis in Ammopitanthus mongolicus [J]. Bulletin of Botanical Research, 2018, 38(1): 141-147.
[7]	LI Tai-Qiang, LIU Xiong-Fang, WAN You-Ming, LI Zheng-Hong, QI Guo-Hai, LI Yu-Ying, LIU Xiu-Xian, HE Rui, MA Yan, MA Hong. Transcriptome Analysis for Rhododendron longipedicellatum (Plant Species with Extremely Small Populations) Based on High Throughput Sequencing [J]. Bulletin of Botanical Research, 2017, 37(6): 825-834.
[8]	TIAN Zun-Zhe, GAO Qing-Bo, CHEN Shi-Long, ZHANG Fa-Qi. Characteristics of SSR and SNP in Gentiana veitchiorum in Qinghai-Tibetan Plateau,by High-throughput Sequencing [J]. Bulletin of Botanical Research, 2016, 36(5): 747-752.
[9]	WANG Yan-Hong1,2；JIA Qing-Bin1；ZHANG Lei1；ZHANG Han-Guo1. Allelic Variation in 4-coumarate:CoA ligase(4CL) Associated with Wood Properties of Larix olgensis* [J]. Bulletin of Botanical Research, 2016, 36(2): 242-251.
[10]	GAO Gang；DENG Jia-Bin；GOU Xue-Mei；WANG Qian；ZHANG Yan；YANG Rui-Wu. Cloning and SNP Analysis of DREB2* Gene in Fourteen Leymus Species [J]. Bulletin of Botanical Research, 2015, 35(1): 47-51.
[11]	XU Yue-Li;ZHANG Han-Guo;SHI Tian-Yuan;YAO Yu;ZHANG Lei;WANG Yan-Hong;LIU Ling. SNP Analysis of 4CL Gene in Larix olgensis* Henry Populations [J]. Bulletin of Botanical Research, 2013, 33(2): 208-213.
[12]	Zhang Gui-yi, Wang Ling-shi, Ma Han-xi. A NEW SUBSPECIES AND A NEW FORM OF Artemisia oxycephala IN CHINA [J]. Bulletin of Botanical Research, 1993, 13(2): 125-126.
[13]	Ling Yeou-ruenn. THE OLD WORLD ARTEMISIA LINN. (COMPOSITAE) [J]. Bulletin of Botanical Research, 1992, 12(1): 1-108.