睡莲品种保罗蓝花器官不同部位的转录组测序分析

为探究热带睡莲花器官不同部位的花香代谢通路及参与萜类香气物质生物合成的差异表达基因,本研究借助转录组测序技术,以热带睡莲品种保罗蓝为研究对象,对其花器官的花瓣（PE,petal）、雄蕊（ST,stamen）和雌蕊（PI,pistil） 3个部位进行转录组测序分析。差异表达基因分析结果显示,花瓣相对于雌蕊（PE-vs-PI）、雄蕊相对于雌蕊（ST-vs-PI）和雄蕊相对于花瓣（ST-vs-PE）的差异表达基因数目分别为7853个、7501个和2526个。GO分类和富集分析显示,3个比较组的差异表达基因主要参与了生物调节、细胞过程、代谢过程和刺激应答的生物学过程;KEGG分类和富集分析显示,PE-vs-PI的差异表达基因显著富集的KEGG通路最多,其次为ST-vs-PI,ST-vs-PE最少。从3个比较组共有的794个差异表达基因中筛选出98个参与萜类物质代谢差异表达基因,富集于4条萜类物质合成通路,且PE-vs-PI和ST-vs-PI的差异表达基因数目均高于ST-vs-PE。已知的金合欢醛和二萜贝壳杉烯合成关键基因HMGR和DXS在花瓣和雄蕊的表达量均高于雌蕊。从98个差异表达基因中随机...     

基于RNA-seq对南蛇藤MADS-box转录因子序列及表达分析

为了分析南蛇藤MADS-box转录因子对其雌花发育和果实形成的调控信息,本文基于南蛇藤的转录组数据,获得15个具有全长开放阅读框的MADS-box转录因子（ColMADS）,并利用生物信息学方法对其性质和结构特性进行了分析。结果表明,基因comp37814_g1和comp41380_g2属于M-type家族,不含有K盒结构,其他均为MIKC家族;除了comp37713_g1之外,其余均属于不稳定的亲水性蛋白;二级结构均含有α-螺旋、扩展链结构、β-转角和无规则卷曲,其中α-螺旋所占比例最高;序列主要包含5种保守基序,基序1和基序2分别含50个氨基酸且属于该基因家族的保守基序,仅基序1含有MADS盒。系统进化分析结果显示南蛇藤ColMADS转录因子被分为10个进化支,分属于不同类型的亚家族及不同的组。另外,利用荧光定量PCR检测方法揭示了这些基因在南蛇藤盛花、落花和幼果不同发育时期的表达情况。结果表明,2个基因在盛花期相对表达值最高,9个基因在落花期相对表达值最高,在幼果形成过程中有4个基因显著上调表达。     

翅果油树种仁脂肪酸合成关键基因筛选

【目的】通过翅果油树不同时期种仁的转录组学分析,挖掘与脂肪酸合成相关的关键基因。【方法】观察记录翅果油树果实发育情况并以成年翅果油树果实开始膨大时期（C6）、果实迅速膨大时期（C7）和果实成熟时期（C8）的新鲜种仁为试验材料,以水解-提取法测定3个时期种仁的脂肪酸组分与含量,利用转录组测序技术分析C6、C7时期样品表达差异,并寻找关键差异表达基因。【结果】3个时期的脂肪酸平均含量随时间推移呈逐渐上升趋势;C6、C7时期的种仁数据中共获得6 375个上调表达基因,8 124个下调表达基因;8 137个差异表达基因在GO数据库中获得注释;4 275个差异表达基因在KEGG数据库中获得注释,主要涉及代谢途径,次生代谢物的生物合成,激素传导,氨基酸合成途径等通路,其中,脂肪酸生物合成、不饱和脂肪酸生物合成等与脂肪酸相关的通路均被富集,并重点挖掘出8个与脂肪酸合成相关的候选基因,发现与果实膨大时期（C6）相比,EmBCCP1、EmSAD1、EmFAD2-1、EmFAD3-1、EmPDCT1基因的上调表达和EmFAbG1、EmPDAT1、EmKASI1基因的下调表达。【结论】EmBCCP1、EmSAD1、EmFAD2-1、EmFAD3-1、EmPDCT1、EmKASI1、EmFAbG1、EmPDAT1基因可能促进了翅果油树种仁油酸、亚油酸的合成积累,可作为翅果油树种仁脂肪酸合成的关键候选基因     

基于转录组挖掘‘云香’水仙萜类合成基因及NaIDI的克隆

为深入挖掘中国水仙转录组中萜类合成相关基因,了解中国水仙萜类代谢途径和分子调控机制,该研究以‘云香’水仙为材料,在其盛花期花瓣与副冠转录组测序的基础上,筛选已注释的萜类合成途径基因并利用NCBI blastn进行再注释,通过对部分候选基因的表达量与生物信息分析进一步筛选出代表基因,采用RT PCR技术克隆了水仙异戊烯基焦磷酸异构酶基因NaIDI,并对其蛋白序列与特异性表达进行分析。结果表明：(1)Blast比对筛选得到52 个与萜类合成上游途径相关的Unigenes,二次筛选获得11个显著差异表达的代表基因。(2)NaIDI基因开放阅读框(ORF)长度为858 bp,编码285个氨基酸,其氨基酸序列与‘金盏银台’水仙相似度97.19%;亚细胞定位预测显示该基因定位在叶绿体;系统进化分析表明其与芦笋亲缘关系比较近。(3)实时荧光定量分析结果表明,NaIDI基因在水仙开花的不同时期和不同组织器官中差异表达显著,且在盛花期时副冠中的表达量最高,与中国水仙挥发性萜类化合物在开花不同时期和不同组织器官中的表达规律一致,表明NaIDI在萜类代谢中发挥着一定作用。     

款冬花不同发育阶段的代谢组学和比较转录组学分析

以不同发育阶段款冬花（发育初期、中期、中后期、解封期及花朵期）为对象,基于核磁共振的代谢组学技术,分析其次生代谢物种类及含量的合成累积规律,并进行高通量转录组学测序,从差异表达的基因中寻找次生代谢物生物合成的关联酶基因。代谢组学分析发现,款冬花不同发育阶段的次生代谢物代谢组成明显不同,苯丙素类成分在发育初期至中后期含量较高,之后逐渐降低;黄酮类成分（芦丁、山奈酚）在发育的各个阶段含量均有波动,但总体变化不大。转录组学分析结果显示,与苯丙素类成分生物合成相关的酶基因（COMT、HCT）,随着花蕾的发育表达量逐渐降低;与黄酮类成分合成相关的酶基因（FLS、F3H、DFR）,其表达量在不同阶段变化不大。转录组测序中的相关酶基因表达量同次生代谢物成分的变化趋势基本一致。本文采用的代谢组学和转录组学结合的方式,对款冬花的次生代谢物累积规律进行分析,为今后款冬花次生代谢物的生物合成调控研究奠定了基础。     

中粒咖啡萜类合成酶基因家族的生物信息学分析

萜类化合物具有重要的生理、生态作用和药用价值,萜类合成酶(TPS)是合成萜类化合物的关键酶。通过整合中粒咖啡(Coffee canephora)的基因组和转录组数据,利用生物信息学方法,鉴定出43个萜类合成酶全长基因,并对这些基因的分子进化、结构、复制、表达及功能分化的机理进行了探究。结果表明,中粒咖啡萜类合成酶基因可以分为5个亚家族(a、b、c、e/f、g),不同亚家族的基因结构差异很大;串联复制是基因家族扩增的主要原因;表达分析结果表明,萜类合成酶基因在不同组织中的表达差异明显;中粒咖啡萜类合成酶基因启动子区的顺式调控元件可能与基因的功能分化相关;不同亚家族之间的功能差异主要由亚家族特异的氨基酸决定。     

冬虫夏草MYB家族基因鉴定及其在生长发育中的表达分析

MYB蛋白是一类广泛存在于真核生物中的转录因子,在真菌的生长发育中发挥调控作用。本研究对冬虫夏草菌中的潜在MYB转录因子家族基因进行了全基因组鉴定和生物信息学分析,并研究了MYB家族成员在冬虫夏草生长发育不同阶段和不同部位的表达模式。结果表明,冬虫夏草MYB转录因子家族包含3个1R-MYB和3个2R-MYB;6个MYB转录因子蛋白均包含近50个保守氨基酸序列,形成螺旋-转角-螺旋的结构。通过分析MYB转录因子基因在不同发育阶段（MYB-1-6）和子实体不同部位（MYB-1、3、6）的相对表达量,发现MYB-1在冬虫夏草整个发育阶段表达较为稳定,MYB-3在子实体成熟阶段（MF）表达量最高,远高于其他阶段,且在MF的顶部可育部分MF-3高表达,表明其可能参与冬虫夏草子实体的有性发育;MYB-6在子座发育初期（YF）的中部YF-2表达量最高,为菌丝阶段（HY）的5倍,表明MYB-6可能参与子座柄部的伸长。     

基于转录组分析高温抑制广东虫草原基形成的调控网络

在大型真菌原基形成和子实体发育的过程中,温度是一个极其关键的因素,高温显著影响多种食用菌原基的形成和子实体的品质。广东虫草是中国华南地区特有的虫草类新食品原料,其子实体富含多种活性成分和营养成分,且能够进行较大规模人工栽培。然而,温度对广东虫草原基形成的影响及调控机制尚不清楚。本研究通过不同温度处理不同时间段,发现29℃处理3d抑制了广东虫草原基的形成;对29℃处理3d处理前后的菌丝阶段（CK）和原基阶段（CK-P和HT-P）样品进行转录组测序分析,高温处理后原基阶段的两个组中发现1 682个差异表达基因,其中1 015个上调表达和667个下调表达。在碳水化合物代谢途径中,多个糖酵解和三羧酸循环及葡聚糖和海藻糖合成相关基因在高温处理后呈现下调表达;在原基样品中,多个热激蛋白基因（Hsp10、Hsp23、DnaJ、Hsp70、Hsp90和Hsp98）和转录因子C2H2转录水平显著上调表达。本研究结果基于分子水平揭示了高温影响原基形成过程中能量代谢和相关基因的差异表达,为后续利用广东虫草抗逆相关基因资源培育新品种奠定了重要的基础。     

光照强度对成熟红颜草莓果实着色和花青素生物合成的影响及可能的分子机制

为了研究光照对红颜草莓（Fragaria×ananassa Duch.‘Benihoppe’）果实着色、花青素含量及花青素合成相关基因表达的影响,本文采用高效液相色谱法（HPLC）和实时荧光定量PCR技术测定了不同遮阴处理下（透光率100%、75%、25%）草莓果实花青素含量及色素相关基因的表达情况。结果显示,与100%透光率下的草莓果实相比,在75%和25%透光率下合成的花青素含量分别下降了41.58%和92.54%;和花青素合成相关基因的表达也都有不同程度的下降,其中二氢黄酮醇4-还原酶基因（FaDFR）、类黄酮-3'羟化酶基因（FaF3'H）和类黄酮3-O-糖基转移酶基因（FaUFGT）的表达与花青素含量相关性达到显著水平;此外,在遮光条件下转录因子FaMYB10、FaMYB1表达也明显下调。可见,光照是影响草莓果实着色的关键环境因素,遮光抑制色素相关功能基因和转录因子的表达,阻碍了果实中花青素的合成,最终导致果实着色差异。     

杜仲雄花芽2个发育时期转录组分析

杜仲（Eucommia ulmoides）雄花富含多种活性成分和营养成分,具有重要的药用和营养价值。为了揭示杜仲雄蕊原基发育相关基因的表达情况,为杜仲雄花芽发育分子调控机制研究提供理论参考。本文以杜仲良种"华仲11号"（ "Huazhong No.11" ）为材料,采用lllumina高通量测序技术,分别对苞叶原基分化期和雄蕊原基分化期的花芽进行转录组测序,通过生物信息学对2个发育时期的转录组进行比较分析,筛选出与雄花芽形态发育相关的差异基因。结果显示,转录组测序共获得40.48 Gb过滤数据,各样品的clean reads与杜仲基因组进行序列比对,比对效率为90.56%~93.01%。在2个发育时期筛选出583个差异表达基因,其中在雄蕊原基发育期上调基因315个,下调基因267个。差异基因GO和KEGG功能分析显示,差异基因富集在与生长发育、光周期途径、激素合成和信号传导、碳代谢等相关的生物过程和代谢通路。结果显示光周期途径是杜仲成花诱导的重要途径,同时雄花芽在形态分化过程中受碳水化合物、植物激素和其他代谢物质调控。此外,MADS-box家族成员FLC、SOC1、AGL3和AGL8参与杜仲雄蕊器官发育。本研究为杜仲花发育基因调控提供了基础数据,也为雄花用杜仲的分子育种提供了参考。     

基于转录组的瑶药半枫荷根和叶的差异分析

为了探究半枫荷(Semiliquidambar cathayensis)根和叶的基因表达差异和关键活性成分合成通路中关键基因的表达规律,本研究对半枫荷的根和叶进行转录组测序和生物信息分析。59378个差异表达基因归类到在GO分类的3个大类中,主要与生物学过程有关(50.59%)。626个差异表达基因注释KOG数据库的24个分类中。81个差异表达基因参与苯丙烷类化合物的生物合成,110个差异表达基因参与黄酮类化合物的生物合成,211个差异表达基因参与萜类化合物的生物合成。本研究获得了半枫荷根和叶的转录组信息特征,为今后半枫荷基因功能鉴定、次生代谢途径解析及调控机制的研究提供依据。     

金针菇菌柄发育的转录组与蛋白组分析

菌柄是金针菇等食用菌的主要商品部位,但其生长机制仍不明确。本研究对金针菇伸长期和成熟期菌柄进行了转录组联合蛋白组分析,结果显示,两样本显著性差异表达基因和蛋白分别为721个和61个,均以上调表达为主。GO（gene ontology）功能聚类分析表明：有72.41%的差异表达基因富集在催化活性（catalytic activity）条目下。细胞组分（cell part）和绑定结合（binding）条目同时富集了较多的差异表达基因和蛋白。KEGG通路富集分析显示：碳水化合物代谢通路（carbohydrate metabolism）和氨基酸代谢通路（amino acid metabolism）富集的差异表达基因较多。差异表达蛋白富集较多的通路是单环菌素生物合成（monobactam biosynthesis,ko00261）、链霉素生物合成（streptomycin biosynthesis,ko00521）和有机含硒化合物代谢（selenocompound metabolism,ko00450）等。内质网蛋白质加工（protein processing in endoplasmic reticulum,ko04141）和MAPK信号通路（MAPK signaling pathway-yeast,ko04011）在转录组和蛋白组的KEGG富集分析中均为差异通路。本研究联合转录组和蛋白组数据筛选了40个金针菇菌柄发育中差异表达基因,为深入研究揭示食用菌菌柄发育过程提供候选基因。     

沙棘刺和芽的组织解剖与转录组分析北大核心CSCD

该研究通过组织切片分析比较了沙棘刺与芽的组织染色差异,通过转录组测序分析比较了刺和芽的基因表达谱差异,通过qRT-PCR验证候选基因的表达水平,为进一步通过基因工程手段调控沙棘刺的发育提供候选基因。结果表明:(1)沙棘芽横切面的组织边界层次分明,而沙棘刺横切面的组织边界模糊,且刺尖中心的红色区域较大,木质素自发荧光结果表明刺中木质化程度显著增强。(2)转录组分析组装得到81560个单基因(unigenes),其中包含9385个差异表达基因(DEGs)。(3)GO和KEGG富集分析表明,富集程度最高的基因主要与细胞壁、表皮、气孔发育有关,尤其是与木质素合成途径相关;芽刺相关DEGs分析结果表明,木质素合成相关DEGs共20个,角质、木栓质和蜡质生物合成相关DEGs共29个,植物表皮发育和气孔复合体发育的相关DEGs相同共14个;结合文献分析推测,部分热头蛋白基因参与控制沙棘刺顶端分生组织消失,木质素合成相关基因则与刺硬化过程密切相关。(4)qRT-PCR分析验证显示,与芽相比,挑选的6个木质素合成相关DEGs中有3个分别上调了9.5倍、41.1倍和13.7倍,3个下调至5%、14%和24%;挑选的4个角质、木栓质和蜡质生物合成相关DEGs中有2个分别上调3.6倍和22.1倍,另外2个下调至3%和6%;挑选的4个表皮和气孔发育相关DEGs均下调至6%、13%、3%和4%;关键基因的相对表达水平与转录组测定的TPM值可相互印证,表明验证结果与转录组分析结果一致。研究发现,沙棘刺和芽具有明显的形态和组织差异,二者间的DEGs主要与木质素合成以及细胞壁、表皮和气孔发育有关。     

PEG模拟干旱胁迫下马铃薯茎段转录组分析

该研究以‘青薯9号’马铃薯无菌苗为材料,采用转录组测序技术分析模拟干旱胁迫下马铃薯茎段的差异表达,探究茎段在干旱胁迫下的分子机制。结果表明:(1)不同程度干旱胁迫下,马铃薯叶片脯氨酸、可溶性糖以及可溶性蛋白含量明显增加;马铃薯茎段差异表达基因下调的数量均多于上调,其中3种处理条件下共有的差异表达基因有657个。(2)GO富集分析表明,马铃薯茎段差异表达基因主要集中在氧化还原过程、激素响应、氧化还原酶活性以及糖基水解酶活性;Pathway富集分析表明,马铃薯茎段差异表达基因主要集中在植物激素信号转导、苯丙酸生物合成、玉米素生物合成、苯丙氨酸代谢、淀粉和蔗糖代谢以及次生代谢产物的生物合成。(3)实时荧光定量PCR验证结果表明,6个差异表达基因在不同程度干旱胁迫中的差异表达与转录组分析的结果基本一致,证明转录组数据的可靠性。该结果对进一步研究马铃薯干旱胁迫响应机制有一定参考价值,也丰富了马铃薯抗旱育种的基因资源。