禽致病性大肠杆菌双组分系统phoP/Q对鞭毛Ⅲ型分泌系统基因的调控预测

本研究旨在探究禽致病性大肠杆菌(APEC)中双组分系统phoP/Q与鞭毛Ⅲ型分泌系统(鞭毛T3SS)蛋白之间的联系。运用软件进行功能分析,预测phoP/Q与鞭毛T3SS的关系,利用基因芯片技术对APEC野生株和ΔphoP/Q缺失突变株进行转录组的差异性比较,筛选鞭毛T3SS基座中变化较大基因,并用荧光定量PCR检测基因表达变化量。结果显示:经STRING预测发现phoP/Q可能通过调控鞭毛T3SS基因座附近的双组分系统来间接调控鞭毛T3SS。基因芯片筛选与荧光定量PCR检测发现phoP/Q缺失能使鞭毛T3SS核心组件的大部分基因表达差异显著,同时也影响了双组分中motA、motB与cheY的表达。本研究对phoP/Q双组分系统与鞭毛T3SS的关系进行了预测与探索,为进一步研究揭示细菌的活动与致病机制提供了新的切入点与研究方向。     

phoP/Q双组分系统对禽致病性大肠杆菌的毒力调控作用

phoP/Q作为一种外在环境感受器,参与调节细菌对外部环境的适应性。本研究拟探讨phoP/Q双组分系统对禽致病性大肠杆菌(APEC)的生长、定植、毒力等生物学特性的影响。利用Red重组系统构建APEC AE 17株的phoP/Q双组分调控系统缺失株,并构建phoP/Q双组分调控系统回复质粒,转化缺失株,构建回复株。通过生长运动性试验、黏附入侵CEF细胞试验、毒力基因转录水平检测以及半数致死量(LD50)等试验比较野生株、缺失株、回复株的生物特性。与野生株相比,phoP-phoQ缺失不改变其生长特性;运动性较野生株下降63%;黏附与入侵CEF细胞能力分别降低69.83%(P0.01)和62.17%(P0.05);LD50显示毒力减弱13.3倍;polA、tsh基因转录水平分别上调41%、54%;sodA、iss分别下调64%和98%。phoP/Q双组分系统参与对APEC致病性的调控,该系统的缺失会降低APEC的致病性。     

禽致病性大肠杆菌phoP/phoQ基因敲除株的构建及其毒力基因转录调控的研究

为研究Pho P/Q二元调控系统对禽致病性大肠杆菌(APEC)致病性的调控作用,本研究采用Red同源重组技术分别构建了APEC AE17株的pho P基因敲除株AE18和pho Q基因敲除株AE19,并通过p STV-28质粒分别构建pho P和pho Q基因回补株CAE18和CAE19。并对这些菌株的运动能力、对鸡胚成纤维细胞(CEF)的黏附和侵入能力以及其13个毒力基因的转录水平进行了检测。结果表明,pho P和pho Q基因敲除株运动性能力分别为AE17的70.92%和83.83%;黏附和侵入CEF细胞的能力比AE17均极显著降低,分别为AE17的63.11%、65.42%和59.83%、57.82%。荧光定量PCR结果显示,AE18的毒力基因(sod A、pol A和iss)和AE19的毒力基因(fim C和iss)比AE17相应的毒力基因m RNA的转录水平均呈极显著下降(p0.01)。此外,回补株CAE18和CAE19株能够部分恢复运动性能力、黏附和侵入CEF细胞的能力及相关毒力基因的转录水平。结果表明pho P和pho Q基因的敲除均能够显著改变APEC AE17株相关的生物学特性,并且能够调控毒力基因的转录水平,使其毒力显著下降。本研究为进一步了解APEC的Pho P/Q二元调控系统对其致病性影响提供实验依据。     

基于非标记定量蛋白质组学技术研究RyhB调控禽致病性大肠杆菌酸性耐受力

为从蛋白水平上探究ryhb对大肠杆菌酸性耐受力相关基因的影响,本研究采用非标记定量(Label-free)蛋白质组学技术对禽致病性大肠杆菌(APEC)野生株和△ryhb基因突变株进行蛋白质组学的差异性比较,并利用荧光定量PCR技术进行验证。本实验共检测到121个差异表达蛋白,其中44个蛋白的表达水平上调,19个蛋白的表达水平下调。为研究ryhb对表达蛋白的影响情况,并鉴定受其调控的相关基因,从上述差异蛋白中筛选出了hdea、hdeb、gada、gadb等与酸性耐受力相关的蛋白,并在基因水平进行验证。本研究表明,ryhb基因的缺失可能与上述4个耐酸基因转录水平上调相关,进而增强APEC耐酸性。本实验为初步探究ryhb调控通路,对APEC致病性影响提供了一定的实验依据。     

禽巴氏杆菌链霉素耐药株转录组测序与分析

为鉴定禽源多杀性巴氏杆菌(Pm)链霉素(Str)敏感株与耐药株(MIC=1024μg/mL)在转录组水平上的差异,本研究采用RNA-seq技术对Str敏感株与耐药株进行双向测序,筛选得到差异表达基因,通过GO、KEGG和ARDB数据库对差异表达基因进行分组和富集性分析,并利用荧光定量PCR对筛选出的差异表达基因进行验证。结果显示,共筛选获得625个差异表达基因,其中包括581个上调基因和44个下调基因。通过GO功能富集分析显示差异表达基因多数发挥细胞膜转运功能;对差异表达基因进行信号通路富集性分析显示,大多数差异表达基因分布在ABC转运蛋白与双组分系统代谢途径中;采用ARDB对差异表达基因进行耐药基因分析,显示其中的耐药基因主要为氨基糖苷转移酶类和外排系统两类。本研究结果表明,禽源Pm对Str的耐药机制可能与细胞膜上多药外排泵的过量表达和氨基糖苷转移酶类灭活抗菌药物有关。同时,也为进一步研究禽源Pm对Str的耐药机制奠定了基础。     

禽巴氏杆菌四环素耐药株转录组测序与分析

为鉴定禽源多杀性巴氏杆菌(avian Pasteurella multocida,Pm)四环素(Tetracycline,Tet)敏感株与耐药株(MIC=32μg/mL)在转录组水平上的差异,本研究采用RNA-seq技术对Tet敏感株与耐药株进行双向测序,通过GO、KEGG和ARDB数据库对表达差异基因进行注释和富集性分析。结果显示,共筛选出442个差异表达基因,其中包括103个上调基因和339个下调基因。通过GO功能富集分析显示差异表达基因多数富集于核糖体结构成分、rRNA结合、跨膜运输和离子转运等生物过程;对差异表达基因进行信号通路富集性分析,大多数差异表达基因分布在ABC转运蛋白与核糖体代谢途径中;对差异表达基因进行耐药基因注释(ARDB),发现其中的耐药基因主要为核糖体保护蛋白和外排系统两类。本研究表明禽源Pm对四环素耐药的耐药机制可能与细胞膜上ABC转运蛋白的过量表达以及核糖体保护蛋白有关。     

捻转血矛线虫阿苯达唑敏感株和耐药株比较转录组学分析

为探明捻转血矛线虫阿苯达唑敏感株和耐药株在转录组水平的差异,挖掘耐药相关基因,本研究采用Illumina HiSeq 4000对敏感虫株和耐药虫株进行转录组测序,构建cDNA文库,通过质控后,对其进行生物信息学分析。结果显示:敏感和耐药虫株相比较共获得218个显著(P≤0.05)差异表达基因,对显著差异表达基因进行GO功能富集,有121、82和102个基因分别注释到生物学过程、细胞组分和分子功能三大类的206个GO terms上。KEGG数据库分析显示,有42个显著差异表达基因参与31条KEGG信号通路,显著富集在药物代谢-细胞色素P450、药物代谢-其他酶等通路。在敏感株和耐药株中,筛选出GST、CYPs和UGT等9个耐药相关基因,这些基因在捻转血矛线虫耐药中发挥重要的作用。本研究利用RNA-Seq预测出敏感虫株和耐药虫株差异表达基因在功能分类和代谢通路等方面的生物学特征,挖掘出耐药相关基因,进一步为捻转血矛线虫药物靶点预测及耐药性早期鉴别诊断的建立提供重要基础。     

基于RNA-Seq筛选禽致病性大肠杆菌损伤雏鸡小肠差异表达基因

为筛选禽致病性大肠杆菌(APEC)感染雏鸡后其小肠损伤的差异表达基因,本研究利用APEC菌株经腿部肌肉途径接种两周龄雏鸡,对照组鸡以相同的途径注射相同体积的生理盐水,选取病变严重的肠道与对照组肠道,利用高通量RNA-Seq技术筛选出雏鸡小肠受APEC损伤后的差异表达基因,对其进行生物信息学GO功能分类和KEGG信号通路分析,并利用荧光定量PCR验证部分差异表达基因。结果显示差异变化在2倍以上的基因共有131个,其中74个上调表达,57个下调表达。GO功能分类结果显示,这些基因主要涉及到蛋白结合、免疫反应、转运活动等功能。KEGG分析表明这些差异表达基因参与PPAR、新陈代谢、细胞因子受体相互作用、Jak-STAT、RIG-I-样受体、吞噬体等信号通路。本研究为进一步研究APEC的致病作用及宿主抗病机制奠定了基础。     

RNA-Seq分析baeSR和acrB对鼠伤寒沙门菌相关毒力基因表达的影响

本研究旨在探索baeSR和acrB基因缺失后对鼠伤寒沙门菌相关毒力基因表达的影响。采用RNA-Seq技术筛选鼠伤寒沙门菌体外诱导环丙沙星耐药株(CR)、acrB基因缺失株(CRΔacrB)、baeSR和acrB基因双缺失株(CRΔbaeSRΔacrB)的差异表达基因,对筛选的基因进行GO富集和KEGG通路富集分析,进一步探究差异表达基因的主要功能及涉及的生物过程,并选出7个差异显著基因进行荧光定量PCR验证。结果表明CRΔacrB与CR相比较,共筛选到1 320个差异表达基因(fold change的绝对值2,Q-value0.005),其中上调426个,下调894个;CRΔbaeSRΔacrB与CR相比较,共筛选到1 377个差异表达基因(fold change的绝对值≥2,Q-value0.005),其中上调405个,下调972个,两个比较组筛选到的毒力相关基因均主要富集在双组分系统、鞭毛组件、沙门菌感染和上皮细胞的细菌入侵等通路中。荧光定量PCR验证结果与RNA-Seq结果基本一致。本研究对baeSR和acrB基因在鼠伤寒沙门菌毒力中的作用进行了初步分析和探索,为进一步研究鼠伤寒沙门菌致病机制奠定基础。     

捻转血矛线虫阿苯达唑耐药株给药前后比较转录组学分析

为探明捻转血矛线虫阿苯达唑耐药株给药前后在转录组水平上的差异,本研究采用Illumina Hiseq4000对给药前后的捻转血矛线虫进行转录组测序,筛选得到差异表达基因并通过GO和KEGG数据库对其进行功能注释及富集性分析,并利用荧光定量PCR验证部分差异表达基因。结果显示,共筛选获得851个显著差异表达基因,其中包括584个上调基因和267个下调基因。通过GO功能富集分析显示,有458、418和367个基因分别注释到生物学过程、分子功能和细胞组分三大类;对差异表达基因进行KEGG富集分析显示,有173个差异表达基因参与到75条KEGG通路中,显著富集在核糖体合成、细胞凋亡、过氧化物酶增殖体激活受体(PPAR)等信号通路。本试验初步筛选了阿苯达唑耐药株给药前后的差异表达基因,为深入探索捻转血矛线虫耐药性分子机制、筛选对耐药性检测的分子标记及耐药性早期鉴别诊断方法的建立提供重要数据。     

双组分系统CpxR影响禽致病性大肠杆菌的耐药性、抗血清杀菌能力及毒力

CpxR是细菌中Cpx双组分系统（two component system,TCS）的反应调控蛋白,通过调控靶基因的转录表达,在细菌细胞膜稳定及毒力方面发挥作用。本研究旨在探究TCS CpxR对禽致病性大肠杆菌（avian pathogenic Escherichia coli,APEC）基本生物学特性、抗血清杀菌能力及致病性的影响。利用Red同源重组系统及互补质粒构建cpxR基因缺失株、互补株,然后比较分析野生株、基因缺失株与互补株的生长曲线、运动性、生物被膜形成能力、药物敏感性、抗血清杀菌能力、动物致病性的差异。结果显示：cpxR基因缺失株与野生株、互补株的生长速度和运动性能无明显差异,且缺失cpxR基因不影响APEC的生物被膜形成能力。然而,缺失CpxR导致APEC对阿米卡星和卡那霉素耐药性降低。血清杀菌试验结果显示,CpxR有助于APEC的抗血清杀菌能力。动物感染试验结果显示,野生株、cpxR基因缺失株和互补株对雏鸭的半数致死量（LD₅₀）分别为7.50×10⁵、7.50×10⁶、1.33×10⁶ CFU,表明CpxR缺失显著降低APEC的毒力。综上表明,TCS CpxR在APEC耐药性、抗血清杀菌能力及毒力方面发挥作用,为阐明APEC的环境适应性、生存能力及致病机制提供参考。     

禽致病性大肠杆菌(APEC)感染对雏鸡脾脏miRNA表达谱的影响

基于miRNA测序技术研究感染禽致病性大肠杆菌(avian pathogenic Escherichia coli,APEC)雏鸡脾脏microRNA(miRNA)的差异表达,为深入了解APEC的致病机制奠定理论基础。将14日龄罗曼鸡随机分为2组,分别腿部肌肉注射APEC分离株AE17和生理盐水后采集脾脏组织,制作病理切片并进行miRNA测序,用miRanda算法预测差异表达miRNA的靶基因,并进行GO和KEGG分析。切片结果显示,与对照组相比,AE17株感染组脾脏组织有病变。通过miRNA测序总共筛选了21个差异表达miRNA的靶基因,其中8个下调和13个上调;对预测靶基因进行GO分析发现其显著富集在细胞、细胞器、膜等功能上;KEGG分析显示其参与Toll样受体信号通路和MAPK信号通路等免疫系统过程。本研究在AE17株感染组脾脏病理切片中发现病理变化,并利用miRNA测序技术成功获取了APEC感染鸡脾脏的miRNA表达谱,发现雏鸡脾脏组织miRNA表达丰富且表达量各异,为探讨APEC致病机制提供新思路。     

鼠伤寒沙门菌伴侣蛋白Hfq调控基因的筛选

为筛选鼠伤寒沙门菌伴侣蛋白Hfq所调控的基因,本实验基于本研究室前期对指数生长期和稳定期的鼠伤寒沙门菌hfq基因缺失株转录组测序结果,从中随机挑选6个在转录水平差异表达的基因,并通过荧光定量PCR技术验证转录组数据的可靠性,然后在转录组数据可靠的前提下将两个转录组中所有显著差异表达(log_2FoldChang1)的基因GO富集细菌致病机制生物学过程,Pathway富集细菌ABC转运通路和群体感应通路,在富集得到的基因中筛选指数生长期和稳定期均表现为转录水平差异表达且趋势相同的基因。荧光定量PCR检测结果显示,6个基因转录水平差异上调或下调的趋势与转录组分析结果中一致,表明转录组数据可靠,可进行后续的研究。GO与Pathway富集结果显示,在鼠伤寒沙门氏菌hfq基因缺失株指数期和稳定期中均表现为转录水平差异上调或下调且差异倍数均大于2的基因共26个。其中,与细菌致病机制生物学过程相关的13个基因均上调;与细菌ABC转运通路相关的基因11个,其中9个基因转录水平均上调2个基因转录水平均下调;与细菌群体感应相关的2个基因转录水平均上调。26个基因中仅rpoE经验证受Hfq直接调控。以上结果表明,Hfq主要参与细菌致病机制生物学过程以及细菌ABC转运通路的负调控。本研究为后续Hfq靶基因的鉴定及其作用机理的研究提供了理论基础,同时丰富了Hfq调控基因的数目。     

江苏及周边地区鸭致病性大肠杆菌血清型、毒力因子及耐药研究

本研究旨在明确江苏及周边地区鸭禽致病性大肠杆菌（avian pathogenic Escherichia coli,APEC）的血清型、毒力基因分布和耐药性之间的相关性,以期为APEC的防控提供依据。从江苏省及周边养鸭场分离了191株APEC,并对其中21株（每个养殖场选取1株）的O抗原血清型、毒力基因分布和耐药性进行检测。对21株APEC的血清型检测结果表明,O65血清型12株,占全部菌株的57.14%,O5、O28、O42、O87、O93、O138、O147血清型均为1株,其他血清型2株;毒力基因检测结果表明,5个毒力基因有较高的分布率,其中fimA基因的阳性率为100%,ECs3737、ECs3703、tsh和irp2基因的阳性率分别为90.5%、85.7%、57.1%和42.9%,含有5个毒力基因的菌株共有6株（28.57%）;药敏试验结果表明,21株APEC均存在多重耐药性,100%的分离菌株对恩诺沙星、强力霉素、万古霉素和红霉素耐药,85.71%的分离株对10种以上抗生素耐药,14.29%的菌株对21种药物都耐药;对血清型、毒力基因分布和耐药性之间的关系分析表明,含有4种以上毒力基因的菌株有13株,其中9株是O65血清型。在13株含有4种以上毒力基因的菌株中,耐15种药物以上的有9株（69.23%）,耐20种以上药物的有3株（23.08%）,表明含有4种以上毒力基因的菌株多重耐药现象严重。研究表明,江苏及周边地区鸭源大肠杆菌血清型复杂,携带多种毒力基因,耐药性严重。     

MarA和SoxS共同调节大肠杆菌K12外排泵AcrAB-TolC的表达

利用Red同源重组技术构建大肠杆菌K12调控基因缺失菌株,探讨耐药发展过程中marRAB和soxRS对大肠杆菌K12外排泵AcrAB-TolC表达的调控作用。在环丙沙星浓度递增的条件下诱导K12及调控基因缺失菌株,分别获得系列突变菌株,测定各菌株对其他药物的MIC,并检测靶位基因突变情况。利用RT-PCR方法检测诱导突变株中外排泵基因、调控基因、膜孔蛋白基因的表达水平。结果表明,K12正向调控基因(marA和soxS)缺失株在环丙沙星选择压力下仅获得环丙沙星敏感性降低的菌株,而负向调控基因(marR和soxR)在环丙沙星选择压力下可获得环丙沙星中介和耐药的菌株,且诱导突变株仅出现与耐药有关的gyrA单突变(Ser83Leu和Asp87His)。正向调控基因marA和soxS缺失后,外排泵基因表达水平变化不明显,而负向调控基因marR和soxR缺失后,marA和soxS的表达水平显著升高,引起外排泵基因表达水平显著升高,最高达到约13倍。当阻遏蛋白SoxR被敲除后,soxS的表达水平没有变化,而marA的表达水平显著升高,外排泵基因acrB表达水平也升高,说明在调控AcrAB-TolC外排泵基因表达方面,marRA与soxRS存在功能上的冗余,即当其中一个调控因子功能缺陷时,另一个调控因子可以发挥互补调控功能。因此,MarA和SoxS对三聚体系统AcrAB-TolC的正向调控作用是通过解除负向调控蛋白MarR和SoxR对其的阻遏作用实现的,MarRA和SoxRS共同调节外排泵AcrAB-TolC的表达水平。     

肠炎沙门氏菌非编码小RNA STnc640对菌毛相关基因的调控分析

STnc640是一种功能未知的沙门氏菌非编码小RNA(sRNA),为研究其对肠炎沙门氏菌(SE)菌毛相关基因的调控作用,本研究通过λ噬菌体Red同源重组系统构建了SE50336株stnc640缺失突变株(50336△stnc640),并利用表达载体p BR322构建了回补株50336△stnc640/pstnc640。利用荧光定量RT-PCR(q RT-PCR)检测STnc640在SE50336株不同生长时期的转录水平,结果显示,细菌生长的稳定后期STnc640转录水平最高,为对数生长早期表达量的16.5倍;同时检测主要菌毛亚单位sefA、safA、stbA、sthA、csgA、csgD、peg 和外膜蛋白基因ompD的表达,结果显示,相比亲本株50336,50336△stnc640中sef A、csgD和ompD基因的表达均显著下调,而50336△stnc640/pstnc640回补株中以上各基因的表达均得到恢复。本实验结果表明STnc640对SE部分菌毛基因的表达具有调控作用。     

外源NO缓解草地早熟禾镉胁迫的转录组分析

为揭示外源一氧化氮对草地早熟禾（Poa pratensis）镉（Cadmium,Cd）胁迫缓解机制,本试验采用高通量Illumina Hiseq测序技术研究了草地早熟禾根施1 000 μmol·L^-1氯化镉和叶面喷施50 μmol·L^-1 NO供体硝普钠24 h后转录组基因的差异表达。对测序数据进行分析结果显示,4个处理组共有22 730个差异表达基因,包括15 332个上调基因和7 398个下调基因。通过富集分析发现外源NO主要通过调控与信号转导、碳水化合物转运与代谢、氨基酸生物合成及苯丙烷生物合成等途径相关的基因缓解镉胁迫。此外,本研究对Cd处理与Cd+NO处理组的差异表达基因和代谢通路进行分析,筛选出了一系列由NO信号介导的草地早熟禾响应镉胁迫关键基因,并将其列为NO缓解草地早熟禾镉胁迫相关的候选基因。本研究为挖掘草地早熟禾耐镉相关调控基因和功能基因提供基础。     

MSTN基因编辑牛肌肉转录组测序及生物信息学分析

为探究肌生长抑制素（myostatin,MSTN）基因对牛肌肉发育的具体调控机制,本研究选取同一牛场健康鲁西黄牛10头,其中通过转基因技术得到的基因编辑牛MSTN^-/-和同种非转基因野生型牛各5头。分别采集两组牛腿臀肌肉样品,利用IIlumina HiSeq高通量测序技术进行转录组测序分析,通过生物信息学方法比较两组样本间的差异表达基因,并进行GO和KEGG富集分析,最后利用实时荧光定量PCR验证转录组测序数据。结果显示,基因编辑型牛和野生型牛之间共检测到18 071个基因。在log₂|FoldChange|≥ 1.48条件下,筛选出406个差异表达基因,其中347个显著上调,59个显著下调。GO功能富集分析显示,MSTN基因编辑后显著影响915个功能类别（P<0.05）,差异基因主要参与结合、生物系统调节、免疫系统等相关功能。KEGG通路富集分析结果共涉及211个通路,差异基因主要富集在细胞黏附分子、趋化因子信号通路、细胞因子互作等信号通路上,进一步从中筛选出可能参与细胞生长、肌肉发育的差异基因（CD14、KIT、CSF1R、FBP1、DUSP4、ULBP21、PRKCB、SPN、CHAD、SRC）。实时荧光定量PCR检测结果显示,所选差异基因表达水平与转录组表达水平一致,证明测序结果的可靠性。本研究结果表明,MSTN基因发挥作用后可以介导多个下游基因表达,从而影响相关信号通路及生物学过程;同时,所筛选出的差异表达基因可作为进一步研究骨骼肌调控机制的候选靶标。     

捻转血矛线虫伊维菌素敏感虫株与耐药虫株差异miRNA的转录组学分析

为探究捻转血矛线虫的伊维菌素(IVM)敏感虫株与耐药虫株显著差异的microRNA (miRNA)转录谱，发掘其耐药相关的miRNA，本实验采用Illumina Hiseq2000平台分别对捻转血矛线虫IVM敏感虫株与耐药虫株测序，构建cDNA文库，通过质控后采用生物信息学软件筛选捻转血矛线虫IVM敏感与耐药虫株的差异mi RNA。结果显示，捻转血矛线虫的IVM敏感与耐药虫株筛选出转录显著差异的miRNA共375个，其中253个转录上调，122个转录下调。利用RNAhybrid (V6.01)、Miranda (V3.3a)、TargetScan (V7.0)软件对筛选到的mi RNA靶基因进行预测，并取交集作为mi RNA靶基因预测结果，将筛选出显著差异的靶基因通过perl脚本与显著差异的miRNA进行关联分析，结果显示共关联到2 106对miRNA-mRNA为负调控关系，其中涉及到477个差异的m RNA和603个差异的miRNA。对显著差异的miRNA靶基因进行GO富集分析，结果显示其中有492个GO terms，包括生物过程(BP) 322个、细胞组分(CC) 47个、分子功能...     

猪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基因的调控及其在卵母细胞成熟过程中的作用提供依据。