草业学报 ›› 2021, Vol. 30 ›› Issue (2): 143-155.DOI: 10.11686/cyxb2020326
周晶1(), 陈思齐2, 史文娇2, 阳伏林2, 林辉1, 林占熺1
收稿日期:
2020-07-13
修回日期:
2020-10-10
出版日期:
2021-02-20
发布日期:
2021-01-19
通讯作者:
周晶
作者简介:
E-mail: zhoujing_lz@hotmail.com基金资助:
Jing ZHOU1(), Si-qi CHEN2, Wen-jiao SHI2, Fu-lin YANG2, Hui LIN1, Zhan-xi LIN1
Received:
2020-07-13
Revised:
2020-10-10
Online:
2021-02-20
Published:
2021-01-19
Contact:
Jing ZHOU
摘要:
巨菌草生物量巨大,根系发达,是一种理想的水土保持、防风固沙禾本科C4植物。为了探讨巨菌草转录组信息,对其幼叶、根样本进行Illumina高通量测序和差异表达基因分析。经组装拼接共获得210806个转录本(transcripts)和150336条单基因簇(unigene),平均长度分别为760和642 bp。将获得的unigene与通用公共数据库进行比对,共有88765条unigene获得了基因注释。通过将巨菌草幼叶与根进行比较,共鉴定到5735个差异表达基因(DEGs),其中上调3435个,下调2300个。对DEGs做GO功能分析,上调DEGs主要集中在电子传递,碳水化合物代谢,金属离子结合等过程,下调DEGs则与代谢过程,催化活性,小分子代谢过程等相关。KEGG pathway富集分析共得到120条代谢通路,其中参与光合生物的固碳作用通路DEGs富集程度最高;参与次生代谢物的生物合成通路的DEGs数目最多。对DEGs转录因子分析显示,bHLH和WRKY是调控叶、根差异表达的主要转录因子。研究结果极大地丰富了巨菌草转录组信息,为今后巨菌草分子生物学研究提供了宝贵的数据资源,同时为巨菌草品种改良提供了理论依据。
周晶, 陈思齐, 史文娇, 阳伏林, 林辉, 林占熺. 巨菌草幼叶及根转录组功能基因测序及分析[J]. 草业学报, 2021, 30(2): 143-155.
Jing ZHOU, Si-qi CHEN, Wen-jiao SHI, Fu-lin YANG, Hui LIN, Zhan-xi LIN. Transcriptome analyses of functional genes in young leaves and roots of Giant Juncao[J]. Acta Prataculturae Sinica, 2021, 30(2): 143-155.
基因 Gene ID | 引物 Primers | 基因 Gene ID | 引物序列 Primers |
---|---|---|---|
c104596_g7 | S: ATGAGACTGGTATAAGGATGTT A: GAAGCAACTTGGTGTTAGAG | c106511_g1 | S: TTGCTGGTTCATTCTCTG A: TTCCTATTGCCATCTCTTG |
c100788_g2 | S: TGGAATGTGGAATGAGAT A: TAGCAGTTGTAGATTACCT | c97722_g1 | S: GTCGCACAACAATCTTCT A: GCATACTATCGGTCTTCATC |
c112282_g1 | S: CAATCTTGTGACTGAATC A: GATGTTGTAGTTCTGAATC | c106537_g6 | S: GACCATCTTACTCTGACT A: CCAACCAATACTTCATCC |
c107530_g1 | S: AGTTAGTACCAGAGGATA A: AGCATCTATGTTCTTGTA | PgACT | S: TGCTCAGTGGAGGGTCTACCAT A: CAGGTGGTGCAATCACTTTAAT |
c99519_g1 | S: TACTTCTCCTCTAATGCTTGTG A: GGTAAGGTGCCTATTGTGAA |
表1 巨菌草幼叶、根8个DEGs和内参基因PgACT的qRT-PCR引物
Table 1 Primers of eight DEGs and reference gene PgACT used for quantitative real-time PCR in leaves and roots of Giant Juncao
基因 Gene ID | 引物 Primers | 基因 Gene ID | 引物序列 Primers |
---|---|---|---|
c104596_g7 | S: ATGAGACTGGTATAAGGATGTT A: GAAGCAACTTGGTGTTAGAG | c106511_g1 | S: TTGCTGGTTCATTCTCTG A: TTCCTATTGCCATCTCTTG |
c100788_g2 | S: TGGAATGTGGAATGAGAT A: TAGCAGTTGTAGATTACCT | c97722_g1 | S: GTCGCACAACAATCTTCT A: GCATACTATCGGTCTTCATC |
c112282_g1 | S: CAATCTTGTGACTGAATC A: GATGTTGTAGTTCTGAATC | c106537_g6 | S: GACCATCTTACTCTGACT A: CCAACCAATACTTCATCC |
c107530_g1 | S: AGTTAGTACCAGAGGATA A: AGCATCTATGTTCTTGTA | PgACT | S: TGCTCAGTGGAGGGTCTACCAT A: CAGGTGGTGCAATCACTTTAAT |
c99519_g1 | S: TACTTCTCCTCTAATGCTTGTG A: GGTAAGGTGCCTATTGTGAA |
核苷酸长度 Nucleotides length | 转录本 Transcripts | 单基因簇 Unigenes (No.) |
---|---|---|
<300 | 67010 | 59354 |
300~500 | 48148 | 37777 |
500~1000 | 45684 | 28059 |
1000~2000 | 34342 | 16704 |
>2000 | 15622 | 8442 |
总数Total | 210806 | 150336 |
最短长度 Shortest length | 201 | 201 |
中间长度 Median length | 444 | 359 |
最长长度 Longest length | 14723 | 14723 |
平均长度 Mean length | 760 | 642 |
不小于总长 50% N50 | 1231 | 995 |
不小于总长90% N90 | 297 | 261 |
表2 巨菌草转录组组装结果
Table 2 Summary of transcriptome assembly for Giant Juncao (bp)
核苷酸长度 Nucleotides length | 转录本 Transcripts | 单基因簇 Unigenes (No.) |
---|---|---|
<300 | 67010 | 59354 |
300~500 | 48148 | 37777 |
500~1000 | 45684 | 28059 |
1000~2000 | 34342 | 16704 |
>2000 | 15622 | 8442 |
总数Total | 210806 | 150336 |
最短长度 Shortest length | 201 | 201 |
中间长度 Median length | 444 | 359 |
最长长度 Longest length | 14723 | 14723 |
平均长度 Mean length | 760 | 642 |
不小于总长 50% N50 | 1231 | 995 |
不小于总长90% N90 | 297 | 261 |
图1 基因注释结果统计(a)及Nr注释的物种分布(b)
Fig.1 Unigene information annotated in different databases (a) and distributed into different species from Nr database (b)
图2 GO功能分布统计X轴代表比对到GO数据库的单基因被分为三大类和41个亚类;Y轴代表每种分类中基因的数量。The X-axis represents BLAST unigenes which were classified into three major categories and 41 sub-categories in GO. The Y-axis represents the number of genes in each categotry.
Fig.2 Functional distribution of GO annotation
GO登录号GO accession | 功能Function | 基因数量Gene number (No.) |
---|---|---|
GO:0007118 | 芽生细胞顶端芽生长 Budding cell apical bud growth | 2 |
GO:0048588 | 细胞生长发育 Developmental cell growth | 52 |
GO:0060560 | 涉及形态发生的生长发育 Developmental growth involved in morphogenesis | 67 |
GO:0080186 | 营养生长发育Developmental vegetative growth | 3 |
GO:0044182 | 单细胞生物种群的丝状生长Filamentous growth of a population of unicellular organisms | 48 |
GO:0035266 | 分生组织增长Meristem growth | 43 |
GO:0045926 | 生长负调节Negative regulation of growth | 4 |
GO:0030307 | 细胞生长的正调控Positive regulation of cell growth | 4 |
GO:0048639 | 发育生长的正调控Positive regulation of developmental growth | 4 |
GO:0090033 | 丝状生长的正调控Positive regulation of filamentous growth | 2 |
GO:0045927 | 生长正调节Positive regulation of growth | 10 |
GO:0007124 | 假菌丝生长Pseudohyphal growth | 3 |
GO:0001558 | 细胞生长调节Regulation of cell growth | 42 |
GO:0048638 | 发育生长调节Regulation of developmental growth | 108 |
GO:0010570 | 丝状生长的调节Regulation of filamentous growth | 8 |
GO:0040009 | 生长速度调节Regulation of growth rate | 4 |
GO:0046620 | 器官生长调节Regulation of organ growth | 10 |
GO:0080112 | 种子生长Seed growth | 1 |
GO:0009826 | 一维细胞生长Unidimensional cell growth | 33 |
表3 与巨菌草生长相关的GO富集分析
Table 3 GO enrichment analysis related to growth annotation in Giant Juncao
GO登录号GO accession | 功能Function | 基因数量Gene number (No.) |
---|---|---|
GO:0007118 | 芽生细胞顶端芽生长 Budding cell apical bud growth | 2 |
GO:0048588 | 细胞生长发育 Developmental cell growth | 52 |
GO:0060560 | 涉及形态发生的生长发育 Developmental growth involved in morphogenesis | 67 |
GO:0080186 | 营养生长发育Developmental vegetative growth | 3 |
GO:0044182 | 单细胞生物种群的丝状生长Filamentous growth of a population of unicellular organisms | 48 |
GO:0035266 | 分生组织增长Meristem growth | 43 |
GO:0045926 | 生长负调节Negative regulation of growth | 4 |
GO:0030307 | 细胞生长的正调控Positive regulation of cell growth | 4 |
GO:0048639 | 发育生长的正调控Positive regulation of developmental growth | 4 |
GO:0090033 | 丝状生长的正调控Positive regulation of filamentous growth | 2 |
GO:0045927 | 生长正调节Positive regulation of growth | 10 |
GO:0007124 | 假菌丝生长Pseudohyphal growth | 3 |
GO:0001558 | 细胞生长调节Regulation of cell growth | 42 |
GO:0048638 | 发育生长调节Regulation of developmental growth | 108 |
GO:0010570 | 丝状生长的调节Regulation of filamentous growth | 8 |
GO:0040009 | 生长速度调节Regulation of growth rate | 4 |
GO:0046620 | 器官生长调节Regulation of organ growth | 10 |
GO:0080112 | 种子生长Seed growth | 1 |
GO:0009826 | 一维细胞生长Unidimensional cell growth | 33 |
GO登录号GO_accession | 功能 Function | 种类 Ontology | P值 P-value | 错误发现率FDR | 差异基因个数 DEGs number (No.) |
---|---|---|---|---|---|
GO:0005507 | 铜离子结合Copper ion binding | MF | 2.49×10-12 | 6.34×10-9 | 35 |
GO:0042221 | 化学刺激反应Response to chemical stimulus | BP | 3.95×10-11 | 5.02×10-8 | 73 |
GO:0016209 | 抗氧化活性Antioxidant activity | MF | 4.07×10-10 | 3.44×10-7 | 42 |
GO:0009415 | 应对水Response to water | BP | 1.20×10-9 | 7.61×10-7 | 25 |
GO:0004601 | 过氧化物酶活性Peroxidase activity | MF | 2.16×10-9 | 9.13×10-7 | 38 |
GO:0016684 | 氧化还原酶活性,作用于过氧化氢为受体Oxidoreductase activity, acting on peroxide as acceptor | MF | 2.16×10-9 | 9.13×10-7 | 38 |
GO:0009269 | 应对干燥Response to desiccation | BP | 6.17×10-9 | 1.96×10-6 | 19 |
GO:0009414 | 对缺水的反应Response to water deprivation | BP | 6.17×10-9 | 1.96×10-6 | 19 |
GO:0043167 | 离子结合Ion binding | MF | 1.34×10-8 | 3.10×10-6 | 152 |
GO:0043169 | 阳离子绑定Cation binding | MF | 1.34×10-8 | 3.10×10-6 | 152 |
GO:0046872 | 金属离子结合Metal ion binding | MF | 1.34×10-8 | 3.10×10-6 | 152 |
GO:0006979 | 氧化应激反应Response to oxidative stress | BP | 5.05×10-8 | 1.00×10-5 | 33 |
GO:0009628 | 对非生物刺激的反应Response to abiotic stimulus | BP | 5.12×10-8 | 1.00×10-5 | 30 |
GO:0005975 | 碳水化合物代谢过程Carbohydrate metabolic process | BP | 1.38×10-7 | 2.51×10-5 | 191 |
GO:0005489 | 电子运输活动Electron transporter activity | MF | 3.03×10-7 | 5.13×10-5 | 56 |
GO:0006118 | 电子传递Electron transport | BP | 4.72×10-6 | 7.50×10-4 | 202 |
GO:0045735 | 水库营养活动Nutrient reservoir activity | MF | 1.51×10-5 | 2.26×10-3 | 18 |
GO:0007155 | 细胞黏附Cell adhesion | BP | 4.67×10-5 | 6.25×10-3 | 24 |
GO:0022610 | 生物黏附Biological adhesion | BP | 4.67×10-5 | 6.25×10-3 | 24 |
GO:0005509 | 钙离子结合Calcium ion binding | MF | 6.44×10-5 | 8.18×10-3 | 43 |
表4 巨菌草幼叶、根上调差异基因GO富集结果
Table 4 Giant Juncao GO enrichment of leaves and roots upregulated DEGs (FDR<0.05)
GO登录号GO_accession | 功能 Function | 种类 Ontology | P值 P-value | 错误发现率FDR | 差异基因个数 DEGs number (No.) |
---|---|---|---|---|---|
GO:0005507 | 铜离子结合Copper ion binding | MF | 2.49×10-12 | 6.34×10-9 | 35 |
GO:0042221 | 化学刺激反应Response to chemical stimulus | BP | 3.95×10-11 | 5.02×10-8 | 73 |
GO:0016209 | 抗氧化活性Antioxidant activity | MF | 4.07×10-10 | 3.44×10-7 | 42 |
GO:0009415 | 应对水Response to water | BP | 1.20×10-9 | 7.61×10-7 | 25 |
GO:0004601 | 过氧化物酶活性Peroxidase activity | MF | 2.16×10-9 | 9.13×10-7 | 38 |
GO:0016684 | 氧化还原酶活性,作用于过氧化氢为受体Oxidoreductase activity, acting on peroxide as acceptor | MF | 2.16×10-9 | 9.13×10-7 | 38 |
GO:0009269 | 应对干燥Response to desiccation | BP | 6.17×10-9 | 1.96×10-6 | 19 |
GO:0009414 | 对缺水的反应Response to water deprivation | BP | 6.17×10-9 | 1.96×10-6 | 19 |
GO:0043167 | 离子结合Ion binding | MF | 1.34×10-8 | 3.10×10-6 | 152 |
GO:0043169 | 阳离子绑定Cation binding | MF | 1.34×10-8 | 3.10×10-6 | 152 |
GO:0046872 | 金属离子结合Metal ion binding | MF | 1.34×10-8 | 3.10×10-6 | 152 |
GO:0006979 | 氧化应激反应Response to oxidative stress | BP | 5.05×10-8 | 1.00×10-5 | 33 |
GO:0009628 | 对非生物刺激的反应Response to abiotic stimulus | BP | 5.12×10-8 | 1.00×10-5 | 30 |
GO:0005975 | 碳水化合物代谢过程Carbohydrate metabolic process | BP | 1.38×10-7 | 2.51×10-5 | 191 |
GO:0005489 | 电子运输活动Electron transporter activity | MF | 3.03×10-7 | 5.13×10-5 | 56 |
GO:0006118 | 电子传递Electron transport | BP | 4.72×10-6 | 7.50×10-4 | 202 |
GO:0045735 | 水库营养活动Nutrient reservoir activity | MF | 1.51×10-5 | 2.26×10-3 | 18 |
GO:0007155 | 细胞黏附Cell adhesion | BP | 4.67×10-5 | 6.25×10-3 | 24 |
GO:0022610 | 生物黏附Biological adhesion | BP | 4.67×10-5 | 6.25×10-3 | 24 |
GO:0005509 | 钙离子结合Calcium ion binding | MF | 6.44×10-5 | 8.18×10-3 | 43 |
GO登录号GO_accession | 功能 Function | 种类 Ontology | P值 P-value | 错误发现率FDR | 差异基因个数DEGs number (No.) |
---|---|---|---|---|---|
GO:0009765 | 光合作用、光收获Photosynthesis, light harvesting | BP | 2.35×10-18 | 2.99×10-15 | 21 |
GO:0019684 | 光合作用、光反应Photosynthesis, light reaction | BP | 2.35×10-18 | 2.99×10-15 | 21 |
GO:0006118 | 电子传递Electron transport | BP | 1.25×10-11 | 1.06×10-8 | 180 |
GO:0015979 | 光合作用Photosynthesis | BP | 6.73×10-11 | 4.28×10-8 | 51 |
GO:0016491 | 氧化还原酶活性Oxidoreductase activity | MF | 1.41×10-8 | 7.18×10-6 | 181 |
GO:0008610 | 脂类生物合成的过程Lipid biosynthetic process | BP | 2.70×10-8 | 1.14×10-5 | 50 |
GO:0046483 | 杂环代谢过程Heterocycle metabolic process | BP | 3.63×10-7 | 1.32×10-4 | 101 |
GO:0009225 | 核苷酸糖代谢过程Nucleotide-sugar metabolic process | BP | 4.20×10-7 | 1.34×10-4 | 34 |
GO:0006629 | 脂质代谢过程Lipid metabolic process | BP | 8.34×10-7 | 2.35×10-4 | 71 |
GO:0003913 | DNA光裂合酶活性DNA photolyase activity | MF | 4.99×10-6 | 1.27×10-3 | 7 |
GO:0016614 | 氧化还原酶活性作用于供体CH-OH基团Oxidoreductase activity, acting on CH-OH group of donors | MF | 6.81×10-6 | 1.57×10-3 | 60 |
GO:0006720 | 类异戊二烯代谢过程Isoprenoid metabolic process | BP | 1.32×10-5 | 2.42×10-3 | 11 |
GO:0008299 | 类异戊二烯生物合成的过程Isoprenoid biosynthetic process | BP | 1.32×10-5 | 2.42×10-3 | 11 |
GO:0016616 | 氧化还原酶活性,作用于供体CH-OH基团,受体NAD或NADP Oxidoreductase activity, acting on the CH-OH group of donors, NAD or NADP as acceptor | MF | 1.34×10-5 | 2.42×10-3 | 57 |
GO:0008462 | 内肽酶Clp活性Endopeptidase clp activity | MF | 1.51×10-5 | 2.55×10-3 | 8 |
GO:0018130 | 杂环生物合成的过程Heterocycle biosynthetic process | BP | 2.38×10-5 | 3.67×10-3 | 81 |
GO:0044271 | 细胞氮化合物生物合成过程Cellular nitrogen compound biosynthetic process | BP | 2.46×10-5 | 3.67×10-3 | 106 |
GO:0003854 | 3-β-羟基-δ5-类固醇脱氢酶活性3-Beta-hydroxy-delta5-steroid dehydrogenase activity | MF | 3.86×10-5 | 4.86×10-3 | 23 |
GO:0016229 | 类固醇脱氢酶活性Steroid dehydrogenase activity | MF | 3.86×10-5 | 4.86×10-3 | 23 |
GO:0033764 | 类固醇脱氢酶活性,作用于供体CH-OH基团,受体NAD或NADP Steroid dehydrogenase activity, acting on the CH-OH group of donors, NAD or NADP as acceptor | MF | 3.86×10-5 | 4.86×10-3 | 23 |
GO:0008152 | 代谢过程Metabolic process | BP | 4.01×10-5 | 4.86×10-3 | 883 |
GO:0006694 | 类固醇生物合成的过程Steroid biosynthetic process | BP | 6.52×10-5 | 7.41×10-3 | 23 |
GO:0016667 | 氧化还原酶活性,作用于硫族供体Oxidoreductase activity, acting on a sulfur group of donors | MF | 6.71×10-5 | 7.41×10-3 | 26 |
GO:0015036 | 二硫氧化还原酶活性Disulfide oxidoreductase activity | MF | 9.85×10-5 | 1.04×10-2 | 25 |
GO:0016564 | 转录抑制因子活性Transcription repressor activity | MF | 1.18×10-4 | 1.12×10-2 | 18 |
GO:0006091 | 前体代谢产物和能量的产生Generation of precursor metabolites and energy | BP | 1.18×10-4 | 1.12×10-2 | 42 |
GO:0006772 | 硫胺素代谢过程Thiamine metabolic process | BP | 1.36×10-4 | 1.12×10-2 | 23 |
GO:0009228 | 硫胺生物合成的过程Thiamine biosynthetic process | BP | 1.36×10-4 | 1.12×10-2 | 23 |
GO:0042723 | 含硫胺的复合代谢过程Thiamine-containing compound metabolic process | BP | 1.36×10-4 | 1.12×10-2 | 23 |
GO:0042724 | 含硫胺复合生物合成工艺Thiamine-containing compound biosynthetic process | BP | 1.36×10-4 | 1.12×10-2 | 23 |
GO:0044281 | 小分子代谢过程Small molecule metabolic process | BP | 1.37×10-4 | 1.12×10-2 | 161 |
GO:0044255 | 细胞脂质代谢过程Cellular lipid metabolic process | BP | 1.41×10-4 | 1.12×10-2 | 30 |
GO:0006808 | 氮素利用调节Regulation of nitrogen utilization | BP | 1.50×10-4 | 1.12×10-2 | 18 |
GO:0019740 | 氮利用率Nitrogen utilization | BP | 1.50×10-4 | 1.12×10-2 | 18 |
GO:0055086 | 含核酸酶的小分子代谢过程Nucleobase-containing small molecule metabolic process | BP | 1.59×10-4 | 1.15×10-2 | 45 |
GO:0004497 | 单氧酶活动Monooxygenase activity | MF | 1.65×10-4 | 1.16×10-2 | 26 |
GO:0005489 | 电子运输活动Electron transporter activity | MF | 1.85×10-4 | 1.27×10-2 | 38 |
GO:0016829 | 裂合酶活性Lyase activity | MF | 2.97×10-4 | 1.97×10-2 | 53 |
GO:0009522 | 光系统I Photosystem I | CC | 3.02×10-4 | 1.97×10-2 | 9 |
GO:0006766 | 维生素代谢过程Vitamin metabolic process | BP | 3.37×10-4 | 1.99×10-2 | 35 |
GO:0006767 | 水溶性维生素代谢过程Water-soluble vitamin metabolic process | BP | 3.37×10-4 | 1.99×10-2 | 35 |
GO:0009110 | 维生素生物合成的过程Vitamin biosynthetic process | BP | 3.37×10-4 | 1.99×10-2 | 35 |
GO:0042364 | 水溶性维生素生物合成过程Water-soluble vitamin biosynthetic process | BP | 3.37×10-4 | 1.99×10-2 | 35 |
GO:0044272 | 硫化合物生物合成工艺Sulfur compound biosynthetic process | BP | 3.80×10-4 | 2.20×10-2 | 26 |
GO:0019438 | 芳香族化合物生物合成工艺Aromatic compound biosynthetic process | BP | 3.90×10-4 | 2.20×10-2 | 34 |
GO:0003824 | 催化活性Catalytic activity | MF | 4.54×10-4 | 2.51×10-2 | 767 |
GO:0044283 | 小分子生物合成过程Small molecule biosynthetic process | BP | 5.38×10-4 | 2.91×10-2 | 71 |
GO:0006790 | 硫化合物代谢过程Sulfur compound metabolic process | BP | 5.81×10-4 | 3.04×10-2 | 31 |
GO:0008202 | 类固醇代谢过程Steroid metabolic process | BP | 5.86×10-4 | 3.04×10-2 | 23 |
GO:0008831 | dTDP-4-脱氢罗汉糖还原酶活性dTDP-4-dehydrorhamnose reductase activity | MF | 6.83×10-4 | 3.34×10-2 | 14 |
GO:0045226 | 胞外多糖生物合成过程Extracellular polysaccharide biosynthetic process | BP | 6.83×10-4 | 3.34×10-2 | 14 |
GO:0046379 | 胞外多糖代谢过程Extracellular polysaccharide metabolic process | BP | 6.83×10-4 | 3.34×10-2 | 14 |
GO:0051186 | 辅因子代谢过程Cofactor metabolic process | BP | 9.29×10-4 | 4.45×10-2 | 56 |
表5 巨菌草幼叶、根下调差异基因GO富集结果
Table 5 Giant Juncao GO enrichment of leaves and roots downregulated DEGs (FDR<0.05)
GO登录号GO_accession | 功能 Function | 种类 Ontology | P值 P-value | 错误发现率FDR | 差异基因个数DEGs number (No.) |
---|---|---|---|---|---|
GO:0009765 | 光合作用、光收获Photosynthesis, light harvesting | BP | 2.35×10-18 | 2.99×10-15 | 21 |
GO:0019684 | 光合作用、光反应Photosynthesis, light reaction | BP | 2.35×10-18 | 2.99×10-15 | 21 |
GO:0006118 | 电子传递Electron transport | BP | 1.25×10-11 | 1.06×10-8 | 180 |
GO:0015979 | 光合作用Photosynthesis | BP | 6.73×10-11 | 4.28×10-8 | 51 |
GO:0016491 | 氧化还原酶活性Oxidoreductase activity | MF | 1.41×10-8 | 7.18×10-6 | 181 |
GO:0008610 | 脂类生物合成的过程Lipid biosynthetic process | BP | 2.70×10-8 | 1.14×10-5 | 50 |
GO:0046483 | 杂环代谢过程Heterocycle metabolic process | BP | 3.63×10-7 | 1.32×10-4 | 101 |
GO:0009225 | 核苷酸糖代谢过程Nucleotide-sugar metabolic process | BP | 4.20×10-7 | 1.34×10-4 | 34 |
GO:0006629 | 脂质代谢过程Lipid metabolic process | BP | 8.34×10-7 | 2.35×10-4 | 71 |
GO:0003913 | DNA光裂合酶活性DNA photolyase activity | MF | 4.99×10-6 | 1.27×10-3 | 7 |
GO:0016614 | 氧化还原酶活性作用于供体CH-OH基团Oxidoreductase activity, acting on CH-OH group of donors | MF | 6.81×10-6 | 1.57×10-3 | 60 |
GO:0006720 | 类异戊二烯代谢过程Isoprenoid metabolic process | BP | 1.32×10-5 | 2.42×10-3 | 11 |
GO:0008299 | 类异戊二烯生物合成的过程Isoprenoid biosynthetic process | BP | 1.32×10-5 | 2.42×10-3 | 11 |
GO:0016616 | 氧化还原酶活性,作用于供体CH-OH基团,受体NAD或NADP Oxidoreductase activity, acting on the CH-OH group of donors, NAD or NADP as acceptor | MF | 1.34×10-5 | 2.42×10-3 | 57 |
GO:0008462 | 内肽酶Clp活性Endopeptidase clp activity | MF | 1.51×10-5 | 2.55×10-3 | 8 |
GO:0018130 | 杂环生物合成的过程Heterocycle biosynthetic process | BP | 2.38×10-5 | 3.67×10-3 | 81 |
GO:0044271 | 细胞氮化合物生物合成过程Cellular nitrogen compound biosynthetic process | BP | 2.46×10-5 | 3.67×10-3 | 106 |
GO:0003854 | 3-β-羟基-δ5-类固醇脱氢酶活性3-Beta-hydroxy-delta5-steroid dehydrogenase activity | MF | 3.86×10-5 | 4.86×10-3 | 23 |
GO:0016229 | 类固醇脱氢酶活性Steroid dehydrogenase activity | MF | 3.86×10-5 | 4.86×10-3 | 23 |
GO:0033764 | 类固醇脱氢酶活性,作用于供体CH-OH基团,受体NAD或NADP Steroid dehydrogenase activity, acting on the CH-OH group of donors, NAD or NADP as acceptor | MF | 3.86×10-5 | 4.86×10-3 | 23 |
GO:0008152 | 代谢过程Metabolic process | BP | 4.01×10-5 | 4.86×10-3 | 883 |
GO:0006694 | 类固醇生物合成的过程Steroid biosynthetic process | BP | 6.52×10-5 | 7.41×10-3 | 23 |
GO:0016667 | 氧化还原酶活性,作用于硫族供体Oxidoreductase activity, acting on a sulfur group of donors | MF | 6.71×10-5 | 7.41×10-3 | 26 |
GO:0015036 | 二硫氧化还原酶活性Disulfide oxidoreductase activity | MF | 9.85×10-5 | 1.04×10-2 | 25 |
GO:0016564 | 转录抑制因子活性Transcription repressor activity | MF | 1.18×10-4 | 1.12×10-2 | 18 |
GO:0006091 | 前体代谢产物和能量的产生Generation of precursor metabolites and energy | BP | 1.18×10-4 | 1.12×10-2 | 42 |
GO:0006772 | 硫胺素代谢过程Thiamine metabolic process | BP | 1.36×10-4 | 1.12×10-2 | 23 |
GO:0009228 | 硫胺生物合成的过程Thiamine biosynthetic process | BP | 1.36×10-4 | 1.12×10-2 | 23 |
GO:0042723 | 含硫胺的复合代谢过程Thiamine-containing compound metabolic process | BP | 1.36×10-4 | 1.12×10-2 | 23 |
GO:0042724 | 含硫胺复合生物合成工艺Thiamine-containing compound biosynthetic process | BP | 1.36×10-4 | 1.12×10-2 | 23 |
GO:0044281 | 小分子代谢过程Small molecule metabolic process | BP | 1.37×10-4 | 1.12×10-2 | 161 |
GO:0044255 | 细胞脂质代谢过程Cellular lipid metabolic process | BP | 1.41×10-4 | 1.12×10-2 | 30 |
GO:0006808 | 氮素利用调节Regulation of nitrogen utilization | BP | 1.50×10-4 | 1.12×10-2 | 18 |
GO:0019740 | 氮利用率Nitrogen utilization | BP | 1.50×10-4 | 1.12×10-2 | 18 |
GO:0055086 | 含核酸酶的小分子代谢过程Nucleobase-containing small molecule metabolic process | BP | 1.59×10-4 | 1.15×10-2 | 45 |
GO:0004497 | 单氧酶活动Monooxygenase activity | MF | 1.65×10-4 | 1.16×10-2 | 26 |
GO:0005489 | 电子运输活动Electron transporter activity | MF | 1.85×10-4 | 1.27×10-2 | 38 |
GO:0016829 | 裂合酶活性Lyase activity | MF | 2.97×10-4 | 1.97×10-2 | 53 |
GO:0009522 | 光系统I Photosystem I | CC | 3.02×10-4 | 1.97×10-2 | 9 |
GO:0006766 | 维生素代谢过程Vitamin metabolic process | BP | 3.37×10-4 | 1.99×10-2 | 35 |
GO:0006767 | 水溶性维生素代谢过程Water-soluble vitamin metabolic process | BP | 3.37×10-4 | 1.99×10-2 | 35 |
GO:0009110 | 维生素生物合成的过程Vitamin biosynthetic process | BP | 3.37×10-4 | 1.99×10-2 | 35 |
GO:0042364 | 水溶性维生素生物合成过程Water-soluble vitamin biosynthetic process | BP | 3.37×10-4 | 1.99×10-2 | 35 |
GO:0044272 | 硫化合物生物合成工艺Sulfur compound biosynthetic process | BP | 3.80×10-4 | 2.20×10-2 | 26 |
GO:0019438 | 芳香族化合物生物合成工艺Aromatic compound biosynthetic process | BP | 3.90×10-4 | 2.20×10-2 | 34 |
GO:0003824 | 催化活性Catalytic activity | MF | 4.54×10-4 | 2.51×10-2 | 767 |
GO:0044283 | 小分子生物合成过程Small molecule biosynthetic process | BP | 5.38×10-4 | 2.91×10-2 | 71 |
GO:0006790 | 硫化合物代谢过程Sulfur compound metabolic process | BP | 5.81×10-4 | 3.04×10-2 | 31 |
GO:0008202 | 类固醇代谢过程Steroid metabolic process | BP | 5.86×10-4 | 3.04×10-2 | 23 |
GO:0008831 | dTDP-4-脱氢罗汉糖还原酶活性dTDP-4-dehydrorhamnose reductase activity | MF | 6.83×10-4 | 3.34×10-2 | 14 |
GO:0045226 | 胞外多糖生物合成过程Extracellular polysaccharide biosynthetic process | BP | 6.83×10-4 | 3.34×10-2 | 14 |
GO:0046379 | 胞外多糖代谢过程Extracellular polysaccharide metabolic process | BP | 6.83×10-4 | 3.34×10-2 | 14 |
GO:0051186 | 辅因子代谢过程Cofactor metabolic process | BP | 9.29×10-4 | 4.45×10-2 | 56 |
通路 Pathway ID | 通路项目 Pathway term | P值 P-value | 错误发现率 FDR | 差异基因个数 DEGs Number (No.) |
---|---|---|---|---|
sita00710 | 光合生物的固碳作用Carbon fixation in photosynthetic organisms | 1.73×10-6 | 2.07×10-4 | 52 |
sita00906 | 类胡萝卜素生物合成Carotenoid biosynthesis | 1.81×10-4 | 1.07×10-2 | 25 |
sita01110 | 次生代谢物的生物合成Biosynthesis of secondary metabolites | 3.53×10-4 | 1.07×10-2 | 396 |
sita00030 | 磷酸戊糖途径Pentose phosphate pathway | 3.56×10-4 | 1.07×10-2 | 35 |
sita01200 | 碳代谢Carbon metabolism | 5.24×10-4 | 1.26×10-2 | 109 |
表6 巨菌草幼叶、根差异基因KEGG富集分析
Table 6 Giant Juncao KEGG enrichment of leaves and roots DEGs (FDR?0.05)
通路 Pathway ID | 通路项目 Pathway term | P值 P-value | 错误发现率 FDR | 差异基因个数 DEGs Number (No.) |
---|---|---|---|---|
sita00710 | 光合生物的固碳作用Carbon fixation in photosynthetic organisms | 1.73×10-6 | 2.07×10-4 | 52 |
sita00906 | 类胡萝卜素生物合成Carotenoid biosynthesis | 1.81×10-4 | 1.07×10-2 | 25 |
sita01110 | 次生代谢物的生物合成Biosynthesis of secondary metabolites | 3.53×10-4 | 1.07×10-2 | 396 |
sita00030 | 磷酸戊糖途径Pentose phosphate pathway | 3.56×10-4 | 1.07×10-2 | 35 |
sita01200 | 碳代谢Carbon metabolism | 5.24×10-4 | 1.26×10-2 | 109 |
1 | Zhou J, Lin X S, Lin H, et al. Advances on Juncao research and application. Journal of Fujian Agriculture and Forestry University (Natural Science Edition), 2020, 49(2): 145-152. |
周晶, 林兴生, 林辉, 等. 菌草研究与应用进展. 福建农林大学学报(自然科学版), 2020, 49(2): 145-152. | |
2 | Liu F S, Lin D M, Lin H, et al. Physiological and photosynthetic responses of Giant Juncao (Pennisetum Giganteum) to drought stress. Fresenius Environmental Bulletin, 2017, 26(6): 3868-3871. |
3 | Zheng H K, Lin X J, Lin H, et al. Research progress on JUJUNCAO (Pennisetum giganteum). Journal of Fujian Agriculture and Forestry University (Natural Science Edition), 2019, 48(6): 681-687. |
郑华坤, 林雄杰, 林辉, 等. 巨菌草 (Pennisetum giganteum)研究进展. 福建农林大学学报 (自然科学版), 2019, 48(6): 681-687. | |
4 | Martin L B B, Fei Z J, Giovannoni J J. Catalyzing plant science research with RNA-seq. Frontiers in Plant Science, 2013, 4: 66. |
5 | Zhang H , He L, Cai L. Transcriptome sequencing: RNA-Seq. Methods in Molecular Biology, 2018, 1754: 15-27. |
6 | Camilios-Neto D, Bonato P, Wassem R, et al. Dual RNA-seq transcriptional analysis of wheat roots colonized by Azospirillum brasilense reveals up-regulation of nutrient acquisition and cell cycle genes. BMC Genomics, 2014, 15: 378. |
7 | Zhang L M, Liu X G, Qu X N, et al. Early transcriptomic adaptation to Na2CO3 stress altered the expression of a quarter of the total genes in the maize genome and exhibited shared and distinctive profiles with NaCl and high pH stresses. Journal of Integrative Plant Biology, 2013, 55(11): 1147-1165. |
8 | Zhou Y, Yang P, Cui F L, et al. Transcriptome analysis of salt stress responsiveness in the seedlings of Dongxiang wild rice (Oryza rufipogon Griff.). PLoS One, 2016, 11(1): e0146242. |
9 | Fracasso A, Trindade L M, Amaducci S, et al. Drought stress tolerance strategies revealed by RNA-Seq in two sorghum genotypes with contrasting WUE. BMC Plant Biology, 2016, 16(1): 115. |
10 | Wang Y Q, Li L,Tang S, et al. Combined small RNA and degradome sequencing to identify miRNAs and their targets in response to drought in foxtail millet. BMC Genetics, 2016, 17(1): 57. |
11 | Liu W X, Zhang Z S, Chen S Y, et al. Global transcriptome profiling analysis reveals insight into saliva-responsive genes in alfalfa. Plant Cell Reports, 2016, 35(3): 561-571. |
12 | Shamila W A, Stephen B, Istvan N, et al. Changes in Lolium perenne transcriptome during cold acclimation in two genotypes adapted to different climatic conditions. BMC Plant Biology, 2015, 15(1): 250-264. |
13 | Ma Q, Bao A K, Chai W W, et al. Transcriptomic analysis of the succulent xerophyte Zygophyllum xanthoxylum in response to salt treatment and osmotic stress. Plant Soil, 2016, 402(1/2): 343-361. |
14 | Sun Y P, Wang F W, Wang N, et al. Transcriptome exploration in Leymus chinensis under saline-alkaline treatment using 454 pyrosequencing. PLoS One, 2013, 8(1): 1-12. |
15 | Leng N, Liu X W, Zhang N, et al. Differential gene analysis of Poa pratensis in response to drought stress. Acta Prataculturae Sinica, 2017, 26(12): 128-137. |
冷暖, 刘晓巍, 张娜, 等. 草地早熟禾干旱胁迫转录组差异性分析. 草业学报, 2017, 26(12): 128-137. | |
16 | Manfred G G, Brian J H, Moran Y, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nature Biotechnology, 2011, 7(29): 644-652. |
17 | Altschul S F, Gish W, Miller W, et al. Basic local alignment search tool. Journal of Molecular Biology, 1990, 3(215): 403-410. |
18 | Kanehisa M, Araki M, Goto S, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Research, 2008, 36(database issue): 480-484. |
19 | Young M D, Wakefield M J, Smyth G K, et al. Method gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biology, 2010, 11: R14. |
20 | Zhao J B, Hou X Y, Wu Z N, et al. Transcriptome analysis of Leymus chinensis under different mowing intensities. Acta Prataculturae Sinica, 2018, 27(2): 105-116. |
赵劲博, 侯向阳, 武自念, 等. 不同刈割强度下羊草转录组研究. 草业学报, 2018, 27(2): 105-116. | |
21 | Qian C, Liu Z W, Zhong X X, et al. Transcriptomic analysis of the self-incompatibility mechansim in Paspalum vaginatum by comparison with an artificial self-compatible mutant. Acta Prataculturae Sinica, 2019, 28(5): 134-144. |
钱晨, 刘智微, 钟小仙, 等. 海滨雀稗自交结实突变体及野生型幼穗组织的转录组分析. 草业学报, 2019, 28(5): 134-144. | |
22 | Chen W, Yao Q, Patil G B, et al. Identification and comparative analysis of differential gene expression in soybean leaf tissue under drought and flooding stress revealed by RNA-Seq. Frontiers in Plant Science, 2016, 7: 1044. |
23 | Jin J, Zhang H, Kongn L, et al. Plant TFDB 3.0: A portal for the functional and evolutionary study of plant transcription factors. Nucleic Acids Research, 2013(D1): 1182-1187. |
24 | Lorenzo C P, Anahit G, Irma R V, et al. Genome-wide classification and evolutionary analysis of the bHLH family of transcription factors in Arabidopsis, poplar, rice, moss, and algae. Plant Physiology, 2010, 153(3): 1398-1412. |
25 | Yang L H, Ying H, Huang Q M, et al. An analysis of the response of the LpWRKY20 gene to abiotic stress and its role in drought resistance. Acta Prataculturae Sinica, 2020, 29(1): 193-202. |
杨柳慧, 尹航, 黄沁梅, 等. 细叶百合LpWRKY20基因对非生物胁迫的响应及抗旱性分析. 草业学报, 2020, 29(1): 193-202. | |
26 | Wei K F, Chen H Q. Comparative functional genomics analysis of bHLH gene family in rice, maize and wheat. BMC Plant Biology, 2018, 18(1): 309. |
27 | Bailey P C, Martin C, Toledo-Ortiz G, et al. Update on the basic helix-loop-helix transcription factor gene family in Arabidopsis thaliana. The Plant Cell, 2003, 15(11): 2497-2501. |
28 | Song X M, Huang Z N, Duan W K, et al. Genome-wide analysis of the bHLH transcription factor family in Chinese cabbage (Brassica rapa ssp.pekinensis). Molecular Genetics and Genomics, 2014, 289(1): 77-91. |
29 | Niu X, Guan Y X, Chen S K, et al. Genome-wide analysis of basic helix-loop-helix (bHLH) transcription factors in Brachypodium distachyon. BMC Genomics, 2017, 18(1): 619. |
30 | Mao K, Dong Q L, Li C, et al. Genome wide identification and characterization of apple bHLH transcription factors and expression analysis in response to drought and salt stress. Frontiers in Plant Science, 2017, 8: 480. |
31 | Ramamoorthy R, Jiang S Y, Kumar N, et al. A comprehensive transcriptional profiling of the WRKY gene family in rice under various abiotic and phytohormone treatments. Plant & Cell Physiology, 2008, 49(6): 865-879. |
32 | Pandey S P, Roccaro M, Schon M, et al. Transcriptional reprogramming regulated by WRKY18 and WRKY40 facilitates powdery mildew infection of Arabidopsis. The Plant Journal, 2010, 64(6): 912-923. |
[1] | 贾雨雷, 廖真, 汪丽芳, 卜建超, 林标声, 林辉, 苏德伟, 鲁国东, 林占熺. 化肥减量配施菌草固氮菌肥对巨菌草生长、营养品质及土壤养分的影响[J]. 草业学报, 2021, 30(3): 215-223. |
[2] | 舒新月, 江波, 马丽, 郑爱萍. 不同侵染时间点稻粒黑粉病菌的转录组分析[J]. 草业学报, 2020, 29(9): 190-202. |
[3] | 滕珂, 张蕊, 檀鹏辉, 岳跃森, 范希峰, 武菊英. 日本结缕草ZjERF1的克隆、转录激活活性、亚细胞定位及表达分析[J]. 草业学报, 2019, 28(6): 56-65. |
[4] | 钱晨, 刘智微, 钟小仙, 吴娟子, 张建丽, 潘玉梅. 海滨雀稗自交结实突变体及野生型幼穗组织的转录组分析[J]. 草业学报, 2019, 28(5): 132-142. |
[5] | 闵学阳, 韦兴燚, 刘文献, 张正社, 金小煜, NDAYAMBAZABoniface, 吴洪林, 李昱, 王彦荣. 箭筈豌豆品种间遗传差异的SSR分析及指纹图谱构建[J]. 草业学报, 2019, 28(4): 116-128. |
[6] | 姜红岩, 滕珂, 檀鹏辉, 尹淑霞. 日本结缕草ZjZFN1基因对拟南芥的转化及其耐旱性分析[J]. 草业学报, 2019, 28(4): 129-138. |
[7] | 宫文龙, 王赞, 赵桂琴, 马琳, 韦宝, 龚攀, 刘希强. 沙打旺EST-SSR分子标记开发及其遗传多样性分析[J]. 草业学报, 2019, 28(11): 147-158. |
[8] | 成启明, 格根图, 撒多文, 王志军, 范文强, 卜振鲲, 司强, 李俊峰, 卢娟, 贾玉山. 不同品种紫花苜蓿转录组分析及营养品质差异的探讨[J]. 草业学报, 2019, 28(10): 199-208. |
[9] | 李孟湛, 尹红菊, 李丁丁, 刘亚琪, 王锁民. 利用CRISPR/Cas9技术敲除拟南芥转录因子MYB 40的两种可变剪接体[J]. 草业学报, 2019, 28(1): 120-127. |
[10] | 吴娟子, 钱晨, 刘智微, 潘玉梅, 钟小仙. 基于转录组测序分析象草木质素合成的研究[J]. 草业学报, 2019, 28(1): 150-161. |
[11] | 柯丹霞, 彭昆鹏, 夏远君, 朱玉莹, 张丹丹. 盐胁迫应答基因GmWRKY6的克隆及转基因百脉根的抗盐分析[J]. 草业学报, 2018, 27(8): 95-106. |
[12] | 宋吉轩, 吕俊, 宗学凤, 何秀娟, 徐宇, 吴潇, 王三根. 干旱胁迫下ALA对羊草生长、内源激素及转录组的影响[J]. 草业学报, 2018, 27(7): 73-83. |
[13] | 朱永群, 彭丹丹, 林超文, 聂刚, 许文志, 黄琳凯, 罗付香, 彭建华, 张新全. 苏丹草转录组SSR分子标记开发及遗传多样性评价[J]. 草业学报, 2018, 27(5): 178-189. |
[14] | 庄黎丽, 王剑, 杨志民. 基于转录组数据库的高羊茅HD-Zip I转录因子的鉴定及表达模式解析[J]. 草业学报, 2018, 27(3): 67-77. |
[15] | 赵劲博, 侯向阳, 武自念, 任卫波, 胡宁宁, 郭丰辉, 马文静. 不同刈割强度下羊草转录组研究[J]. 草业学报, 2018, 27(2): 105-116. |
阅读次数 | ||||||
全文 |
|
|||||
摘要 |
|
|||||