盐胁迫下两个甜瓜品种转录因子的转录组分析

利用新一代高通量测序手段——转录组测序（RNA-Seq）技术研究在300 mmol.L-1NaCl胁迫下两个甜瓜（Cucumis melo L.）品种‘玉露’和‘冰雪脆’的转录因子基因表达变化。这两个甜瓜品种在叶绿素荧光参数上的差异表明其在盐胁迫下有不同的生理反应。转录组测序结果表明,盐胁迫与对照相比,‘玉露’共有属于19个转录因子家族的56个转录因子基因表达发生变化（在转录水平,属于7个转录因子家族的22个转录因子上调表达,属于14个转录因子家族的34个转录因子下调表达）。‘冰雪脆’有属于20个转录因子家族的47个转录因子基因表达发生变化（在转录水平上,属于5个转录因子家族的17个转录因子上调表达,属于17个转录因子家族的30个转录因子下调表达）。盐胁迫下,两个甜瓜品种差异表达的转录因子既表现特异性,也存在部分重叠。‘玉露’有29个转录因子特异响应,‘冰雪脆’有20个特异响应。盐胁迫响应重叠的转录因子有27个,其中9个上调表达,18个下调表达。采用实时荧光定量PCR对几个转录因子进行了盐胁迫下的表达检测,其趋势与转录组分析结果基本一致。     

低温胁迫下两种不同形态群体拟穴青蟹转录组分析

为研究拟穴青蟹(Scylla paramamosain)不同形态群体对温度变化的适应差异,以Sp1群体(腕节外缘中部都不具有发达的刺)和Sp2群体(至少一只螯足腕节外缘中部具有发达的刺)为研究对象,选用低温胁迫组(8℃)和对照组(20℃)拟穴青蟹肝胰腺为实验材料,进行转录组测序。结果显示,测序共获得81853个Unigene,平均长度为420 bp, N₅₀为1460 bp,其中22.33%为已知基因。差异表达基因筛选发现Sp1群体在低温胁迫下共有15773个Unigene显著差异表达,其中15628个Unigene表达下调, 145个Unigene表达上调,这些Unigene显著富集于细胞质基质等生物学过程和RNA结合、核糖体结构组分等分子功能,并进一步富集在核糖体、剪接体和内质网蛋白质加工等信号通路。Sp2群体在低温胁迫下共有323个Unigene显著差异表达,其中上调114个,下调209个。这些Unigene显著富集于DNA整合等生物学过程和甾体羟化酶活性、碳酸氢盐跨膜转运体活性等分子功能,并进一步富集在赖氨酸降解、甘氨酸、丝氨酸和苏氨酸的代谢、精氨酸和脯氨...     

鲤在低温胁迫下肝胰腺转录组测序分析

检测低温耐受与低温敏感鲤品种低温胁迫下的差异表达基因,探究鲤低温适应的分子应答过程。低温胁迫后对松浦镜鲤及荷包红鲤肝胰腺进行转录组测序,并对表达差异基因进行功能、通路富集分析。通过差异表达分析,共获得10 521个差异表达基因,其中5 246个基因仅在低温耐受品种发生表达变化。功能、通路富集分析将差异基因显著富集在包括糖酵解/糖质新生代谢途径在内的144个代谢途径。选取5个基因进行qRT-PCR验证,定量结果与RNA-seq结果相关性达0.89。研究结果表明低温胁迫会显著改变鲤肝胰腺的转录活动,编码糖酵解/糖质新生途径的多个关键酶基因在低温耐受品种的表达量显著高于低温敏感品种,推测这些基因功能有助于其长期低温适应。     

萱草叶片响应低温胁迫的转录组分析

低温胁迫是萱草（Hemerocallis fulva）生长过程中经常会遭遇的一种非生物胁迫。比较了萱草叶片在低温处理（10、5、0 ℃）下转录组与对照（15 ℃）数据的差异,共筛选出差异表达基因2 457个,其中上调基因1 253个,下调基因1 204个。差异表达基因主要富集在细胞过程、代谢过程和催化活性等49个GO过程,代谢途径、次生代谢产物的生物合成、植物激素信号转导等42条KEGG代谢途径中。其中参与植物激素信号转导通路的差异表达基因发生了不同程度的变化,GH3.10基因上调至对照组的13.624倍,IAA1基因下调0.120倍;参与可溶性糖合成通路的差异基因发生了0.076~28.114倍不同程度的变化。随后对3个低温处理组共有的29个差异表达基因进行热图和网络调控分析,基于基因在网络调控中的位置,对ABCF5、OFPs和SWEETs等基因在冷应答的作用进行了分析。本研究结果为进一步挖掘萱草低温响应的关键基因及耐寒萱草种质开发、分子育种提供了一定的理论支撑。     

肺形侧耳低温胁迫时期的转录组分析

本文以肺形侧耳栽培菌株X57为研究对象,利用高通量RNA测序技术对4℃低温处理0h、6h和12h后的肺形侧耳进行基因表达分析来探讨肺形侧耳低温应答机制。分析结果筛选出低温处理6h差异表达基因742个,其中上调表达基因374个,下调表达基因368个;低温处理12h差异表达基因1 489个,其中上调表达基因占53%,下调表达基因占47%。Gene Ontology（GO）功能聚类分析表明,差异表达基因主要富集在结合、催化分子功能组和代谢过程、细胞过程生物学过程中。KEGG功能富集分析结果显示,差异基因主要富集在氨基酸、核糖体和类固醇生物合成,及氮类物质代谢的通路上,富集到与低温胁迫相关通路MAPK signaling pathway-yeast上的基因表达随低温处理时间的增加呈上调趋势。已报道HOG-MAPK通路是MAPK途径研究较为明确、信号传递单一的真菌低温胁迫中重要的通路,本文运用生物信息学软件构建肺形侧耳的HOG-MAPK通路,并利用荧光定量PCR对相关基因表达进行验证,结果显示以低温处理0h为参照,随着低温胁迫时间增加通路上大部分基因的表达量持续上升,且差异显著,与RNA‐Seq分析结果一致。本文通过全面分析肺形侧耳低温胁迫时期基因表达情况,为进一步研究肺形侧耳低温胁迫相关重要基因的功能鉴定与信号通路调控机理提供基础。     

冰菜盐胁迫下的转录组分析

为了解冰菜(Mesembryanthemum crystallinum)叶片抗盐相关基因组学,利用Illumina Hi-seq TM2500高通量测序技术研究冰菜叶片在400 mmol L~(–1) NaCl胁迫下转录组基因的差异表达。结果表明,从400 mmol L~(–1) NaCl胁迫和对照的冰菜叶片中共获得13.01 Gb Clean data,Q30碱基均大于90.08%。共获得123个差异表达基因(DEGs),包括73个上调基因,50个下调基因,其中功能注释的基因有96个。根据Unigene库序列进行GO、COG和KEGG注释,筛选出8个与抗盐性相关差异表达基因,植物激素代谢相关基因,脱落酸8'-羟基化酶、吲哚-3-乙酰酸酰胺合成酶和茉莉酮酸酯ZIM结构域蛋白基因均下调表达,生长素响应蛋白、细胞分裂素合酶基因则上调表达,糖代谢相关基因棉子糖合成酶基因上调表达,质膜H+-ATPase基因上调表达,脱水蛋白基因下调表达。这为冰菜耐盐基因组学和分子生物学的研究奠定基础。     

东乡野生稻苗期响应低温胁迫的转录组分析

为了解东乡野生稻（Oryza rufipogon）对低温胁迫的响应机制,对苗期的RNA-seq转录表达谱进行了研究。结果表明,与对照相比,共检测到10 200个差异表达基因（DEGs）,其中5 201个上调表达,4 999个下调表达,其中有426个DEGs位于已报道的水稻耐冷QTL区间,且37个为耐冷调控相关的家族基因。GO功能分类和KEGG代谢路径分析表明,核酸结合转录因子活性、氨基酸生物合成以及光合作用代谢等均参与响应低温胁迫过程。实时荧光定量分析表明,ABA响应蛋白基因、MYB转录因子和40S核糖体蛋白SA基因等12个可能与低温胁迫响应相关的DEGs表达模式与RNA-seq的一致。可见,植物激素传导途径和转录因子相关调控基因在东乡野生稻苗期响应低温胁迫过程中起重要作用。     

多年生黑麦草温度胁迫差异表达转录因子的比较转录组分析

多年生黑麦草是禾本科冷季型草种,温度胁迫严重影响其分布和产量。转录因子调控基因表达在植物响应逆境胁迫中发挥重要作用。该研究以多年生黑麦草品种雅晴为材料,采用高通量RNA-seq技术,对热(40℃)、冷冻(-10℃)和对照(22℃)处理的样本进行转录因子应答分析,以探索多年生黑麦草转录因子对温度胁迫的响应规律,并筛选抗逆转录因子候选基因。结果显示:(1)共鉴定了694个转录因子unigenes,分属于AP2/ERF、GTF、HSF、MYB、NAC、WRKY、bHLH和bZIP等32个家族。(2)在热和冷冻胁迫下,ERF(AP2/ERF)、MYB、NAC和bZIP家族基因均以上调为主,而WRKY家族则多数下调。HSF、GTF和DREB(AP2/ERF)家族成员大多数受热诱导上调,而多数bHLH家族成员受冷冻胁迫上调。(3)功能富集结果表明,多年生黑麦草差异表达的转录因子基因主要参与植物激素信号转导、昼夜节律和病原体互作等通路。研究表明,多年生黑麦草中与胁迫适应相关的转录因子基因普遍上调表达,而与生长和抗病相关的基因多数下调。该研究筛选到大量抗逆转录因子候选基因,为分子育种提供抗逆基因资源。     

准噶尔小胸鳖甲短时低温胁迫响应的转录组分析

【目的】拟步甲科昆虫小胸鳖甲Microdera punctipennis是分布于中国新疆古尔班通古特沙漠的特有物种,具有很强的耐寒性,但其低温响应的分子机制尚不明确。本研究旨在利用转录组测序技术丰富小胸鳖甲的已知基因信息,分析在短时低温胁迫下,上调表达基因的主要类群及参与的主要代谢通路。【方法】利用Trinity软件对分别在4℃低温和25℃(对照)下处理3 h的小胸鳖甲成虫RNA-Seq数据进行从头组装,利用Blast2go软件进行基因注释。通过DESeq和GFOLD软件分析差异表达基因,并对上调表达基因进行GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)代谢途径富集分析。【结果】测序过滤后得到68 165 001个有效读序,51 712个平均长度为984 bp的非冗余基因,其中29 925个基因(约57.87%)有同源序列(E-value10-5),与拟步甲科的赤拟谷盗Tribolium castaneum相应基因相似性最高(核苷酸序列一致性为72.75%)。低温响应上调表达基因有514个,富集到35个GO类群,18个KEGG途径。其中胁迫响应类群、生物调节类群和免疫系统过程类群都占有重要比例,嘌呤代谢、硫胺素代谢、糖酵解/糖异生等代谢通路显著富集。在胁迫响应类群中,热激蛋白Hsp90基因、超氧化物歧化酶(SOD)基因和钙联接蛋白Calnexin基因等8个非生物胁迫相关基因显著上调表达。【结论】荒漠昆虫小胸鳖甲低温转录组揭示,在4℃冷驯化过程中,代谢过程、胁迫响应、生物调节和免疫系统等生物学过程相关基因表达显著上调,其中几个非生物胁迫响应基因提示小胸鳖甲可能从分子伴侣、细胞周期阻滞、线粒体稳定性、抗氧化胁迫等多方面应对低温胁迫。KEGG富集分析所揭示的小胸鳖甲在低温下的转录组代谢通路与其他物种有较大差异。研究结果为后续生物信息学分析及小胸鳖甲耐寒关键基因的发掘及耐寒机制的揭示提供了基础数据。     

基于转录组测序的花烟草低温胁迫响应转录因子挖掘

为了鉴定参与花烟草低温胁迫的转录因子,对花烟草幼苗进行了4℃低温处理,并在处理后12 h采集其幼苗样本,提取总RNA后,采用高通量测序技术,进行了转录组测序。在对差异表达基因进行GO及KEGG分析的基础上,对参与其中的转录因子进行了挖掘,并采用qRT-PCR的方法,对转录组测序的结果进行了验证。结果表明,低温处理后,花烟草基因表达量变化在2倍以上的基因有8388个(P<0.01),其中,上调表达4229个,下调表达4159个。这些差异表达基因的功能归类于生物过程、细胞组分及分子功能3大类69个GO条目,并显著富集在40条KEGG代谢通路中。同时,在低温胁迫下,花烟草有118个转录因子的表达发生了显著改变,其中,上调表达82个,下调表达36个。这些转录因子属于28个家族,其中,数量最多的为NAC家族19个,其次为ERF家族16个、MYB家族15个、WRKY家族15个。本研究的结果为花烟草低温响应分子机制的研究提供了借鉴。     

基于转录组的不同火龙果品种抗性差异分析

不同的品种抗性不同,为进一步探究不同火龙果品种之间的抗性差异,为后续火龙果抗性育种提供参考,该研究利用Illumina HiSeq 2000测序平台对'普通白肉'(BR)和'厄瓜多尔黄龙'(EY)两个品种进行转录组测序分析,并参考GO Ontology、KEGG等公共数据库对差异表达基因进行功能分类与富集分析.结果表明...     

低温胁迫下不同甘蔗品种(系)光合特性的变化及其与耐寒性的关系

以广西农科院甘蔗研究所自育的7个新材料和2个生产上的主栽品种为研究对象,在甘蔗苗期进行低温胁迫处理,研究了各品种(系)甘蔗形态特征的冷害指数、叶绿素含量及光合特性相关指标的变化及其光合特性相关指标与甘蔗抗寒性间的相关性。结果表明:随着低温胁迫处理时间的延长,冷害指数不断增大,但变化的大小与快慢因品种(系)不同表现不一样。各甘蔗品种(系)叶片叶绿素含量均随时间延长而降低。叶片净光合速率、气孔导度在低温处理与常温处理间具有显著差异。低温胁迫处理显著降低了各甘蔗品种(系)最大光化学效率(Fv/Fm)、PSⅡ实际光能转化效率ΦPSⅡ、光适应下PSⅡ反应中心的最大光能转化效率Fv′/Fm′、光化学猝灭系数qP、电子传递速率ETR,而显著提高了初始荧光Fo、稳态荧光Fs、非光化学猝灭系数qNP。相关性分析表明整个测定时期各指标间相关显著,Fv/Fm、Fv′/Fm′、ΦPSⅡ与冷害指数I之间的相关系数在0.800以上,Fv/Fm、Fv′/Fm′、ΦPSⅡ可以作为甘蔗品种(系)抗寒性鉴定的重要参考指标。     

不同水分条件下杜鹃花转录因子的转录组分析

为了解杜鹃花(Rhododendron pulchurum)转录因子在干旱胁迫下的表达模式,利用高通量测序技术,对不同水分条件下"白凤4号"叶的转录因子表达谱进行分析。结果表明,在不同水分处理间有34~161个差异表达的转录因子,杜鹃花响应不同水分的转录调控主要是通过ERF、b HLH和MYB基因的表达协同完成的;干旱胁迫时,特异调动了NAC的差异表达和WRKY、b ZIP、PLATZ的上调表达,干旱后复水时特异调动了GATA表达来调控。RT-q PCR验证结果表明,基因表达趋势与测序结果一致,证明了测序结果的有效性。这为明晰杜鹃花抗旱的分子机理及分子育种奠定了理论基础。     

NaCl胁迫下黑果枸杞转录组测序分析

研究黑果枸杞在不同浓度盐胁迫下基因表达谱变化情况,为进一步研究黑果枸杞抗盐分子机制奠定研究基础。对0(CK)、50、250 mmol/L NaCl溶液胁迫的黑果枸杞组培苗的根和叶在胁迫时间为0、1、12 h时分别取样,采用转录组测序(RNASeq)技术进行测序分析。结果表明,转录组测序共产生222.49 Gb原始数据,拼接出Unigenes 86 037条,注释到7大功能数据库(GO、KEGG、KOG、NR、Pfam、Swiss-Prot和egg NOG)上的Unigenes总数为46 594个,占总Unigenes的54.76%,还有38 929个Unigenes在这些数据库中没有得到注释。通过GO分类和KEGG Pathway富集性分析,分别归于51个GO类别和211条代谢途径。差异表达基因分析显示,黑果枸杞叶片和根的上调基因和下调基因数随着NaCl浓度的增大和处理时间的延长均呈增加趋势,叶片中的上调基因数(7 514)小于下调基因数(9 032),根中的上调基因数(12 347)大于下调基因数(11 559)。在黑果枸杞盐胁迫下转录组中发现28 325个SSR位点,最多的为单核苷...     

镉胁迫下麻疯树转录组测序分析

麻疯树是一种能适应多种恶劣环境条件的能源植物,目前关于其抵抗重金属胁迫的分子调控机理尚不清楚。从组学水平整体分析其基因表达模式对于筛选关键基因、解析镉胁迫响应调控网络和促进分子育种具有重要意义。利用Illumina测序技术对水培条件下培养的处理组(Cd 100)和对照组(CK)麻疯树幼苗叶的转录组进行高通量测序,两个数据库的测序数据经de novo组装得到50448条高质量的Unigene。两个样品中发现了2551条差异表达基因,其中539条上调基因,2012条下调基因。根据不同数据库的注释信息,发现麻疯树镉胁迫引起叶片中多种代谢途径的变化,包括碳代谢,光合作用,植物激素信号转导以及植物病原响应途径。DAVID分析显示镉胁迫引起了麻疯树叶中与离子转运相关基因的变化,导致叶片中Na离子和铁离子稳态的变化。转录因子分析发现WRKY和ZIP在镉胁迫中发挥重要作用。用qRT-PCR技术对随机挑选的5个基因进行荧光定量验证,结果与测序数据一致,证实了差异表达基因数据的有效性。深入探讨了麻疯树镉胁迫的分子机理,为进一步应用于基因工程和植物修复提供基础。     

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

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

短期盐胁迫下盐穗木的转录组分析

盐穗木（Halostachys caspica）是荒漠盐碱地广泛分布的盐生植物,具有极强的耐盐性。为揭示盐胁迫下盐穗木基因组层面的基因表达变化特性,通过对300和500mmol·L^-1 NaCl胁迫3h的盐穗木同化枝进行了转录组测序。有效序列组装共得到153298条平均长度为643bp的unigenes,进行GO和KEGG功能聚类,分别获得47个GO功能小类和118个KEGG通路。差异表达基因分析显示,短期低盐（300mmol·L^-1）响应基因有4432个,高盐（500mmol·L^-1）响应基因有2580个,两个胁迫的共差异基因有1245个,主要富集在细胞过程、代谢过程和响应刺激等类别中。从短期盐胁迫下盐穗木转录组筛选出渗透调节和活性氧清除的相关基因,大多为上调基因。说明盐穗木能够通过促进渗透调节和增强活性氧清除提高短期的盐胁迫适应能力。     

甘薯响应冷胁迫的转录组分析

甘薯是世界第七大粮食作物,虽然甘薯对盐、旱等胁迫的抗逆性较强,但因其起源热带而对低温耐受性差,目前低温已经成为影响甘薯产量及限制甘薯种植面积的重要非生物胁迫。本研究以耐冷品种辽薯36为试验材料,在低温处理0 h、6 h、24 h和48 h后,利用转录组测序技术筛选与冷胁迫紧密相关的代谢通路和差异表达基因。结果共获得甘薯转录组数据75.01 Gb,与参考基因组比对效率达76.33%。差异表达基因分析结果表明,甘薯响应冷胁迫的过程中有6282个基因上调表达,3881个基因下调表达,其中283个共有基因在低温处理6 h、24 h和48 h上调表达,26个基因在3个处理中均下调表达。GO功能富集分析发现,差异表达基因被大量富集到蛋白激酶活性、过氧化物酶活性、氧化还原酶活性、蛋白质磷酸化、氨基酸跨膜转运、蛋白质生色团连锁及果胶分解代谢分子功能和生物学过程等分类条目中。KEGG代谢通路富集分析表明,苯丙烷生物合成、淀粉和蔗糖代谢、植物激素信号转导、植物-病原相互作用、谷胱甘肽代谢这5个代谢通路富集了大量的差异表达基因。本研究为甘薯耐冷基因克隆及耐冷机制解析提供了理论依据。     

不同甘蔗品种叶片气孔对水分胁迫的响应

干旱是甘蔗面临最主要的环境胁迫之一,为了解不同甘蔗品种在干旱胁迫时的气孔响应,该研究以F172、GT21、YT93/159和 YL6四个抗旱性有显著差异的甘蔗品种为材料,采用桶栽,在伸长期进行四种不同程度的干旱胁迫(不浇水)处理：土壤持水量在①65％~70％为轻度干旱;②45％~50％为中度干旱;③25％~30％为重度干旱;④以土壤含水量为75％为对照(CK).检测不同品种不同处理甘蔗的叶片相对持水量变化,并利用扫描电镜技术观察甘蔗叶片下表皮气孔特性.结果表明：在干旱胁迫下,四个甘蔗品种叶片气孔导度急剧下降,重度干旱时耐旱性强的 F172和 GT21的气孔导度低于耐旱性弱的 YT93/159和 YL6的;复水后3 d,F172和 GT21的气孔导度上升至82．07和88．85 mmol·m-2·s-1,而 YT93/159和 YL6的仅有18．88和33．08 mmol·m-2·s-1.干旱还导致气孔下陷、闭合,气孔器的长、宽明显减小,且品种间气孔器长度变化差异显著;干旱胁迫下气孔密度增大,尤以耐旱性最强的 F172在重度干旱时达到显著差异.重度干旱时 F172与GT21的气孔闭合百分比是 YT93/159和 YL6近3~4倍.在水分胁迫下,叶片相对含水量降低,但 F172和GT21在重度干旱时仍可以保持相对较高的含水量,其它两个品种相对较低,尤以 YT93/159的最低.在复水后叶片含水量都有所恢复.这些研究结果表明不同甘蔗品种抗旱能力与叶片气孔特性和含水量密切相关.     

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

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