玉米种胚内质网胁迫相关基因对人工老化处理的响应

【目的】在植物中,内质网胁迫（endoplasmic reticulum stress,ERS）和未折叠蛋白应答（unfolded protein response,UPR）参与环境胁迫响应过程,然而,玉米种子老化过程中内质网胁迫相关基因表达情况尚未见报道。文章利用基因数字表达谱技术探究玉米种子老化过程中内质网胁迫相关基因表达规律,以期为揭示种子衰老的分子机制提供理论依据。【方法】以玉米杂交种郑单958种子为材料,采用高温（45℃）高湿（相对湿度100%）的方法进行人工老化处理。分别提取未老化处理（对照）和老化处理3 d的玉米种胚总RNA,利用Illumina HiSeqTM 2000平台进行高通量测序。去除原始数据中的接头序列、包含模糊碱基的序列以及低质量序列,获得Clean reads,利用短序列比对软件SOAPaligner/ SOAP2将Clean Reads分别比对到玉米参考基因组和参考基因序列,采用RPKM（reads per kb per million reads）方法计算基因的表达量,根据FDR（false discovery rate）＜0.001和|log2 ratio(T/CK)|≥1的标准筛选差异表达的基因,对获得的差异表达基因（differentially expressed genes,DEGs）进行KEGG（kyoto encyclopedia of genes and genomes）数据库功能注释分析,筛选出响应人工老化的内质网胁迫相关差异表达基因。利用qRT-PCR技术定量分析内质网胁迫相关基因在不同人工老化时间内的表达特性。【结果】基因数字表达谱鉴定结果表明,有104个差异表达基因在人工老化过程中参与内质网蛋白质加工（protein processing in endoplasmic reticulum）通路,其中内质网胁迫相关基因有97个（81个上调表达,16个下调表达）。对差异表达基因功能注释分析表明,内质网胁迫的标志性蛋白基因BiP以及分子伴侣蛋白基因CRT、CNT和GRP94等显著上调表达。参与内质网相关性降解（endoplasmic reticulum-associated degradation,ERAD）途径的有83个差异表达基因（70个上调,13个下调）,其中启动ERAD途径的关键酶基因EDEM （ER degradation enhancing mannosidase I-like protein）下调,参与蛋白泛素化的E2泛素结合酶基因UbcH5、E3泛素连接酶基因Hrd1和Doa10等也发生显著的表达变化。qRT-PCR结果表明,内质网胁迫相关基因在不同人工老化时间内表现表达多样性和复杂性。【结论】人工老化处理能造成玉米种胚细胞发生内质网胁迫。细胞通过上调分子伴侣基因表达和诱导ERAD途径响应内质网胁迫,但ERAD途径受阻可能引起错误折叠蛋白聚集,从而进一步加剧细胞损伤,最终导致种子活力降低甚至丧失。     

低氧—复氧胁迫下脊尾白虾鳃组织差异表达基因趋势分析

【目的】探究脊尾白虾（Exopalaemon carinicauda）鳃组织在渐变式低氧—复氧胁迫下的分子调控机制,为今后开展脊尾白虾耐低氧品系（种）的选育提供理论指导。【方法】通过模拟自然低氧环境的形成过程,分别于低氧处理0（对照）、3和6 h及复氧后1和8 h采集脊尾白虾鳃组织,利用Illumina HiSeqTM 4000测序平台进行转录组测序分析,经过滤和Trinity组装获得Unigenes,选取Nr、Swiss-Prot、KEGG和COG/KOG等数据库进行注释分析,在Omicsmart平台上完成差异表达基因筛选及其表达趋势分析,然后进行GO功能注释分析和KEGG信号通路富集分析,并随机选取5个差异表达基因进行实时荧光定量PCR验证。【结果】脊尾白虾鳃组织样本转录组测序数据经过滤后的Cleanreads进行Trinity组装共获得93227条Unigenes,其长度范围在201~35402 bp,平均长度为834 bp,N50长度为1352 bp。通过组间两两比较分析,共鉴定出4750个差异表达基因,其中上调差异表达基因3557个、下调差异表达基因2829个;超过50%的差异表达基因被显著富集到6种基因表达趋势模式中（P<0.01）,具体表现为:Profile 0模式富集到415个基因,Profile 5模式富集到201个基因,Profile 11模式富集到371个基因,Profile 13模式富集到841个基因,Profile 17模式富集到387个基因,Profile 18模式富集到411个基因。6种基因表达趋势模式中的差异表达基因被注释到代谢进程、细胞进程、单一有机体进程、细胞、细胞零件、大分子复合物及催化活性等GO功能条目上;而KEGG信号通路富集分析结果显示,以Profile 13模式中的差异表达基因富集到最多信号通路（86条）,其中呈显著富集的有8条,分别为核糖体、碳代谢、氧化磷酸化、氨基酸生物合成、内质网蛋白质加工、糖酵解/糖异生、谷胱甘肽代谢和蛋白输出。【结论】脊尾白虾鳃组织在受低氧胁迫早期通过合成蛋白质及提高代谢能力来抵御低氧环境,随着低氧胁迫时间的延长,物质合成和能量代谢活动均显著下降;但在复氧后随着复氧时间的延长,其蛋白质合成和能量代谢水平又逐渐升高恢复至常氧水平。     

1-MCP处理采后红阳猕猴桃果实生理效应的研究

【目的】为探究1-甲基环丙烯（1-MCP）处理结合冷藏（1～4℃）贮藏采后‘红阳’猕猴桃果实条件下保鲜的分子调控机制。【方法】以‘红阳’猕猴桃果实为材料，采用1-MCP配合冷藏进行红阳猕猴桃保鲜，研究贮藏过程中果实生理效应的变化。【结果】1-MCP处理结合冷藏保藏可以减缓猕猴桃果实软化速度，抑制TSS、可滴定酸、VC及花色苷的降解，减缓MDA的积累，维持了SOD的活力，降低果实中乙烯释放量。利用RNA-Seq技术进行转录组学比较分析，共筛选获得2 380个差异表达基因，其中上调表达基因有1 503个，下调表达基因877个。GO富集分析表明，差异表达基因显著富集在细胞组分、分子功能和生物过程中的肽生物合成与代谢、酰胺生物合成、蛋白质运输、防御反应等功能类别。KEGG富集分析显示，上调基因主要富集于植物激素信号转导通路、半胱氨酸和蛋氨酸代谢通路以及苯丙氨酸代谢、谷胱甘肽代谢等抗氧化等相关通路，推测蛋氨酸代谢表达的增加是减少乙烯合成的主要原因之一。下调基因主要富集在丙酮酸代谢通路、叶绿素代谢、烟酸和烟酰胺代谢以及蛋白酶体代谢通路等影响呼吸作用强度的通路。【结论】1-MCP结合冷藏处理‘红阳’...     

鸭胚胎发育中后期胸肌发育阻滞的RNA-seq分析

【目的】选择两个中国地方品种高邮鸭和金定鸭,对鸭胚胎发育中后期胸肌进行转录组分析研究,旨在探明胸肌发育阻滞的分子变化机制,为鸭骨骼肌调控机理研究打下基础。【方法】在21胚龄和27胚龄两个时间点,分别解剖高邮鸭、金定鸭各3只,采集胸大肌,提取总RNA构建文库,利用Illumina的HiseqTM2000进行高通量测序,并利用生物信息学方法进行差异表达基因挖掘、基因功能注释等分析,探讨21胚龄和27胚龄两个时间点之间胸肌发育阻滞的分子机制。【结果】高邮鸭、金定鸭21胚龄和27胚龄胸大肌组织RNA-seq质量Q20均在94%以上,Q30均在89%以上,测序得到的结果可靠,可用于后续分析。RNA水平相关性检查和基因mRNA表达量聚类图结果都表明,21胚龄(27胚龄)高邮鸭和金定鸭之间表达模式的相关性高于高邮鸭(金定鸭)21胚龄和27胚龄之间的表达模式。不同品种内时间点之间的差异基因数量(高邮鸭6 128个,金定鸭6 452个)远多于同一时间点不同品种间的差异基因数量(21胚龄522个,27胚龄299个)。qRT-PCR验证试验结果与RNA-seq分析结果相关性较强。通过GO和KEGG富集分析发现,高邮鸭、金定鸭胸肌在21胚龄到27胚龄阶段,能量代谢相关基因(主要为辅酶Q相关基因、ATP酶合成相关基因和细胞色素C相关基因)均显著上调,DNA复制和细胞周期相关基因(主要为微型染色体维持蛋白(MCM)相关基因、复制因子C(RFC)相关基因)均显著下调。相关基因表达的变化可能与此阶段成肌细胞增殖速度减慢,逐渐退出细胞周期开始准备下一阶段融合成多核肌管并形成肌纤维有关。对肌肉生长发育相关的关键基因分析发现,促进肌肉生长的IGF1和诱导成肌细胞末端分化的MyoG显著下调,促进肌纤维分化融合的MUSTN1基因、诱导肌祖细胞向成肌细胞转化的MyoD1基因显著上调。【结论】鸭胚胎中后期胸肌发育过程中大量基因差异表达。其中能量代谢相关基因的上调和DNA复制、细胞周期相关基因的下调,以及肌肉发育相关基因MUSTN1显著上调,IGF1、MyoG等显著下调,可能与鸭胚胎中后期胸肌发育阻滞现象密切相关。     

百子莲胚性愈伤组织在超低温过程中的 蛋白质组学研究

【目的】从蛋白质组学层面揭示植物胚性细胞对超低温保存复合逆境的响应调控机制,筛选出百子莲胚性愈伤组织在超低温保存中的差异表达蛋白,进一步丰富观赏植物超低温保存的分子生物学理论。【方法】利用iTRAQ技术鉴定百子莲胚性愈伤组织在超低温保存四个关键步骤(预培养、玻璃化脱水处理、洗涤和恢复培养)中的差异表达蛋白,应用生物信息学分析方法构建细胞响应超低温保存复合逆境的蛋白调控网络,筛选出植物超低温逆境保护类蛋白。【结果】共鉴定到2730种蛋白,筛选出821个差异表达蛋白。GO功能注释分析结果表明差异表达蛋白主要参与代谢、细胞过程和刺激应答等生物学过程,而催化活性和分子结合是最主要的分子功能类别。KEGG通路分析结果表明差异表达蛋白主要富集于能量代谢、蛋白的合成与降解、信号转导及细胞质膜修复等通路。通过STRING软件构建差异表达蛋白互作调控网络,发现上调表达蛋白主要参与糖代谢、氧化还原过程和核糖体翻译过程,下调表达蛋白主要参与ROS代谢、DNA修复、抑制细胞程序性死亡过程和细胞骨架。此外,筛选到96个与逆境保护相关的差异表达蛋白,分别参与能量代谢、氧化还原平衡、蛋白质合成与降解和信号转导。【结论】在超低温保存过程中,百子莲胚性愈伤组织的能量代谢、氧化应激与蛋白降解及PCD有所加剧,而细胞的ROS清除和胁迫防御过程明显下降,即百子莲胚性愈伤组织在超低温复合逆境下形成以糖代谢和抗氧化系统为主的抗逆保护机制。     

基于转录组测序对紫花苜蓿细胞质雄性不育系相关代谢通路的鉴定

【目的】研究紫花苜蓿(Medicago sative L.)细胞质雄性不育(CMS)的机理,为紫花苜蓿利用不育系的育种工作提供理论基础。【方法】本研究以紫花苜蓿细胞质雄性不育系MSGN1A及其保持系MSGN1B为材料,进行转录组测序(RNA-Seq)以及基因功能注释,筛选与不育系机理相关的差异表达基因及代谢通路。【结果】Illumina测序共得到紫花苜蓿95 679条功能基因,经过差异表达分析筛选出187个显著差异表达的基因(DEGs),其中包括18个上调基因,169个下调基因。采用GO、KEGG、COG注释库分别对187个DEGs进行功能注释,其主要涉及催化活性、代谢活动、信号传导机制、合成、膜功能、碳水化合物运输和代谢、能量代谢等相关通路。【结论】筛选出与紫花苜蓿细胞质雄性不育性相关的代谢通路,初步探讨了紫花苜蓿细胞质雄性不育系的分子机理。     

东方蜜蜂微孢子虫孢子中长链非编码RNA的竞争性内源RNA调控网络及潜在功能

【背景】 长链非编码RNA（long non-coding RNA,lncRNA）是一类长度>200 nt且不具有编码能力的转录本,能够在真核生物的转录水平和转录后水平通过顺式（cis）、反式（trans）或竞争性内源RNA（competing endogenous RNA,ceRNA）机制调控基因的表达,已被证明在生物体的生长和发育等生物学过程中发挥重要的调控作用。【目的】 结合测序得到的东方蜜蜂微孢子虫（Nosema ceranae）纯化孢子的small RNA（sRNA）组学数据和前期获得的lncRNA组学数据对lncRNA的调控方式进行深入细致的分析和探讨,揭示lncRNA在东方蜜蜂微孢子虫孢子中的潜在功能。【方法】 利用small RNA-seq（sRNA-seq）技术对东方蜜蜂微孢子虫纯化孢子进行测序,通过相关的生物信息学软件对测序数据进行质控。根据lncRNA与mRNA的位置关系预测lncRNA的上下游基因;使用Blast软件比对GO和KEGG数据库,对上下游基因进行功能和代谢通路注释;利用Target Finder软件预测lncRNA-miRNA及lncRNA-miRNA-mRNA调控网络,进而通过Cytoscape v3.7.2软件对调控网络进行可视化;利用RT-PCR对调控网络中的部分lncRNA、miRNA和mRNA的表达进行验证。【结果】 东方蜜蜂微孢子虫纯化孢子样品的sRNA-seq分别得到16 597 883、15 451 791和12 248 316条raw reads,经质控分别得到15 608 370、14 249 255和11 440 684条clean reads,各样品的Q30≥98.58%。共预测出lncRNA的上下游基因310个,它们可注释到代谢进程、细胞进程、催化活性、结合和细胞等相关的35个功能条目;此外还可注释到56条代谢通路,涉及嘌呤代谢、碳代谢和丙酮酸代谢等物质代谢通路,以及甲烷代谢、氧化磷酸化和糖酵解/糖异生等能量代谢通路。上述结果表明相应的lncRNA可能通过cis作用调节上下游基因的表达,从而参与调控东方蜜蜂微孢子虫孢子中的物质和能量代谢以及细胞生命活动。进一步构建lncRNA-miRNA调控网络,分析结果显示MSTRG.3636.1、MSTRG.4498.1和MSTRG.4883.1可靶向结合4个miRNA（nce-miR-7729、nce-miR-7502、nce-miR-8639和nce-miR-8565）,说明此3个lncRNA可作为ceRNA在东方蜜蜂微孢子虫孢子中发挥潜在作用。LncRNA-miRNA-mRNA调控网络分析结果显示,三者之间形成较为复杂的调控网络关系,miRNA处于调控网络的中心并联系lncRNA和mRNA,nce-miR-7502和nce-miR-8639结合的mRNA最多,达到28个;与MSTRG.3636.1、MSTRG.4498.1和MSTRG.4883.1存在靶向关系的miRNA结合多个mRNA,表明上述3个lncRNA可能通过ceRNA机制发挥作用,从而影响东方蜜蜂微孢子虫孢子中的新陈代谢和生命活动。【结论】 lncRNA可能通过cis作用调控上下游基因的表达,以及作为ceRNA通过吸附miRNA间接影响靶基因的表达,从而参与调节东方蜜蜂微孢子虫孢子中的物质和能量代谢等基本生命活动。     

谷子萌发吸水期关键代谢途径的筛选与分析

【目的】谷子适应性强,抗旱耐瘠,是起源于中国的重要作物。通过转录组测序技术分析谷子萌发不同吸水期的转录组差异,以期获得谷子萌发过程中的差异表达基因,寻找调控谷子萌发的重要代谢途径和代谢物。【方法】以晋谷20为材料,构建谷子萌发过程中开始快速吸水期、滞缓吸水期和重新大量吸水期的cDNA文库,进行转录组分析;采用K-Means开展基因表达聚类分析;利用DESeq筛选差异表达基因;通过COG、GO、KEGG等对差异表达基因进行功能注释;利用KEGG富集挖掘不同吸水期调控种子萌发的关键代谢途径和关键基因;并采用qRT-PCR验证其可靠性;用HPLC分析关键代谢物含量。【结果】转录物测序分析获得谷子萌发开始快速吸水期、滞缓吸水期和重新大量吸水期覆盖整个基因组的基因表达谱,共获得33 643个基因,识别9个具有不同表达模式的共表达基因簇。比较种子萌发的开始快速吸水期与滞缓吸水期、滞缓吸水期与重新大量吸水期、开始快速吸水期与重新大量吸水期,分别筛选出3 893、4 612和8 472个差异表达基因。KEGG富集分析表明,3个比较的差异表达基因都显著富集到phenylpropanoid biosynthesis、phenylalanine metabolism、starch and sucrose metabolism代谢途径;开始快速吸水期与滞缓吸水期、开始快速吸水期与重新大量吸水期的差异表达基因还显著富集到plant hormone signal transduction途径。并且3个比较中富集到phenylpropanoid biosynthesis和phenylalanine metabolism代谢途径的差异表达基因数都最多,其中过氧化物酶基因（peroxidase）比例最高。通过qRT-PCR对4个苯丙烷生物合成途径相关基因的分析表明,其表达趋势与转录组分析结果基本一致,其中,4-香豆酸-CoA连接酶3（4-coumarate-CoA ligase 3）在谷子种子中存在已形成mRNA,萌发吸水过程中呈先下调后上调再下调的表达趋势。苯丙烷类相关代谢物含量分析显示,芥子酸在种子中大量储备,在萌发过程中呈下调趋势;阿魏酸、对香豆酸和咖啡酸呈先上调后下调趋势。【结论】谷子萌发过程中,不同吸水期的差异表达基因显著与苯丙烷生物合成途径和苯丙氨酸代谢途径相关;其上游基因4-香豆酰-辅酶A连接酶和下游基因过氧化物酶家族成员在谷子萌发响应水分过程中发挥调控作用;中间产物芥子酸可能参与种子的休眠与萌发。     

谷子萌发期响应干旱胁迫的基因表达谱分析

【目的】谷子是一种耐旱作物,通过二代测序技术获得大量的谷子萌发期响应干旱胁迫的差异基因,进而挖掘谷子萌发期抵御干旱的关键基因及其相关的分子机制。【方法】以晋谷45为材料,谷子萌发期分别用18%PEG-6000干旱胁迫（处理组）和蒸馏水（对照组）处理种子并测定1、10和18 h种子的SOD、POD和CAT活性。SOD活性用氮蓝四唑（NBT）法测定,POD活性用愈创木酚法测定,CAT活性用比色法测定;对萌发10和18 h种子的对照组和处理组构建cDNA文库并进行差异基因表达谱分析;利用Bowtie将reads比对到参考基因组,采用RSEM对bowtie的比对结果进行表达量估计;使用DESeq进行差异表达基因分析;利用NR、Swiss-Prot、KEGG、COG和GO在线数据库对差异基因进行功能注释,挖掘调控谷子萌发的关键基因;利用qRT-PCR验证测序结果的可靠性。【结果】处理组的SOD活性整体比对照组高,而POD活性和CAT活性与之相反;随着萌发时间的变化,SOD活性在不断地增加,但CAT和POD活性逐渐减小。基因表达谱序列与所选参考基因组序列高度一致,基因表达呈现出高度不均一性。通过高通量测序最后获得35 470个基因,以RPKM≥0.01为筛选标准,对照样本中分别筛选出24 030和24 486个表达基因,PEG干旱胁迫处理10和18 h的样本分别筛选出24 019和23 877个表达基因;差异表达基因分析表明,谷子萌发10和18 h分别筛选出456和545个差异基因,其中87和267个上调表达基因,369和278个下调表达基因;GO功能显著性富集分析表明,差异基因主要涉及代谢过程,细胞进程和响应刺激;KEGG富集分析表明,差异基因参与到苯丙烷代谢和植物激素信号转导过程;通过qRT-PCR对5个差异基因在干旱胁迫下种子萌发时的表达分析表明,其表达趋势与表达谱分析结果基本一致。【结论】差异表达基因广泛涉及到糖、蛋白质、核酸等生物大分子代谢、次生代谢和能量代谢等过程;SnRK2和PAL可能在干旱胁迫下调节种子的萌发。     

静原鸡肌肉组织肌苷酸特异性沉积相关LNC_003828- gga-miR-107-3p-MINPP1的关联分析

【目的】探究静原鸡肌肉组织肌苷酸沉积过程中关键调控因子的调节作用,利用lncRNA-miRNA-mRNA关联分析鉴定与肌苷酸特异性沉积相关的LNC_003828、gga-miR-107-3p和MINPP1,其作为肉质研究的候选基因,为分子辅助育种提高肌肉品质提供理论基础。【方法】测定15只静原鸡胸肌和腿肌的肌苷酸含量,筛选高肌苷酸含量的胸肌和低肌苷酸含量的腿肌各3个样本提取总RNA,质量检测合格后构建cDNA文库、PCR扩增,利用Agilent 2100对文库质量进行评价,库检合格后送Illumina-Hiseq平台进行转录组测序。利用生物信息学方法筛选出静原鸡肌肉组织不同部位差异表达的MINPP1、gga-miR-107-3p和LNC_003828,进行GO注释和蛋白互作网络分析MINPP1的功能。采用qRT-PCR方法检测LNC_003828、gga-miR-107-3p和MINPP1在静原鸡胸肌和腿肌组织中的表达情况,并分析其与肌苷酸含量的相关性。【结果】测序样品间基因表达水平相关性R²>0.9,即试验样本之间基因表达可用于后续的差异基因分析。参与肌苷酸合成和代谢的糖酵解/糖异生途径中检测出3个差异表达基因MINPP1、PKM和ALDH9A1。互作分析发现lncRNA-miRNA-mRNA网络图中共有17个miRNA（9个上调,8个下调）、44个mRNA（16个上调,28个下调）和155个lncRNA（68个上调、87个下调）,核心节点gga-miR-107-3p互作的靶基因有MINPP1、靶lncRNA有LNC_003828。GO富集分析发现MINPP1基因具有磷酸酶活性、双磷酸甘油酸酯磷酸酶活性等功能;蛋白互作网络中MINPP1基因与参与糖酵解/糖异生和氨基酸生物合成通路中的PGAM1、ENO1、BPGM基因均有互作关系。qRT-PCR结果表明,静原鸡胸肌LNC_003828和gga-miR-107-3p的相对表达量低于腿肌,但差异不显著;胸肌MINPP1的相对表达量显著低于腿肌（P<0.05）。静原鸡胸肌和腿肌组织中gga-miR-107-3p的表达量与LNC_003828表达量均呈正相关,与MINPP1的表达量均呈负相关。胸肌和腿肌组织中LNC_003828、gga-miR-107-3p的表达量与肌苷酸含量均呈正相关,且差异均不显著;胸肌MINPP1表达量与肌苷酸含量呈负相关,腿肌MINPP1表达量与肌苷酸含量呈显著负相关（P<0.05）。综上所述,推测静原鸡肌肉组织中gga-miR-107-3p作为核心调节因子吸附LNC_003828,影响MINPP1基因调控肌肉肌苷酸特异性沉积,从而改善肉质。【结论】筛选出LNC_003828、gga-miR-107-3p和MINPP1为影响肌苷酸特异性沉积的候选调控因子。     

意大利蜜蜂响应东方蜜蜂微孢子虫胁迫的免疫应答

【目的】通过对意大利蜜蜂（Apis mellifera ligustica,简称意蜂）工蜂中肠响应东方蜜蜂微孢子虫（Nosema ceranae）胁迫的差异表达基因（differentially expressed gene,DEG）及免疫通路进行深入细致的分析,揭示宿主的细胞和体液免疫应答,为关键免疫应答基因的筛选及功能研究打下基础。【方法】基于前期获得的正常及东方蜜蜂微孢子虫胁迫的意蜂7日龄和10日龄工蜂中肠（Am7CK、Am7T、Am10CK、Am10T）转录组数据,按照FDR≤1,P≤0.05和|log₂ fold change|≥1的标准筛选出各组的DEG,利用相关生物信息学软件对DEG进行Pearson相关性分析、Venn分析、GO分类和KEGG代谢通路富集分析,进而对免疫通路富集基因进行统计和分析,利用实时荧光定量PCR（real-time quantitative PCR,RT-qPCR）验证转录组数据的可靠性。【结果】 差异分析结果显示,Am7CK vs Am7T比较组包含472个上调基因和385个下调基因;Am10CK vs Am10T比较组包含611个上调基因和360个下调基因。Venn分析结果显示,两个比较组特有的DEG分别为739和853个,共有的DEG为118个。GO分类结果显示,Am7CK vs Am7T比较组中上调和下调基因分别涉及23和29个功能条目,其中富集上调基因数最多的前5位分别为结合、催化活性、代谢进程、细胞进程和单组织进程;富集下调基因数最多的前5位分别为代谢进程、单组织进程、催化活性、细胞进程和结合。Am10CK vs Am10T比较组中上调和下调基因分别涉及36和26个功能条目,其中富集上调基因数最多的前5位分别为单组织进程、结合、细胞进程、催化活性和代谢活性;富集下调基因数最多的前5位分别为结合、细胞进程、催化活性、代谢进程和单组织进程。KEGG代谢通路富集分析结果显示,Am7CK vs Am7T比较组中上调和下调基因分别富集在38和33条代谢通路,富集上调基因数最多的前5条分别为胆汁分泌、内质网蛋白加工、泛素介导的蛋白水解、PI3K-Akt信号通路和神经营养因子信号通路;富集下调基因数最多的前5条分别为胞质DNA传感通路、嘌呤代谢、嘧啶代谢、RNA聚合酶和核糖体;涉及泛素介导的蛋白水解等3条细胞免疫通路,以及PI3K-Akt信号通路等7条体液免疫通路。Am10CK vs Am10T比较组中上调和下调基因分别富集在54和43条代谢通路,富集上调基因数最多的前5条分别为Hippo信号通路、药物代谢-细胞色素P450、P450对外源物质的代谢、泛素介导的蛋白水解和鞘脂类代谢;富集下调基因数最多的前5条分别为mRNA监测、鞘脂类信号通路、果糖和甘露糖代谢、半乳糖代谢和鞘脂类代谢;涉及泛素介导的蛋白水解等7条细胞免疫通路,以及NF-κB信号通路等2条体液免疫通路。RT-qPCR验证结果显示6个随机挑选的DEG的表达量变化趋势与测序结果一致,证实了本研究中测序数据的可靠性。进一步分析发现意蜂7日龄和10日龄工蜂中肠的NF-κB信号通路均被东方蜜蜂微孢子虫激活,随即启动了3种抗菌肽apidaecin、defensin-1和hymenoptaecin的合成,体现了它们在宿主抵御东方蜜蜂微孢子虫入侵中的重要性。【结论】在转录组水平解析了意蜂工蜂中肠对东方蜜蜂微孢子虫的免疫应答,揭示宿主在胁迫前期同时作出细胞和体液免疫应答,前者可能在抵御病原入侵方面发挥主要作用;宿主的细胞免疫在胁迫后期持续增强,但体液免疫较大程度地减弱;泛素介导的蛋白水解通路及富集DEG,NF-κB信号通路及富集DEG,以及抗菌肽编码基因apidaecin、defensin-1和hymenoptaecin可能在意蜂工蜂对东方蜜蜂微孢子虫的免疫应答中起到关键作用。     

利用RNA-Seq发掘玉米叶片形态建成相关的调控基因

【目的】叶片宽度和长度等叶形特性是决定植株形态,进而影响种植密度的重要农艺性状,通过转录组测序技术筛选并挖掘玉米叶片形态建成相关的代谢路径及调控基因,为深入认识叶片发育的分子机理和鉴定叶宽、叶长候选基因奠定基础。【方法】以极端窄叶自交系NL409和宽叶自交系WB665为材料,利用RNA-Seq技术鉴定7叶期第七片叶近基部的差异表达基因（DEGs）,通过生物信息学分析,筛选与叶片发育密切相关的代谢通路,利用qRT-PCR验证不同激素路径叶形相关基因的表达结果,并结合启动子区域的序列差异挖掘叶形功能基因。【结果】分析对照（WB665）和样品（NL409）高通量测序结果,在叶宽形成关键部位共筛选出5 199个DEGs,其中,2 264（43.55%）个基因表达上调,2 935（56.45%）个基因下调表达,下调基因明显多于上调基因;GO功能富集分析表明,差异基因主要富集在细胞膜相关的细胞组分中,涉及代谢过程和细胞响应刺激;KEGG富集分析表明,差异基因主要参与到核糖体、植物激素信号转导、苯丙烷类代谢、乙醛酸和二羧酸代谢等过程,其中核糖体、植物激素信号转导、鞘脂类代谢下调表达基因较多的路径与叶片发育密切相关。核糖体路径富集到多个PRS（PRESSED FLOWER）基因,分析发现PRS13（PFL2）可能在调控窄叶发育过程中发挥重要作用。鞘脂代谢路径富集的基因几乎全部下调表达,引起抑制叶片发育的AP1（APETALA1）类和MAPK（Mitogen-Activated Protein Kinase）类基因上调,以及促进叶片发育的LFY（LEAFY）下调,与窄叶发育受抑制的表型一致。植物激素信号转导路径富集到的油菜素内酯（BR）响应基因和赤霉素（GA）代谢基因下调,细胞分裂素（CTK）和大部分生长素（Auxin）响应基因上调,与窄叶中DELLA蛋白基因上调表达,抑制GA并促进CTK基因表达的作用模式一致。通过qRT-PCR对18个叶片发育相关基因进行分析,结果表明,其表达趋势与转录组结果一致,分析发现BR相关的ROT3、Auxin相关的NAL7-like、AGO7-like以及TCP类转录因子CYC/TB1等基因与窄叶的形成密切相关。【结论】明确了一些与玉米叶片发育密切相关的代谢路径,还发现植物激素间的动态平衡对叶片发育有着重要影响,尤其是生长素与油菜素内酯、细胞分裂素与赤霉素之间的相互作用对调控叶片形态可能发挥重要作用。     

拟轮枝镰孢与玉米籽粒互作的差异基因筛选及代谢通路分析

【目的】拟轮枝镰孢(Fusarium verticillioides)引起的玉米穗腐病是我国玉米产区发生严重的病害之一，论文旨在了解病原菌与玉米籽粒互作过程中的基因表达差异，为揭示病原菌的致病机制和玉米抗病机制提供依据。【方法】对拟轮枝镰孢侵染玉米籽粒0、4、12和72 h的样品进行转录组测序，之后采用生物信息学分析，分别以玉米和拟轮枝镰孢基因组为参考，以|log₂FC|≥1,P-adjust<0.05为阈值筛选互作过程中玉米和拟轮枝镰孢的差异表达基因，利用GO和KEGG对其进行功能注释及富集分析。Goatools软件分析植物-病原互作、MAPK途径和植物激素信号转导通路相关差异基因的表达变化，采用实时荧光定量PCR(qRT-PCR)方法对测序差异基因进行验证。【结果】在互作4、12和72 h后拟轮枝镰孢分别有140、400和1 945个基因上调表达，有9、302和1 784个基因下调表达；玉米分别有293、692和1 426个基因上调表达，320、482和153个基因下调表达。GO和KEGG富集分析显示，侵染早期拟轮枝镰孢在细胞间隙生长，差异基因主要富集在...     

利用RNA-Seq鉴定甘蓝型油菜叶片干旱胁迫应答基因

【目的】利用RNA Sequencing（RNA-Seq）技术比较2种不同生长条件下甘蓝型油菜苗期叶片转录组,鉴定油菜叶片干旱胁迫应答相关基因,从转录组水平揭示油菜适应干旱胁迫环境的分子机制。【方法】提取正常生长（ZY）和自然失水处理（ZY8D）的六叶期甘蓝型油菜中油821的叶片总RNA,以Illumina Hiseq 2000平台进行RNA-Seq分析。利用NGSQCTookit v2.3.3去除低质量和包含模糊碱基的reads。以甘蓝型油菜亲本物种白菜染色体v1.5和甘蓝Scaffold v1.0为参考序列,采用TopHat2-Cufflinks-Cuffmerge-Cuffdiff标准流程进行差异表达基因（differential expressed genes,DEGs）筛选。对上调和下调DEGs分别采用Cytoscape v3.1.0中的BiNGO和KOBAS2.0进行基因本体（gene ontology,GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes,KEGG）代谢途径富集分析。选择上调和下调DEGs各3个,以实时荧光定量PCR（quantitative real-time PCR,qRT-PCR）验证RNA-Seq结果的可靠性。【结果】过滤低质量reads后,ZY和ZY8D分别保留了26 192 312和28 378 899对高质量reads用于DEGs筛选,其中86.6%和85.8%的reads能准确比对到参考序列上,说明RNA-Seq结果和参考序列可靠。DEGs鉴定结果表明3 657个基因受干旱胁迫诱导差异表达,其中上调表达基因1 431个,下调表达基因2 226个。GO富集分析发现上调表达基因主要与非生物胁迫响应和化学刺激响应相关,其中,参与水分胁迫响应和脱落酸（abscisic acid,ABA）刺激响应的基因分别有127和141个,而下调表达基因与植物病原菌防御、蛋白激酶活性和水杨酸（salicylic acid,SA）刺激相关。KEGG富集分析表明上调表达基因主要富集于苯丙烷和类胡萝卜素的生物合成及淀粉与蔗糖代谢途径,而下调表达基因主要富集于植物-病原菌互作和植物激素ABA、SA和茉莉酸（jasmonic acid,JA）信号转导途径。qRT-PCR检测6个DEGs的表达模式与RNA-Seq分析结果一致,证实了RNA-Seq结果的可靠性。【结论】RNA-Seq分析鉴定出3 657个甘蓝型油菜叶片干旱胁迫应答基因。GO和KEGG代谢途径分析明确了差异表达基因富集的分子功能与代谢途径。     

玉米pTAC2影响苗期叶片叶绿素合成的转录组分析

【目的】叶绿素是参与光合途径最为重要的光合色素。叶绿体的发育及叶绿素的合成在很大程度上依赖于质体基因组与核基因组之间的双向信号传导来精确协调基因表达。通过对白化表型的CRISPR/Cas9-ZmpTAC2转基因阳性纯合突变材料进行RNA-seq研究,筛选和鉴定参与叶绿素合成的相关基因,为明确叶绿素的合成途径奠定基础。【方法】以CRISPR/Cas9-ZmpTAC2玉米转基因编辑纯合突变株系为研究材料,使用透射电镜观察叶绿体超微结构和分光光度法测定叶片叶绿素含量,确定叶绿体发育状态及叶绿素合成情况。对转基因阴性材料（CK）和CRISPR/Cas9-ZmpTAC2转基因纯合编辑材料（zmptac2）苗期叶片取样进行转录组测序。通过生物信息学分析,寻找CK与zmptac2间差异表达的基因;qRT-PCR对差异表达基因进行验证。通过酵母双杂交筛选与玉米pTAC2互作的蛋白质。【结果】共获得15株T₀转基因植株,包括绿色植株（7株）和白色植株（8株）。绿色幼苗中3株为转基因阴性材料,4株为转基因阳性（2株为未编辑,2株为杂合编辑突变）,白色植株(8株)均为转基因阳性纯合编辑。与CK相比,突变体（zmptac2）叶绿体发异常,叶绿素含量显著降低。RNA-seq的结果显示,CK与zmptac2之间共检测到1 367个基因差异表达,其中618个基因上调表达（zmptac2/CK）,749个基因下调表达（zmptac2/CK）。GO富集分析显示,下调基因显著富集到叶绿体和质体中。KEGG分析表明下调表达基因显著富集在苯丙氨酸代谢、酪氨酸代谢和异喹啉生物碱生物合成等途径。选取的15个差异基因表达模式均与测序数据相一致,表明测序结果是可靠的。与CK相比,zmptac2中依赖PEP（plastid-encoded RNA polymerase）转录的基因表达量显著降低,而依赖NEP（nuclear gene-encoded RNA polymerase）转录的基因表达量则显著上升。通过对玉米cDNA文库筛选和互作验证,鉴定出ZmpTAC3与ZmpTAC2存在互作。【结论】ZmpTAC2突变会导致叶绿体早期生物合成受阻,该基因参与叶绿体发育及叶绿素合成,且该种作用是由ZmpTAC2调控PEP相关基因表达而实现的。     

枯草芽孢杆菌NCD-2对棉花根系分泌物L-脯氨酸响应的转录-蛋白质组学联合分析

【目的】棉花根系分泌物L-脯氨酸是影响生防菌株定殖的关键因素之一,前期研究发现L-脯氨酸能够提高枯草芽孢杆菌（Bacillus subtilis）NCD-2生物膜形成能力。通过高通量测序技术分析L-脯氨酸调控枯草芽孢杆菌NCD-2生物膜形成和生防潜力相关的基因,为深入认识棉花根系分泌物与生防菌株的分子互作关系打下基础。【方法】以枯草芽孢杆菌NCD-2菌株为供试材料,外源添加浓度为10 mg·mL^-1的L-脯氨酸共培养24 h后,分别进行转录组（RNA-seq）和同位素标记相对定量蛋白质组技术（iTRAQ）分析,对所获得的转录组和蛋白质组数据进行生物信息学分析,利用实时荧光定量PCR（RT-qPCR）验证不同代谢通路中部分差异基因的表达情况。 【结果】转录组分析发现,L-脯氨酸和NCD-2菌株共培养后,获得1 071个差异表达基因（DEG）,其中602个基因上调,469个基因下调。GO分析发现在生物过程、细胞组分和分子功能方面分别有49、14和30个功能条目显著富集。KEGG代谢通路主要富集在化合物代谢、鞭毛组装、细菌运动和趋化作用。蛋白质组学分析发现,筛选到211个差异表达蛋白（DEP）,其中118个蛋白上调,93个蛋白下调。GO分析发现在生物过程和分子功能方面分别有13和8个功能条目显著富集。KEGG代谢通路主要富集在氨基酸代谢、碳水化合物类代谢、鞭毛组装和ABC转运蛋白。进一步转录-蛋白质组学联合分析发现,在L-脯氨酸作用下,检测到共有的显著差异表达基因（或蛋白）112个,其中38个基因（或蛋白）下调,74个基因（或蛋白）上调。GO功能显著富集主要集中在营养库活性、催化活性、细胞膜、定位、细胞脂质代谢过程、氧化还原过程、sigma因子活性、转运活性、芽孢形成9个方面。KEGG代谢通路主要富集在能量代谢、ABC转运蛋白、抗生素生物合成、鞭毛组装、运动或趋化作用和双组分系统方面。RT­qPCR验证了26个差异表达显著基因,结果发现在表达量上存在一定的差异,但表达趋势与RNA­seq和iTRAQ组学分析结果基本一致。【结论】棉花根系分泌物L-脯氨酸与枯草芽孢杆菌NCD-2之间的互作存在一个复杂的生物过程,依赖于不同代谢通路网络中的多个基因。双组分系统、抗生素生物合成、能量代谢、运动或趋化作用、鞭毛组装和ABC转运蛋白通路中的差异基因（或蛋白）可能在棉花根系分泌物与枯草芽孢杆菌互作过程方面发挥着重要作用。     

红橘响应褐斑病菌侵染的转录组学分析

【目的】 明确红橘（Citrus reticulata Blanco, cv. Hongjv）接种褐斑病致病菌——链格孢菌橘致病型（Alternaria alternata tangerine pathotype）后基因种类和表达量在转录水平的变化规律,确定红橘响应该致病型侵染的关键基因。【方法】 采用链格孢菌橘致病型接种红橘离体叶片,28 h后选取感病叶片和未接种叶片提取RNA,进行转录组高通量测序,然后利用生物信息学分析,以甜橙基因组为参考,以|log₂ fold change|≥1,q-value≤0.01为阈值选取感病和健康红橘叶片转录组的差异表达基因,应用GO数据库对差异表达基因（differentially expressed gene,DEG）进行功能分类,KEGG分析代谢途径,MapMan软件分析生物胁迫信号通路相关基因的表达变化。采用qRT-PCR方法对测序结果进行验证。【结果】 红橘接种链格孢菌橘致病型28 h后产生大量与胁迫相关的差异表达基因,获得上调差异基因5 173个,下调差异基因6 555个。GO功能分类显示差异基因主要与蛋白结合、膜、氧化还原过程等相关。通过KEGG富集和MapMan软件分析发现,红橘在受链格孢菌橘致病型胁迫的过程中基础代谢被严重破坏。乙烯、水杨酸和生长素等寄主防御相关的植物激素信号转导途径多个基因表达出现差异,其中乙烯起主导作用,乙烯受体ETR 3个成员被不同程度激活,下游激酶和乙烯响应因子均上调,生长素大部分关键信号基因、绝大部分生长素响应因子ARF和水杨酸合成途径的基因均下调表达。同时,黄酮醇、花青素、萜类化合物和生物碱合成相关基因受该菌诱导显著变化,萜类合成中大部分基因下调,而黄酮类合成相关上调基因数量和表达趋势均强于表达下调基因,有抗虫和抑菌作用的硫代葡萄糖苷基因呈上调趋势。进一步研究发现,大量参与抗逆过程的转录因子如WRKY、bZIP、ERF、MYB、NAC被诱导激活,其中大部分WRKY和bZIP转录因子受该菌正向调控,超过50%的ERF家族基因表达上调;在转录因子调控下,PTI及ETI响应基因如受体激酶、R蛋白、NBS抗病蛋白等大量表达,多个PR家族抗菌蛋白基因上调表达,22个抗氧化保护酶系统POD成员基因受到活性氧信号激发大量表达。以上结果表明链格孢菌橘致病型侵染对寄主内部生理状态产生显著影响。选取了19个与植物抗病相关基因进行qRT-PCR分析,其基因表达趋势与测序数据一致。【结论】 获得了红橘响应链格孢菌橘致病型侵染的差异表达基因及显著上调表达基因,其主要富集于代谢过程、应激反应及转录调控等条目中,这些基因的相互协同调控是红橘对该致病型产生防御反应的重要机制。