球囊菌胁迫中华蜜蜂幼虫肠道过程中病原的转录组学研究

【目的】本研究利用RNA-seq技术对球囊菌胁迫的中华蜜蜂(中蜂)幼虫肠道进行深度测序,经趋势分析得到差异表达基因(DEGs)的显著表达模式,进而对胁迫过程中的球囊菌进行转录组学分析。【方法】利用Illumina HiSeq 2500平台对球囊菌胁迫的中蜂幼虫肠道进行深度测序,并利用相关软件进行了深入分析。最后,通过RT-qPCR对RNA-seq数据进行了验证。【结果】本研究共得到球囊菌的41133932条高质量clean reads。22865个DEGs共聚类为8个基因表达模式,其中,16769个DEGs聚类为2个显著上调趋势与2个显著下调趋势。GO富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于40与37个GO term,基因富集数最多的为细胞进程(2486 unigenes)。KEGG代谢通路(pathway)富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于119和112个pathway,基因富集数最多的分别是氨基酸生物合成(127 unigenes)与核糖体(98 unigenes)。进一步分析表明球囊菌在胁迫中蜂幼虫肠道的过程中通过提高物质合成促进其增殖,而宿主通过抑制球囊菌的蛋白合成抵御病原入侵。富集在MAPK信号通路的11个DEGs的表达水平随着胁迫时间的延长而逐渐下降,推测中蜂幼虫通过抑制该通路而阻遏球囊菌增殖。【结论】本研究不仅为揭示白垩病过程中的球囊菌-中蜂幼虫互作提供了重要信息,也为阐明不同抗性蜂种的球囊菌抗性差异奠定了基础。     

侵染中华蜜蜂6日龄幼虫的蜜蜂球囊菌的微小RNA差异表达谱及调控网络

【目的】蜜蜂球囊菌(Ascosphaeraapis,简称球囊菌)是专性侵染蜜蜂幼虫的致死性真菌病原。MicroRNA(miRNA)作为一类重要的基因表达调控因子,能够广泛参与真菌及其宿主的相互作用过程。本研究通过比较分析球囊菌孢子(AaCK)和侵染中华蜜蜂(Apis cerana cerana,简称中蜂) 6日龄幼虫肠道内的球囊菌(AaT)的smallRNA(sRNA)组学数据对球囊菌的差异表达miRNA(differentiallyexpressed miRNA,DEmiRNA)、靶mRNA及二者间的调控网络进行全面解析,旨在揭示miRNA介导的球囊菌对中蜂幼虫的侵染机制。【方法】对于球囊菌侵染的中蜂6日龄幼虫肠道的small RNA-seq (sRNA-seq)数据,利用BLAST工具连续比对东方蜜蜂(Apiscerana)和球囊菌的参考基因组筛滤得到AaT的sRNA组学数据。分别将AaCK和AaT的sRNA组学数据比对miRBase数据库,对球囊菌侵染宿主前后miRNA的数量和结构特征进行分析。联用RNAhybrid+svm_light、Miranda和Tar...     

参与调控意大利蜜蜂工蜂中肠基因表达的东方蜜蜂微孢子虫miRNA的组学解析及其调控网络

[目的]本研究结合前期已获得的miRNA和mRNA组学数据对东方蜜蜂微孢子虫Nosema ceranae 的差异表达 miRNA(differentially expressed miRNA,DEmiRNA)靶向意大利蜜蜂 Apis mellifera ligustica 工蜂中肠的 mRNA 和差异表达 mRNA(d...     

蜜蜂球囊菌的microRNA鉴定及其调控网络分析

【目的】本研究利用small RNA-seq技术对球囊菌的纯培养进行测序,对球囊菌的micro RNAs miRNAs)进行预测、鉴定和分析,进而构建miRNAs-mRNAs的调控网络。【方法】利用Illumina Hiseq Xten平台对球囊菌菌丝与孢子进行测序,通过相关生物信息学软件对球囊菌的miRNAs进行预测和分析,通过茎环(Stem-loop)PCR对部分miRNAs进行鉴定,利用Cytoskype软件构建miRNAs-mRNAs的调控网络。【结果】本研究共获得48268696条clean reads,预测出118个球囊菌的miRNAs,它们的长度分布介于18–25 nt之间,不同长度的mi RNA的首位碱基偏好性差异明显。Stem-loop PCR验证结果显示共有10个miRNAs能够扩增出符合预期的目的片段,说明多数miRNAs可能真实存在。共预测出6529个球囊菌miRNAs的靶基因,其中5725个能够注释到Nr、Swissprot、KOG、GO和KEGG数据库。进一步分析结果显示有24个靶基因注释在MAPK信号通路。Cytoskype软件分析结果显示球囊菌的miRNAs与mRNAs之间存在复杂的调控网络,绝大多数的miRNAs处于调控网络的内部且同时结合多个mRNAs。【结论】本研究率先对球囊菌的miRNAs及miRNAs-mRNAs调控网络进行全面分析,研究结果丰富了对球囊菌miRNAs的认识,为其基础生物学信息提供了有益补充,也为阐明球囊菌致病的分子机理打下了一定基础。     

施氮对水稻土N2O释放及反硝化功能基因(narG/nosZ)丰度的影响

以紫潮泥和红黄泥两种不同质地的水稻土壤作为研究对象,通过室内培养试验,分析施用硝态氮肥对N2O释放和反硝化基因(narG/nosZ)丰度的影响,并探讨反硝化基因丰度与N2O释放之间的关系。结果表明,施用硝态氮显著增加两种水稻土的N2O释放量。在72h培养过程中,施氮改变了紫潮泥反硝化基因(narG/nosZ)的丰度,但并未明显影响红黄泥反硝化基因(narG/nosZ)丰度。通过双变量相关分析发现,除了紫潮泥narG基因外,其它的反硝化基因丰度和N2O释放之间并没有显著相关性。     

东方蜜蜂微孢子虫的基因结构优化及新基因鉴定

东方蜜蜂微孢子虫Nosema ceranae是一种寄生于蜜蜂中肠上皮细胞的单细胞真菌,对蜜蜂的健康危害严重,给世界各国的养蜂业造成较大损失。本研究基于前期获得的N.ceranae孢子的转录组数据对其已注释基因进行结构优化,并对未注释基因进行预测和分析。通过将测序得到的clean reads比对参考基因组和转录本重构,共对10个N.ceranae的已注释基因的5'端或3'端进行了延长。利用Cuffcompare软件将重构转录本与参考基因组进行比对,共鉴定出27个新基因,随机挑选9个新基因进行RT-PCR验证,均能扩增出符合预期的目的片段,表明预测出的新基因真实存在。有6个新基因能够注释到GO数据库和6个基因注释到KEGG数据库。进一步分析结果显示上述新基因注释到细胞等10个GO条目上,它们可能在N.ceranae的生命活动中具有重要功能。研究结果为N.ceranae的基因结构和功能注释信息的完善提供了有益补充,也为新基因的功能研究打下了基础。     

禾本科植物微卫星序列的特征分析和比较

微卫星或简单序列重复(simple sequence repeats,SSRs)广泛存在于真核生物的基因组中,是目前最有用的分子标记之一.本研究首次利用全基因序列对4个禾本科植物种,包括二穗短柄草(Brachypodium distachyon)、水稻(Oryza sativa)、高粱(Sorghum bicolor)和玉米(Zea mays)的SSR的类型、丰度(数量)、相对丰度、频率、长度和偏好性等进行了分析和比较.结果表明,禾本科植物中SSR的丰度与基因组大小成正比,而相对丰度与基因组大小相关不显著;三核苷酸和六核苷酸基序类型比其它类型更为丰富;不同物种皆呈现出明显的基序偏好性.总的来讲,除个别例外,禾本科植物的基因组中更倾向于富含A/T的基序.本研究结果为研究禾本科植物的基因组进化和SSR标记的开发提供了有价值的信息.     

意大利蜜蜂工蜂中肠响应Nosema ceranae胁迫的高表达基因分析

东方蜜蜂微孢子虫Nosema ceranae是专性寄生蜜蜂中肠上皮细胞的单细胞真菌,对意大利蜜蜂(意蜂)具有较强的侵染性。本研究利用RNA-seq技术对正常10 d (Apis mellifera ligustica control group,Am CK)和N. ceranae胁迫10 d的意蜂工蜂中肠(Apis mellifera ligustica treatment group,Am T)进行测序,共得到160 847 237 100条原始读段,过滤后得到1 066 955 298条有效读段。主成分分析结果显示Am CK与Am T测序样品的组内重复性较好,组间的基因表达模式差异明显。GO分类结果显示,Am CK与Am T的前100位高表达基因(HEGs)分别分布于32和33个GO terms,基因富集数最多的均为代谢进程。KEGG代谢通路(pathway)富集分析结果显示,Am CK的前100位HEGs富集在21个pathways,基因富集数最多的是内吞作用、信号通路和嘌呤代谢; Am T的前100位HEGs富集在26个pathways,基因富集数最多的是内吞作用、RNA转运和蛋白质的内质网加工。前100位HEGs的Venn分析结果显示Am CK与Am T的共有HEGs为87个,特有HEGs分别均为13个。研究结果揭示了N. ceranae胁迫意蜂工蜂中肠过程中的基因表达谱信息,也为解析N. ceranae致病的分子机理提供了有益信息和线索。     

ATPase旋转催化运动的随机跃迁动力学研究

ATP合酶既可在跨膜质子势的推动下催化合成ATP,也可以利用水解ATP释放的化学能而充当质子泵,把质子从线粒体基质中输送到内膜外侧,其能量转化效率却高得惊人,几乎达到100%。在旋转分子马达ATP合酶结构为基础上,结合随机主方程方法,提出了描述旋转分子马达ATPase合酶四态随机跃迁不等距旋转催化运动的理论模型;得到其角速度、扩散系数与ATP浓度之间的变化关系,并且得出了符合旋转分子马达生物机理的结果,定性半定量地解释了其动力学行为。     

转录组分析揭示东方蜜蜂微孢子虫侵染意大利蜜蜂的分子机制

【目的】本研究旨在通过差异表达基因(differentially expressed gene, DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制。【方法】基于前期已获得高质量的东方蜜蜂微孢子虫纯化孢子(NcCK)及侵染意大利蜜蜂工蜂7和10 d的东方蜜蜂微孢子虫(分别为NcT1和NcT2)转录组数据,根据P≤0.05且|log₂(Fold change)|≥1的标准,通过比较分析筛选出NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组的DEG。通过相关生物信息学软件对上述DEG进行Venn分析、GO分类和KEGG代谢通路富集分析。根据Nr和KEGG数据库注释信息和相关文献进行对东方蜜蜂微孢子虫的毒力因子和侵染相关因子的统计和分析。通过RT-qPCR验证转录组数据及DEG表达趋势。【结果】从NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组分别鉴定出1 397, 1 ...