湖羊卵巢不同发育阶段的miRNA鉴定与分析

旨在分析不同发育阶段绵羊卵巢组织中miRNA的表达水平,了解miRNA表达及其调控机理。本研究采用转录组测序技术鉴定1和8月龄湖羊卵巢组织中的miRNA,运用生物信息学方法对差异表达的miRNA进行靶基因预测,并对靶基因进行GO和KEGG功能注释,筛选参与卵巢发育的miRNAs。结果显示,6个样品中检测到的已知miRNA个数分别为2 186、2 091、2 217、2 136、2 176、2 088,新miRNA的个数分别为613、567、668、640、662、450。筛选到701个差异表达novel miRNAs,38个差异表达已知miRNAs。对差异表达miRNAs的靶基因进行GO和KEGG分析,发现主要富集在血管内皮生长因子通路和卵巢类固醇生成等过程。对3个差异miRNAs进行qRT-PCR验证,与测序结果一致。oar-let-7b、oar-miR-29a、oar-miR-134-3p、oar-miR-541-3p等的靶基因富集在TGF-beta通路、卵巢卵泡生长等与卵巢发育相关的生物过程。oar-miR-148a、oar-miR-136的靶基因调控颗粒细胞增殖和类固醇生成,可能参与调控绵羊卵巢发育。     

猪GV/MⅡ期卵母细胞miRNAs表达谱及Npm2基因相关miRNAs筛选

旨在探索体外成熟前后猪卵母细胞miRNAs表达谱，并筛选参与调控Npm2基因表达的miRNAs。本试验回收屠宰场猪卵巢，收集GV期卵母细胞，经体外成熟培养后获得MⅡ期卵母细胞，分别提取约130个GV期和MⅡ期卵母细胞总RNA进行Illumina HiSeqTM 2500测序，获取miRNA测序数据，每个处理重复3次，筛选差异表达microRNAs并进行靶基因GO和KEGG聚类分析。结果表明，本研究成功构建了GV期和MⅡ期猪卵母细胞miRNAs表达谱，GV期和MⅡ期卵母细胞差异表达miRNAs有95个，具有相似表达模式的miRNA聚为一类，对95个差异表达miRNAs进行靶基因预测，共得到3 967个靶基因，经GO和KEGG富集分析表明，这些靶基因被富集于5 194个GO条目和212个信号通路，其中参与卵母细胞减数分裂成熟的信号通路有2个。在差异表达miRNAs中筛选到5个与Npm2基因相关的miRNAs。采用qRT-PCR对其中2个miRNAs进行验证，证实其表达趋势与测序结果一致。本研究筛选了GV期和MⅡ期卵母细胞差异表达miRNAs，推测差异表达的miRNAs可能通过代谢、卵母细胞减数分裂相关通路、孕激素介导的卵母细胞成熟通路等途径在卵母细胞体外成熟过程中发挥作用，并筛选出调控Npm2基因表达的miRNAs，研究结果可为进一步阐明miRNA对Npm2基因的调控及其在卵母细胞成熟过程中的作用提供依据。     

乏情和发情初产母猪下丘脑-垂体轴中miRNA表达谱比较分析

【目的】从神经内分泌系统揭示微小RNA(microRNA,miRNA)调控初产母猪情期活动的遗传机制。【方法】利用RNA-Seq技术对乏情和发情初产母猪下丘脑-垂体轴中的small RNA进行测序,获得了母猪下丘脑和垂体miRNA表达谱,并对筛选的差异表达miRNA进行GO和KEGG等功能分析。【结果】从所构建的4个文库中共鉴定出1 096个miRNAs,其中已知miRNAs 364个,新预测miRNAs 732个。以发情母猪为对照,在乏情母猪的下丘脑有16个miRNAs存在差异表达,其中6个上调,10个下调;在乏情母猪垂体有9个miRNAs存在差异表达,其中2个上调,7个下调。下丘脑中16个差异表达miRNAs共预测到5 212个靶基因,垂体中9个差异表达miRNAs共预测到6 897个靶基因。下丘脑中靶基因主要在细胞投射组装、质膜结合的细胞投射组装、神经元投射发育的负调控、蛋白丝氨酸/苏氨酸激酶活性等GO条目及Ras信号通路、神经营养因子信号通路、ErbB信号通路、mTOR信号通路等KEGG信号通路中显著富集;垂体中靶基因主要在mRNA加工调控、Ras GTP酶结合、泛素样蛋白连接...     

营养诱导舍饲哈萨克羊非繁殖季节卵巢组织中miRNAs差异表达分析

旨在利用Solexa测序技术分析、挖掘营养诱导非繁殖季节绵羊发情期和乏情期卵巢组织micro RNAs(miRNA),为进一步研究特定miRNA参与非繁殖季节绵羊发情的分子调控机制提供理论基础。以不同营养水平饲喂的哈萨克绵羊作为研究对象,在非繁殖季节前后分别屠宰已鉴定的3只发情(EN)和乏情(AN)母羊,分别采集卵巢(O)组织,构建两个小片段RNA文库(OEN和OAN),进行高通量测序和生物信息学分析,实时荧光定量PCR(qRT-PCR)验证差异表达的miRNAs在哈萨克羊卵巢中表达水平。结果显示,两文库分别获得8 665 978(发情期)和9 071 102条(乏情期)clean reads,且绝大多数小RNA序列长度为21~23nt。两文库相比较,发现有9个已知和104个未知miRNAs表达差异显著(P0.01)。利用qRT-PCR随机选择的6个显著差异表达的miRNAs进行验证,其表达水平和测序分析结果一致。KEGG通路分析表明,这些差异miRNAs的靶基因主要参与了甘氨酸、丝氨酸和苏氨酸代谢,氨基糖和核苷酸糖代谢,精氨酸和鸟氨酸代谢及甘油磷脂代谢等营养调控通路。结合靶基因预测及通路富集分析,推测差异表达的miRNAs是通过营养代谢通路对绵羊的非季节性发情起重要作用。     

BMP2对猪前体脂肪细胞差异表达miRNA1的影响

采用高通量测序技术构建经骨形态发生蛋白-2(BMP2)处理前后的猪前体脂肪细胞差异miRNA表达谱,以明确BMP2影响的功能性miRNA。结果显示:BMP2刺激猪前体脂肪细胞前后共有55个miRNAs存在较大的表达差异,包括30个上调表达的miRNAs和25个下调表达的miRNAs;实时荧光定量PCR(qRT-PCR)检测表明,高通量测序数据结果正确可靠;经miRnada和RNAhybrid这2个软件对差异表达的miRNAs进行靶基因预测及靶基因进行GO和KEGG Pathway分析,发现靶基因主要富集在Notch信号通路、胰岛素信号通路、甘油磷酯代谢、MAPK信号通路、半乳糖代谢、类固醇激素的生物合成等通路。结果表明:BMP2作为前体脂肪细胞分化重要调控因子,可能通过介导miRNA及其调控的某些关键基因的表达,从而影响前体脂肪细胞外基质与受体的结合、胰岛素利用、脂类和糖的合成与代谢等生物学过程,最终作用于前体脂肪细胞的分化和脂质沉积。     

感染产肠毒素大肠杆菌的猪小肠上皮细胞microRNA表达谱分析

试验旨在探索产肠毒素大肠杆菌（enterotoxigenic Escherichia coli，ETEC）感染猪小肠上皮细胞（IPEC-J2）诱导的microRNA （miRNA）表达谱变化，为解析宿主miRNA在ETEC感染过程中的调控作用提供理论基础。利用Illumina 6000 Novoseq SE50测序平台分别对ETEC感染前后的IPEC-J2进行高通量测序，用Bowtie与参考基因组比对，用DESeq R Package进行miRNA差异性分析。通过miRanda和RNAhybid共同预测差异表达miRNA的靶基因，对差异表达miRNA靶基因进行GO功能和KEGG通路分析。随机选取5个miRNAs，对测序结果进行实时荧光定量PCR验证。结果显示，IPEC-J2在感染前后的sRNA文库经过滤分别得到12 889 260和11 203 056条clean reads。感染前后文库中，miRNA所占比例最高，分别为73.16%和54.10%；分别有97.98%和69.83%长度为18~40 nt的sRNA可比对到参考基因组，表明测序质控良好。长度在22~24 nt的序列大部分首位碱基偏向U，2~8位点出现频率最高的碱基分别为AGCUUAU。共发现311个已知miRNAs，128个新miRNAs。在2个文库中，长度为23 nt的miRNA序列占比最高，分别为41.42%和23.56%。感染后共筛选到140个差异表达miRNAs，其中74个表达上调，66个表达下调。GO分析表明，miRNA靶基因显著富集于代谢过程、正向调节代谢过程、细胞成分或生物合成、免疫系统、细胞内部分和细胞器等功能。KEGG分析表明，差异表达miRNA靶基因显著富集于赖氨酸降解、生产IgA的肠道免疫网络、NF-κB信号通路和T细胞受体信号通路等。实时荧光定量PCR验证结果表明，随机选取的5个miRNAs表达趋势与测序结果一致，表明测序准确可靠。综上所述，IPEC-J2的miRNAs参与了ETEC感染过程，为进一步揭示调控ETEC感染的关键miRNA及其作用机制提供科学依据。     

鸡胚胎期肌纤维形成相关miRNAs的筛选与分析

研究旨在挖掘鸡胚胎期肌纤维形成相关miRNAs,为进一步探究miRNA在肌肉发育过程中的调控机制提供理论基础。研究采集玫瑰冠鸡、科宝鸡胚胎期第8天（E8）和第14天（E14）腿肌组织进行HE染色和small RNA测序,通过生物信息学分析差异表达miRNAs,随机选择14个差异表达的miRNAs进行荧光定量PCR(qRT-PCR)验证结果的准确性。结果显示,玫瑰冠鸡、科宝鸡E8的肌纤维直径和横截面积极显著小于E14(P<0.01),E8的肌纤维密度极显著高于E14(P<0.01)。通过生物信息学分析,在12个测序文库中共鉴定出2 455个已知miRNAs和329个新的miRNAs,其中在E8和E14(C8T vs C14T＿R8T vs R14T)中筛选出499个差异表达miRNAs。GO富集结果显示,差异miRNA靶基因在细胞迁移、代谢过程、细胞增殖、细胞发育和小分子结合等GO条目中显著富集。KEGG富集分析发现,大量的靶基因富集到与细胞周转和代谢相关的信号通路中,同时富集到TGF-β、Notch、Hedgehog、Wnt、FoxO等与肌肉发育相关的信号通路。miRNA-...     

鸡白痢沙门菌感染雏鸡的骨髓miRNA表达谱分析

旨在探究鸡白痢沙门菌(Salmonella enterica serovar Pullorum,S.Pullorum)感染雏鸡后骨髓miRNA的差异表达特征,为深入了解鸡白痢沙门菌的致病机制提供理论基础。将7日龄SPF雏鸡随机分为两组,分别口服鸡白痢沙门菌和PBS,于感染24 h后采集骨髓进行miRNA测序,筛选差异倍数≥2且P值≤0.05的差异表达miRNAs进行靶基因预测以及GO、KEGG富集分析,随机选取6个miRNAs进行qRT-PCR验证,构建与免疫过程相关的miRNA-mRNA靶点网络。通过miRNA测序,共获得20个已知的差异表达miRNAs,其中11个上调,9个下调。qRT-PCR结果表明,miRNA变化趋势与测序结果一致。GO分析结果表明,差异表达基因主要富集在膜运输、信号转导、免疫系统、碳水化合物代谢、糖类的生物合成和代谢等,KEGG的信号通路主要富集在Notch信号通路、Hedgehog信号通路、PPAR信号通路、AMPK信号通路、Hippo信号通路等,miRNA-mRNA网络互作发现,gga-miR-1466和gga-miR-6643-5p可能是参与免疫过程相关的关键miRNA。本研究解析了鸡白痢沙门菌感染雏鸡的骨髓miRNA表达谱特征,为了解鸡白痢沙门菌和鸡之间相互作用的复杂分子致病机制提供了新的见解,为防控鸡白痢沙门菌感染提供了新策略。     

外源皮质酮对雏鸡法氏囊microRNA表达的影响

MicroRNA(miRNA)是在真核生物中发现的一类内源性的具有调控功能的非编码小RNA,在鸡的许多生物学功能中发挥重要作用,但在鸡的应激性免疫抑制的分子调控机制中尚不完全清楚。本试验通过在日粮中添加外源皮质酮激素,构建雏鸡应激性免疫抑制模型,高通量Solexa测序技术建立雏鸡法氏囊差异miRNA表达谱,共鉴定出77个差异表达的miRNAs,34个上调,43个下调,生物学分析预测差异表达miRNA的靶基因有282个,对这些预测的靶基因进行GO和KEGG通路富集分析,结果显示这些通路主要富集到心肌收缩、精氨酸和脯氨酸代谢、RNA聚合酶、肾素-血管紧张素系统。试验结果为揭示家禽应激性免疫抑制转录后分子调控机制提供了依据。     

猪卵泡发育过程中差异表达miRNA的鉴定

本研究旨在探讨miRNAs对猪卵泡发育过程的影响。以卵泡期第4天杜洛克猪中等M1卵泡(直径3.0~4.9 mm)与M2卵泡(直径5.0~6.9 mm)为试验材料,通过miRNA芯片杂交检测筛选到在杜洛克母猪中等M1卵泡与M2卵泡之间的差异表达miRNAs,随机抽选若干个检测到的miRNA用Q-PCR的方法对芯片结果进行验证,并预测差异表达miRNA的靶基因,进行GO及KEGG分析。结果表明:芯片检测出347个猪miRNA的表达,筛选得到表达量高且差异大的34个miRNA;随机筛选出6个差异表达的miRNA经Q-PCR验证,与基因芯片结果一致;分析得到的靶基因涉及最多的为Jak-STAT、细胞因子及其受体的相互作用信号通路,且两个信号通路均与卵泡发育相关。结果提示:检测筛选得到的差异表达miRNA可能对猪的卵泡发育过程具有调控作用。     

多浪羊与小尾寒羊皮下脂肪组织转录组分析

旨在探索多浪羊（D组）与小尾寒羊（X组）皮下脂肪组织中的特异性表达基因,并研究其潜在的作用,为理解绵羊脂肪组织发育规律以及对脂代谢相关疾病的预防和治疗研究提供依据。本研究选取脂肪沉积能力存在差异、健康无病、体况良好、种内个体体重相近（约50 kg）的雌性成年多浪羊和小尾寒羊为试验材料,分为D组（试验组）和X组（对照组）,每组3个重复,采集位于背最长肌的皮下脂肪组织,应用RNA-Seq技术和生物信息学方法进行转录组测序并对结果进行分析。以|Fold change|≥2,P adjust≤0.05为标准筛选差异表达基因,通过对差异表达基因进行GO功能注释和KEGG通路富集分析,得到与脂肪沉积和脂代谢有关的差异基因。为了验证测序数据的可靠性,本研究随机选取了6个差异表达基因进行qRT-PCR验证。结果显示,在6个样本中共检测到38 672个已知的mRNAs,新的mRNAs为1 606个,在两组中共有839个差异表达基因,其中有320个差异基因在多浪羊组中上调表达,有519个差异基因在多浪羊组中下调表达。通过GO功能注释分析发现,差异表达基因主要参与脂质分解代谢过程、脂质生物合成过程、脂质分解代谢负调控过程、MAPK级联反应调控、对甘油三酯的反应等生物学过程。KEGG通路富集结果显示,差异表达基因显著富集到了PI3K-Akt、MAPK、胰岛素以及PPAR等信号通路中。qRT-PCR结果与测序结果一致,表明测序结果可靠。通过对多浪羊和小尾寒羊皮下脂肪组织进行转录组测序以及生物信息学分析,筛选到与脂肪沉积和脂代谢相关的差异表达基因,这些基因主要参与脂质生物合成、脂质代谢等过程,其中COL1A1、AKT2、SCD、LPL、PCK1与PPP2R5A可能在多浪羊与小尾寒羊的皮下脂肪组织的沉积与代谢中发挥重要作用。     

藏鸡不同发育阶段腿部肌肉组织转录组及microRNA联合分析

旨在从转录组和miRNA角度探讨藏鸡肌肉发育的调控机制,了解藏鸡肌肉发育的特殊性。本研究对120和150日龄藏鸡的肌肉组织进行转录组和small RNA测序,筛选出两个日龄阶段差异表达的基因和miRNA,并利用qRT-PCR技术对测序结果进行验证。结果,共筛选出1 691个差异表达基因,其中上调基因330个,下调基因1 361个。差异表达miRNA共有22个,其中9个上调,13个下调。随机选择的5个基因和miRNAs的qPCR验证结果表明表达趋势与测序结果一致。GO富集分析显示,在富集前10的条目中,与免疫相关的条目占较大比例。KEGG分析结果显示,在19个显著富集通路中,与免疫相关的通路占较大比例,且83个miRNA的靶基因出现在显著富集通路上。转录组和small RNA测序数据联合分析表明,343个miRNA-mRNA对为负相关调控模式。本研究从多组学层面为进一步理解藏鸡的肌肉发育提供了理论基础。     

不同繁殖季节绵羊卵巢组织中piRNAs序列特征及表达差异分析

为了分析绵羊卵巢组织中piRNAs的序列特征及不同繁殖季节绵羊卵巢组织中piRNAs的表达差异,本试验采用Solexa测序技术对休情季节（休情期）和发情季节（卵泡期和黄体期）滩羊卵巢组织进行高通量测序,采用生物信息学方法筛选绵羊卵巢中的piRNAs,对piRNAs的长度分布、首位碱基偏好、基因组来源、染色体分布、链特异性、ping-pong循环特征、piRNAs比对重复序列类型进行分析,并对发情季节与休情季节绵羊卵巢之间差异表达的piRNAs进行靶基因预测及其靶基因的GO功能和KEGG Pathway富集分析。结果显示,绵羊卵巢中不同长度的piRNAs首位碱基偏好性不一致;3个时期（休情期、黄体期和卵泡期）的piRNAs均主要比对到基因内含子、编码区及lncRNA区域,而比对到重复序列的piRNAs相对较少;绵羊卵巢piRNAs在染色体上呈不均匀分布;绵羊piRNAs簇存在很强的链特异性,但首位碱基偏好性不明显;3个时期的piRNAs均无明显的ping-pong循环特征。KEGG富集分析发现,差异piRNA靶基因主要富集在RNA运输通路（RNA transport pathway）、嘌呤代谢通路（purine metabolism pathway）、黏着斑通路（focal adhesion）、PI3K-Akt信号通路（PI3K-Akt signaling pathway）、Hippo信号通路（Hippo signaling pathway）等与绵羊繁殖相关的通路。本研究揭示了不同繁殖季节绵羊卵巢组织中的piRNAs序列特征,暗示绵羊卵巢组织中的piRNAs可能通过BAD等基因和黏着斑通路等调控绵羊季节性繁殖,为深入解析绵羊季节性发情分子机制提供一定参考。     

Integrated analysis of miRNA–mRNA interaction in ovaries of Turpan Black Sheep during follicular and luteal phases

To investigate the regulatory mechanism of the follicular–luteal phase transition in Turpan black sheep (Ovis aries), the genome-wide expression patterns of microRNAs (miRNAs) and genes were investigated in ovaries of six sheep (3 years and single lamb with 3 consecutive births) during follicular and luteal phases of the oestrous cycle. Bioinformatic analysis was used to screen potential miRNAs and genes related to Turpan black sheep ovarian function. RT-qPCR was used to validate the sequencing results. In total, we identified 139 known and 71 novel miRNAs in the two phases with miRNA-seq, and a total of 19 miRNAs were significantly differentially expressed, of which 7 were up-regulated and 12 were down-regulated in the follicular phase compared with luteal phase. A total of 150 genes were significantly differentially expressed, including 63 up-regulated and 87 down-regulated in the follicular phase compared with the luteal phase by RNA-seq data analysis. Those DEGs were significantly enriched in 103 GO terms and several KEGG pathways, including metabolic pathway, ovarian steroidogenesis, steroid hormone biosynthesis and oestrogen signalling pathway. In addition, we created a miRNA–mRNA regulatory network to further elucidate the mechanism of follicular–luteal transition. Finally, we identified key miRNAs and genes including miR-143, miR-99a, miR-150, miR-27a, miR-125b, STAR, STAT1, which might play crucial roles in reproductive hormone biosynthesis and follicular development. The miRNA–mRNA interactive network clearly illustrates molecular basis involving in follicular–luteal transition.     

野生盘羊×巴什拜羊杂交二代尾部脂肪差异表达lncRNA的鉴定与分析

【目的】 分析不同脂尾型绵羊尾部脂肪组织中的长链非编码RNA(long non-coding RNA,lncRNA)的表达谱,探究lncRNA与绵羊尾部脂肪沉积机制的关系,为揭示绵羊尾脂沉积机制提供理论依据。【方法】 选择新疆地方绵羊巴什拜羊(肥尾型)和野生盘羊×巴什拜羊杂交二代(小尾型)作为研究对象,采集2个群体尾部脂肪组织,利用转录组测序技术筛选差异表达lncRNA并预测靶基因,对其进行GO功能和KEGG通路富集分析;利用实时荧光定量PCR技术验证转录组测序的表达结果。【结果】 通过差异表达分析共筛选出728个差异表达lncRNAs,其中270个表达下调,458个表达上调;通过差异表达lncRNA靶基因的GO功能和KEGG通路富集分析,筛选到634个靶基因参与疾病、免疫、蛋白修饰、细胞代谢等相关功能分类,共涉及76条信号通路,部分lncRNA介导的靶基因SCD、GPAM、THRSP、FASN等与绵羊尾部脂肪沉积相关。实时荧光定量PCR与转录组测序结果趋势相同。【结论】 lncRNA在绵羊脂尾进化中对尾部脂肪沉积具有重要的调控作用,为从lncRNA的角度分析绵羊脂肪沉积提供了理论依据。     

我国蒙古羊系统主要地方绵羊品种遗传分化的研究

应用中心产区典型群简单随机抽样的方法,检测控制湖羊、同羊、滩羊、小尾寒羊和洼地羊的血液酶和其他蛋白质变异的5个结构基因座位(Tf、X-p、Ary-Es、Hb-β、Ke)上的遗传变异,引用国内外学者以相同实验方法获得的其他分布于我国周边国家和地区的11个绵羊品种作为分析背景,对16个群体进行遗传距离分析及亲缘关系聚类。结果表明:中亚以东南绵羊群体可以划分为蒙古羊、藏羊和南亚-东南亚羊三大集团,湖羊、同羊、滩羊、小尾寒羊和洼地羊均属于"蒙古羊"系统。     

发情周期不同阶段绵羊卵巢lncRNA的鉴定和功能分析

旨在分析不同繁殖周期绵羊卵巢组织中长链非编码RNA (long non-coding RNA,lncRNA)的表达谱,了解lncRNA表达及其调控机理,为绵羊繁育研究提供理论依据。本研究以湖羊为研究对象,选取年龄在1.5～2.5岁的黄体期和卵泡期母羊各3只,通过高通量测序筛选出卵巢组织中的lncRNA,运用生物信息学分析对差异表达的lncRNA进行靶基因预测,通过GO和KEGG富集分析找出与绵羊繁殖相关的通路。结果显示,本研究共获得1 379个差异表达的lncRNAs,其中1 158个表达上调,221个表达下调。GO和KEGG富集分析表明,差异表达的lncRNAs及其靶基因主要参与卵泡发育、排卵周期过程、钙离子信号通路、卵母细胞减数分裂、催产素信号通路、MAPK信号通路、甲状腺激素合成通路、雌激素信号通路。关键的lncRNAs可能通过调控参与这些信号通路和生物学过程的相关基因来调控生殖。其中,LNC_011239、LNC_012847、LNC_003902、LNC_003906、LNC_003907等靶向的MAPK1、ADCY1、ADCY5、PPP3CA和CDC23可能发挥关键的调控作用。qRT-PCR验证证明,随机选取的5个差异lncRNAs定量结果与测序结果基本一致。本研究利用RNA-Seq技术筛选出黄体期和卵泡期卵巢组织中的lncRNAs,并进行差异分析,为揭示绵羊繁殖能力的分子机制提供依据。     

胚胎移植前供体牛与受体牛血浆外泌体miRNA差异表达分析

【目的】 探索胚胎移植前供体牛与受体牛血浆外泌体miRNA的表达差异, 以明确血浆外泌体miRNA在牛早期妊娠中的作用及其调控机制。【方法】 以3~6岁、体重480~600 kg的夏南牛作为研究对象, 选取10头供体牛进行同期发情、超数排卵和人工授精, 23头受体牛只做同期发情处理。在人工授精后第7天, 冲洗供体牛子宫以获取囊胚, 选取3头囊胚数相近的供体牛及3头与供体牛体重和年龄均相近的受体牛, 颈静脉采血, 进行血浆外泌体的分离与鉴定; 然后提取血浆外泌体miRNA, 并检测其表达量; 采用R语言中的DESeq差异算法计算P值, 并筛选出P<0.05的miRNA, 对差异表达的miRNA进行靶基因预测、GO功能富集分析和KEGG信号通路分析。【结果】 6个样本的囊泡粒径均在135 nm左右, 符合外泌体的特征。与供体牛相比, 受体牛中有9个miRNAs表达显著上调(P<0.05), 13个miRNAs显著下调(P<0.05);22个差异表达的miRNAs中, 有15个miRNAs预测出无重复的靶基因2 990个。GO功能富集分析和KEGG信号通路分析的结果表明, 这些靶基因主要富集在与生物黏附(biological adhesion)、定位(localization)、细胞连接(cell junction)功能有关的通路上, 显著富集的信号通路与黏着斑(focal adhesion)、黏着连接(adherens junction)有关, 提示血浆外泌体miRNA可能参与调控胚胎着床。【结论】 研究结果可为筛选和探究影响胚胎着床的血浆外泌体miRNA提供参考, 并为进一步阐明血浆外泌体miRNA在母牛早期妊娠调控中的作用提供依据。