灰树花菌丝体与原基转录组差异表达分析

• 花匠小妙招
• 2024-12-25 06:51

灰树花(Grifola frondosa)又名栗子蘑、莲花菌、贝叶多孔菌, 在分类系统中属于担子菌亚门、层菌纲、非褶菌目、多孔菌科[1]。其子实体肉质柔嫩、口味鲜美, 富含氨基酸、多种维生素、微量元素、多糖与生物活性物质, 是近年来开发的珍奇药食两用菌之一[2]。研究表明, 从灰树花中提取的多糖具有抗肿瘤、抗氧化、抗HIV 病毒、增强免疫功能等多种生物活性[3]。

灰树花原基的形成是其生活史中重要的生理过程, 对灰树花的品质及出菇有重要影响[4]。原基的形成不仅需要一定的环境条件和营养条件, 更是由一系列基因协调作用的结果[5]。研究灰树花原基形成时的差异表达基因, 对阐明灰树花子实体形成的分子机理具有重要的意义。

通过对灰树花菌丝体、原基2个不同生长发育期转录组的比较研究, 筛选出2个时期差异表达的基因, 并对这些差异表达的基因进行生物信息学分析, 旨在发现控制灰树花原基形成和子实体生长的功能基因, 为通过基因工程培育优质、高产的灰树花新品种提供理论依据。

1 材料与方法

1.1 试材及取样

供试菌株为灰树花小黑汀, 引自山东泰安。菌丝长满菌包后移入出菇大棚培养7~10 d, 待菌丝扭结形成原基后, 收集菌丝体与原基样品。

1.2 RNA提取

用Trizol法分别提取菌丝体与原基总RNA, 并用RNAeasy plant mini kit对提取的总RNA进行纯化, 70 ℃变性2 min后, NanoDrop ND-2000检测其浓度、琼脂糖凝胶电泳分析RNA的完整性[6]。检测合格的RNA用于mRNA的富集及cDNA的合成。

1.3 测序

用Ultra RNAlibrary prep kit for illumina进行文库构建, 纯化后Agilent high sensitivity DNA kit检测文库插入片段大小, 定量后Illumina HiSeq TM 2000对建好的测序文库进行测序[7]。

1.4 参考基因组比对

使用Trinity(版本r20131110, 默认参数)对RNA-seq的原始reads数据进行拼接, 并使用测序序列(reads)比对软件Tophat v2.0.5将clean reads分别比对到NR数据库[8]。后续分析都基于上述比对结果, 然后统计clean reads比对到NR数据库的分布情况。

1.5 基因表达量分析

基因表达量的计算使用FPKM数值[9], FPKM代表每百万个比对上的reads中每kb外显子的片段数, 用于衡量基因的表达水平, 判断差异表达显著的方式为满足条件:表达量变化倍数>2且FPKM差值>20, 或表达量变化倍数>4且FPKM差值>2。

1.6 Go功能富集分析

对不同样本中差异表达的基因进行基因本体(gene ontology, GO)分析能确定差异表达基因的功能[10]。卡方检验P值≤0.05 条件的GO 条目被认为是显著富集的GO条目。

1.7 Pathway 显著性富集分析

使用http://www.genome.jp/tools/kaas/对灰树花原基和菌丝体差异表达基因进行Pathway显著性富集分析。卡方检验P值≤0.05的Pathway确定为不同样本间差异表达的基因所参与的最主要代谢途径和信号转导途径。

2 结果与分析

2.1 RNA提取

灰树花菌丝体和原基RNA经过NanoDrop定量后, 样本浓度(表1)完整性及28S∶18S(图1)均符合转录组测序质量要求, 可用于开展下一步试验。

2.2 两个样本的转录本分析

通过高通量测序, 在菌丝体和原基样本库中分别得到了35913116个和32873980个原始测序标签(raw reads)。原始数据去除低质量的标签后, 在菌丝体和原基这2个文库得到高质量测序标签(clean reads)总数分别为35788532和32755254。将所有高质量测序标签与参考基因序列进行比对, 在菌丝体和原基文库中, 在参考序列上有唯一比对位置的测序序列分别为22085081和21547780, 占高质量标签种数的61.7%和65.78%。整段比对到外显子的测序序列总数为22061957和21537230(表2)。

表1   灰树花菌丝体、原基RNA浓度

Tab.1   Concentration of RNA from mycelium and primordium of Grifola frondosa

样品
Sampels浓度
Concentration单位
UnitOD260/280OD260/2301-JST266.4ng·μL-12.141.922-YJ697.8ng·μL-12.162.21

Note: 1-JST mean mycelium, 2-YJ mean primordium.注1-JST表示菌丝体, 2-YJ表示原基。

新窗口打开

图1   灰树花菌丝体、原基RNA电泳检测1-JST表示菌丝体, 2-YJ表示原基

Fig.1   Electropherogram of RNAs from mycelium and primordium of Grifola frondosa1-JST mean mycelium, 2-YJ mean primordium

2.3 两样本间差异表达基因分析

菌丝体共有41858个表达基因, 原基中有34255个表达基因, 最终在菌丝体和原基数据库中, 共发现显著差异表达的基因4094个, 其中在原基中上调1886个, 下调2208个, 只在菌丝阶段表达的基因915个, 只在原基中表达的基因284个。

为进一步了解菌丝体阶段与原基阶段差异表达基因的情况, 对两者具有明显差异表达的4094个基因进行分析, 对同源比对(blastnr)得到的有确切功能的基因列于表3(明显差异表达前20个基因)。从表中可知, 在菌丝体及原基阶段差异表达的基因主要与酸性蛋白酶类、凝集素、细胞色素、NADPH-P450还原酶、酯酶、胺氧化酶、克拉维胺合成酶、糖苷水解酶家族相关。其中凝集素能凝集细胞并参与糖结合作用; 酯酶协助参与降解纤维素、半纤维素和木质素; 胺氧化酶参与分解氨基酸; 糖苷水解酶能作用各种糖苷或寡糖使糖苷键水解, 参与糖的代谢。上述基因表达量的变化可能在原基的形成中起到一定的作用。

表2   两个样本中标签分布一览表

Tab.2   Distributions of tags from two libraries

数据产出及比对
Data output and comparison菌丝体
Mycelium原基
Primordium原始序列标签
Rawreads3591311632873980高质量测序标签
Cleanreads3578853232755254未比对上的标签
Unmapped reads90281397574132非唯一比对位置标签
Mapq<mapq_cut(non-unique)46658943623751唯一比对位置标签
Mapq>=apq_cut(unique)2208508121547780第一端测序标签
Read-11101004010756337第二端测序标签
Read-21107504110791443比对到基因组上正列
Reads map to‘+’1098187310704107比对到基因组上负列
Reads map to‘-’1110320810843673整段比对到外显子的序列
Non-splice reads2206195721537230分段比对到两个外显子上的序列
Splice reads2312410550双端测序序列定位到基因组上相对
的距离符合测序片段的长度分布
Eads mapped in proper pairs1959784019326088

新窗口打开

2.4 原基特异性表达基因分析

为了进一步研究原基阶段特异基因的表达情况, 对原基中特异表达的284个基因进行了同源比对(blastnr)。284个基因中在NCBI 中比对到的同源基因绝大多数是功能未知的假定蛋白(hypothetical protein)或预测蛋白(predicted protein)。比对得到有确切功能的基因与糖代谢、脂类代谢、核酸代谢及细胞膜、叶绿体膜有关(表4)。其中, E1脱氢酶及转酮醇酶在戊糖磷酸循环中起着重要作用, HAD水解酶在核酸代谢过程中起重要作用。研究该类基因可能有助于揭示原基的形成。

2.5 GO功能富集分析

在富集到的前20个具有明显差异的Go Term中, 有13个显著富集的Go Term富集在参与分子功能(molecular function)本体上, 有5个显著富集的Go Term富集在参与生物过程(biological process)本体上, 有2个显著富集的Go Term富集在参与细胞组成(cellular component)本体上。其中, 在细胞组成本体中, 与线粒体膜成分相关的基因均上调表达; 在生物过程本体中, 与谷氨酰胺代谢、脂肪酸生物合成相关的基因也均上调表达(表5)。

表3   菌丝体、原基中差异表达基因的同源比对

Tab.3   Blastnr of differently expressed genes between mycelium and primordium libraries

基因ID
Gene ID菌丝中表达值
Expression value of
mycelium原基中表达值
Expression value of
primordiumQ值
Q value同源基因及其功能
Homologous gene and its functionunigene23873483.46.00.023酸性蛋白酶(gi|392558359|)(gb|EIW51548.1|)
Acid proteaseunigene1392911.65860.50.037凝集素(gi|18478668|)(gb|AAL73235.1|)unigene30943186.41.40.037DUF1479结构域蛋白(gi|392564554|)(gb|EIW57732.1|)
DUF1479 domain containing proteinunigene122108.22412.50.042凝集素(gi|50980806|)(gb|AAT91249.1|)
lectinunigene418732.381.10.042细胞色素P450(gi|395327404|)(gb|EJF59804.1|)
Cytochrome P450unigene459371.767.00.042NADPH-P450还原酶(gi|392559250|)(gb|EIW52435.1|)
NADPH-P450 reductaseunigene27346223.55.90.044IlvN结构域蛋白(gi|395334382|)(gb|EJF66758.1|)
IlvN domain containing proteinunigene283233.5162.90.044酯酶(gi|395328576|)(gb|EJF60967.1|)
esteraseunigene418742.168.90.044细胞色素P450(gi|395327404|)(gb|EJF59804.1|)
Cytochrome P450unigene459362.273.60.044NADPH-P450还原酶(gi|392559250|)(gb|EIW52435.1|)
NADPH-P450 reductaseunigene418753.0170.30.048细胞色素P450(gi|395327404|)(gb|EJF59804.1|)
Cytochrome P450unigene146353.5118.40.052真菌疏水蛋白结构(gi|390605204|)(gb|EIN14595.1|)
Fungal hydrophobinunigene47574248.01.30.056胺氧化酶(gi|392563541|)(gb|EIW56720.1|)
Amine oxidaseunigene1987022.8653.10.058真菌疏水蛋白结构(gi|395323171|)(gb|EJF55662.1|)
Fungal hydrophobinunigene360891.834.70.058细胞色素P450(gi|390599619|)(gb|EIN09015.1|)
Cytochrome P450unigene305250.813.60.06克拉维胺合成酶(gi|392569544|)(gb|EIW62717.1|)
Clavaminate synthase-like proteinunigene360881.527.10.06细胞色素P450(gi|390599619|)(gb|EIN09015.1|)
Cytochrome P450unigene3962714.61.00.061类CDF转运蛋白(gi|395334914|)(gb|EJF67290.1|)
CDF-like metal transporterunigene305230.514.70.063克拉维胺合成酶(gi|392569544|)(gb|EIW62717.1|)
Clavaminate synthase-like proteinunigene298469.60.30.065糖苷水解酶家族(gi|392564918|)(gb|EIW58095.1|)
Glycoside hydrolase family protein

新窗口打开

2.6 Pathway 显著性富集分析

菌丝体和原基中Pathway 显著性富集分析见表6, 在富集到的前11个具有明显差异的pathway中, 核糖体蛋白代谢途径中的全部基因均上调表达, 即核糖体数量在增多, 核糖体是蛋白质的合成场所, 核糖体数量在原基阶段增加, 表明原基形成时蛋白质合成增加。另外, 萜类生物合成途径、苯丙氨酸、酪氨酸、色氨酸生物合成途径均上调表达, 氧化磷酸化途径、花生四烯酸代谢途径、脂肪酸代谢相关途径基因基本上调表达, 表明在原基形成时期细胞代谢旺盛。

表4   原基中特异表达基因的同源比对

Tab.4   Blastnr comparison of genes only expressed in the primordium

同源基因功能
Functions of
homologous genes同源基因Homologous gene糖代谢相关
Glucose metabolismE1脱氢酶及转酮醇酶结构域蛋白(gi|395325950|)(gb|EJF58365.1|),
dehydrogenase E1 and transketolase domain-containing protein
醛脱氢酶(gi|392561496|)(gb|EIW54677.1|),
Aldehyde dehydrogenase
MFS通用底物转运蛋白(gi|392561630|)(gb|EIW54811.1|)
MFS general substrate transporter脂类代谢相关
Lipid metabolism酮脂酰合成酶结构域蛋白(gi|395326786|)(gb|EJF59192.1|)
ketoacyl-synt-domain-containing protein核酸代谢相关
Nucleotide etabolismHAD水解酶(gi|395328072|)(gb|EJF60467.1|)
HAD hydrolase细胞膜蛋白
Cell membrane proteinDUF6结构域蛋白(gi|392559753|)(gb|EIW52937.1|),
DUF6-domain-containing protein
NAD(P)结合蛋白(gi|392588861|)(gb|EIW78192.1|)
NAD(P)-binding protein叶绿体膜蛋白
Chloroplast membrane protein磷酸丙糖转运蛋白家族(emb|AJ437264.1)
phosphate translocator-like protein其他
Others细胞色素(gi|395325039|)(gb|EJF57468.1|),
cytochrome
叶绿素A-B 结合蛋白PF00504(gb|EF576336.1|),
Chlorophyll A-B binding protein PF00504
类克拉维胺合成酶(gi|395330466|)(gb|EJF62849.1|),
Clavaminate synthase-like protein
乌头酸水合酶(gi|4029334|)(emb|CAA76360.1|)
Aconitate hydratase

新窗口打开

3 结论

食用菌子实体的形成和发育, 是食用菌育种研究的主要领域[1]。子实体由原基诱导分化而来, 对菌丝体与原基形成过程中基因表达差异的研究有助于进一步了解原基形成、分化机制[12]。灰树花转录组序列不含内含子及其他非编码序列, 可识别表达丰度较低的RNA, 对灰树花菌丝体、原基进行全转录组测序, 组建并比较两者完整的基因表达谱, 有助于我们更深入研究灰树花原基形成机制, 为食用菌灰树花的分子育种奠定基础。

采用新一代Illumina高通量测序技术, 对灰树花菌丝体和原基进行了全转录组测序, 共发现两者具有显著差异表达的基因4094 个, 其中有1886个基因在原基中上调表达, 有2208个基因在原基中下调表达, 这些差异表达的基因主要涉及酸性蛋白酶类、凝集素、细胞色素、NADPH-P450还原酶、酯酶、胺氧化酶、克拉维胺合成酶、糖苷水解酶家族等。同时分析得到284个只在原基中表达的基因, 这些基因分别与糖代谢、脂类代谢、核酸代谢及细胞膜、叶绿体膜结构相关。

通过GO功能富集分析和Pathyway显著性分析, 发现灰树花原基形成时细胞代谢、蛋白质合成均升高, 这表明灰树花原基形成是一个复杂的生理过程, 在多种能量代谢的调控下合成了新蛋白质。本研究得到的灰树花菌丝体及原基形成时差异表达基因, 将为今后进一步探索其原基形成的分子机制奠定基础。

表5   菌丝体与原基GO功能富集分析

Tab.5   Significantly enriched GO terms of mycelium and primordium

基因本体条目
Gene ontology entries分类
ClassificationP值
P value上调基因数
Up-regulated genes下调基因数
Down-regulated genes结合ATP
ATP binding分子功能
Molecular function9.9E-067430酸酐水解酶活性
Hydrolase activity分子功能
Molecular function2.1E-07465胞内成分
Intracellular part细胞组成
Cellular component3.1E-079124结合过渡金属离子
Transition metal ion binding分子功能
Molecular function9.8E-0910252线粒体膜成分
Mitochondrial membrane part细胞组成
Cellular component9.9E-0860结合金属离子
Metal ion binding分子功能
Molecular function1.6E-1111760结合核酸
Nucleic acid binding分子功能
Molecular function8.0E-1311524谷氨酰胺家族氨基酸代谢过程
Glutamine family amino acid metabolic process生物过程
Biological process1.4E-0580肽酶活性
Peptidase activity分子功能
Molecular function6.6E-112725结合阳离子
Cation binding分子功能
Molecular function2.9E-1111760脂肪酸生物合成过程
Fatty acid biosynthetic process生物过程
Biological process4.8E-0640单糖分解代谢过程
Monosaccharide catabolic process生物过程
Biological process7.7E-0625胞内蛋白代谢过程
Cellular protein metabolic process生物过程
Biological process2.2E-06412小分子代谢过程
Small molecule metabolic process生物过程
Biological process6.7E-062110金属内肽酶活性
Metalloendopeptidase activity分子功能
Molecular function1.1E-0691水解酶活性
Hydrolase activity分子功能
Molecular function1.5E-0912062结合嘌呤核苷三磷酸
Purine ribonucleoside triphosphate binding分子功能
Molecular function3.2E-079532氧化还原酶活性
Oxidoreductase activity分子功能
Molecular function1.1E-1614159依赖NTP嘌呤解旋酶活性
Purine NTP-dependent helicase activity分子功能
Molecular function2.6E-07173单价无机阳离子跨膜运输活性
Monovalent inorganic cation transmembrane
transporter activity分子功能
Molecular function1.1E-0565

新窗口打开

表6   菌丝体与原基pathway 显着富集分析列表

Tab.6   List of pathway enrichment analysis of mycelium and primordium

代谢途径
Metabolic pathwayP值
P value上调基因数
Up-regulated genes下调基因数
Downk-regulated genes氧化磷酸化途径
Oxidative phosphorylation7.2E-05121萜类生物合成途径
Terpenoid backbone biosynthesis0.00011860核糖体
Ribosome1.8E-06210花生四烯酸代谢途径
Arachidonic acid metabolism6.0E-0541莨菪、哌啶、吡啶类生物碱合成
Tropane, piperidine and pyridine alkaloid biosynthesis5.7E-0552mRNA检测途径
mRNA surveillance pathway0.000155140次级代谢产物生物合成途径
Biosynthesis of secondary metabolites2.4E-075015色氨酸代谢
Tryptophan metabolism0.000225145丙氨酸、天冬酰胺、谷氨酸代谢相关
Alanine, aspartate and glutamate metabolism0.00019461脂肪酸代谢相关
Fatty acid metabolism3.4E-05131苯丙氨酸、酪氨酸、色氨酸生物合成途径
Phenylalanine, tyrosine and tryptophan biosynthesis0.00023880

新窗口打开

The authors have declared that no competing interests exist.

参考文献

[1]汤庆莉, 张洪, 张勇, 等.

灰树花胞外多糖提取工艺优化

[J].核农学报, 2013, 27(11):1704-1708.

Magsci     [本文引用: 2]    摘要

<p>通过对乙醇浓度、醇沉时间、醇沉温度3个单因素试验对灰树花胞外多糖提取过程中的变化规律进行研究,并且通过正交试验法,以灰树花EPS提取量为考察指标,对灰树花胞外多糖提取工艺条件进行优化。结果表明:灰树花胞外多糖的最佳提取条件为乙醇浓度95%、醇沉时间10h、醇沉温度4℃。在此条件下灰树花胞外多糖的提取量为321.94mg&middot;L^-1,较未优化的灰树花胞外多糖提取条件提高了57.79%。</p>

[2]杨海, 耿传信, 周学锋.

灰树花活性多糖药理研究综述

[J].中国执业药师, 2012, 9(3):30-33.

DOI:10.3969/j.issn.1672-5433.2012.03.008    URL     [本文引用: 1]    摘要

灰树花是一种名贵的食用、药用真菌,具有增强机体免疫、抗肿瘤、调血脂及保肝等多种生物活性,本文就近年来灰树花多糖的药理研究进展进行了综述。

[3]张桂春, 刘玉静, 解卫海, 等.

响应面法优化灰树花发酵全液多糖提取工艺及其抗肿瘤活性

[J].食品科学, 2013, 34(24):30-35.

URL     [本文引用: 1]    

[4]许占伍, 张引芳, 金力, 等.灰树花工厂化袋栽原基形成条件研究[J].食用菌, 2012(2):7-9.

DOI:10.3969/j.issn.1000-8357.2012.02.004    URL     [本文引用: 1]    摘要

灰树花为食用菌工厂化生产研究的新品种。对袋栽灰树花后熟期、环 境条件等对原基形成的影响研究。结果表明：菌丝发满袋后熟10～15d，袋壁割口搔菌诱导出菇，菌丝恢复期和现蕾期温度控制 为：20～22℃；16～18℃，光照控制为：0h／24h，空气相对湿度控制为：97％-99％；89％～91％，通风量控制 为：4min／10min；10min／10min，有利于原基的形成。

[5]潘辉, 李正鹏, 王瑞娟, 等.

灰树花子实体发育过程研究

[J].食用菌, 2010, 32(5):9-11.

DOI:10.3969/j.issn.1000-8357.2010.05.005    URL     [本文引用: 1]    摘要

在工厂化栽培设施条件下,使用 工厂化栽培菌株"灰M云"进行袋栽试验,观察记录灰树花原基形成及分化全过程。从形态学角度出发,将灰树花原基形成过程划分为:菌丝恢复期、现蕾期、原基 愈合期,其中关键时期为菌丝恢复期;灰树花原基分化过程分为:脑状体期、蜂窝期、珊瑚期及成熟期,其中关键时期为珊瑚期。

[6]Tomas Maseka, Vaclav Vopalenskya, Petra Suchomelovab, et al.

Denaturing RNA electrophoresis in TAE agarose gels

[J].Analytical Biochemistry, 2005, 336(1):46-50.

DOI:10.1016/j.ab.2004.09.010    URL    PMID: 15582557     [本文引用: 1]    摘要

Abstract Current methods of analytical RNA electrophoresis are based on the utilization of either complicated laboratory instrumentation or toxic, carcinogenic, or expensive chemicals. We suggest here the use of classical Tris-acetate-ethylenediamine tetraacetic acid (TAE) agarose gels combined with prior denaturation of RNA samples in hot formamide for the electrophoretic separation of RNA species. We present a brief comparison of the proposed TAE/formamide method with the most common 3-(N-morpholino)propanesulfonic acid/formaldehyde agarose gel protocol and show that both methods produce comparable results for size determination of RNA molecules and subsequent Northern blotting of gels. In addition to purified RNA samples, the robustness of the TAE/formamide protocol is demonstrated by its suitability for the analysis of RNA quality in crude yeast cell lysates containing large amounts of proteins, DNA, and other contaminating molecules. We therefore propose the TAE/formamide agarose electrophoresis as a rapid, simple, and cheaper alternative to current methods of RNA electrophoresis. Additionally, another benefit is the reduced exposure of laboratory personnel to hazardous chemicals.

[7]Shen R, Fan J B, Campbell D, et al.

High-throughput SNP genotyping on universal bead arrays

[J].Mutat Res, 2005, 573(1/2):70-82.

DOI:10.1016/j.mrfmmm.2004.07.022    URL    PMID: 15829238     [本文引用: 1]    摘要

We have developed a flexible, accurate and highly multiplexed SNP genotyping assay for high-throughput genetic analysis of large populations on a bead array platform. The novel genotyping system combines high assay conversion rate and data quality with >1500 multiplexing, and Array of Arrays鈩 formats. Genotyping assay oligos corresponding to specific SNP sequences are each linked to a unique sequence (address) that can hybridize to its complementary strand on universal arrays. The arrays are made of beads located in microwells of optical fiber bundles (Sentrix 庐; Array Matrix) or silicon slides (Sentrix BeadChip). The optical fiber bundles are further organized into a matrix that matches a 96-well microtiter plate. The arrays on the silicon slides are multi-channel pipette compatible for loading multiple samples onto a single silicon slide. These formats allow many samples to be processed in parallel. This genotyping system enables investigators to generate approximately 300,000 genotypes per day with minimal equipment requirements and greater than 1.6 million genotypes per day in a robotics-assisted process. With a streamlined and comprehensive assay, this system brings a new level of flexibility, throughput, and affordability to genetic research.

[8]Haas B J, Papanicolaou A, Yassour M, et al.

De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis

[J].Nature Protocols, 2013, 8(8):1494-1512.

[本文引用: 1]    

[9]Mortazavi A, Williams B A,

McCue K, et al.Mapping and quantifying mammalian transcriptomes by rna-seq

[J].Nat Methods, 2008, 5(7):621-628.

DOI:10.1038/nmeth.1226    URL    PMID: 18516045     [本文引用: 1]    摘要

We have mapped and quantified mouse transcriptomes by deeply sequencing them and recording how frequently each gene is represented in the sequence sample (RNA-Seq). This provides a digital measure of the presence and prevalence of transcripts from known and previously unknown genes. We report reference measurements composed of 4109恪52 million mapped 25-base-pair reads for poly(A)-selected RNA from adult mouse brain, liver and skeletal muscle tissues. We used RNA standards to quantify transcript prevalence and to test the linear range of transcript detection, which spanned five orders of magnitude. Although >90% of uniquely mapped reads fell within known exons, the remaining data suggest new and revised gene models, including changed or additional promoters, exons and 309銇05 untranscribed regions, as well as new candidate microRNA precursors. RNA splice events, which are not readily measured by standard gene expression microarray or serial analysis of gene expression methods, were detected directly by mapping splice-crossing sequence reads. We observed 1.45 01— 105 distinct splices, and alternative splices were prominent, with 3,500 different genes expressing one or more alternate internal splices.

[10]Ye J, Fang L, Zheng H, et al.

WEGO:a web tool for plotting GO annotations

[J].Nucleic Acids Research, 2006, 34(2):293-297.

DOI:10.1093/nar/gkl031    URL    PMID: 16845012     [本文引用: 1]    摘要

Abstract Unified, structured vocabularies and classifications freely provided by the Gene Ontology (GO) Consortium are widely accepted in most of the large scale gene annotation projects. Consequently, many tools have been created for use with the GO ontologies. WEGO (Web Gene Ontology Annotation Plot) is a simple but useful tool for visualizing, comparing and plotting GO annotation results. Different from other commercial software for creating chart, WEGO is designed to deal with the directed acyclic graph structure of GO to facilitate histogram creation of GO annotation results. WEGO has been used widely in many important biological research projects, such as the rice genome project and the silkworm genome project. It has become one of the daily tools for downstream gene annotation analysis, especially when performing comparative genomics tasks. WEGO, along with the two other tools, namely External to GO Query and GO Archive Query, are freely available for all users at http://wego.genomics.org.cn. There are two available mirror sites at http://wego2.genomics.org.cn and http://wego.genomics.com.cn. Any suggestions are welcome at wego@genomics.org.cn.

[11]刘朋虎, 邓优锦, 江玉姬, 等.

草菇菌丝体与原基差异表达基因分析

[J].基因组学与应用生物学, 2012, 31(2):147-153.

DOI:10.5376/gb.cn.2012.01.0003    URL    摘要

采用Solexa测序技术,对草菇菌丝体和原基进行了数字基因表达谱（DGE）测序,在菌丝 体和原基文库中分别得到5701781个和5659262个高质量测序标签（cleantags）,对应的标签种数（distinct clean tags）分别为85626和95363。将所有高质量测序标签与参考基因库进行比对,在菌丝体和原基文库中,占标签种数的43.32%和52.57%的 标签可以唯一定位（map）到参考序列上,占标签种数的21.65%和21.47%的标签可以被定位到基因组序列上。最终,被菌丝体和原基标签唯一定位的 基因数（unambiguous tag-mapped genes）分别为14794和15534。差异基因分析显示,两个文库中共有显著性差异表达的基因4163个,其中在原基中上调、下调的基因数分别为 2486和1677,只在原基中表达的基因321个。经过Blastnr比对,在原基中特异表达的基因,涉及蛋白质（氨基酸）合成与代谢、糖代谢、脂类代 谢和抗逆反应等多个代谢途径。GO功能富集分析结果表明,葡萄糖、己糖和乙醇等代谢途径大部分基因下调表达。Pathway功能富集分析结果表明,合成核 糖体蛋白的基因均下调表达,表明原基形成时细胞代谢减弱,蛋白质合成量减小。

[12]刘芳, 王威, 谢宝贵.金针菇菌丝与原基差异表达基因分析[J].食用菌学报, 2014(1):1.

DOI:10.3969/j.issn.1005-9873.2014.01.001    URL     [本文引用: 1]    摘要

以金针菇(Flammulina velutipes)1123菌株的单孢 W23基因组为参考完成其菌丝和原基转录组测序与数据分析，研究两样本间差异基因，并对差异基因进行 GO功能和 KEGG Pathway 显著性富集分析。结果表明：两个样本中共有显著性差异表达的基因3310个，其中在原基中上调、下调的基因数分别为1686、1624个，只在原基中表达 的基因有26个。GO功能分析结果表明，在原基中膜封闭腔、内膜系统、细胞器腔、核糖核蛋白复合体、翻译调节器、发育过程、免疫系统过程、多细胞生物过程 这8个GO基本单元中的差异基因全部呈现上调表达。Pathway 功能富集分析结果表明，核糖体与DNA复制中的差异基因全部呈现上调表达，糖酵解途径中从D-葡萄糖转化成丙酮酸过程相关的11个差异基因全部呈下调表 达，磷酸戊糖途径中相关的13个差异基因中有11个基因呈下调表达。

相关知识

真菌灰树花菌丝体转录组测序及分析.pdf
‘月月粉’响应灰霉菌侵染的转录组与代谢组分析
岷江百合响应灰霉菌侵染的转录组分析及抗灰霉病相关基因的挖掘
不同碳源条件下香蕉枯萎病菌2个小种转录组差异分析
基于转录组分析玉米弯孢霉叶斑病菌Clf敲除突变体差异表达基因
不同品种灰毡毛忍冬ACS3基因克隆、表达及生物信息学分析
不同品种茶树花发育过程中的转录组分析及激素响应
红鳍笛鲷皮肤转录组分析及体色相关基因克隆、表达与进化研究
菊花‘神马’花芽分化期的数字基因表达谱及关键基因的克隆和表达分析
国内首篇！10x Genomics 植物空间转录组文章已上线！

网址: 灰树花菌丝体与原基转录组差异表达分析 https://m.huajiangbk.com/newsview1278765.html

上一篇: 「美罗牌灵芝菌丝体松茸灰树花胶囊

下一篇: 那些你需要知道的关于食用菌菌丝体