乏情和发情初产母猪卵巢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酶结合、泛素样蛋白连接...     

营养诱导舍饲哈萨克羊非繁殖季节卵巢组织中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是通过营养代谢通路对绵羊的非季节性发情起重要作用。     

香猪卵巢circ_CSPP1的鉴定及功能预测分析

研究旨在揭示香猪卵巢组织中circ_CSPP1的潜在功能。利用CIRI-full拼接circ_CSPP1的序列，circBase数据库在线blat分析circ_CSPP1的保守性，采用RT-PCR和Sanger测序技术从香猪卵巢分离和鉴定circ_CSPP1，采用miRanda、RNAhybrid和PITA软件预测和分析circ_CSPP1结合的miRNA分子及miRNA的靶基因，对靶基因进行GO和KEGG富集分析。结果显示，分离到的circ_CSPP1是保守的，由CSPP1基因的外显子8~12组成，可作为ssc-miR-324、ssc-miR-10a-5p、ssc-miR-9847-3p、ssc-miR-365-5p、ssc-miR-92b-5p、ssc-miR-885-3p、ssc-miR-7139-3p、ssc-miR-9851-3p和ssc-miR-33b-3p等14个miRNAs的分子海绵，靶向755个蛋白编码基因。GO富集分析表明，circ_CSPP1的靶基因显著富集到153个GO条目，其中细胞成熟、细胞分裂位点、细胞迁移的负调控、Wnt信号通路、细胞周期蛋白结合、子宫发育、MAPK活性的激活、黏着斑和卵裂沟等条目与卵巢功能相关。KEGG富集分析表明，circ_CSPP1的靶基因富集到260个信号通路，其中显著富集的有29个信号通路，主要富集于癌症通路、细胞周期、雌激素信号通路、PI3K-AKT信号通路、卵巢类固醇生成、催产素信号通路和细胞凋亡等信号通路。综上，circ_CSPP1可作为许多miRNA的分子海绵调控香猪卵巢颗粒细胞的增殖与凋亡、卵泡成熟、排卵过程、卵巢类固醇生成和卵巢早衰等生物学过程。     

生产单一性别羔羊的奶山羊卵巢组织差异表达microRNAs筛选

运用Solexa/Illumina高通量测序技术分析生产单一性别羔羊的奶山羊卵巢miRNAs表达情况,探讨miRNAs在繁殖调控中的作用。采集4~6岁始终生产单一性别羔羊的奶山羊母羊双侧卵巢,提取总RNA,建立cDNA文库,并进行克隆、测序及生物信息学分析。结果表明,在单一产公羔奶山羊(ASGO1)卵巢中有744 782个reads,包含127 205条miRNAs。单一产母羔奶山羊(ASGO2)中有1 602 751个reads,包含454 911条miRNAs,ASGO1和ASGO2差异表达的miRNA有630种,有143种仅在单一产母羔奶山羊卵巢组织中表达,有70种在单一产公羔奶山羊中高表达,其中新miRNA有31种,且都是单一产母羔奶山羊的表达量明显高于单一产公羔奶山羊,单一产母羔奶山羊中有2个miRNA(ssc-miR-451_R-1和bta-miR-486_R-1)与单一产公羔奶山羊相比差异极显著,bta-miR-486_R-1差异最大,达到了800 000个拷贝数的高表达量水平。结果提示,ssc-miR-451_R-1和bta-miR-486_R-1可能是影响奶山羊选择性授精的重要miRNA。     

饲粮中添加共轭亚油酸对猪肌肉组织miRNA表达谱的影响

旨在探讨共轭亚油酸(CLA)通过影响猪肌肉miRNA表达调控肌肉代谢的分子机制。本研究选用体重相近的初产荣昌母猪12头,随机分为对照组和CLA组,CLA组母猪饲粮从妊娠初期开始添加1.5%CLA,持续到仔猪28日龄断奶;断奶后分别从对照组和CLA组挑选仔猪各30头,原CLA组仔猪饲粮中继续添加1.5%CLA。待仔猪体重约30kg时,每组各选择6头进行屠宰,采集背最长肌和腿肌样品。提取RNA后将对照组和CLA组的背肌和腿肌进行建库;通过Solexa高通量测序技术检测CLA对猪肌肉组织miRNA表达谱的影响,利用生物信息学进行差异显著表达miRNA的靶基因预测和功能分析。结果显示:1)背肌和腿肌分别获得了44 869 982条和45 105 806条clean reads,大多数序列长度在20~23nt。2)背肌和腿肌分别鉴定到306和304个已知miRNAs,有295个miRNAs共表达。3)添加CLA后能改变猪肌肉中miRNA的表达,分析表明,CLA对腿肌miRNA表达的影响要强于背肌;背肌和腿肌中分别发现5个和12个差异显著表达的miRNAs(P0.05),其中ssc-miR-224在两种组织中都差异表达,因此共有16个差异表达miRNAs。4)对16个差异表达miRNAs进行靶基因预测和KEGG通路富集分析,发现它们的靶基因主要富集于292个通路,其中24个通路显著富集(P0.05),包括与肌肉代谢密切相关的MAPK信号通路和Notch信号通路。5)利用荧光定量PCR对随机选择的6个差异表达miRNAs进行验证,证实其表达水平和测序结果基本一致。以上结果说明,CLA可能通过影响肌肉miRNA表达,调控重要肌肉代谢通路。     

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

鸡性成熟启动前后下丘脑microRNA表达谱分析

本研究旨在探讨miRNA在鸡性成熟启动调控中的作用。利用Solexa高通量测序技术检测了鸡性成熟启动前(Before puberty onset,BPO)和性成熟启动后(After puberty onset,APO)下丘脑中miRNA表达谱,通过生物信息学方法进行靶基因预测和功能分析。结果表明,在鸡下丘脑中鉴定出374个已知miRNAs、46个新miRNAs,144个已知miRNAs(reads10)的表达水平在性成熟启动前后发生了显著改变(P0.05)。5个差异表达miRNAs被qRT-PCR方法进一步验证。15个差异最显著miRNA(|log2(fold-change)|2.0,P0.01)预测到2 013个靶基因,GO富集和KEGG调节通路分析结果表明,这些差异表达miRNAs在性成熟启动转录调节和信号转导通路中发挥了重要作用。综上表明,miRNA是参与鸡性成熟启动新的调节因子。鉴于miRNA功能的保守性,该研究结果将为深入开展动物(包括人类)性成熟的分子调控机制研究提供新的理论依据。     

乌头驴与三粉驴骨骼肌中miRNAs的表达鉴定与分析

试验旨在分析microRNAs （miRNAs）在德州驴骨骼肌中的表达情况,探讨小RNA （small RNA）在驴骨骼肌发育中的调控机制。采集德州驴2个品系（乌头驴和三粉驴）的背最长肌组织,构建乌头驴背最长肌（WB）和三粉驴背最长肌（SB）的small RNA测序文库并进行测序,运用生物信息学技术对WB和SB之间差异表达的miRNAs进行分析,预测和解析可能影响驴产肉性状的miRNAs,从中挑选了6个表达量较高的miRNAs,利用实时荧光定量PCR技术来验证测序结果的准确性。结果表明,在WB和SB文库中分别鉴定出保守miRNAs有982和1 149个,新miRNAs有106和143个。在WB与SB比较组中基于条件|log₂（fold_change）|≥2、P≤0.05获得17个差异表达的miRNAs,其中表达量上调的有5个,下调的有12个。通过GO注释,发现差异表达的miRNAs显著富集于调控细胞迁移分化、细胞分裂和骨骼肌细胞收缩等生物学过程（P<0.05）,并且挖掘到与驴骨骼肌发育相关的bta-miR-378d_L+1R-1_1ss22CA、bta-miR-378d_R-2,其靶向NPY1R、GPR152、BEST1、KCNH8、SPAG9等多个基因。KEGG富集结果表明,miRNAs富集于Wnt、MAPK、泛素蛋白水解系统等与骨骼肌发育相关的信号通路中,其中在Wnt信号通路中共涉及到17个差异表达miRNAs的207个靶基因,尤其以PC-5p-84483_31、hsa-miR-186-3p_R-3和cfa-miR-329b_1ss19CT的靶基因最多。实时荧光定量PCR结果与miRNAs测序结果一致,并发现除eca-miR-181b外,其他5个miRNAs在SB中的表达水平均高于WB。本研究首次对驴骨骼肌miRNAs进行转录组学分析,揭示了乌头驴和三粉驴miRNAs的表达差异,研究结果为阐明德州驴骨骼肌生长发育的分子调控机制提供参考。     

罗伊氏乳杆菌对仔猪空肠miRNAs表达谱的影响

试验旨在研究罗伊氏乳杆菌对仔猪空肠miRNAs表达的影响。选用1日龄健康约克×荣昌仔猪6窝,随机分为两组,每组3窝（共27头）。即对照组（每天每头灌喂0.1%的无菌蛋白胨溶液）,LR组（每天每头灌喂活菌总数为1.0×10¹⁰CFU的罗伊氏乳杆菌）,分别在10和20日龄,每组选取6头进行屠宰,采集空肠样品。通过Solexa高通量测序技术检测空肠miRNA的表达谱,利用生物信息学技术进行靶基因预测和功能分析。结果显示：1）试验构建了8个文库,共获得187 103 840条纯净reads序列,其中22 nt长度序列占的比例最大;2）10和20日龄分别鉴定出354和352个已知的miRNAs,其中有329个miRNAs共表达;3）灌喂罗伊氏乳杆菌组仔猪空肠miRNA差异表达,10和20日龄分别获得7和9个差异显著的miRNAs,ssc-miR-218在两个日龄中都差异表达;4）对差异表达的miRNAs进行靶基因预测和KEGG通路分析,发现有10 221个靶基因富集到293条通路上,有66条通路显著富集（P<0.05）,包括与肠道发育及免疫调控相关的MAPK和PI3K-Akt信号通路;5）随机挑选6个差异表达miRNAs进行荧光定量PCR验证,验证结果与测序结果基本一致。本研究表明,罗伊氏乳杆菌可能通过影响空肠miRNAs表达来调控有关肠道健康的信号通路。     

宫内生长受限对胎猪肝脏microRNA表达谱的影响

本试验旨在分析宫内生长受限(IUGR)和正常胎猪肝脏microRNA(miRNA)表达的差异,以揭示IUGR对胎猪肝脏miRNA表达谱的影响。试验选用6头达兰(Dalland)母猪,在孕110 d左右致死,剖腹并取出整个子宫。在每窝中选取1头IUGR胎猪和1头正常胎猪,分别取出肝脏备用。利用双荧光标记miRNA芯片技术对IUGR和正常胎猪肝脏进行分析。结果表明:IUGR对胎猪肝脏miRNA的表达产生显著影响,共有41种miRNA存在表达差异,其中27种在IUGR胎猪肝脏中表达上调,14种表达下调。其所导致的胎猪葡萄糖吸收代谢紊乱、脂肪代谢紊乱、蛋白质合成受阻;肝脏解毒能力下降以及肝脏细胞增殖抑制,可能是造成IUGR仔猪出生后营养利用率低、患病率和死亡率高以及生长发育不良的重要原因。     

猪繁殖与呼吸综合征病毒相关猪microRNA的初步筛选

本试验旨在初步筛选出与猪繁殖与呼吸综合征病毒(PRRSV)感染相关的猪microRNA(miRNA)。采用生物信息学方法预测可能与PRRSV 3′非编码区(3′utr)产生相互作用的猪miRNA,构建PRRSV 3′utr的双荧光素酶报告基因载体psi-CHECK-utr,同时合成预测到的miRNA(ssc-miR-323、ssc-miR-105-1)及其相应的抑制物,共转染猪脐静脉血管内皮细胞系(SUVEC),检测双荧光素酶活性。结果发现,ssc-miR-323与psiCHECK-utr共转染后,荧光比率略有升高,ssc-miR-105-1与psiCHECK-utr共转染后,荧光比率降低。初步判定ssc-miR-323对PRRSV 3′utr没有抑制作用,而ssc-miR-105-1对PRRSV3′utr有一定的抑制作用。     

Effect of Lactobacillus reuteri on the Expression Profile of miRNAs in the Jejunum of Piglets

This study was conducted to investigate the effect of Lactobacillus reuteri on the expression of miRNAs in the jejunum of piglets. Six litters of healthy 1-day-old York×Rongchang piglets were selected and randomly divided into the two group,3 litters per group(n=27). The two groups were control group (daily orally adminitrated with 0.1% sterile peptone solution per piglet),and the LR group (daily orally administrated with 1.0×10¹⁰ CFU of L. reuteri per piglet), respectively. At 10 and 20 days of age, six piglets per group were selected for slaughter, respectively, and the jejunum samples were collected. After RNA extraction, the library was constructed and the miRNA expression profile of the jejunum was detected by Solexa high-throughput sequencing technology. Bioinformatics technology was used to conduct target genes prediction and functional analysis. The results showed that: 1) A total of 8 libraries were constructed and a total of 187 103 840 clean reads were obtained, of which the 22 nt length sequence accounted for the largest proportion; 2) Three hundred and fifty-four and 352 known miRNAs were respectively identified at 10 and 20 days of age, of which 329 miRNAs were co-expressed; 3) Differentially expressed (DE) miRNAs in the jejunum were found after L. reuteri administration, 7 and 9 DE miRNAs were obtained at 10 and 20 days of age,respectively. Notably, ssc-miR-218 was expressed differentially at two ages. 4) KEGG pathway analysis and target gene prediction on the differentially expressed miRNAs revealed that 10 221 target genes were enriched in 293 biological pathways, of which 66 were significantly enriched (P<0.05), including the MAPK signaling pathway and PI3K-Akt signaling pathway, which were known to regulate intestinal development and immune. 5) Six miRNAs were selected to validate the sequencing results by qRT-PCR, the qRT-PCR expression results corresponded well with those from the sequencing. The results indicates that L. reuteri may regulate intestinal health related pathways by influencing the expression of jejunal miRNAs.     

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

发情和乏情初产母猪下丘脑-垂体-卵巢轴差异表达基因筛选及分析

研究旨在分析发情和乏情初产母猪下丘脑-垂体-卵巢轴基因表达的差异,探讨母猪乏情调控的分子机制。选用6头健康的断奶初产母猪,分为发情组（3头）和乏情组（3头）,屠宰后分别采集下丘脑、垂体、卵巢进行转录组测序（RNA-Seq）、GO功能富集、KEGG通路及相关蛋白的Western blotting分析。结果发现,发情组下丘脑、垂体、卵巢分别获得52 432 800、52 573 730和52 209 252条clean reads,与猪参考基因组（Sus scrofa 11.1）的比对率分别为78.69%、81.49%和79.26%;乏情组下丘脑、垂体、卵巢分别获得52 516 724,52 476 820和52 195 962条clean reads,与猪参考基因组的比对率分别为82.38%、83.05%和80.20%。与发情组母猪相比,乏情组母猪下丘脑中共有710个差异表达基因,其中上调基因392个,下调基因318个;垂体中共有707个差异表达基因,其中上调基因283个,下调基因424个;卵巢中共有956个差异表达基因,其中上调基因635个,下调基因321个。3种组织中的共有差异表达基因为36个。与母猪情期调控相关的亲吻素-1（KISS1）、G-蛋白偶联受体54（GPR54）、速激肽3（TAC3）、速激肽3受体（TACR3）、17-α-羟化酶/17,20碳链裂解酶（CYP17A1）、芳香化酶（CYP19A1）、类固醇激素合成急性调节蛋白（STAR）、促性腺激素释放激素受体（GnRHR）和雌激素受体1（ESR1）基因在乏情组母猪体内显著下调（P<0.05;P<0.01）。Western blotting分析结果显示,TAC3、TAC3R、KISS1和GPR54相关蛋白在乏情组母猪下丘脑中也显著下调（P<0.05）。GO和KEGG通路分析发现,差异表达基因显著富集于正调控胆固醇酯化反应、卵巢类固醇生成、GnRH信号通路、FoxO信号等一些与类固醇激素生成和卵泡发育相关的信号通路。本研究结果表明,发情和乏情初产母猪下丘脑-垂体-卵巢轴上与情期调控相关的基因存在差异表达,尤其是Kisspeptin/GPR54和TAC3/TACR3两大系统的mRNA及蛋白表达水平在乏情组母猪下丘脑中均显著下调,推测可能是Kisspeptin/GPR54和TAC3/TACR表达受损导致了下丘脑调节功能障碍,进而导致了初产母猪乏情。该研究为揭示初产母猪乏情分子调控机制提供了重要支持,为今后利用分子技术改善母猪乏情问题提供了重要理论依据。     

感染肺炎支原体姜曲海猪肺组织差异表达circRNAs的筛选与分析

为探明姜曲海猪感染猪肺炎支原体（Mycoplasma hyopneumoniae,Mhp）后肺组织环状RNA （circular RNA,circRNA）差异性表达谱及其在抗Mhp感染中的作用,试验以姜曲海猪为研究对象,分为感染组和对照组,人工感染Mhp 28 d后,解剖采集肺组织,采用高通量测序和生物信息学软件分析circRNA的表达情况。测序数据比对参考猪基因组序列,共鉴定到23 632个circRNAs,差异表达circRNAs为213个,其中97个上调,116个下调。随机选择4个差异表达circRNAs进行实时荧光定量PCR验证,检测结果与测序结果基本一致。差异表达circRNA来源基因可注释到包括抗原加工递呈、溶酶体、白细胞跨内皮迁移等免疫应答信号通路。circRNA-miRNA-mRNA靶标关系分析显示,筛选到海绵结合miRNA数量最多的3个差异表达circRNAs,分别为：circRNA-17284（21个）、circRNA-04848（19个）和circRNA-17270（19个）;预测到6个与免疫调控相关的靶向miRNAs,分别为：ssc-miR-4331、ssc-miR-370、ssc-miR-328、ssc-miR-30c-3p、ssc-miR-122和ssc-miR-125b。本研究测定了感染Mhp的猪肺组织circRNA表达谱,筛选到与免疫调控相关的差异表达circRNA,这有助于阐明姜曲海猪对Mhp的易感机制,为抗病育种研究提供了参考依据。     

Identification of a novel miRNA from the ovine ovary by a combinatorial approach of bioinformatics and experiments

MicroRNAs (miRNAs) are a class of short endogenous, single-stranded, non-coding small RNA molecules, about 19–25 nucleotides in length that regulate gene expression at the translation level and influence many physiological process, such apoptosis, metabolism, signal transduction, and occurrence and development of diseases. In this study, we constructed a library from the ovine luteal phase ovary by using next-generation sequencing technology (Solexa high-throughput sequencing technique) and identified 267 novel miRNAs by bioinformatics. One of the novel miRNAs (ovis_aries_ovary-m0033_3p), which expressed in the sheep ovary and testis, was confirmed by real time PCR and northern blot. Ovis_aries_ovary-m0033_3p was 21 nucleotides in length and located on chromosome 12, and it had 100% similarity to hsa-miR-214-3p, mmu-miR-214-3p, dre-miR-214and ssc-miR-214. Meanwhile, the pre-miRNA was 82 nucleotides in length and had a standard hairpin stem-loop structure. From the consistency of the sequence and structure, we speculated that ovis_aries_ovary-m0033_3p had a function similar to hsa-miR-214-3p, which is involved in the fine regulation of cell survival, embryonic development, breeding activities and resistance to ovarian cancer, so we defined it as oar-miR-214-3p. These experimental results will enrich the miRNA database for ovis aries and provide the basis for researching the regulation mechanism of miRNA in relation to breeding activities of seasonal breeding animals.     

MicroRNA fingerprints in serum and whole blood of sarcoid‐affected horses as potential non‐invasive diagnostic biomarkers

Serum and whole blood microRNA (miRNA) fingerprints have been proposed as a new class of non‐invasive human cancer biomarkers. In this study, we compared equine sarcoid (ES) disease‐specific serum and whole blood miRNA fingerprints and correlated them to miRNA expression in sarcoid tissue. After high throughput sequencing, miRNA differential expression analysis between six ES‐affected and five control horses was carried out in serum and whole blood using a DESeq algorithm, accounting for the influence of hemolysis and the white blood cell count. Target gene, pathway prediction and enrichment analyses were conducted using TarBase, mirPath and GeneCodis. After exclusion of 4 hemolyzed out of a total of 11 serum samples, 9 miRNAs were found to be differentially expressed in serum of ES vs control horses. In whole blood, all 11 samples showed normal white blood cell counts and 19 miRNAs were found to be differentially expressed. A total of 2/9 serum and 7/19 whole blood differentially expressed miRNAs were also highly expressed at the tissue level and their predicted target genes were associated with cancer pathways. Serum and whole blood miRNA expression allowed discrimination between ES and control horses and merits further validation in a larger study cohort. The use of whole blood might be superior because it has higher miRNA content and is less influenced by pre‐analytical variables compared to serum. Concurrent dysregulation of single miRNAs in tissue and blood suggests a possible biological function of circulating miRNAs.