利用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个基因与植物不亲和性密切相关,这些蛋白与基因可能在远缘杂交不亲和方面发挥着重要作用。     

授粉方式对牡丹和芍药远缘杂交受精前障碍的影响

【目的】探寻克服牡丹与芍药远缘杂交受精前障碍的授粉方式,为克服牡丹、芍药远缘杂交不亲和提供理论依据。【方法】以芍药品种粉玉奴(Paeonia lactiflora‘Fenyunu’)为母本,牡丹品种凤丹白(Paeonia ostii‘Fengdanbai’)、洛阳红(Paeonia suffruticos a‘Luoyanghong’)为父本分别进行远缘杂交,采用常规授粉、延迟授粉、重复授粉、切割柱头授粉和蒙导授粉5种不同授粉方法,对授粉后花粉萌发及花粉管伸长过程进行荧光显微观察,研究母本心皮形态学特征,统计杂交组合结实率。【结果】杂交组合粉玉奴×凤丹白、粉玉奴×洛阳红采用重复授粉、蒙导授粉、延迟授粉的效果均优于常规授粉,在花粉萌发量、花粉管伸长及单花结实率方面表现为延迟授粉重复授粉蒙导授粉常规授粉切割柱头授粉,说明延迟授粉对于克服受精前障碍效果最佳,为最佳授粉方式。切割柱头授粉虽效果不佳,但对克服受精前障碍有一定影响。【结论】在牡丹、芍药远缘杂交育种中,可以采取延迟授粉的方法,减缓甚至克服母本柱头或花柱对异源花粉萌发及花粉管生长的抑制,克服受精前障碍,提高育种效率。     

牡丹与芍药的授粉亲和性表现及其生理机制分析

【目的】探究花粉-柱头相互作用及主要保护酶、内源激素与牡丹芍药远缘杂交不亲和性的关系,为阐明不亲和的生理机制奠定基础。【方法】对芍药与牡丹杂交授粉及芍药自交授粉过程中,花粉在雌蕊上不同时间(1,3,5,8,12,24,36,48,72,96h)的生长动态进行荧光显微观察,并对此过程中雌蕊内的保护酶(SOD、POD)活性和内源激素(ZR、ABA、IAA、GA3)含量变化进行比较,分析其形态变化的内在生理机制。【结果】芍药自交授粉后花粉大量萌发,花粉管伸长速度先慢后快,并通过花柱基部;芍药与牡丹杂交后,花粉大部分不能萌发,并在柱头上出现扭曲、肿胀等现象,且产生强烈的胼胝质反应。自交亲和及杂交不亲和授粉雌蕊的POD和SOD活性在授粉过程中发生明显变化,且2种雌蕊的变化趋势存在明显差异;当POD和SOD处于高活性时,花粉管迅速向胚囊伸长并进入胚囊(授粉后3h)。高含量的ZR、IAA和GA3有利于花粉-柱头识别、萌发和受精作用(授粉后3h);整个受精过程中,自交亲和授粉雌蕊的IAA和GA3含量均显著高于杂交不亲和授粉雌蕊,而杂交授粉雌蕊中的ABA含量高于自交授粉,说明高含量的ABA与杂交不亲和性相关。【结论】芍药与牡丹远缘杂交过程中的不亲和性与内部酶活性变化及内源激素含量的动态变化有关。     

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

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

甘蓝型油菜授粉功能缺失突变体及其 野生型柱头差异蛋白研究

【目的】基于蛋白质组学水平分析甘蓝型油菜野生型和授粉功能缺失突变体的柱头差异蛋白。【方法】运用双向电泳技术研究甘蓝型油菜柱头授粉功能缺失的突变体FS-M1及其野生型（宁油10号）柱头差异蛋白的表达,结合质谱（MALDI-TOF-TOF/MS）分析与蛋白质数据检索技术,鉴定甘蓝型油菜野生型柱头上调表达蛋白的性质和功能。【结果】获得了甘蓝型油菜野生型柱头显著上调且可鉴定的蛋白有33个,包括抗胁迫/防御、氧化还原反应平衡调节、碳水化合物代谢、蛋白质水解、蛋白折叠、氨基酸和氮代谢、核苷酸代谢等7类功能蛋白。【结论】甘蓝型油菜野生型柱头中上调表达蛋白功能广泛,这些表达蛋白可能与甘蓝型油菜柱头授粉功能的调控有关。     

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

【目的】通过对意大利蜜蜂（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可能在意蜂工蜂对东方蜜蜂微孢子虫的免疫应答中起到关键作用。     

5-氮杂胞苷对苦瓜果实成熟及转录表达谱的影响

【目的】研究5-氮杂胞苷（5-azaC）对苦瓜果实成熟和基因表达模式的影响,为解析苦瓜果实成熟的分子调控机制提供理论依据。【方法】以苦瓜（Momordica charantia L.）高世代自交系E12201-e1为材料,利用50 mmol/L 5-azaC处理绿熟期（授粉后18 d）苦瓜果实,以ddH₂O处理为对照,观察处理组（TR）和对照组（CK）果实的表型变化,并利用RNA-seq技术比较TR和CK苦瓜果肉组织中的转录水平差异。【结果】50 mmol/L 5-azaC处理能明显延迟苦瓜果实成熟。RNA-seq获得的数据质量较高,可满足后续分析需求。通过筛选共获得1 773个差异表达基因（DEGs）,其中1 007个上调基因,766个下调基因。GO和KEGG分析结果表明,DEGs主要富集在细胞组分和分子功能模块,217个DEGs共富集到93条代谢通路中,其中30个DEGs显著富集到光合作用和苯丙烷类生物合成通路。8个DEGs参与细胞壁代谢过程,其中6个基因在5-azaC处理中下调表达。共获得681个差异表达转录因子（TFs）,分属于220个TF家族。12个DEGs的qRT-PCR结果与RNA-seq结果变化趋势一致。【结论】5-azaC处理抑制苦瓜果实成熟,并影响果实成熟相关基因的表达模式。     

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₂诱导营养物质代谢、能量代谢、细胞壁代谢途径相关基因的表达,该过程受到多种转录因子调控,最终导致葡萄浆果脱落。     

持续感染口蹄疫病毒VP1蛋白对细胞转录组的影响

【目的】通过转录组测序（RNA-Seq）方法分析持续感染口蹄疫病毒结构蛋白VP1诱导的PK-15细胞基因表达谱的变化,探索VP1蛋白对宿主细胞信号通路的影响。为后续开展VP1影响FMDV持续感染相关机制研究提供基础数据和理论支持。【方法】PK-15细胞分别转染pCAGGS-HA（对照载体）、pCAGGS-HA-Pi-VP1质粒,转染36 h后,Western blot检测VP1蛋白表达,同时提取两组细胞总RNA进行转录组测序（每组3个重复,2组6个样本）。将测序所得的序列进行过滤比对,获得差异表达基因后,对差异表达基因进行了生物信息学分析。采用qRT-PCR验证关键差异表达基因的表达。【结果】RNA-Seq共鉴定出3 376个差异表达基因（DEGs）|log2(Fold Change)|>0和padj<0.05,其中上调基因1 744个,下调基因1 632个。许多参与免疫反应和炎症的效应基因差异表达。GO和KEGG富集分析显示,差异基因GO条目绝大部分与核糖体有关,而KEGG富集中涉及了与天然免疫相关的信号通路,如TNF信号通路、Toll受体信号通路。【结论】pCAGGS-...     

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

【目的】拟轮枝镰孢(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富集分析显示，侵染早期拟轮枝镰孢在细胞间隙生长，差异基因主要富集在...     

芍药花茎强度与褪黑素含量的关系分析

花茎挺直程度是影响芍药切花品质的重要指标,而花茎强度是衡量茎秆挺直程度的物理学标准。为了明确芍药花茎强度与褪黑素含量的关系,以高、低花茎强度的2组芍药品种为材料,测定了花茎强度、植株形态指标、植株光合参数、花茎木质素含量和褪黑素含量。结果表明:与低花茎强度芍药品种相比,高花茎强度芍药品种具有较高的茎粗、茎质量、光合能力、木质素含量和褪黑素含量,并且花茎强度与褪黑素含量之间呈正相关,两者均与木质素含量达显著正相关,表明褪黑素可能通过调控木质素合成来间接影响花茎强度,这为芍药花茎强度的调控奠定了理论基础。     

福清山羊与努比亚黑山羊发情期卵巢组织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在山羊繁殖过程中可能发挥重要作用,为进一步探索山羊繁殖性状相关机理提供参考依据。     

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

【目的】通过对桃品种‘仓方早生’及其早熟芽变不同发育时期的果实进行转录组分析，挖掘参与调控桃果实成熟的关键因子，为深入研究桃果实成熟调控机理提供理论依据。【方法】以桃品种‘仓方早生’及其早熟芽变为试材，每个品种分别选择长势一致的样品树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个。其...     

匍匐翦股颖接种立枯丝核菌后基因表达变化的转录组学分析

【目的】确定匍匐翦股颖(Agrostis stolonifera)接种立枯丝核菌(Rhizoctonia solani)后基因种类和表达量在转录水平的变化规律,明确草坪草病原菌侵染响应的关键基因。【方法】匍匐翦股颖生长14 d后采用麦粒培养物接种立枯丝核菌,接种3 d后选取感病叶片和未接种叶片提取RNA,进行转录组高通量测序,然后利用生物信息学分析,用Trinity组装匍匐翦股颖转录组,以组装子为参考,以|log2(fold change)|1,q-value0.005为阈值选取感病和健康匍匐翦股颖叶片转录组的差异表达基因,并用i TAK软件分析其转录因子家族及表达变化,与植物R基因库进行blast分析R蛋白分类、利用Mapman软件分析生物胁迫信号通路相关基因的表达变化。【结果】高通量测序得到125 253 092条高质量待分析reads。经Trinity从头组装后,得到466 761条转录本。过半数的转录本长度为700 bp以上,组装结果 N50=1 100 bp。使用CD-HIT选择334 212条转录本(所有转录本的71.60%)作为Unigene,平均长度573 bp,N50=791 bp。接种后植物比接种前植物基因有7 937个上调表达,1 570个下调表达。上调基因中296个,下调基因有142个都可被定义为转录因子,分布在58个转录因子家族中,其中锌指蛋白包含转录因子C2H2最多,有54个,C3H次之,为22个。差异基因中451个可定义为植物R蛋白表达基因,可分为33类,其中包含NBS-LRR结构域的抗病蛋白、LRR受体蛋白激酶、ABC-2类型的转运蛋白、U-box结构域蛋白激酶和热激蛋白这5类基因变化最显著。差异基因中大量上调表达基因可富集在病原识别、活性氧消除、信号传导、细胞凋亡、病程相关蛋白等生物胁迫相关基因类别,下调基因显著富集在植物生长发育相关类别和通路。q RT-PCR验证了随机挑选差异基因的表达量变化,均与RNA-seq分析结果一致,其中包括12个C2H2转录因子基因,10个C3H转录因子基因,以及12个R蛋白基因。【结论】病原菌侵染后,引起匍匐翦股颖大量基因表达变化,其中转录因子、R蛋白以及抗性相关基因多为上调表达,作用为抑制病原菌扩散,而生长发育相关基因表达下调,这些进程共同使匍匐翦股颖产生对立枯丝核菌的先天基础抗性。     

西方蜜蜂工蜂不同虫态发育的转录组学分析

【目的】基于转录组学对西方蜜蜂工蜂不同虫态间的差异表达基因（DEGs）进行筛选和功能注释分析,揭示与工蜂生长发育相关的信号通路,为深入解析工蜂生长发育的分子调控机理提供基础数据。【方法】以西方蜜蜂工蜂的3日龄幼虫、1日龄蛹和1日龄羽化工蜂3个虫态为研究对象,利用llumina NovaSeq 6000平台进行转录组测序,采用DESeq2筛选不同虫态样品间的表达差异基因,然后分别进行GO功能注释分析及KEGG信号通路富集分析,并通过实时荧光定量PCR进行验证。【结果】经转录组测序,在西方蜜蜂工蜂3日龄幼虫与1日龄蛹间筛选出4823个差异表达基因（51.86%上调,48.14%下调）,在1日龄蛹与1日龄羽化工蜂间筛选出3295个差异表达基因（57.51%上调,42.49%下调）,在3日龄幼虫与1日龄羽化工蜂间筛选出5267个差异表达基因（52.95%上调,47.05%下调）。GO功能注释分析结果显示,3日龄幼虫与1日龄蛹间的差异表达基因注释到43个GO功能条目,1日龄蛹与1日龄羽化工蜂间的差异表达基因注释到45个GO功能条目,3日龄幼虫与1日龄羽化工蜂间的差异表达基因注释到44个GO功能条目,主要涉及细胞过程、细胞部分及结合等。KEGG信号通路富集分析发现,3日龄幼虫与1日龄蛹间有2905个差异表达基因富集到332条KEGG信号通路上,其中17条KEGG信号通路呈显著富集,涉及核糖体、氧化磷酸化和昆虫激素生物合成等;1日龄蛹与1日龄羽化工蜂间有1644个差异表达基因富集到331条KEGG信号通路上,其中45条KEGG信号通路呈显著富集,涉及氧化磷酸化、生热作用和胰岛素分泌等;3日龄幼虫与1日龄羽化工蜂间有2958个差异表达基因富集到337条KEGG信号通路上,其中14条KEGG信号通路呈显著富集,涉及核糖体、蛋白酶体和胰岛素分泌等。6个随机挑选差异表达基因的实时荧光定量PCR检测结果与转录组测序结果相符,即转录组测序结果可靠。【结论】昆虫激素生物合成通路相关差异表达基因调控与西方蜜蜂工蜂各虫态JH滴度变化规律一致,氧化磷酸化信号通路则与各虫态的营养摄入和活动行为相关,而胰岛素分泌通路涉及各虫态的营养调控、脂肪体合成及细胞凋亡。可见,昆虫激素生物合成、胰岛素分泌和氧化磷酸化3种信号通路在西方蜜蜂工蜂幼虫、蛹和成虫的发育调控中发挥着重要作用。