热应激奶牛血清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表达谱及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的调控功能提供参考,并为树木低氮营养高效利用育种提供重要的候选基因。方法根据miRNA的保守性及与靶基因的严谨互补配对关系,以杨树miRNA为探针利用靶基因预测软件psRNATarget,通过与毛白杨转录组的基因序列进行比对鉴定靶基因,进一步开展毛白杨受低氮胁迫后靶基因的差异表达分析及功能注释。结果获得了131个miRNA家族的242个miRNA成员对应的3 024个靶基因,分别参与了植物激素信号转导、次生代谢产物的生物合成、氨基酸合成代谢、碳代谢和RNA运输等通路。57个靶基因在低氮胁迫处理后发生显著变化,其中受到诱导（29个）和抑制（28个）的基因数目相当。14个低氮胁迫响应的miRNA,其对应的11个靶基因也发生了显著的差异表达变化,其中miRNA和靶基因表达量发生相反变化的有8个miRNA。本研究发现参与植物激素信号转导的靶基因（2个）及参与代谢途径的靶基因（6个）发生了差异表达。miR162的靶基因编码ABC转运蛋白,miR393运用于靶基因KAT2调节Na+和K+动态平衡,miR399的靶基因PIF3编码光敏色素互作因子PIFs蛋白,这些miRNA及靶基因可能在杨树响应低氮胁迫中发挥重要作用。结论本文鉴定到了毛白杨中一批低氮胁迫响应miRNA的靶基因,可调控杨树对氮逆境胁迫信号的反应。这些miRNA及靶基因为进一步揭示miRNA及靶基因在低氮胁迫下的调控功能提供了研究线索,为树木氮营养的高效利用改良提供了重要候选基因。      

家蚕miRNA对核型多角体病毒靶基因的预测

了解家蚕miRNA对家蚕核型多角体病毒(BmNPV)的调控作用,为家蚕病毒的防控研究提供参考,利用已发布的563条家蚕成熟miRNA预测其在家蚕核型多角体病毒(BmNPV)中的靶基因,并利用在线靶基因预测软件RNAhybrid进行预测。结果表明:预测得到55个可能的靶基因,功能涉及家蚕体内病毒的复制增殖、病毒早期及晚期的转录和表达、抗家蚕细胞凋亡以及与病毒的结构蛋白相关等。miRNA与病毒的相互作用存在于病毒的侵染到释放的整个过程中,可能是一种重要的调控方式。     

热应激对羊成肌细胞部分miRNA的影响及其生物信息学分析与靶基因预测

[目的]本文旨在探究热应激对羊成肌细胞部分miRNA的影响,并对miR-23a、miR-24、miR-27a-3p和miR-30a-5p进行生物信息学分析与靶基因预测,揭示热应激下miRNA影响骨骼肌发育的调控机制。[方法]通过RT-qPCR技术研究热应激对羊成肌细胞miRNA表达量的影响;利用miRBase数据库获取不同物种中miR-23a、miR-24、miR-27a-3p和miR-30a-5p的前体序列和成熟序列,并利用MEGA X软件构建系统进化树进行相似性分析;通过TargetScan、miRDB、RNA22、miRWalk对上述miRNA进行靶基因预测,使用DAVID对预测出的共同靶基因进行GO功能与KEGG通路富集分析,最后进行miRNA、共同靶基因与KEGG通路的关联分析。[结果]热应激显著下调羊成肌细胞在增殖阶段与分化阶段miRNA的表达。miR-23a、miR-24、miR-27a和miR-30a在不同物种间的保守性较高,除了牛的miR-24以外,成熟序列基本相同,种子序列完全一致。预测出的共同靶基因主要富集在细胞增殖、分化、周期、凋亡、迁移等GO条目,以及mTOR...     

意大利蜜蜂幼虫肠道在球囊菌侵染前期的 差异表达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信号通路和内吞作用的调控,可作为白垩病治疗的潜在分子靶标。     

奶牛淋巴白血病中细胞凋亡与细胞增殖病理研究

利用细胞凋亡试剂盒（POD）及小鼠抗大鼠增殖细胞核抗原（PCNA）单克隆抗体将淋巴白血病中的凋亡细胞与增殖细胞准确地染色显示出来,并进行对比观察分析.结果表明,奶牛淋巴白血病中有凋亡细胞与增殖细胞,细胞凋亡与细胞周期调控是肿瘤发生发展的两方面.单克隆抗体免疫细胞化学方法可以很好地检测PCNA和细胞凋亡的表达,PCNA与细胞凋亡可作为奶牛淋巴白血病细胞增殖活性的指标,其表达可能是急性白血病预后不良的一个指标.可以认为牛淋巴白血病乃是细胞死亡被抑制,而细胞的增殖率相对升高的病理过程.     

藏绵羊不同毛色皮肤组织miRNA表达谱及靶基因分析

【目的】藏绵羊又称藏系绵羊,是中国三大粗毛羊品种之一,主要分布在西藏及其毗邻的高寒牧区,如青海、甘肃、四川、云南等地。藏绵羊由于其被毛纤维粗长,毛质细腻且保暖性好,是制造藏式地毯的上好原料。而毛色作为一种重要的经济性状,同时制约羊毛的经济价值。然而,目前对于绵羊毛色的分子调控机制尚不明确,以藏绵羊为研究对象,对其不同颜色皮肤组织（黑色和白色）进行转录组测序,旨在探讨miRNA在藏绵羊不同毛色皮肤组织转录后水平中发挥的作用及其可能参与的调控通路。【方法】选取2只具有黑白花毛色的健康藏绵羊,减去体侧被毛使皮肤完全裸露并用75%酒精消毒处理,屠宰后快速切取皮肤组织并保存在组织样品保存液中防止RNA降解,提取RNA后通过miRNA测序和生物信息学分析,获得藏绵羊不同毛色（黑色和白色）皮肤组织miRNA表达谱,筛选组织差异表达miRNA并预测相关靶基因。通过靶基因分析,获得与调控藏绵羊毛色性状可能相关的信号通路。【结果】黑色和白色皮肤组织共获得85.76兆条原始读长,经过滤后得到85.08兆条clean reads;共鉴定出334个已知的miRNA和59个新的miRNA;从中获得23个差异表达miRNA,其中14个差异表达miRNA在白色皮肤组织中表达显著上调,9个表达下调。miR-2284b和miR-744仅在藏绵羊白色皮肤中表达,而miR-23b、miR-411a-5p、miR-30c、miR-423-3p和miR-324-5p则仅在藏绵羊黑色皮肤中表达,但表达量相对较低。表达量最高的miRNA为miR-10a,而倍数差异最大的为miR-23b。其中,miR-10a可参与调节多种信号通路,其作用涉及增殖、分化、凋亡、细胞粘附等各个方面,目前并无调控绵羊毛色方面的相关报道,但其功能及作用机制的研究仍在不断扩展。miR-23b则与Wnt、Notch等信号通路有关,这些信号通路中的部分基因也和黑色素的生成有密切关系。结合倍数差异和miRNA表达量分析,miR-411a-5p、miR103和miR-200b可能也和毛色调控密切相关。本研究共预测出981个靶基因可能参与毛色调控,多数基因和WNT信号通路、MAPK通路、EDNRB信号通路及黑素原生成通路（melanogenesis pathway）相关,其中黑素原生成通路显示和WNT信号通路、KIT信号通路及EDNRB信号通路相互关联。黑素原生成通路显示藏绵羊不同毛色可能最终通过MITF基因调控下游的TYR基因（酪氨酸酶基因）影响黑色素的合成。同时,筛选7个差异表达miRNA进行定量验证,结果表明,5个miRNA在黑色皮肤组织中上调表达,2个miRNA在白色组织中上调表达,定量结果与测序结果一致。【结论】获得了藏绵羊皮肤组织miRNA表达谱,并通过生物信息学分析得到可能和调控毛色性状相关的miRNA和信号通路。这些结果表明,藏绵羊毛色性状可能受到多种miRNA的调控并涉及多个信号通路,这有助于更进一步了解毛色转录后水平的调控过程,并对进一步的功能验证奠定基础。     

北京油鸡和来航鸡脾脏差异表达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淋巴细胞分化、增殖、激活等过程,进而影响脾脏的重量、组织结构及其免疫应答能力。     

淫羊藿苷对结肠癌增殖、侵袭抑制作用及机制研究

淫羊藿苷(Ica)梯度处理人结肠癌细胞HT29,采用MTT法、流式细胞术、Transwell小室检测Ica对HT29细胞增殖、周期、凋亡、侵袭和迁移的影响;蛋白组织学技术寻找Ica作用前后的差异蛋白, Western blot对部分蛋白进行验证。结果表明:Ica能够抑制HT29增殖;阻滞细胞由G_0/G_1期向S期的发展;诱导HT29细胞凋亡;抑制细胞侵袭和迁移; Ica干预24 h后,乙醛脱氢酶2 (ALDH2)和钙离子依赖性磷脂结合蛋白1 (CPNE1)上调, Western blot验证结果一致。这一研究提示Ica可显著抑制HT29细胞的增殖、侵袭和迁移,诱导细胞凋亡,其机制可能与上调ALDH2和CPNE1有关。     

单双羔绵羊卵巢组织差异表达microRNAs筛选

【目的】全面了解绵羊卵巢组织miRNAs表达情况,分析产单羔和产双羔绵羊卵巢miRNAs表达差异,从而为探讨miRNAs在繁殖力调控中的作用提供依据。【方法】应用多物种miRNA芯片对绵羊卵巢组织miRNA进行表达分析。首先选择经产单羔羊和双羔羊,发情后采集双侧卵巢提取总RNA,分离小片段 RNA与miRNA芯片杂交,然后进行数据分析获得绵羊卵巢组织miRNAs表达谱;利用生物信息学软件,以q-value ≤ 5%,且Fold Change ≥2 或≤0.5的标准筛选单羔羊和双羔羊卵巢差异表达miRNAs,并利用q-PCR技术验证芯片结果;分别采用microT和miRDB两种方法预测差异表达miRNAs的靶基因,然后合并两种方法结果数据取二者的交集,用在线软件对靶基因进行GO功能注释。【结果】在检测的来自所有物种的miRNAs中,有5 448个miRNAs在单羔和双羔羊卵巢组织中共同表达,22个在单羔羊卵巢中特异性表达,15个在双羔羊卵巢中特异性表达;对绵羊中已报道的103个miRNAs,比较其在两组母羊卵巢组织中的表达变化,共获得11个差异表达miRNAs,其中4个上调表达,7个下调表达;随机选择一个上调和一个下调表达的miRNAs进行q-PCR验证,定量结果与芯片结果一致,说明芯片结果准确、可信;预测获得7个miRNAs的靶基因,分别是：oar-miR-370-5p、oar-miR-376b-5p、oar-miR-381-5p、oar-miR-412-5p、oar-miR-541-3p、oar-miR-544-5p和oar-miR-1185-5p,每个miRNA获得的靶基因个数分别为：115、71、1、5、8、135和23个。GO注释结果显示,miR-376b-5p和miR-1185-5p的靶基因主要参与形成细胞内组分和细胞器,在分子功能分类中,绝大部分基因为连接分子类和催化活性分子类,在生物学过程分类中,主要参与细胞增殖、分化、凋亡过程和代谢过程;同时miR-376b-5p的靶基因IGF-1、DAZL、MTOR、MET、NEDD4和miR-376b-5p的靶基因AHR参与生殖过程的调控。【结论】成功构建了绵羊卵巢miRNAs表达谱,获得产单羔母羊和产双羔母羊卵巢组织差异表达miRNAs,这些miRNAs可能与绵羊卵泡发育和产羔数多少有关。     

靶向猪内质网应激通路的microRNAs预测与验证

【目的】预测并验证靶向猪内质网应激通路中关键基因的microRNAs（miRNAs）,为进一步研究miRNAs对猪内质网应激信号通路调控提供理论基础。【方法】首先利用伪狂犬病毒（PRV）感染猪肾上皮（PK15）细胞,高通量测序检测差异表达的miRNAs。然后通过TargetScan预测靶向内质网应激通路关键基因ATF6、IRE1、PERK、GRP78、XBP1的miRNAs。构建含有候选miRNAs作用位点的双荧光素酶报告基因重组载体,并分别与miRNA-mimics共转染幼仓鼠肾（BHK-21）细胞,通过测定荧光素酶活性来验证内质网应激通路关键基因与候选miRNAs的靶标关系。然后在PK15细胞中分别过表达候选miRNAs,利用qRT-PCR和Western blot检测候选miRNAs对内质网应激通路关键基因mRNA和蛋白表达的影响。【结果】miRNA测序结果显示,PRV感染后共引起35条miRNAs差异表达。TargetScan预测显示,靶向ATF6、IRE1、PERK、GRP78、XBP1的交集miRNAs为miR-142-5p、miR-145-5p、miR-150和miR-199a-5p,并将这些交集miRNAs作为候选miRNAs。随后成功构建psiCHECK-2-ATF6-m142-3′UTR、psiCHECK-2-ATF6-m145-3′UTR、psiCHECK-2-ATF6-m150-3′UTR、psiCHECK-2-ATF6-m199-3′UTR、psiCHECK-2-IRE1- m150-3′UTR、psiCHECK-2-IRE1-m142/145/199-3′UTR、psiCHECK-2-PERK-m145/150-3′UTR、psiCHECK- 2-XBP1-m142/ 145/150/199-3′UTR、psiCHECK-2-GRP78-m145/199-3′UTR双荧光素酶报告基因载体。双荧光素酶检测结果显示,miR-142-5p显著抑制psiCHECK-2-ATF6-m142-3′UTR荧光素酶活性。psiCHECK-2-IRE1-m142/145/199-3′UTR分别与miR-142-5p mimics、miR-145-5p mimics、miR-199a-5p mimics共转染,以及psiCHECK-2-XBP1-m142/145/ 150/199-3′UTR分别与miR-142-5p mimics、miR-199a-5p mimics共转染,过表达组的荧光素酶活性均极显著低于阴性对照组。同时miR-145-5p能够显著抑制psiCHECK-2-PERK-m145/ 150-3′UTR荧光素酶活性。这些结果表明miR-142-5p、miR-145-5p和miR-199a-5p均有可能分别靶向ATF6、IRE1、XBP1,而其中miR-142-5p可能同时靶向这三个关键基因调控内质网信号通路。通过qRT- PCR和Western blot分析发现,过表达miR-142-5p后显著抑制ATF6的mRNA和蛋白表达,表明miR- 142-5p靶向ATF6参与调控内质网应激信号通路。【结论】验证了靶向内质网应激通路关键基因ATF6的miRNA- miR-142-5p,为进一步研究miR-142-5p通过调控ATF6的表达而影响内质网应激信号通路奠定了基础。     

外泌体microRNA在猪成熟和闭锁卵泡中的表达差异及功能分析

【目的】通过分析成熟卵泡液外泌体（mature follicular fiuid Exosomes, mffEXs）和闭锁卵泡液外泌体（atretic follicular fiuid Exosomes, affEXs）miRNA的表达差异,探索卵泡液外泌体（EXs）miRNA在卵泡发育和闭锁过程中的调控作用。【方法】本研究通过抽提4—6 mm猪成熟发育和闭锁卵泡的卵泡液分离外泌体,进行粒径分析及Western Blot检测对EXs进行鉴定,接下来对特征性EXs携带的miRNA测序和功能富集分析,筛选关键信号通路和差异基因。最后,将mffEXs和affEXs作为添加剂进行颗粒细胞培养,利用Q-PCR检测技术分析关键基因的表达,验证两类卵泡液内EXs miRNA在卵泡发育中的调控功能。【结果】成功分离了mffEXs和affEXs,对比mffEXs测序结果,affEXs中有90个miRNA上调表达,220个miRNA下调表达,表明了卵泡液中的miRNA表达水平与调控卵泡发育有关;KEGG富集分析结果显示两类卵泡的差异信号通路主要集中在Ras、cAMP、P53和MAPK等信号通路,涉及调控卵母细胞发育、减数分裂以及颗粒细胞细胞周期等生物学功能。在闭锁卵泡中,上调表达的ssc-let-7a和ssc-miR-133a-3p分别潜在靶向调控细胞周期蛋白依赖性激酶（CDK1）和胰岛素生长因子（IGF1）,抑制了G1和G2/M期的运转和类固醇激素代谢,促使颗粒细胞周期运转受阻和颗粒细胞凋亡,引起卵泡闭锁的发生;下调的ssc-miR-21-5p潜在靶向肿瘤抑癌基因（P53）,抑制细胞周期运转,促使颗粒细胞凋亡。在体外培养的颗粒细胞中分别添加mffEXs和affEXs,Q-PCR结果显示CDK1在mffEXs中显著上调表达,而P53显著下调表达,表明了测序分析结果的可靠性。这些结果均显示了affEXs中miRNA表达水平的变化促使颗粒细胞凋亡和细胞周期阻滞,引起卵泡闭锁。【结论】猪affEXs携带miRNA增加了对CDK1、IGF1和P53的表达调控,抑制颗粒细胞细胞周期运转和类固醇激素代谢等信号通路,引起颗粒细胞凋亡,导致卵泡闭锁。     

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

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

高温条件下miRNA-24对奶牛乳腺上皮细胞增殖与凋亡的影响

【目的】MicroRNA(miRNA)是一类长度约21-23个核苷酸的非编码小RNA分子,本研究旨在分析高温条件下miRNA-24对乳腺上皮细胞生长、凋亡的作用及其初步机制。【方法】体外培养奶牛乳腺上皮细胞,高温(41℃)处理所培养的细胞;应用miRNA基因沉默技术,抑制高温处理过的乳腺上皮细胞中miRNA-24的表达;采用细胞计数法分析细胞增殖情况;Hoechst33342/PI荧光染色、流式细胞术检测miRNA-24对乳腺上皮细胞凋亡的作用,Western blot检测凋亡相关因子表达变化。【结果】抑制miRNA-24的表达,高温条件下的乳腺上皮细胞增殖能力显著增强(P0.05),凋亡细胞数显著减少(P0.05),促凋亡因子caspases-8和caspases-3表达量下调。【结论】内源性miRNA-24对乳腺组织发育具有重要的调控作用,抑制miRNA-24的表达可以缓解高温诱导下奶牛乳腺上皮细胞凋亡发生率、促进细胞的生长发育。     

基于基因组的绒毛状烟草和林烟草microRNA及其靶基因分析

【目的】填补绒毛状烟草（Nicotiana tomentosiformis）和林烟草（Nicotiana sylvestris）在miRNA相关领域的研究空白,揭示普通烟草（Nicotiana tabacum）的生长发育调控机理。【方法】在绒毛状烟草和林烟草全基因组中预测并分析miRNA及其靶基因,通过同源比对及miRNA前体二级结构特征进行预测：参考序列在绒毛状烟草和林烟草基因组的比对中,允许最多1-2个错配;miRNA二级结构为经典茎环结构,其中MEF最大值为-25,MEFI最小值为0.85,预测的miRNA与已知的同一家族的miRNA位于发夹结构的同一条臂上;去除E值小于等于1e-6的编码蛋白的序列。【结果】在绒毛状烟草中得到39个家族的162条miRNA,包括14对正/反义miRNA和5个基因簇。在林烟草中得到40个家族的169条miRNA,包括13对正/反义miRNA和3个基因簇。2个野生烟草在保守度高的miRNA家族中,其成员分布相似,且成员数相近。在保守度相对较低的家族中,2个野生烟草其成员分布差异较为明显,其中,miR5021、miR5203等9个家族仅在绒毛状烟草中有成员,miR1446、miR1509等10个家族仅在林烟草中有成员。2种野生烟草的正义miRNA与反义miRNA都有着1-4个碱基差异,这些差异位点在不同的家族中呈现出偏好性,而且在2种野生烟草中偏好性相似：miR164家族的9、12、13个碱基处,miR172家族的1、21个碱基处,miR396家族的2、17个碱基处,miR399家族的15、20个碱基处。2种野生烟草的基因簇主要是由miR156、miR169家族组成,其前体的间距小于350 nt,同时在绒毛状烟草中首次发现miR6019/miR6020基因簇。以普通烟草的unigene数据库作为靶基因集进行预测与分析,在绒毛状烟草122条miRNA中得到749个靶基因,去掉重复基因得到非冗余靶基因206条,其中89条（43%）得到GO功能注释;在林烟草中117条miRNA得到650个靶基因,去掉重复基因得到非冗余靶基因169条,其中78条（46%）得到GO功能注释。在分子功能方面,大多数靶基因具有结合等活性。在生物学过程中,靶基因主要参与了发育过程、生殖过程、多细胞器官发育过程、胁迫应答过程等。【结论】控制发育和多细胞器官发育过程的靶基因数方面以林烟草居多,而胁迫应答的靶基因数以绒毛状烟草较多。     

水稻根系盐胁迫响应miRNA和tRF的鉴定

【目的】水稻（Oryza sativa L.）是中国最重要的粮食作物,也是盐胁迫敏感作物。研究水稻盐胁迫响应基因表达,发掘水稻耐盐基因对水稻抵御盐胁迫的分子机制及解析其调控网络,为培育耐盐水稻品种奠定基础。【方法】以水稻品种日本晴种子为试验材料,在1/2MS培养基培养条件下对其进行盐处理（150 mmol·L^-1 NaCl）,对盐处理和非盐处理21 d的根系进行小分子量RNA组学测序。通过生物信息学分析,以log₂ Fold Change（log₂FC）>1或＜-1为筛选条件,寻找在盐胁迫和非盐胁迫条件下差异表达的miRNA和tRF（tRNA-derived RNA fragments）,分析其靶基因,并通过实时荧光定量PCR对测序结果及靶基因进行验证。【结果】以至少一组数据RPM（Reads Per Million reads）>500和log₂FC>1或＜-1为筛选条件,共得到31个差异表达miRNA,其中8个为盐胁迫诱导miRNA,23个为盐胁迫抑制miRNA,这些差异表达miRNA属于12个miRNA家族,其中8个miRNA家族在拟南芥、玉米和小麦等物种中也被报道为盐胁迫信号响应miRNA,包括盐胁迫抑制的osa-miR397、osa-miR396、osa-miR156、osa-miR167、osa-miR1432和盐胁迫诱导的osa-miR159、osa-miR168、osa-miR164。其余4个miRNA家族osa-miR1882、osa-miR1876、osa-miR1423和osa-miR5077尚未见与盐胁迫相关的报道。通过靶基因预测,得到这31个差异表达miRNA的靶基因共162个。此外,盐处理后,水稻根系产生的34—38nt tRF数量显著多于非盐处理材料,说明tRF的产生并非随机,而是通过某种特定机制响应高盐信号,诱发tRNA特异性加工而产生。以RPM>50和log₂FC>1或＜-1为筛选条件,检测到盐胁迫诱导的5'端tRF 3个,3'端tRF 3个,这些tRF由6个tRNA加工产生。推测这些差异tRF是潜在的水稻盐胁迫响应tRF。【结论】共检测到12种水稻根系中盐胁迫响应miRNA,其靶基因多为转录因子编码基因,推测其通过对其靶基因转录因子的转录后调控参与了盐胁迫响应的表达调控。其中8个水稻盐胁迫响应miRNA家族是不同物种间保守的通用盐胁迫响应miRNA。另外,从转录组水平挖掘出水稻盐胁迫响应tRF,并鉴定了6个盐胁迫诱导表达的tRF。     

山羊卵巢富集chi-miR-100f-5p的鉴定与表达分析

【目的】探讨山羊卵巢中富集的miR-100及预测的靶基因与山羊卵巢功能及卵泡发育的潜在关系,为山羊繁殖机理研究提供数据。【方法】 以山羊卵巢组织为材料,通过miRNA的分离、克隆及测序得到山羊卵巢中表达的miR-100,采用RT-qPCR进行表达检测,用生物信息学方法预测miR-100的靶基因及其功能。【结果】① 山羊卵巢的测序结果中,一段成熟体为21nt 的序列与果蝇dan-miR-100高度同源,将其命名为chi-miR-100f-5p,其前体位于基因组负（-）链上,长度为59 bp,定位于山羊CM000899.1:33 997 894—33 997 914序列;② chi-miR-100f-5p 在山羊卵巢、垂体及心脏等组织中均有表达,说明此miRNA不是卵巢组织所特有;心脏和垂体组织中均有较高表达,其次是卵巢,而在肝脏、脾脏、肺及肾脏中表达显著低于卵巢（P ＜0.01）;在单产羔内蒙古绒山羊卵巢中chi-miR-100f-5p的表达极显著高于高繁殖力大足黑山羊卵巢（P＜0.01）;③ chi-miR-100f-5p主要靶作用于整合蛋白β1亚基基因及垂体特异性转录因子1基因,可能参与整合蛋白信号通路,促进垂体内分泌细胞的生长与分化,调节下丘脑-垂体-性腺轴激素分泌。【结论】山羊卵巢中新发现的chi-miR-100f-5p 可能是与繁殖调控有关的重要候选miRNAs 之一。