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

【目的】通过对意大利蜜蜂（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介导意大利蜜蜂工蜂对东方蜜蜂微孢子虫的跨界调控

【目的】东方蜜蜂微孢子虫（Nosema ceranae）感染意大利蜜蜂（Apis mellifera ligustica,简称意蜂）导致蜜蜂微孢子虫病。本研究结合前期已获得的miRNA和mRNA组学数据,通过生物信息学方法对意蜂工蜂中肠的差异表达miRNA（differentially expressed miRNA,DEmiRNA）靶向结合的东方蜜蜂微孢子虫的mRNA和差异表达mRNA（DEmRNA）进行预测、数据库注释和调控网络分析,以期在组学水平解析miRNA介导意蜂工蜂对东方蜜蜂微孢子虫的跨界调控机制。【方法】通过比较东方蜜蜂微孢子虫侵染7 d和10 d的意蜂工蜂中肠（AmT1、AmT2）和未受侵染的工蜂中肠（AmCK1、AmCK2）的miRNA组学数据筛选出宿主的显著性DEmiRNA,通过比较侵染意蜂工蜂中肠的东方蜜蜂微孢子虫（NcT1、NcT2）和东方蜜蜂微孢子虫纯净孢子（NcCK）的mRNA数据筛选出病原的DEmRNA。利用TargetFinder软件预测宿主显著性DEmiRNA靶向结合的病原mRNA和DEmRNA。利用相关生物信息学工具对上述靶DEmRNA进行GO和KEGG数据库注释。结合前期研究结果筛选出孢壁蛋白、极管蛋白、蓖麻毒素B凝集素、ABC转运蛋白、ATP/ADP移位酶和糖酵解/糖异生途径等毒力因子和能量代谢通路相关的病原DEmRNA及与其存在靶向结合关系的宿主显著性DEmiRNA,并构建和分析二者的调控网络。【结果】AmCK1 vs AmT1比较组中宿主的48条显著上调miRNA和36条显著下调miRNA分别靶向病原的1 345和1 046条mRNA;进一步分析发现,宿主的47条显著上调miRNA和34条显著下调miRNA可分别靶向NcCK vs NcT1比较组中病原的584条显著下调mRNA和265条显著上调mRNA,它们可分别注释到19和22个功能条目以及66和64条通路。AmCK2 vs AmT2比较组中宿主的56条显著上调miRNA和51条显著下调miRNA分别靶向病原的1 260和1 317条mRNA;进一步分析发现,宿主的52条显著上调miRNA和49条显著下调miRNA可分别靶向NcCK vs NcT2比较组中病原的587条显著下调mRNA和336条显著上调mRNA,它们可分别注释到20和23个功能条目以及64和65条通路。AmCK1 vs AmT1和AmCK2 vs AmT2比较组的8条共同显著上调miRNA和1条共同显著下调miRNA分别靶向NcCK vs NcT1和NcCK vs NcT2比较组中的144条共同显著下调和10条共同显著上调mRNA,可分别注释到18和13个功能条目以及38和7条通路。此外,AmCK1 vs AmT1和AmCK2 vs AmT2比较组中宿主的显著上调miRNA可靶向结合NcCK vs NcT1和NcCK vs NcT2比较组中与RNAi途径,孢壁蛋白和蓖麻毒素B凝集素等毒力因子,糖酵解/糖异生途径以及MAPK信号通路相关的病原下调表达mRNA。【结论】在东方蜜蜂微孢子虫的侵染过程中,意蜂工蜂中肠的DEmiRNA与病原的DEmRNA之间存在复杂的靶向结合关系以及潜在的跨界调控关系;宿主的DEmiRNA可能通过抑制或降解病原的RNAi途径、毒力因子、糖酵解/糖异生通路、ATP/ADP移位酶、ABC转运蛋白及MAPK信号通路相关靶DEmRNA影响病原的侵染和增殖。     

微小RNA介导东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂的分子机制

【目的】通过对东方蜜蜂微孢子虫（Nosema ceranae）纯化孢子与侵染意大利蜜蜂（Apis mellifera ligustica,简称意蜂）工蜂的东方蜜蜂微孢子虫的差异表达miRNA（DEmiRNA）及其靶mRNA进行系统分析,筛选、分析和探讨病原毒力因子和侵染因子相关的DEmiRNA及调控网络,在miRNA组学层面揭示东方蜜蜂微孢子虫对意蜂的侵染机制。【方法】利用small RNA-seq（sRNA-seq）技术对东方蜜蜂微孢子虫感染7 d和10 d的意蜂工蜂中肠和东方蜜蜂微孢子虫纯化孢子（NcCK）进行深度测序,通过连续比对rRNA数据库、西方蜜蜂（Apis mellifera）基因组和东方蜜蜂微孢子虫基因组筛滤出处于侵染过程的东方蜜蜂微孢子虫（NcT1和NcT2）数据和东方蜜蜂微孢子虫孢子的测序数据。根据P≤0.05,|log₂ fold change|≥1的标准,通过比较分析筛选出各比较组中的差异表达miRNA（differentially expressed miRNA,DEmiRNA）。通过相关生物信息学软件对DEmiRNA进行表达谱分析,靶mRNA预测及功能和代谢通路注释,以及调控网络的构建与分析。通过Stem-loop RT-qPCR验证DEmiRNA的差异表达趋势及测序数据的可靠性。【结果】NcCK vs NcT1、NcCK vs NcT2和NcT1 vs NcT2比较组分别包含164、122和60个DEmiRNA。Venn分析结果显示,3个比较组共有的上调和下调miRNA分别为5和6个。上述DEmiRNA分别预测出1 885、1 733和1 524个靶mRNA。这些靶mRNA分别注释到27、25和26个功能条目,其中注释数量最多的是新陈代谢进程、催化活性、细胞进程、结合和细胞。上述靶mRNA可分别注释到84、84和84条代谢通路,其中注释数量最多的是代谢途径、核糖体和次级代谢产物生物合成。此外,对于NcCK vs NcT1、NcCK vs NcT2和NcT1 vs NcT2中的DEmiRNA,分别有35、26和12个靶向结合MAPK信号通路相关靶mRNA,分别有49、40和17个DEmiRNA靶向结合糖酵解/糖异生通路相关靶mRNA。进一步分析发现,东方蜜蜂微孢子虫的DEmiRNA参与调控蓖麻毒素B凝集素、细胞凋亡抑制因子、极管蛋白和孢壁蛋白等病原毒力因子的基因表达,以及己糖激酶、ATP/ADP移位酶、ABC转运蛋白和转录因子ste12等侵染因子的基因表达。【结论】通过对东方蜜蜂微孢子虫纯化孢子与侵染意蜂工蜂的东方蜜蜂微孢子虫进行深入细致的miRNA组学分析和探讨,解析了病原侵染过程的miRNA差异表达谱,揭示了东方蜜蜂微孢子虫可能通过调节相应miRNA的表达水平对蓖麻毒素B凝集素、细胞凋亡抑制因子、极管蛋白和孢壁蛋白等毒力因子及己糖激酶、ATP/ADP移位酶、ABC转运蛋白和转录因子ste12等侵染因子基因表达进行调控,从而适应宿主细胞内的环境并促进自身的增殖与侵染。     

利用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等基因与窄叶的形成密切相关。【结论】明确了一些与玉米叶片发育密切相关的代谢路径,还发现植物激素间的动态平衡对叶片发育有着重要影响,尤其是生长素与油菜素内酯、细胞分裂素与赤霉素之间的相互作用对调控叶片形态可能发挥重要作用。     

人工劣变处理对玉米种胚差异基因表达的影响

【目的】利用数字基因表达谱技术（digital gene expression tag profiling,DGE）探索人工控制劣变（controlled deterioration treatment,CDT）条件下玉米种胚差异基因的表达变化,为揭示种子劣变的分子机理提供依据。【方法】以玉米杂交种郑单958种子为材料,采用高温（45℃）高湿（相对湿度100%）方法处理72 h,以未处理材料为对照,利用DGE分别进行高通量测序,测序结果与参考基因组和参考基因数据库比对获得表达基因,利用RPKM（Reads Per Kb per Million reads）方法计算基因的表达量,根据FDR（false discovery rate)<0.001和|log2 ratio(T/CK)|≥1的标准筛选差异表达的基因,对获得的差异表达基因（digital expression genes,DEGs）进行GO（gene ontology）和KEGG（kyoto encyclopedia of genes and genomes）数据库功能注释,并对注释结果进行富集分析处理。【结果】对照玉米种胚中检测到32 000多种mRNAs。CDT处理后差异表达基因有4 713个。其中上调表达2 874个,下调表达1 839个。GO富集分析表明,这些基因涉及细胞组分、分子功能和生物学过程方面,其编码的蛋白主要分布在细胞器和细胞膜上,参与能量代谢、信号转导、刺激响应和衰老死亡等过程,具有结合、催化和抗氧化等功能。这些基因中,能被KEGG数据库注释的有2 470个,参与了288条代谢通路,其中有16条通路存在着显著性富集。这些通路中,参与能量代谢的基因共有113个,分别参与糖酵解/糖异生过程（59个）、磷酸戊糖途径（31个）和丙酮酸代谢过程（50个）;其中,调节糖酵解/糖异生过程的烯醇化酶基因、3-磷酸甘油醛脱氢酶基因等、丙酮酸代谢中的丙酮酸激酶基因等和磷酸戊糖途径中的α-L-岩藻糖苷酶基因等上调幅度最大。还发现,调节NADH代谢相关的基因有25个（9个上调,16个下调）;NADPH代谢的有10个（4个上调,6个下调）。这些基因的表达变化能够调节活性氧的产生和积累。【结论】DGE可以作为研究种子劣变和活力丧失的有效手段。CDT能够影响玉米种胚差异基因的表达,进而影响细胞能量代谢相关过程,主要表现在糖酵解途径受到抑制,导致ROS产生和积累,从而加速玉米种胚细胞的衰老及死亡,最终造成种子的劣变和活力丧失。在这个过程中,差异表达基因可能扮演重要角色。     

东方蜜蜂微孢子虫孢子中长链非编码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间接影响靶基因的表达,从而参与调节东方蜜蜂微孢子虫孢子中的物质和能量代谢等基本生命活动。     

高原雨点鸽与詹森鸽胸肌转录组差异表达基因筛选与功能分析

筛选高原雨点鸽与詹森鸽胸肌飞行能力差异的关键基因,为赛鸽选育奠定基础.选择性别相同、体况与日龄相近的高原雨点鸽、詹森鸽各3只,屠宰后取胸肌,进行转录组测序,筛选差异表达基因,对差异基因进行Gene Ontology(GO)分析、Kyoto Encyclopedia of Genes and Genomes(KEGG)分...     

SO2引起巨峰葡萄采后落粒的共表达网络和转录调控分析

【目的】 二氧化硫（SO₂）处理可有效防治葡萄灰霉病和采后腐烂,但会导致浆果脱落,研究SO₂引起葡萄浆果脱落的分子机制,明确关键基因和转录调控机制。【方法】 使用SO₂处理‘巨峰’葡萄,分别在2、4和6 d进行采样并统计对照组（CK）和SO₂处理组的葡萄浆果脱落率。通过高通量测序技术对CK对照组和SO₂处理组‘巨峰’葡萄采后2、4和6 d的样品测序,以葡萄基因组作为参考基因组进行序列比对,利用TPM算法计算基因表达量,使用基因集富集分析（GSEA）、基因共表达网络（GCN）以及转录调控网络预测方法对转录组数据进行系统的分析,并利用实时荧光定量PCR（qRT-PCR）技术对其表达量进行验证。【结果】 SO₂处理能够显著诱导葡萄浆果脱落,2 d时,SO₂处理的葡萄浆果脱落率为9.88%,4 d时达到19.24%,且均显著高于对照组,6 d时脱落率达到38.25%,而对照组脱落率仅11.85%。通过GSEA分析,发现在CK组富集的GO生物过程主要与氧化应激反应、细胞壁代谢以及苯丙氨酸代谢通路相关,其中4 d时CK组富集到植物细胞壁组织、果胶代谢、细胞壁修饰、聚半乳糖醛酸等通路。在SO₂组富集的GO生物过程主要与能量代谢通路相关,其中2 d时在SO₂组富集到光合作用、四吡咯代谢、前体代谢物和能量的产生、葡萄糖代谢等通路。CK组富集的KEGG代谢通路主要有戊糖和葡萄糖醛酸的相互转换、半乳糖代谢、植物激素信号转导、柠檬酸循环（TCA循环）等,SO₂组有光合作用、柠檬酸循环（TCA循环）、糖酵解等。GCN将GSEA分析中的领头基因分为12个水平,水平4—9均富集到了能量代谢、糖代谢相关通路,其中水平4富集到激素反应、氧化应激反应。对GCN关键水平中的基因启动子序列进行转录调控预测分析,结果显示95个转录因子（TFs）与269个靶基因存在987对调控关系。WRKY14、WOX8、KUA1在SO₂处理下持续下调,MYB60、MYB73、ANL2、ERF2、DOF3.6、GATA25、WRKY57、KAN2、ATHB6在SO₂处理后持续上调。此外,MYB15、WRKY11、WRKY33、WRKY40、WRK75先上调后下调。ERF2、MYB60和WRKY40的转录调控网络发现调控的靶基因涉及细胞壁代谢、糖代谢等相关途径。qRT-PCR结果显示PME36和ERF2上调表达趋势相似,GAUT7、MYB60、UGE3上调表达趋势相似,此外,WRKY40在SO₂处理的2和4 d被诱导上调,PPME1、COMT1表达持续下调,SO₂处理4 d时LAC15被显著诱导上调。【结论】 SO₂诱导营养物质代谢、能量代谢、细胞壁代谢途径相关基因的表达,该过程受到多种转录因子调控,最终导致葡萄浆果脱落。     

福清山羊与努比亚黑山羊发情期卵巢组织RNA-Seq分析

【目的】 研究福清山羊与努比亚黑山羊发情期卵巢组织转录组差异表达水平,一方面为山羊繁殖性状形成的分子机制提供理论依据,另一方面为利用努比亚黑山羊杂交改良福清山羊提供可加快遗传进展的分子标记。【方法】 利用转录组测序方法对福清山羊和努比亚黑山羊发情期卵巢组织进行研究,筛选品种间的差异表达基因（differentially expressed genes,DEGs）,并对DEGs功能进行注释和若干基因荧光定量PCR（quantitative real-time PCR,qRT-PCR）验证;同时通过与参考基因组比对,分析和筛选测序样品中存在的SNP/InDel。【结果】 6个样品共得到46.68Gb Clean Data,DESeq分析发现了福清山羊和努比亚黑山羊发情期卵巢组织的DEGs149个（包含25个新转录本）,其中表达上调53个,表达下调96个;初步认为输卵管素基因（oviductin,OVN）、类固醇合成快速调节蛋白基因(steroidogenic acute regulatory protein,STAR)、早期生长应答1基因(early growth response 1,EGR1)可作为福清山羊繁殖性能的候选基因。149个DEGs中的108个基因能被GO（gene ontology）数据库注释,30个DEGs能够被COG（Cluster of Orthologous Groups of proteins）数据库注释,91个DEGs能够被KEGG（kyoto encyclopediaof genes and genomes）数据库注释。KEGG通路分析表明DEGs共富集到21条信号通路中,6条通路显著富集。经过进一步的与参考基因组序列比对分析,共发掘1 506个新转录本。经qRT-PCR验证,所选基因（转录本）表达变化模式与转录组测序结果一致,表明测序结果可靠。【结论】 在转录组水平上筛选出了福清山羊和努比亚黑山羊发情期卵巢组织的149个DEGs,发掘了1 506个新转录本,初步揭示了OVN、STAR和EGR1在山羊繁殖过程中可能发挥重要作用,为进一步探索山羊繁殖性状相关机理提供参考依据。     

利用iTRAQ技术和转录组筛选芍药属远缘杂交不亲和基因

【目的】远缘杂交育种是目前牡丹、芍药品种改良和育种的主要方法,而远缘杂交不亲和一直是制约其快速发展的主要因素。本研究从牡丹、芍药远缘杂交授粉后不亲和应答相关的柱头差异蛋白与转录组方面深入研究,揭示牡丹、芍药远缘杂交不亲和的分子机理,为杂交育种提供理论依据。【方法】以芍药‘粉玉奴’自交、芍药‘粉玉奴’与牡丹‘凤丹白’杂交为供试材料,在授粉后24 h采取柱头,分别进行同位素标记相对定量（iTRAQ）和转录组技术分析。对所获得的蛋白和转录组数据进行生物信息学分析,并对其中可能与远缘杂交不亲和相关的基因进行定量PCR验证。【结果】利用iTRAQ技术分析牡丹、芍药远缘杂交后柱头中蛋白质的表达差异,共鉴定到685个差异蛋白,富集到了188条通路,其中显著富集的Pathway有18条。与不亲和授粉相关代谢通路有RNA降解、钙信号途径、丝裂原活化蛋白激酶（mitogen-activated protein kinase signaling pathway,MAPK）信号途径、磷脂酰肌醇信号系统。在RNA降解代谢通路中,烯醇酶（Enolase）、热休克蛋白DnaK（HSP70）及病菌抗原（GroEL）均表达下调。在钙信号途径中,钙调蛋白（CALM）表达下调,腺苷酸转运酶（adenine nucleotide translocase,ANT）表达量增加,表达上调。MAPK信号途径中,乙二醛酶Ⅰ（GloI）表达下调。磷脂酰肌醇信号系统中的钙调蛋白（CALM）表达下调。随机选取与差异蛋白相关的6个基因进行qRT-PCR验证,结果显示,6个基因的表达与蛋白质水平趋势相一致,均表达下调。通过转录组测序,共获得了52 998个有注释信息的Unigene,占所有Unigene的40.37%。基于6组样品的RPKM（Reads Per Kilobase per Million）值,共筛选到16 224个差异基因。其中上调基因13 361,下调基因2 863个。对差异基因进行Pathway显著富集分析,杂交与自交相比,不亲和差异表达的基因主要富集在氧化磷酸化代谢、ABC转运蛋白、次级代谢产物等通路。与远缘杂交不亲和相关且发生显著变化的基因有CalS-5、CalS-12（胼胝质酶）和SPL（squamosa promoter binding protein-like）表达上调,ABCF（ABC transporter family protein）表达下调。【结论】在转录组和蛋白数据共注释到6个蛋白、4个基因与植物不亲和性密切相关,这些蛋白与基因可能在远缘杂交不亲和方面发挥着重要作用。     

枯草芽孢杆菌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转运蛋白通路中的差异基因（或蛋白）可能在棉花根系分泌物与枯草芽孢杆菌互作过程方面发挥着重要作用。     

‘仓方早生’桃及其早熟芽变不同发育时期果实的转录组分析

【目的】通过对桃品种‘仓方早生’及其早熟芽变不同发育时期的果实进行转录组分析，挖掘参与调控桃果实成熟的关键因子，为深入研究桃果实成熟调控机理提供理论依据。【方法】以桃品种‘仓方早生’及其早熟芽变为试材，每个品种分别选择长势一致的样品树5株，分别于花后30 d（对应‘仓方早生’c1、早熟芽变y1）、45 d（对应c2、y2）、59 d（对应c3、y3）、71 d（对应c4、y4）及89 d（对应c5）对不同发育时期的桃去皮果肉进行取样和转录组测序，并利用实时荧光定量PCR(qRT-PCR)对筛选的差异表达基因进行定量验证；利用GO和KEGG对‘仓方早生’及其早熟芽变的差异表达基因进行分析；基于差异表达基因构建加权基因共表达网络分析（weighted gene co-expression network analysis,WGCNA），从中鉴定出与果实成熟密切相关的枢纽模块和枢纽基因。【结果】将处于果实相同发育时期的转录组数据进行比较，得到y1与c1、y2与c2、y3与c4和y4与c5四组对比数据，共筛选出差异表达基因4 395个，其中上调表达基因2 212个，下调表达基因2 183个。其...     

基于转录-代谢联合分析哈茨木霉ACCC32527 对NaCl胁迫的分子调节

【目的】 通过NaCl胁迫下哈茨木霉（Trichoderma harzianum）ACCC32527差异转录组和代谢组数据分析,挖掘关键功能型差异表达基因（DEG）及次级代谢产物。【方法】 采用RNA-seq和GC-TOF-MS技术分别完成0（T1）、0.4（T2）、0.6（T3）mol?L ^-1 NaCl胁迫下哈茨木霉ACCC32527的差异表达基因和次级代谢产物检测,利用相关软件及数据库（GO、COG、KEGG pathway等）完成对差异表达基因（|log₂fold change)|>1 & FDR<0.01）和次级代谢产物（P-value≤0.05 & VIP>1）的筛选、注释和分类,并进行RT-qPCR验证。 【结果】 NaCl胁迫处理后,ACCC32527的生长受到抑制,繁殖速率明显降低,并在该条件下分别获得T1 vs T2和T1 vs T3比较组637、1 570个DEG,获得DEG的16种表达模式,包括9种上调表达模式（950 gene）、3种下调表达模式（662 gene）和4种不规则表达模式（309 gene）,共有的GO功能分类为催化活性、结合、细胞器、细胞膜部分、细胞膜、细胞部分、细胞、单组织过程和代谢过程,且占主要比例,约为61%—94%,其中处理前后基因比例差异较大的分类为催化活性。KEGG代谢通路富集分析发现,不同NaCl浓度处理下的DEG富集到的代谢通路种类及数量存在较大差异,分别有225和535个差异基因富集于20条KEGG通路,包括代谢通路、抗生素的合成和次级代谢产物合成, 其中T2和T3处理下核糖体的生物合成途径分别有12个和18个相关基因受到抑制。从NaCl胁迫下差异转录组中筛选出活性氧清除、离子转运和细胞壁及胞外结构等共73个相关基因,并且多数为上调表达基因。另外,差异代谢组学分析表明,T3处理下,胞内小分子代谢物出现明显变化,累积量增加的代谢产物有30种,包括氨基酸及其衍生物、糖类及其衍生物和醇类,而减少的代谢物有53种,涉及脂肪酸和有机酸等。其中,甘油是已知的真菌重要渗透保护物,T3处理下ACCC32527胞内甘油水平提高,可参与细胞的渗透调节,维持渗透压平衡。RT-qPCR验证了7个差异表达显著基因,虽然在表达量上存在一定的差异,但表达趋势与RNA-seq分析结果基本一致,表明转录组测序结果的可靠性。【结论】 NaCl胁迫引起哈茨木霉ACCC32527大量基因及次级代谢产物变化,获得了与NaCl胁迫调节相关基因和代谢产物,这些机制共同作用减少盐离子对细胞的毒害作用,研究结果可为研究木霉菌的耐盐机理提供重要信息。     

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

【目的】拟轮枝镰孢(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及其介导的竞争性内源RNA调控网络在意大利蜜蜂工蜂中肠发育过程中的潜在作用

【目的】微小RNA（microRNA,miRNA）能够通过靶向结合导致mRNA的抑制或降解,从而对基因表达进行负调控。MiRNA在昆虫的生长、发育、免疫和细胞生命活动调控方面扮演关键角色。本研究旨在对意大利蜜蜂（Apis mellifera ligustica,简称意蜂）工蜂中肠的差异表达miRNA（differentially expressed miRNA,DEmiRNA)及其调控网络进行深入分析,并结合前期在mRNA、长链非编码RNA(long non-coding RNA,lncRNA)和环状RNA(circular RNA,circRNA)组学层面的相关研究结果,系统解析中肠发育过程的miRNA差异表达谱及竞争性内源RNA（competing endogenous RNA,ceRNA）调控网络介导的中肠发育机理。【方法】利用small RNA-seq技术对正常的意蜂7日龄和10日龄工蜂中肠（Am7和Am10）进行测序。将质控后的测序数据比对西方蜜蜂基因组,再将比对上的序列标签（tags）比对到miRBase数据库以鉴定已知miRNA。利用每百万标签序列（tags per million,TPM）算法对miRNA表达量进行计算和归一化处理。根据|log₂fold change|≥1且P≤0.05的标准筛选得到显著性DEmiRNA;利用相关软件进行靶mRNA的预测及GO和KEGG数据库注释。根据靶向结合关系及前期研究结果,利用Cytoscape软件对AMPK、PI3K-Akt、Wnt、cAMP、Hippo、mTOR、Toll/Imd、TGF-beta和MAPK 9条信号通路相关的DElncRNA/DEcircRNA-DEmiRNA-DEmRNA网络进行可视化;构建以miR-342-y为核心的ceRNA调控网络,进一步对网络中的靶mRNA进行代谢通路注释。利用茎环实时荧光定量PCR（Stem-loop RT-qPCR）验证测序数据及miRNA差异表达的可靠性。【结果】Am7 vs Am10比较组中共有112个显著性DEmiRNA,包括38个显著上调miRNA和74个显著下调miRNA,分别靶向结合7 434和9 559个mRNA,它们可注释到生物学进程相关的21和23个功能条目,细胞组分相关的16和17个功能条目,以及分子功能相关的10和11个功能条目;此外还能注释到果糖和甘露糖代谢、嘌呤代谢、甘氨酸、丝氨酸和苏氨酸代谢等物质代谢相关的83和86条通路;内吞作用、泛素蛋白水解、黑色素生成等细胞免疫相关的10和10条通路;MAPK、Jak-STAT、NF-κB等体液免疫相关的5和5条通路;以及Hippo、FoxO、Notch等涉及生长发育的13和11条信号通路。进一步构建上述9条信号通路相关的ceRNA调控网络,分析发现意蜂中肠发育过程中DEmiRNA与DElncRNA、DEcircRNA和DEmRNA之间存在复杂的调控关系;DEmiRNA位于网络中心,而DElncRNA、DEcircRNA和DEmRNA位于网络的外周。miR-342-y在意蜂工蜂中肠和幼虫肠道发育过程中均为显著性下调表达,可靶向结合3个DEcircRNA、4个DElncRNA和327个mRNA。Stem-loop RT-qPCR结果显示4个DEmiRNA的差异表达趋势与测序结果一致,说明本研究中测序数据及miRNA差异表达趋势真实可靠。【结论】DEmiRNA可能通过参与调控物质和能量代谢通路、Hippo和Wnt等信号通路以及细胞和体液免疫通路相关基因的表达影响意蜂中肠的生长和发育;miR-182-x、miR-291-y、miR-342-y、ame-miR-6001-3p等关键DEmiRNA介导的ceRNA调控网络可能在意蜂中肠发育过程发挥重要的调控作用。     

不同温度下TMV侵染枯斑三生烟的LncRNA差异表达

【目的】 筛选不同温度下烟草花叶病毒（Tobacco mosaic virus,TMV）侵染后枯斑三生烟（Nicotiana tabacum var. Samsun NN）差异表达的长链非编码RNA（long non-coding RNA,lncRNA）,研究lncRNA在枯斑三生烟抗性反应中的作用。【方法】 N基因的温度敏感性使枯斑三生烟在25℃时具备对TMV的抗性、在31℃抗性丧失,在这两个温度条件下对枯斑三生烟接种TMV和磷酸盐缓冲盐水（phosphate buffered saline,PBS）,48 h后提取系统叶总RNA,构建链特异性文库后进行深度测序。对测序结果进行过滤后利用HTSeq将有效数据与近缘品种TN90（N. tabacum var. TN90）基因组比对,筛选得到lncRNA后利用FPKM法估计lncRNA的表达水平。通过edgeR筛选差异表达lncRNA（differentially expressed lncRNA,DElncRNA）,并利用qRT-PCR技术对这一结果进行验证。通过共定位及共表达分析预测DElncRNA的靶基因,通过参考基因组注释、GO和KEGG富集分析研究靶基因的功能。【结果】 4个处理共12个样本经lncRNA-seq各测得约8 000万条clean reads,共获得4 737条已知lncRNA、40 169条新lncRNA。其中64个lncRNA在不同温度条件下TMV侵染后存在差异表达,qRT-PCR测定结果显示这些lncRNA的测序正确率在80%左右,表明本研究所得测序数据具备较高的可信度。对DElncRNA进行靶基因预测,发现一些基因同时被25℃下调和31℃上调的DElncRNA靶向。靶基因注释功能丰富,主要参与植物抗病、激素和代谢等生理过程。部分可能与激素通路相关的lncRNA,在25℃下TMV侵染时呈现下调趋势,而在31℃下TMV侵染则呈现上调趋势。GO富集分析显示靶基因主要参与构成膜、囊泡等组分,具备钙、钾离子通道抑制剂活性等分子功能,使相应离子得以转运引发随后的反应,同时也参与发病、抗原加工和呈现、细胞分裂素代谢等生理过程。KEGG分析发现靶基因显著富集在植物激素信号转导通路,25℃下调和31℃上调的DElncRNA靶基因同时富集在激素信号传导、ABC运输蛋白、苯丙烷类生物合成等通路。【结论】 不同温度（25℃和31℃）条件下TMV侵染枯斑三生烟后,长链非编码RNA差异表达,DElncRNA通过作用于激素信号传导、物质转运等过程参与寄主系统获得性抗性反应。研究结果可为揭示植物系统获得性抗性中lncRNA的调控功能以及新型抗病毒技术开发提供依据。     

2种不同生长规格墨瑞鳕肌肉组织转录组测序分析

【目的】掌握不同生长规格墨瑞鳕个体间差异表达基因的表达特点,为其功能相关基因深度挖掘及分子遗传育种提供科学依据。【方法】挑选同一养殖条件下极大个体和极小个体的墨瑞鳕,构建肌肉组织cDNA文库后,采用Illumina HiSeqTM 4000测序平台对存在生长差异的墨瑞鳕肌肉组织进行转录组测序分析,获得的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等数据库中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行GO功能注释及KEGG信号通路富集分析,并采用MISA进行SSR鉴定分析。【结果】从墨瑞鳕肌肉组织中共测序获得39749条Unigenes,其长度范围在301~55230 bp,平均长度为1705 bp。注释到Nt、Nr、Swiss-Prot、Pfam数据库的Unigenes分别有27046、20824、18268和17772条,在7个数据库中均得到注释的Unigenes共计6742条,占Unigenes总数的16.96%。根据差异表达基因筛选条件P<0.05且|log₂ Fold Change|>1,共筛选出722个差异表达基因,其中上调基因308个、下调基因414个。差异表达基因GO功能注释分析结果表明,注释基因数目较多的GO功能条目包括细胞过程、代谢过程、膜、细胞器及结合等;KEGG信号通路富集分析发现,差异表达基因被成功富集到234条信号通路上,主要涉及磷脂酰肌醇3激酶/蛋白激酶信号通路、MAPK信号通路、胰岛素信号通路及FoxO信号通路等。在39749条Unigenes中鉴定筛选出22120个SSRs,占Unigene总数的55.65%,SSR的平均间距为3063 bp。【结论】基于转录组测序分析获得的墨瑞鳕肌肉组织差异表达基因以发挥结合、细胞过程及代谢过程等功能为主,且主要富集在PI3K-Akt信号通路、核糖体信号通路、FoxO信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。