西门塔尔牛肌内和皮下脂肪miRNA表达谱及miR-27b靶基因分析

【目的】通过分析西门塔尔牛肌内脂肪和皮下脂肪差异表达miRNA,及对重要的候选miRNA-27b进行靶基因预测及相关生物信息学分析,探索miRNA在动物肌内脂肪生长、发育过程中的作用及其调控机制。【方法】利用miRNA微阵列芯片技术对肌内脂肪和皮下脂肪miRNA表达谱进行检测和分析以筛选差异表达miRNA;qRT-PCR方法验证4个在肌内脂肪显著高表达的候选miRNAs;采用Targetscan生物信息学计算方法预测目标miR-27b（对肌内脂肪沉积可能起重要调节作用）的靶基因,并对靶基因进行GO注释描述、富集分析及KEGG富集。【结果】miRNA芯片分析发现肌内脂肪和皮下脂肪共有88个差异表达极显著miRNA（P＜0.01）。候选的4个在肌内脂肪高表达miRNAs（miR-140、miR-145、miR-143、miR-27b）的qRT-PCR检测结果与芯片分析结果具有很好的一致性。KEGG富集显示,目标miR-27b的预测靶基因在MAPK、Wnt、Hedgehog、TGF-beta、GnRH等信号通路中显著富集,其中在MAPK信号通路中富集度最高。【结论】牛肌内脂肪和皮下脂肪组织中确实存在一些差异表达的miRNAs,某些重要的在差异高表达miRNA（如miR-27b）可能以独特的通路（如MAPK信号）来调节牛肌内脂肪的形成过程。     

家兔肌肉组织中保守miRNA的鉴定及进化分析

【目的】旨在鉴定家兔肌肉组织中的保守miRNA，获取其表达谱并分析其进化特征，进而为研究保守性miRNA在家兔生长发育过程中的生物学功能提供理论依据。【方法】采用Solexa技术对家兔背最长肌组织中miRNA进行了测序和分析。【结果】通过序列比对鉴定出521个保守的miRNA序列，分属于238个不同的家族。其中miR-1、miR-133、miR-206和miR-499等肌肉特异性表达的miRNA在肌肉中均有表达，且除miR-499外的其它3个miRNA在所测定样本中均高表达。10个动物保守性miRNA家族在S86专门化父系肌肉基因组中出现了丢失。靶基因显著富集于13条信号通路，其中Focal adhesion和PI3K-Akt信号通路与肌肉发育密切相关。【结论】研究鉴定的影响肌肉发育的miRNA及其靶基因富集通路，将有助于家兔生长调控机制的研究，加快可用于标记辅助选择的分子标记的开发。     

1月龄五指山猪与长白猪骨骼肌miRNA转录组比较

为研究microRNA(miRNA)在五指山猪和长白猪肌肉生长发育中的差异和探究差异miRNA对骨骼肌发育转录后调控的影响,以1月龄的五指山猪和长白猪背最长肌为材料,通过转录组测序和生物信息学分析筛选与骨骼肌发育相关的差异miRNA。结果显示,从2个文库中共鉴定获得318种己知的猪miRNA,五指山猪和长白猪分别鉴定出312个和306个已知miRNA,同时分别发现了83个和88个新的miRNA。获得17个差异显著的miRNA,其中16个miRNA表达上调,1个miRNA表达下调。靶基因预测和生物学功能预测发现17个差异miRNA,共预测到535个靶基因,靶向作用于mTOR信号通路和肌动蛋白细胞骨架调控的ssc-miR-486,ssc-miR-204,ssc-miR-338和ssc-miR-885-3p 4个miRNA可能参与骨骼肌生长发育过程。     

热应激奶牛血清miRNA表达谱及miR-181a靶基因分析

通过分析热应激奶牛和正常奶牛血清差异表达miRNA,及对重要的候选miRNA进行靶基因预测及相关生物信息学分析,探索miRNA在荷斯坦奶牛发生热应激过程中的作用及其调控机制。利用miRNA高通量测序技术对奶牛血清中miRNA表达谱进行检测和分析以筛选差异表达miRNA;用RT-qPCR方法验证4个在热应激奶牛血清与正常奶牛血清中显著差异表达的候选miRNAs;采用Targetscan生物信息学计算方法预测目标miR-181a(对应激以及免疫反应可能起重要调节作用)的靶基因,并对靶基因进行GO注释描述、富集分析及KEGG富集。结果表明:miRNA高通量测序分析发现共有52个差异表达极显著miRNA(P0.01)。候选的4个miRNAs(miR-19a、miR-19b、miR-181a、miR-1246)的RT-qPCR检测结果与测序分析结果具有很好的一致性。GO和KEGG分析结果显示,miR-181a的靶基因与应激反应以及免疫功能密切相关,并且靶基因在B细胞和T细胞受体信号通路中显著富集。结论:在热应激奶牛与正常奶牛血清中存在差异表达的miRNAs,某些重要的miRNA(如miR-181a)可能以独特的通路(如B细胞和T细胞受体信号通路)来调节奶牛热应激的发生过程。     

两品种鸭下丘脑miRNA的鉴定及其表达验证分析

为了鉴定与鸭性成熟启动相关的miRNA,根据鸭卵巢组织的发育情况,利用Solexa技术分别对105日龄连城白鸭和樱桃谷鸭的下丘脑组织进行小RNA测序分析。两个小RNA文库共检测到939个成熟miRNA,共有miRNA数量806个。以樱桃谷鸭为对照,两个文库之间总共检测到290个差异表达miRNA,其中上调miRNA 85个,下调miRNA 205个。靶基因预测分析结果表明,较多的靶基因被注释到神经活性配体-受体相互作用、黏着斑、MAPK信号通路、胞吞作用等通路中。定量PCR结果表明,2个差异表达miRNA靶基因的mRNA表达量在两品种鸭下丘脑组织中存在着极显著差异。     

育成期高能饲喂下开产与未开产蛋鸡肝脏miRNA表达谱差异分析

【目的】对育成期高能饲喂下开产与未开产蛋鸡的肝脏进行miRNA高通量测序分析,探明高能饮食状态下蛋鸡肝脏中影响其开产的miRNA,为提高产蛋性能打下基础。【方法】以代谢能水平为12.14 MJ/kg的高能饲粮饲喂育成期蛋鸡,通过Illumina NextSeq500高通量测序平台对开产与未开产蛋鸡肝脏进行small RNA测序,使用DESeq分析miRNA表达量,并以实时荧光定量PCR进行验证;采用miRanda对差异表达miRNA进行靶基因和靶位点预测,同时以超几何分布对差异表达的靶基因进行GO功能注释分析和KEGG信号通路富集分析。【结果】 6个样本（3个开产蛋鸡样本,3个未开产蛋鸡样本）能注释到的miRNA均超过300个,约占miRBase中已鉴定鸡miRNA的30.00%,且前体鉴定结果显示各样本注释到的miRNA有部分定位在同一前体上。与未开产组样本相比,开产组样本有9个miRNA表达上调、3个miRNA表达下调。其中,gga-miR-132c-5p、gga-miR-132b-5p、gga-miR-2184a-5p、gga-miR-132c-3p、gga-miR-132b-3p及gga-miR-2184a-3p属于同一miRNA基因簇,且gga-miR-1682、gga-miR-132b-3p和gga-miR-2184a-3p等3个上调miRNA在未开产蛋鸡样本中不表达。通过miRanda共预测得到129个差异表达潜在靶基因,以miR-203a预测得到的靶基因最多,为32个;其次是gga-miR-2184a-5p、gga-miR-148a-5p和gga-miR-12211-5p,分别预测到30、25和23个靶基因;而miR-132c-5p预测得到的靶基因最少,仅为9个。129个靶基因显著注释到8条GO功能条目上,生物学过程（Biological process）主要注释到脂质代谢过程、细胞脂质代谢过程、脂质生物合成过程、磷脂生物合成过程、磷脂代谢过程、甘油磷脂生物合成过程和解剖学结构发育,分子功能（Molecular funcion）仅注释到1条GO功能条目,即利钠肽受体活性,未发现涉及细胞组分（Cellular component）的GO功能条目;在注释到的8条GO功能条目中有6条与脂质代谢相关,涉及的靶基因有22个,占总潜在靶基因的17.10%。KEGG信号通路富集分析结果显示共显著富集到7条KEGG信号通路,其中脂肪酸降解、脂肪酸生物合成、酮体合成与降解及PPAR信号通路等4条信号通路与脂质代谢相关。【结论】育成期高能饲喂下开产与未开产蛋鸡中共存在12个差异表达miRNA,涉及129个差异表达潜在靶基因,且主要富集在肝脏脂质代谢相关过程和信号通路,说明miRNA是通过调控脂质代谢及其相关基因表达而影响蛋鸡开产。     

禽致病性大肠杆菌phoP/Q基因缺失对雏鸡脾脏microRNAs功能的影响

从microRNA(miRNA)水平探究禽致病性大肠杆菌(APEC)phoP/Q基因缺失对雏鸡脾脏功能表达的影响,为APEC的防制提供参考资料。用禽致病性大肠杆菌和其phoP/Q基因的缺失株分别攻毒14日龄雏鸡,采集其脾脏组织为材料进行高通量测序获得miRNA表达谱,选择差异倍数大于8的miRNA进行Real-time PCR,对差异表达的miRNA进行靶基因预测以及GO、KEGG富集分析。测序结果获得10个差异表达的miRNA,其中1个miRNA表达量下调,9个miRNA表达量上调,Real-time PCR结果表明,miRNA变化趋势与测序结果一致,GO分析结果表明,差异表达基因主要富集在能量代谢、免疫系统、细胞生长与凋亡、糖合成与代谢等生物学过程。KEGG分析结果显示预测靶基因参与MAPK信号通路、糖代谢通路、Wnt信号通路等重要通路。本试验分析禽致病性大肠杆菌phoP/Q基因缺失后感染雏鸡脾脏miRNA的表达差异,从miRNA水平探讨了phoP/Q基因缺失对APEC致病性的影响,为APEC的防治提供参考资料。     

低致病性A/Anhui/1/2013(H_7N_9)流感病毒感染小鼠肺组织的miRNA表达谱分析

[目的]本试验旨在利用基因芯片测序技术分析宿主miRNA在调控流感病毒感染小鼠肺组织过程中的作用。[方法]应用低致病性H_7N_9亚型A/Anhui/1/2013流感病毒及PBS感染BALB/c小鼠,感染后3 d取小鼠肺组织分别利用转录组测序技术和miRNA芯片测序技术筛选差异表达mRNA和miRNA,并利用RT-qPCR验证差异miRNA。之后分别采用Targetscan、PITA及microRNAorg软件预测靶基因并结合转录组mRNA信息筛选候选miRNA的假定靶基因,最后利用Gene Ontology(GO)和Kyoto Encyclopedia of Gene and Genomes(KEGG)分析差异miRNA的生物学功能及其可能调控的信号通路。[结果]与PBS组相比,低致病性H_7N_9病毒感染组共筛选到265个差异表达miRNA,其中143个miRNA显著上调,122个miRNA显著下调。差异表达miRNA经RT-qPCR验证,10条候选miRNA的RT-qPCR结果与芯片结果有非常好的一致性。进一步KEGG分析表明,这些差异表达miRNA主要富集在Rap1、PI3K-Akt、Hippo、MAPK、Wnt、黏着斑、自噬等免疫相关信号通路中。结合差异mRNA信息进行miRNA-mRNA调控网络分析,显示主要有15条差异表达miRNA和31个差异靶基因富集到这些信号通路中。[结论]成功筛选到了低致病性H_7N_9感染小鼠肺组织后引起的差异表达miRNA,且RT-qPCR证明其与基因芯片测序结果表达趋势具有很好的一致性,为深入研究宿主miRNA在调控低致病性H_7N_9流感病毒与宿主相互作用的分子机制奠定了基础。     

捻转血矛线虫丙硫咪唑敏感虫株和耐药虫株miRNA-mRNA关联比较分析

旨在比较分析捻转血矛线虫敏感虫株和耐药虫株miRNA表达谱,探讨miRNA与捻转血矛线虫丙硫咪唑耐药性相关基因之间的调控机制。运用Illumina Hiseq2000平台进行测序并进行cDNA文库的构建,利用RNAhybrid、Miranda和TargetScan对测序得到的数据进行靶基因预测,并取交集作为miRNA的靶基因预测结果,使用Bowtie、RepeatMasker、MIREAP、Rfam、miRBase和DAVID等生物信息学分析软件和权威数据库系统,筛选出差异的miRNA以及与捻转血矛线虫耐药相关的miRNA并进行靶基因预测,同时对KEGG pathway和GO功能富集到的靶基因和可能参与调控的通路进行预测分析。结果表明,敏感虫株和耐药虫株共筛选显著差异表达的miRNA共294个,其中113个上调,181个下调。对miRNA和其靶向调节的mRNA进行关联分析,共关联到1 770个miRNA-mRNA对为负调控关系,其中涉及到274个差异的miRNA和603个差异的mRNA。显著差异表达的miRNA靶基因被注释到了1 430条GO terms和327条KEGG pathway,富集到FoxO信号通路(FoxO signaling pathway),mTOR信号通路(mTOR signaling pathway),PI3K-Akt信号通路(PI3K-Akt signaling pathway)等与耐药相关的一些抗性通路。综上,通过对miRNA表达谱分析以及miRNA-mRNA靶向分析,发现捻转血矛线虫敏感虫株和耐药虫株显著差异表达的miRNA和其对应的靶基因,并对部分耐药相关通路中存在的靶基因及显著差异的靶基因所在的通路进行筛选和分析,为后续进一步通过转录组学深层次挖掘捻转血矛线虫产生耐药性的关键基因及相关调控分子提供科学依据。     

中国美利奴羊羊毛弯曲相关microRNA的筛选

【目的】 研究miRNA在不同羊毛弯曲中国美利奴羊皮肤组织中的表达差异及其可能参与的调控通路,筛选出与中国美利奴羊羊毛弯曲相关的miRNA,为进一步探究细毛羊羊毛弯曲性状的调控机理提供理论依据。【方法】 运用高通量测序技术对羊毛弯曲频率有极端差异表型值(大弯曲组与小弯曲组)细毛羊个体皮肤的miRNA组进行测序与比较分析。【结果】 共鉴定出425个miRNAs,其中包括143个已知miRNAs和282个新发现的miRNAs。在大弯曲组与小弯曲组共筛选到8个上调和17个下调的miRNAs。对差异miRNA的预测靶基因进行简单的Gene Ontology与KEGG Pathway富集分析。差异表达的miRNA的靶基因主要参与跨膜信号受体活动、信号转导活动、分子转导活动、膜的组成等过程。15 761个靶基因注释到252个信号通路。【结论】 筛选出了25个与细毛羊羊毛弯曲性状相关的候选miRNAs。     

意大利蜜蜂幼虫肠道在球囊菌侵染前期的 差异表达microRNA及其调控网络

【目的】 微小RNA（microRNA,miRNA）是一类重要的基因表达调控因子,可影响宿主与病原间的互作过程。蜜蜂球囊菌（Ascosphaera apis）是一种特异性侵染蜜蜂幼虫的致死性真菌病原。本研究旨在对意大利蜜蜂（Apis mellifera ligustica,简称意蜂）幼虫肠道在球囊菌胁迫前期的差异表达miRNA（differentially expressed miRNA,DEmiRNA）及其靶基因进行深入分析,在miRNA组学水平探究意蜂幼虫在球囊菌侵染前期的胁迫应答,并通过构建显著DEmiRNA的调控网络筛选出与宿主应答相关的关键miRNA。【方法】 利用small RNA-seq（sRNA-seq）技术对正常及球囊菌侵染的意蜂4日龄幼虫肠道（AmCK和AmT）进行高通量测序,首先对原始数据进行质控和评估,随后将过滤后的数据与西方蜜蜂（Apis mellifera）参考基因组进行比对;将比对上的序列标签（tags）注释到miRBase数据库,得出已知miRNA的表达量;通过TPM（tags per million）算法对各样本中miRNA的表达量进行归一化处理,以|log₂ fold change|≥1且P≤0.05作为标准筛选得到显著DEmiRNA;利用TargetFinder 软件预测显著DEmiRNA的靶基因,并对其进行GO和KEGG代谢通路（pathway）富集分析。利用Cytoscape软件对miRNA-mRNA调控网络进行可视化。最后,利用茎环反转录PCR（Stem-loop RT-PCR）和荧光定量PCR（qPCR）验证测序数据的可靠性。【结果】 AmCK和AmT样品的测序分别得到13 553 302和10 777 534条原始读段（raw reads）,经严格过滤后得到的有效读段（clean reads）数分别为13 186 921和10 480 913条。各样品的生物学重复间的Pearson相关性系数分别在0.9822和0.9508以上。共有10个显著DEmiRNA,包括4个上调miRNA和6个下调miRNA。显著DEmiRNA在AmT的整体表达水平低于AmCK。10个显著DEmiRNA可靶向结合3 788个靶基因,其中上调miRNA的1 240个靶基因可注释到GO数据库中的39个GO条目,主要富集在结合、细胞进程、代谢进程和应激反应等;下调miRNA的749个靶基因可注释到34个GO条目,主要富集在细胞进程、结合、代谢进程和应激反应等。KEGG数据库注释结果显示,上调miRNA和下调miRNA的靶基因分别注释到95和66条代谢通路,富集基因数最多的分别是Wnt信号通路、Hippo信号通路、光传导以及内吞作用、磷脂酰肌醇信号系统、嘌呤代谢。对于上调和下调miRNA,分别有31和52个靶基因注释到内吞作用,15和7个靶基因注释到泛素介导的蛋白水解,11和5个靶基因注释到Jak-STAT信号通路,1和3个靶基因注释到MAPK信号通路。显著DEmiRNA与靶mRNA之间形成复杂的调控网络,7个显著DEmiRNA靶向结合96个与Wnt信号通路相关的mRNA,8个显著DEmiRNA靶向结合55个与内吞作用相关的mRNA。Stem-loop RT-PCR和qPCR结果验证了测序数据的可靠性。【结论】 对意蜂幼虫肠道在球囊菌侵染前期的DEmiRNA及其靶基因进行预测和分析,并构建和分析了DEmiRNA-mRNA调控网络,研究结果提供了宿主miRNA的表达谱和差异表达信息,揭示了DEmiRNA通过调控细胞生命活动、新陈代谢以及部分细胞和体液免疫等生物学过程参与宿主的胁迫应答。miR-4331-y、miR-4968-y、miR-8440-y、novel-m0023-5p和novel-m0025-3p共同参与了宿主的Wnt信号通路和内吞作用的调控,可作为白垩病治疗的潜在分子靶标。     

捻转血矛线虫阿苯达唑耐药虫株给药前后miRNA差异表达谱及网络调控分析

旨在通过miRNA测序技术对捻转血矛线虫阿苯达唑给药前后的耐药虫株差异表达miRNA及网络调控分析,获取其转录本中全miRNA表达谱及耐药相关差异表达基因和miRNA,探究miRNA在捻转血矛线虫阿苯达唑耐药虫株给药前后的差异。采用Illumina Hiseq2 000平台对捻转血矛线虫阿苯达唑耐药虫株给药前和给药后的虫体进行Small RNA-seq测序,将Raw Data与参考基因组和miRBase数据库进行比对分析,使用RNAhybrid、Miranda和TargetScan对miRNA靶基因进行预测,筛选捻转血矛线虫耐药虫株给药前后表达量显著差异的miRNA并通过Stem-loop qRT-PCR验证,同时对差异显著的miRNA-mRNA调控关系采用Cytoscape软件进行可视化分析。将靶基因向KEGG和GO两大权威数据库映射,预测和分析靶基因功能注释情况和参与调控的通路情况。结果显示：1）捻转血矛线虫阿苯达唑耐药虫株给药前和给药后共筛选出显著差异表达的miRNA 188个,其中125个上调,63个下调,且qRT-PCR结果与转录组测序结果表达情况一致。2）显著差异表达的14...     

NDV基因Ⅶ型与基因Ⅱ型感染CEF后microRNA表达谱分析

本研究旨在探究新城疫病毒（Newcastle disease virus,NDV）流行株HB株（基因Ⅶ型）与疫苗株La Sota株（基因Ⅱ型）在感染鸡胚成纤维细胞（Chicken embryo fibroblasts,CEF）36 h后宿主细胞microRNA(miRNA)表达谱变化情况。通过高通量测序技术分别对La Sota株与HB株感染CEF后36 h后的样本与参考鸡源基因组进行系统比对，筛选差异表达miRNA并对其靶基因进行GO功能注释以及KEGG通路富集。最后随机筛选5条miRNA通过实时荧光定量PCR(Real-time quantification PCR,RT-qPCR)进行验证。结果表明，La Sota与HB感染组中已知miRNA分别为56.9%和47.8%;La Sota与HB感染组通过差异性比较后共筛选出10条差异表达miRNA；通过对差异表达基因的靶基因进行GO与KEGG富集发现，其主要在细胞代谢、生物合成以及应激信号传导等通路上出现显著富集；RT-qPCR结果与RNA-seq结果趋势一致，证明了测序结果的可靠性。为进一步揭示不同NDV毒株感染CEF以及在感染过程中miRNA发挥的调控作用提供了研究基础。     

基于高通量测序分析肺炎支原体感染姜曲海猪肺组织miRNA表达谱

为探明姜曲海猪感染肺炎支原体后肺组织miRNA(microRNA)表达谱和分子机制,选取50日龄姜曲海猪为实验猪,随机分成感染组和对照组,人工感染肺炎支原体28 d后,采集肺部组织,进行高通量miRNA测序,采用生物信息学软件进行miRNA鉴定和靶基因预测。结果显示:感染组和对照组的肺组织分别筛选到14 265 786条和14 000 588条小RNA纯净序列(clean reads)。与对照组相比,感染组中筛选到73个显著差异表达的miRNAs,其中39个表达上调,34个表达下调,从中随机选取4个miRNAs进行定量PCR(quantitative real time PCR,qRT-PCR)验证,感染组和对照组的表达水平与测序结果基本一致。73个差异表达miRNAs预测到1 685个靶基因和4 220个靶位点,靶基因预测筛选到8个与免疫调控相关的miRNAs。靶基因GO(gene ontology)分析显示,miRNA广泛参与生物过程、细胞组成和分子功能的调控。靶基因KEGG(Kyoto encyclopedia of genes and genomes)分析显示,miRNA参与调控细胞凋亡、ECM受体相互作用、粘附斑通路等信号通路。本研究通过高通量测序获得姜曲海猪感染肺炎支原体后肺组织miRNA表达谱,通过生物信息学分析筛选到与免疫调控相关的miRNAs和信号通路,为进一步阐明姜曲海猪的肺炎支原体感染机制奠定了基础。     

北京油鸡和来航鸡脾脏差异表达microRNA的鉴定与分析

【目的】明确北京油鸡和来航鸡脾脏重量、组织结构及其miRNA表达谱的差异,筛选与两个鸡种免疫应答能力差异相关的候选基因,为家禽抗病育种研究奠定基础。【方法】选取1日龄同期孵化的北京油鸡和来航鸡母雏各45只作为试验材料,在相同营养水平和环境条件下饲养,42日龄测量体重和脾脏重。每个品种随机选取3只个体（体重接近群体均值）,取其一半脾脏组织制作石蜡切片,利用H.E.染色,显微镜下观察脾脏组织结构、脾小体的数量和动脉周围淋巴鞘厚度。提取另一半脾脏组织中的小RNA,等量混合组成2个RNA池,构建cDNA文库,利用Solexa平台进行高通量测序。测序结果与参考基因组数据库比对获得表达基因。利用Mireap软件预测新miRNAs基因。利用RPKM（reads per kb per million reads）方法计算基因的表达量,根据FDR（false discovery rate）＜0.001和|log₂ ratio(T/CK)|≥1的标准筛选差异表达的基因。利用TargetScan和Pictar软件预测差异表达miRNA的靶基因,并对其进行GO和KEGG数据库功能注释。利用DAVID软件对靶基因蛋白进行富集分析。将预测靶基因与已知脾脏重和脾脏指数相关QTL区域注解的基因取交集。【结果】在42日龄,北京油鸡的体重和脾脏重均显著高于来航鸡（P＜0.01）。组织学观察显示,两个鸡种的脾脏组织结构发育完整,可见清晰的白髓、红髓、脾小结和动脉周围淋巴细胞鞘。与来航鸡相比,北京油鸡脾小体的数量更多、动脉周围淋巴鞘更厚,鞘内淋巴细胞也更为密集。高通量测序在两组样本共鉴定出486种已知miRNA 、130个候选miRNA 和216个共表达miRNA 。两组间显著差异表达的miRNA共计35种,其中在北京油鸡中有8个表达上调,27个表达下调,它们的变化倍数在1.002-10.41之间。差异基因表达丰度在前5位的miRNA是miR-2954、miR-6606-5p、miR-146b-5p、miR-21-3p和miR-21-5p,对其2 767个靶基因的生物信息学分析结果显示,它们参与了T和B淋巴细胞的分化、免疫器官的发育、蛋白质磷酸化和细胞形态发生等生物学过程,并在泛素介导的蛋白水解、凋亡和磷脂酰肌醇信号系统等14个通路中显著富集。miR-6606-5p、miR-146b-5p的共同靶基因BLM和miR-2954的靶基因IRF8均定位于脾重或脾指数相关的QTL区域。【结论】（1）42日龄北京油鸡和来航鸡脾重、组织结构及其miRNA表达谱的确存在一些差异。（2）某些重要的差异高丰度表达miRNA(如miR-2954、miR-6606-5p、miR-146b-5p 、miR-21-3p和miR-21-5p ) 可能通过对靶基因调控,来参与调节T、B淋巴细胞分化、增殖、激活等过程,进而影响脾脏的重量、组织结构及其免疫应答能力。     

新疆野苹果幼苗叶片响应NaCl胁迫的miRNAs及靶基因筛选

【目的】 筛选新疆野苹果响应NaCl胁迫相关差异miRNA及靶基因。【方法】 以新疆野苹果组培苗为材料,对150 mmol/L NaCl处理6和48 h时的新疆野苹果幼苗叶片进行small RNA测序。【结果】 NaCl处理6 h时3个miRNA差异表达,其中2个上调表达,1个下调表达,miRNA对应的靶基因数目为64个;NaCl处理48 h时13个miRNA差异表达,其中4个上调表达,9个下调表达,miRNA对应的靶基因数目为108个。对差异靶基因进行GO和KEGG分类注释,NaCl处理6 h时,GO主要富集包括:肌醇磷酸2激酶活性、寡糖基转移酶活性、蛋白丝氨酸/苏氨酸磷酸酶抑制剂活性、木聚糖生物合成过程等;NaCl处理48 h时,GO主要富集包括:氧氧化还原酶活性、DNA结合、质外体、次生代谢过程等。KEGG共同富集在N-glycan生物合成通路。【结论】 mdm-miR168b、mdm-miR159c、mdm-miR827、mdm-miR390f、mdm-miR171i、mdm-miR399j基因在NaCl处理6和48 h时表达量存在差异,其中mdm-miR390f及其靶基因HF10976、mdm-miR171i及其靶基因HF33844、mdm-miR399j及其靶基因HF11095表达量呈负相关,3个miRNA参与了新疆野苹果对NaCl胁迫的响应过程。     

中华蜜蜂幼虫肠道中微小RNA的鉴定及分析

[目的]利用small RNA-seq(sRNA-seq)技术和生物信息学方法对中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫肠道的微小RNA(microRNA,miRNA)进行全转录组鉴定和分析,旨在丰富中蜂的miRNA信息,并为深入研究miRNA调控中蜂幼虫肠道发育的分子机理提供依据.[方法]利用s...     

蜜蜂球囊菌中长链非编码RNA的调控作用

【目的】长链非编码RNA（long non-coding RNA,lncRNA）是一类二、三级结构高度保守,长度>200 nt且不具蛋白编码能力的RNA,在转录和转录后水平广泛参与调控剂量补偿、细胞分化和生长发育等生命活动。本研究基于前期获得的蜜蜂球囊菌（Ascosphaera apis,简称球囊菌）菌丝和孢子混合样品的高质量lncRNA组学数据进行球囊菌lncRNA的顺式（cis）作用、反义lncRNA（antisense lncRNA）作用和竞争性内源RNA（competing endogenous RNA,ceRNA）作用的分析和探讨,以期揭示lncRNA在球囊菌中的潜在功能。【方法】基于lncRNA基因在染色体上的位置,预测lncRNA上下游100 kb以内的蛋白编码基因;使用Blast软件将上下游基因比对到GO和KEGG数据库,以获得功能和通路注释。利用LncTar软件对反义lncRNA的靶mRNA进行预测,并使用Blast软件将上述靶mRNA比对到KEGG和eggNOG数据库。利用TargetFinder软件预测lncRNA靶向结合的miRNA及miRNA靶向结合的mRNA,根据靶向结合关系建立lncRNA-miRNA和lncRNA-miRNA-mRNA调控网络,进而通过Cytoscape v3.7.1软件进行调控网络的可视化。利用RT-PCR对调控网络中的lncRNA、靶miRNA和靶mRNA进行表达验证。【结果】共预测出371个lncRNA的5 852个上下游基因,可注释到细胞进程、代谢进程和应激反应等48个功能条目,以及新陈代谢途径、次生代谢产物的生物合成和抗生素的生物合成等121条通路,表明部分lncRNA可通过顺式作用调节上下游基因的表达,从而参与调控球囊菌的生长发育及物质能量代谢等基础生命活动。球囊菌的7个lncRNA与7个靶mRNA存在序列互补关系,其中5个mRNA在eggNOG数据库中仅注释为假定蛋白,仅gene3444在KEGG数据库注释为核孔复合体蛋白An-Nup120和假定蛋白,表明上述1个反义lncRNA可能参与调控核孔复合体蛋白的生物合成等生物学过程。此外,共预测出227个lncRNA与73个miRNA之间存在靶向结合关系,其中多数lncRNA(79.02%)仅能结合1—2个miRNA,部分miRNA可被多个lncRNA靶向结合;进一步构建氧化磷酸化通路和MAPK信号通路相关的lncRNA-miRNA-mRNA调控网络,分析结果显示氧化磷酸化通路相关的222个lncRNA靶向78个miRNA及50个mRNA,MAPK信号通路相关的222个lncRNA靶向76个miRNA及46个mRNA,表明部分lncRNA通过ceRNA作用调控此两条通路,从而影响球囊菌的能量合成、环境适应以及生长发育等过程。【结论】部分lncRNA可能通过顺式作用和ceRNA作用调节球囊菌的生长发育、物质能量代谢以及环境适应等生物学过程;MSTRG.5393.1可能作为反义lncRNA调控球囊菌的核孔复合体的蛋白合成。