青藏高原地区3种蒿属植物转录组比较分析

doi:10.7525/j.issn.1673-5102.2022.02.005

植物研究 ›› 2022, Vol. 42 ›› Issue (2): 200-210.doi: 10.7525/j.issn.1673-5102.2022.02.005

青藏高原地区3种蒿属植物转录组比较分析

余静雅¹^,², 夏铭泽¹^,², 徐浩¹^,², 张发起¹()

^1.中国科学院西北高原生物研究所高原生物适应与进化重点实验室，西宁 810001
^2.中国科学院大学生命科学学院，北京 100049

收稿日期:2020-10-08 出版日期:2022-03-20 发布日期:2022-02-22
通讯作者: 张发起 E-mail:fqzhang@nwipb.cas.cn
作者简介:余静雅（1995—），女，博士研究生，主要从事植物系统与进化研究。
基金资助:
第二次青藏高原综合科学考察研究(2019QZKK0502);青海省基础研究计划项目(2019-ZJ-7042)

Comparative Transcriptome Analysis of Three Artemisia Species in Qinghai Tibet Plateau

Jingya Yu¹^,², Mingze Xia¹^,², Hao Xu¹^,², Faqi Zhang¹()

^1.Key Laboratory of Adaptation and Evolution of Plateau Biota，Northwest Institute of Plateau Biology，Chinese Academy of Sciences，Xining 810001
^2.College of Life Sciences，University of Chinese Academy of Sciences，Beijing 100049

Received:2020-10-08 Online:2022-03-20 Published:2022-02-22
Contact: Faqi Zhang E-mail:fqzhang@nwipb.cas.cn
About author:Yu Jingya（1995—），female，doctoral student，major in plant systematic and evolution.
Supported by:
Second Tibetan Plateau Scientific Expedition and Research（STEP） Program(2019QZKK0502);Basic Research Project of Qinghai Province(2019-ZJ-7042)

摘要/Abstract

摘要：

米蒿（Artemisia dalai-lamae）、冷蒿（Artemisia frigida）和臭蒿（Artemisia hedinii）作为西北干旱半干旱地区常见植物，具有防风固沙、阻止草原退化、杀菌消炎和治疗多种疾病的功效。通过Illumina测序技术分别对3种蒿属植物进行转录组测序和分析，分别获得54 268 322，46 434 864和43 971 646条Clean reads，并对测序结果进行了SNP和可变剪接分析。在米蒿和冷蒿对比组（ADL vs AF）中筛选到6 107个差异基因，在米蒿和臭蒿对比组（ADL vs AH）中筛选到4 822个差异基因，在冷蒿和臭蒿对比组（AF vs AH）中筛选到3 755个差异基因。GO富集将差异基因注释到生物过程、细胞组成和分子功能这3个类别共47个条目。KEGG代谢通路分析将3个对比组的差异基因分别富集到198、198和197条通路。对每个对比组的前10个高表达差异基因进行了分析，并筛选出25个生物碱、萜类和类黄酮的活性成分可能相关的差异基因。本研究为蒿属植物的物种鉴定、抗逆性研究和资源利用提供一定的科学依据。

关键词: 蒿属, 转录组比较分析, SNP, 可变剪切

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

中图分类号:

Q949.783.5

余静雅, 夏铭泽, 徐浩, 张发起. 青藏高原地区3种蒿属植物转录组比较分析[J]. 植物研究, 2022, 42(2): 200-210.

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.

图/表 8

表1

图1

表2

图2

图3

图4

表3

高表达基因功能分析

对比组 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	—	—

表3

表4

活性成分相关基因分类

基因编号 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

表4

参考文献 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]	张昊楠, 陈珊珊, 徐建民, 罗萍, 王晓萍, 许志茹, 范春节. 巨桉EgrWAT1基因克隆和功能初步分析[J]. 植物研究, 2023, 43(4): 601-611.
[2]	尹一卜, 李吉祥, 郭樱杰, 芦子廷, 肖英, 刘华领, 詹亚光, 曾凡锁. 白蜡属种间杂交子代木质素含量变异及FmPAL核苷酸多态性关联分析[J]. 植物研究, 2022, 42(2): 191-199.
[3]	李义良, 赵奋成, 刘阳, 钟岁英, 林昌明, 谭志强, 李福明, 廖仿炎, 吴惠姗, 郭文冰, 王哲. 湿加松亲本间SNP遗传距离与生长性状杂种优势关系研究[J]. 植物研究, 2021, 41(5): 738-743.
[4]	张序, 张秀姣, 马永鹏, 李正红, 万友名, 马宏. 基于GBS简化基因组技术的宽杯杜鹃遗传多样性分析[J]. 植物研究, 2021, 41(3): 429-435.
[5]	张珊珊, 康洪梅, 杨文忠. 基于简化基因组技术的云南蓝果树群体遗传分析[J]. 植物研究, 2019, 39(6): 899-907.
[6]	杨赟, 陈梦娇, 杜庆鑫, 朱景乐, 杜红岩, 杨绍彬. 基于全基因组重测序技术的‘红叶’杜仲SNP位点开发[J]. 植物研究, 2019, 39(6): 947-954.
[7]	段义忠, 王建武, 杜忠毓, 亢福仁. 基于SLAF-seq简化基因组技术的沙冬青SNP位点开发及遗传分析[J]. 植物研究, 2018, 38(1): 141-147.
[8]	王艳红1,2；贾庆斌1；张磊1；张含国1*. 长白落叶松4CL基因单核苷酸多态性及其与木材材性的关联分析[J]. 植物研究, 2016, 36(2): 242-251.
[9]	高刚；邓家彬；苟学梅；王茜；张艳；杨瑞武. 14个赖草属植物DREB2*基因的克隆与SNP分析[J]. 植物研究, 2015, 35(1): 47-51.
[10]	王铁娟;韩国栋. 中国蒿属沙生地理替代种亲缘关系的RAPD分析[J]. 植物研究, 2007, 27(4): 449-454.
[11]	陆毛珍;廉永善*. 甘肃绿绒蒿属（罂粟科）一新变种—光果红花绿绒蒿[J]. 植物研究, 2006, 26(1): 8-9.
[12]	李冀云. 中国东北蒿属光砂蒿的两个新变种[J]. 植物研究, 2003, 23(4): 388-389.
[13]	朱宗元, 梁存柱, 王炜, 刘钟龄. 紊蒿属一新种和对该属分类及演化的讨论[J]. 植物研究, 2003, 23(2): 147-153.
[14]	王臣, 路芳, 关旸, 张贵一. 东北蒿属艾组12种植物核型研究[J]. 植物研究, 2001, 21(2): 215-221.
[15]	王凌诗. 中国东北蒿属龙蒿组植物核型研究[J]. 植物研究, 2000, 20(4): 402-410.

青藏高原地区3种蒿属植物转录组比较分析

Comparative Transcriptome Analysis of Three Artemisia Species in Qinghai Tibet Plateau

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 8

参考文献 33

相关文章 15

编辑推荐

Metrics

本文评价