不同繁殖季节绵羊卵巢组织中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等基因和黏着斑通路等调控绵羊季节性繁殖,为深入解析绵羊季节性发情分子机制提供一定参考。     

乏情和发情初产母猪卵巢microRNA差异表达分析

为了探索卵巢microRNA(miRNA)在初产母猪生殖调控中的作用,本研究利用Illumina高通量测序技术检测了乏情和发情初产母猪卵巢miRNA的表达谱,并对表达量较高且差异显著的miRNA进行生物信息学分析。结果显示,本研究构建的两个small RNA文库共鉴定出503个miRNAs,其中已知的303个,新预测的200个。在已知的miRNA中,ssc-miR-10b在两个文库中的表达量最高,其次为ssc-miR-143-3p和ssc-miR-26a;在新预测的miRNA中,chr13_2637_mature在两个文库中的表达量最高,其次为chr8_9994_mature。与发情母猪相比,共有145个miRNAs发生显著变化(read counts10,∣log2(fold-change)∣1),上调114个,下调31个。在进一步筛选的31个表达量较高且差异显著的miRNAs(read counts1 000,∣log2(fold-change)∣1)中,新预测的chr13_2585_mature上调倍数最高,且只在乏情母猪卵巢中表达。31个miRNAs共预测到7 388个靶基因,KEGG信号通路分析显示,有2 788个靶基因注释到了297个KEGG通路,前20个最富集的通路部分与生殖过程或生殖活动调控相关,表明这31个miRNAs参与了初产母猪的生殖调控。本研究结果丰富了猪miRNA数据资源,为进一步深入研究初产母猪的繁殖性能提供了理论依据。     

乏情和发情初产母猪下丘脑-垂体轴中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酶结合、泛素样蛋白连接...     

小尾寒羊间情期和发情期microRNAs差异表达分析

本研究通过构建和分析小尾寒羊间情期和发情期卵巢组织microRNA(miRNA)表达谱,筛选两个时期差异表达的microRNAs,为研究microRNA调控小尾寒羊繁殖过程提供相应的理论基础。利用活体手术法在发情期(第1天)和间情期(第13天)分别采集一侧卵巢。从卵巢中提取总RNA,利用illumina Hiseq2000测序平台获取RNA数据,对表达谱和差异表达microRNAs进行生物信息学分析。实时荧光定量PCR(qRT-PCR)验证差异表达的microRNAs在小尾寒羊卵巢中表达水平。结果,成功地构建出小尾寒羊间情期和发情期microRNA的表达谱,oar-miR-99a和oar-miR-143分别是间情期和发情期表达量最高的microRNA,并筛选出在两个时期间3个显著差异表达的microRNAs,分别为oar-miR-200a、oar-miR-200b和oar-miR-200c。利用qRT-PCR对随机选择的2个显著差异表达的microRNAs进行验证,其表达水平和RNA-Seq分析结果一致。microRNA表达谱为后续的绵羊卵巢microRNA研究提供更详尽的信息。结合靶基因预测及通路富集分析,推测差异表达的microRNAs是通过代谢和免疫途径调控卵巢周期性活动。     

乏情和发情初产母猪卵巢microRNA差异表达分析

为了探索卵巢microRNA （miRNA）在初产母猪生殖调控中的作用,本研究利用Illumina高通量测序技术检测了乏情和发情初产母猪卵巢miRNA的表达谱,并对表达量较高且差异显著的miRNA进行生物信息学分析。结果显示,本研究构建的两个small RNA文库共鉴定出503个miRNAs,其中已知的303个,新预测的200个。在已知的miRNA中,ssc-miR-10b在两个文库中的表达量最高,其次为ssc-miR-143-3p和ssc-miR-26a;在新预测的miRNA中,chr13_2637_mature在两个文库中的表达量最高,其次为chr8_9994_mature。与发情母猪相比,共有145个miRNAs发生显著变化（read counts>10,∣log₂（fold-change）∣>1）,上调114个,下调31个。在进一步筛选的31个表达量较高且差异显著的miRNAs（read counts>1 000,∣log₂（fold-change）∣>1）中,新预测的chr13_2585_mature上调倍数最高,且只在乏情母猪卵巢中表达。31个miRNAs共预测到7 388个靶基因,KEGG信号通路分析显示,有2 788个靶基因注释到了297个KEGG通路,前20个最富集的通路部分与生殖过程或生殖活动调控相关,表明这31个miRNAs参与了初产母猪的生殖调控。本研究结果丰富了猪miRNA数据资源,为进一步深入研究初产母猪的繁殖性能提供了理论依据。     

不同繁殖季节绵羊卵巢组织中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等基因和黏着斑通路等调控绵羊季节性繁殖,为深入解析绵羊季节性发情分子机制提供一定参考。     

湖羊卵巢不同发育阶段的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的靶基因调控颗粒细胞增殖和类固醇生成,可能参与调控绵羊卵巢发育。     

Smad1基因在敖汉细毛羊发情不同时期卵巢中的表达规律研究

为了探索Smad1 mRNA及其蛋白在敖汉细毛羊不同发情时期卵巢中表达量的变化规律,试验采用实时荧光定量PCR、免疫组化方法分析绵羊不同繁殖状态下(乏情期、发情间期、发情前期、发情期)Smad1基因在卵巢中mRNA表达量的变化规律及卵巢中该基因蛋白的表达,并进行图像分析。结果表明:发情期Smad1基因表达量最高,极显著高于其他三个时期(P0.01);发情前期Smad1基因表达量显著高于乏情期和发情间期(P0.05)。在绵羊不同发情时期Smad1蛋白主要表达于乏情期原始卵泡、次级卵泡的颗粒细胞、发情间期次级卵泡的卵泡膜细胞,在发情前期与发情期卵泡中未见表达。卵巢中Smad1基因可能促进黄体溶解,对卵泡发育、卵子的成熟以及释放起作用,与绵羊发情的启动有一定关系。     

发情期和乏情期哈萨克羊垂体差异表达基因筛选及功能预测

垂体是哺乳动物重要的内分泌器官之一,主要通过调控激素合成来调节动物发情。为了研究绵羊垂体系统在发情状态和乏情状态之间的基因表达差异,本研究使用RNA-seq技术分析了发情期和乏情期哈萨克羊垂体前叶的转录组数据。通过比较转录组数据发现:在发情期和乏情期的哈萨克羊垂体中共筛选出了3 211个差异表达的基因,其中包含298个上调的和2 913个下调的基因。利用GO和KEGG富集分析这些差异表达的基因发现:它们被显著地富集在卵母细胞减数分裂和孕酮介导的卵母细胞成熟等调控发情的信号通路上,提示这些基因在调控哈萨克羊发情状态中发挥着不可或缺的作用。本研究丰富了绵羊的转录组资源,为今后深入研究绵羊季节发情机制、提高绵羊繁殖率打下了良好的基础。     

BAD基因在敖汉细毛羊发情不同时期卵巢表达规律及对生殖激素的影响

本实验研究旨在探究BAD基因在敖汉细毛羊中的不同发情时期卵巢中表达量的变化规律以及对生殖激素的影响,为探讨其与绵羊发情的关系奠定基础。首先利用放射免疫方法测定敖汉细毛羊血液中促卵泡激素(FSH)、促黄体生成素(LH)、雌二醇(E_2)和孕激素(P_4)的浓度,然后利用qRT-PCR检测乏情期、发情前期、发情间期和发情期卵巢BAD基因的表达量变化。结果表明:卵巢中BAD基因在发情前期的表达量极显著高于乏情期、发情期和发情间期(P0.01),在发情期的表达量极显著高于乏情期和发情间期(P0.01),乏情期和发情间期的表达量基本相同(P0.05);E_2和LH浓度在BAD基因表达量最高的发情前期维持较高水平,P_4的浓度维持较低水平。综上可知,BAD基因的表达可能促进了E_2和LH浓度升高,抑制了P_4产生,进而影响细毛羊发情,可为绵羊发情的研究提供借鉴。     

Smad2基因在敖汉细毛羊不同发情时期卵巢中的表达分析

Smad2基因是Smads蛋白家族的成员之一,Smad家族中的信号分子对绵羊的繁殖、卵泡发育有重要影响。本实验选取24只3~4岁健康空怀敖汉细毛羊,采用荧光定量PCR、蛋白免疫印迹和免疫组织化学技术,测定Smad2基因在绵羊不同发情时期(乏情期、发情间期、发情前期、发情期)的表达量并分析其在卵泡中表达规律。结果表明:Smad2基因m RNA和蛋白均在绵羊发情前期表达量最高,m RNA水平极显著高于乏情期、发情间期和发情期(P0.01);而发情前期蛋白水平极显著高于乏情期和发情期(P0.01),虽高于发情间期,但差异不显著(P0.05);Smad2蛋白主要在卵巢颗粒细胞、卵泡膜细胞以及卵母细胞表达。据此,Smad2基因可能会促进敖汉细毛羊的发情启动。     

成年滩羊和小尾寒羊皮肤毛囊差异表达miRNA的筛选

为比较miRNA在小尾寒羊和滩羊皮肤毛囊中的表达差异,本研究利用高通量测序技术分析了成年滩羊(TY_1)和小尾寒羊(XWHY_1)皮肤毛囊组织中miRNAs的表达谱,在2个品种绵羊皮肤毛囊组织中共鉴定出561个miRNAs,其中包括138个已知和423个新发现的miRNAs。鉴定出的miRNAs进行表达量差异分析发现,在TY_1与XWHY_1共筛选到63个上调和16个下调的miRNAs。对差异表达miRNA靶基因预测后与基因本体数据库(GO)和京都基因与基因组百科全书(KEGG)比对,分别获得靶基因的注释信息为3 886个和4 449个。GO统计发现,差异表达miRNA的靶基因主要参与代谢过程、催化活性、细胞进程和细胞组分等;而KEGG通路分析表明,4 449个靶基因富集到113个信号通路上,其中在嘌呤代谢、内吞作用和糖酵解/糖异生等信号通路上富集显著。综上,在小尾寒羊和滩羊皮肤毛囊中筛选到的差异表达miRNA可能通过调控其靶基因最终参与了绵羊皮肤毛囊的发育。     

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及其调控的某些关键基因的表达,从而影响前体脂肪细胞外基质与受体的结合、胰岛素利用、脂类和糖的合成与代谢等生物学过程,最终作用于前体脂肪细胞的分化和脂质沉积。     

小尾寒羊高繁殖力候选基因ESR在同期发情处理过程中的表达分析

试验选择高繁殖力小尾寒羊为试验动物,采用公认的CIDR结合外源激素处理使其同期发情,与自然发情及乏情的小尾寒羊作对比,利用反转录荧光定量RT-PCR和Western blotting技术研究雌激素受体(estrogen receptor, ESR)基因在卵巢细胞中的表达差异,同时对其启动子区的甲基化岛进行了分析。结果发现,自然发情绵羊和外源激素处理发情的羊在发情表现和时间上无明显差别;在转录水平的检测结果上也无显著性差异(P＞0.05);但是它们显著高于在乏情期绵羊卵巢细胞内的表达(P＜0.05)。尽管激素处理的同期发情羊卵巢细胞内的ESR蛋白表达量显著高于乏情期绵羊,但是仍然显著低于自然发情期绵羊卵巢细胞的表达(P＜0.05)。此外,自然发情和激素处理的绵羊ESR基因启动子区CpG岛甲基化水平分别为0和10.77%,而发情期绵羊ESR基因启动子处于超甲基化状态（92.30%）。以上结果说明,外源激素处理的同期发情小尾寒羊尽管发情表现、时间等与自然发情的绵羊类似,但是ESR基因在蛋白水平上差异显著,且其启动子区CpG岛甲基化水平也存在显著性差异,这也许是目前胚胎移植后着床率低下的原因之一。     

调控香猪繁殖候选miRNA筛选

研究旨在从香猪卵巢small RNA测序数据中挖掘调控香猪繁殖的microRNA （miRNA）,解析miRNA调控香猪卵巢功能及繁殖性状的分子机制,对于猪的遗传选育具有重要意义。研究选取发情期和间情期的香猪各4头,屠宰后取其卵巢组织,提取总RNA进行small RNA测序,利用生物信息学方法检测miRNA表达谱,筛选差异表达的miRNA,预测差异表达miRNA靶基因,并对靶基因进行GO功能富集和KEGG通路富集分析。结果显示,香猪卵巢组织中miRNA在染色体上呈不均匀分布,主要分布于1号染色体和X染色体上。香猪卵巢中共有627个已知猪miRNAs表达,其中有34个差异表达miRNAs,表达量前五的分别是miR-23、let-7i-5p、miR-103、miR-30e-5p和miR-1271-5p。GO功能分析结果显示,靶基因主要参与的生物学过程是细胞过程（cellular process）,主要分布于细胞（cell）,主要分子功能是结合（binding）;KEGG显著富集通路中,促性腺激素释放激素受体通路（gonadotropin-releasing hormone receptor pathway）、胰岛素样生长因子通路（IGF pathway）和表皮生长因子受体信号通路（EGF receptor signaling pathway）与卵母细胞的发育成熟相关,因此推测miR-23b、let-7i-5p、miR-103、miR-30e-5p和miR-1271-5p可能参与了香猪繁殖调控。本研究初步筛选出5个可能调控香猪繁殖性能的miRNAs,可为从分子水平提高香猪产仔数提供理论基础。     

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

非繁殖季节补饲精料诱导哈萨克羊发情及相关生殖激素变化规律的探究

为了研究补饲精料诱导非繁殖季节哈萨克母羊发情,解析绵羊非繁殖期发情的生殖激素变化规律,为利用营养因素调控非繁殖季节绵羊发情提供激素方面的理论依据,研究选用成年哈萨克母羊36只,随机分成2组,试验组在2—6月份补饲精料,在非繁殖季节(4—6月份)对发情绵羊及不发情母羊(对照)颈静脉采血,采用酶联免疫法(ELISA)测定外周血液4种重要生殖激素雌二醇(E2)、孕酮(P4)、促卵泡激素(FSH)和促黄体素(LH)的浓度。结果表明:血液中E2和P4呈现波动式变化,FSH和LH呈现规律性脉冲式变化。说明E2和FSH与绵羊的非繁殖季节发情有关,高浓度的P4和LH对开启或维持绵羊的发情有重要作用。     

鸡胚胎期肌纤维形成相关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-...     

感染产肠毒素大肠杆菌的猪小肠上皮细胞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及其作用机制提供科学依据。     

鸡白痢沙门菌感染雏鸡的骨髓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表达谱特征,为了解鸡白痢沙门菌和鸡之间相互作用的复杂分子致病机制提供了新的见解,为防控鸡白痢沙门菌感染提供了新策略。