单双羔绵羊卵巢组织差异表达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可能与绵羊卵泡发育和产羔数多少有关。     

鸡白痢沙门氏菌感染对雏鸡脾脏miRNA表达谱的影响

【目的】鉴定感染鸡白痢沙门氏菌后雏鸡脾脏差异表达的miRNA，为阐明miRNA在鸡白痢沙门氏菌感染中的调控机制提供理论基础。【方法】采集感染和未感染鸡白痢沙门氏菌的SPF雏鸡脾脏，提取脾脏RNA，构建6个miRNA文库；利用Illumina Hiseq 2500测序技术和生物信息学分析筛选差异表达miRNA，通过TargetScan和miRanda算法预测差异表达miRNA的靶基因，并进行GO功能富集和KEGG通路富集分析，最后利用RT-qPCR方法验证测序结果。【结果】鉴定到29个差异表达miRNAs，其中15个显著上调、14个显著下调。GO分析表明：脂多糖介导的信号通路正调控以及NIK/NF-κB介导的信号通路负调控等免疫相关生物学过程显著富集；KEGG分析显示：差异表达的miRNA主要参与溶酶体和内吞作用等免疫相关通路。高通量测序结果与RT-qPCR结果一致。【结论】成功鉴定到鸡白痢沙门氏菌感染雏鸡脾脏miRNA表达谱，获得29个鸡白痢沙门氏菌感染潜在相关miRNAs。     

基于RNA-seq技术筛选山羊卵泡发育相关下调基因

【目的】对山羊优势卵泡(DF)和从属卵泡(SF)颗粒细胞进行高通量测序,筛选并研究与山羊卵泡发育相关的下调基因。【方法】以1岁龄贵州白山羊为供试对象,在其发情3 d后,采集第一卵泡波中的优势卵泡(DF)和从属卵泡(SF),分别分离其颗粒细胞,提取RNA并进行文库构建、Illumina测序,将获得的序列与山羊的RefSeq数据库进行比对后,用DESeq2软件对获得的RNA进行差异表达分析,并对表达下调的基因进行GO分析和KEGG信号通路分析,最后通过Real-time PCR对筛选出的表达下调基因进行验证。【结果】将测序结果与山羊的RefSeq数据库比对后,共获得了33 896个基因。经DESeq2软件对获得的mRNA进行差异表达分析,共获得695个差异表达m-RNA,其中462个表达显著下调,233个表达则显著上调;对462个表达下调基因进行GO分析,共分为2大类39组,其中参与生物学过程(BP)的基因占48.7%;与细胞组分(CC)有关的基因占51.3%。KEGG信号通路分析发现20条通路,其中参与癌症通路基因富集最为显著。从下调基因中筛选出6个可能会抑制山羊卵泡发育的基因,分别为PTX3、LDLR、RGN、NID1、PDLIM5、PRLR。Real-time PCR结果显示,RGN、NID1、PDLIM5、PRLR在SF与DF中的表达趋势与高通量测序结果一致,且RGN在SF颗粒细胞中的表达量极显著高于DF(P0.01),NID1、PDLIM5在SF颗粒细胞中的表达量显著高于DF(P0.05);PTX3、LDLR在SF与DF中表达量与高通量测序结果相反。【结论】获得的6个表达下调基因在山羊卵泡发育过程中可能会抑制卵泡的发育,最终导致卵泡闭锁。     

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

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

外泌体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的表达调控,抑制颗粒细胞细胞周期运转和类固醇激素代谢等信号通路,引起颗粒细胞凋亡,导致卵泡闭锁。     

藏黄牛与宣汉黄牛心脏miRNA表达谱比较

【目的】miRNA作为一类非编码RNA,广泛参与机体多种生命活动,研究旨在挖掘miRNA在藏黄牛和宣汉黄牛心脏组织中的差异miRNA,为进一步研究藏黄牛低氧适应的分子机制提供基础数据。【方法】随机选取健康的藏黄牛和宣汉黄牛各3头,迅速采集其心脏组织,使用Trizol法提取RNA,琼脂糖凝胶电泳切胶选取18—30nt的片段,连接3'端和5'端,反转录后进行扩增,凝胶电泳切胶纯化后建立藏黄牛和黄牛各3个文库,利用Illumina HiSeq4000测序平台进行高通量测序;将测序得到的序列进行过滤,比对GenBank和Rfam数据库,筛选出藏黄牛和宣汉黄牛的差异miRNA,进行功能注释及信号通路富集分析;随机选择8个miRNAs,利用实时荧光定量PCR检测其在心脏组织的表达量,以验证测序数据的准确性。【结果】藏黄牛和宣汉黄牛心脏组织的高质量读值的序列分别为17 463 446条和13 662 812条,干净读值为16 552 296条和12 055 304条,且在藏黄牛和宣汉黄牛中高质量核酸序列长度富集最多的均是21nt,分别为37.5%和32.1%;且共筛选出219个差异miRNAs,其中48个显著上调,171个显著下调;GO功能注释得到差异miRNA靶基因分子功能中显著富集的条目有22条,如,GO:0005488(结合)、GO:0005515(蛋白质结合)和GO:0043167(离子结合);细胞组分中显著富集的条目有20条,如,GO:0005623(细胞)、GO:0044464(细胞组分)和GO:0005622(细胞内);生物过程中显著富集的条目有13条,如,GO:0035556(细胞内信号转导)、GO:0032774(RNA生物合成过程)和GO:0006351(转录,DNA模板化),KEGG信号通路分析得到差异表达miRNAs靶基因显著富集到胰岛素信号通路(139个靶基因)、mTOR信号通路(38个靶基因)和HIF-1信号通路(92个靶基因)等232个信号通路中,其中有12个靶基因共同作用于这3个信号通路,说明miRNAs可能通过这3个信号通路,共同参与藏黄牛低氧适应性的调控;随机选取8个miRNAs进行荧光定量验证,其表达趋势与测序结果表达趋势基本一致,表明高通量测序数据可信度较高。【结论】得到了miRNA在藏黄牛与宣汉黄牛心脏组织中的表达谱,为揭示藏黄牛低氧适应性的调控机制研究奠定了基础。     

深度测序鉴定玉米雄穗花器官分化期响应干旱胁迫的miRNA和其靶基因

【目的】 运用现代分子生物学及高通量测序技术,研究不同品种玉米雄穗花器官分化期响应干旱胁迫的差异miRNA和其靶基因,挖掘和鉴定与玉米发育相关的基因,分析其性状差异的信号通路及分子调控网络。【方法】 以耐旱自交系"PHBA6"和干旱敏感自交系“吉63”为研究对象,应用深度测序技术进行miRNA文库构建,鉴定其差异表达的 miRNA,对差异表达mi RNA及其靶基因进行功能注释、聚类分析及通路富集。【结果】 鉴定了总共337种前体miRNA,其中包含289种已知的miRNA和48种新的miRNA。在3个文库,两个分组中共有155种差异表达miRNA。基于GO功能分类及遗传和基因组(KEGG)的功能富集显示,这些miRNA可能通过靶向一系列与胁迫相关的基因而在干旱胁迫中发挥作用。至少55个预测的靶基因进一步被60个miRNA调控。NAC,MYB和MAPK基因家族在干旱胁迫下评分最高,在植物抗旱中重要作用。一系列miRNA与这些排名靠前的基因相关,包括miR164、miR172、miR1520、miR6158、ghr-n24、ghr-n56等。【结论】 miRNAs可能在玉米雄穗花器官分化期耐旱中发挥重要作用,miRNAs的筛选为分子辅助育种和转基因育种将提供新的靶标。     

藏黄牛与宣汉黄牛心脏miRNA表达谱比较

【目的】miRNA作为一类非编码RNA,广泛参与机体多种生命活动,研究旨在挖掘miRNA在藏黄牛和宣汉黄牛心脏组织中的差异miRNA,为进一步研究藏黄牛低氧适应的分子机制提供基础数据。【方法】随机选取健康的藏黄牛和宣汉黄牛各3头,迅速采集其心脏组织,使用Trizol法提取RNA,琼脂糖凝胶电泳切胶选取18—30nt的片段,连接3'端和5'端,反转录后进行扩增,凝胶电泳切胶纯化后建立藏黄牛和黄牛各3个文库,利用Illumina HiSeq4000测序平台进行高通量测序;将测序得到的序列进行过滤,比对GenBank和Rfam数据库,筛选出藏黄牛和宣汉黄牛的差异miRNA,进行功能注释及信号通路富集分析;随机选择8个miRNAs,利用实时荧光定量PCR检测其在心脏组织的表达量,以验证测序数据的准确性。【结果】藏黄牛和宣汉黄牛心脏组织的高质量读值的序列分别为17 463 446条和13 662 812条,干净读值为16 552 296条和12 055 304条,且在藏黄牛和宣汉黄牛中高质量核酸序列长度富集最多的均是21nt,分别为37.5%和32.1%;且共筛选出219个差异miRNAs,其中48个显著上调,171个显著下调;GO功能注释得到差异miRNA靶基因分子功能中显著富集的条目有22条,如,GO:0005488（结合）、GO:0005515（蛋白质结合）和GO:0043167（离子结合）;细胞组分中显著富集的条目有20条,如,GO:0005623（细胞）、GO:0044464（细胞组分）和GO:0005622（细胞内）;生物过程中显著富集的条目有13条,如,GO:0035556（细胞内信号转导）、GO:0032774（RNA生物合成过程）和GO:0006351（转录,DNA模板化）,KEGG信号通路分析得到差异表达miRNAs靶基因显著富集到胰岛素信号通路（139个靶基因）、mTOR信号通路（38个靶基因）和HIF-1 信号通路（92个靶基因）等232个信号通路中,其中有12个靶基因共同作用于这3个信号通路,说明miRNAs可能通过这3个信号通路,共同参与藏黄牛低氧适应性的调控;随机选取8个miRNAs进行荧光定量验证,其表达趋势与测序结果表达趋势基本一致,表明高通量测序数据可信度较高。【结论】得到了miRNA在藏黄牛与宣汉黄牛心脏组织中的表达谱,为揭示藏黄牛低氧适应性的调控机制研究奠定了基础。     

山羊环状RNA circ_ZCCHC24的鉴定和表达分析

【目的】本文对山羊卵泡差异表达的环状RNA circ_ZCCHC24进行了成环鉴定和表达定位分析。【方法】分别采用生物信息学技术、Sanger测序技术、实时荧光定量PCR和荧光原位杂交技术,研究了circ_ZCCHC24的序列和结构组成、接头序列中的环化位点和表达定位情况。【结果】circ_ZCCHC24是由山羊ZCCHC24基因的部分第1内含子和第2外显子环化而成,sanger测序结果也证明了这一结果;circ_ZCCHC24在麻城黑山羊卵泡中的表达量显著高于波尔山羊卵泡(P0.01)。circ_ZCCHC24表达具有组织特异性,在波尔山羊肾脏高表达,卵泡中低表达;circ_ZCCHC24定位于山羊卵泡颗粒细胞的胞质中,发挥"miRNA吸附海绵"的功能。【结论】本研究为进一步研究circ_ZCCHC24调控母羊卵泡发育的分子机理,挖掘出高产母羊的特有调控通路奠定了基础,也为高繁殖力山羊品种的培育提供新靶点。     

猪心肌和骨骼肌miRNA转录组比较分析

【目的】研究microRNA(miRNA)在猪不同类型肌肉生长发育中的作用。【方法】以猪心肌以及2种骨骼肌(背最长肌和腰大肌)为材料,采用Illumina高通量测序技术鉴定其miRNA转录组,并进行差异表达分析和靶基因富集分析。【结果】从猪心肌、背最长肌和腰大肌组织的3个小RNA测序文库中共获得约56M的序列,通过生物信息分析共鉴定出907个非冗余miRNA,其中441个miRNA在3个文库共同表达,123个miRNA在文库间表现为显著的差异表达(P10-6)。对3个文库差异表达miRNA的相关性分析发现,背最长肌与腰大肌的相关性最高(r=0.94),证明2种骨骼肌的表达模式十分相近,而心肌与2种骨骼肌的表达模式不同。通过对各文库高表达的miRNA分析发现,表达量前十的miRNA中有6个miRNA的表达量在心肌中极显著上调。对在心肌内上调的6个miRNA进行的靶基因预测和生物学功能预测发现,这些差异miRNA主要靶向作用于钙信号通路和胰岛素信号通路。荧光定量PCR(qRT-PCR)验证结果表明高通量测序结果准确可靠。【结论】从miRNA层面揭示了猪心肌和骨骼肌的差异,这种差异在一定程度上印证了心肌和骨骼肌之间迥异的生理学特征。     

新疆野苹果幼苗叶片响应NaCl胁迫的miRNAs及靶基因筛选

【目的】 筛选新疆野苹果响应NaCl胁迫相关差异miRNA及靶基因。【方法】 以新疆野苹果组培苗为材料,对150 mmol/L NaCl处理6和48 h时的新疆野苹果幼苗叶片进行small RNA测序。【结果】 NaCl处理6 h时3个miRNA差异表达,其中2个上调表达,1个下调表达,miRNA对应的靶基因数目为64个;NaCl处理48 h时13个miRNA差异表达,其中4个上调表达,9个下调表达,miRNA对应的靶基因数目为108个。对差异靶基因进行GO和KEGG分类注释,NaCl处理6 h时,GO主要富集包括:肌醇磷酸2激酶活性、寡糖基转移酶活性、蛋白丝氨酸/苏氨酸磷酸酶抑制剂活性、木聚糖生物合成过程等;NaCl处理48 h时,GO主要富集包括:氧氧化还原酶活性、DNA结合、质外体、次生代谢过程等。KEGG共同富集在N-glycan生物合成通路。【结论】 mdm-miR168b、mdm-miR159c、mdm-miR827、mdm-miR390f、mdm-miR171i、mdm-miR399j基因在NaCl处理6和48 h时表达量存在差异,其中mdm-miR390f及其靶基因HF10976、mdm-miR171i及其靶基因HF33844、mdm-miR399j及其靶基因HF11095表达量呈负相关,3个miRNA参与了新疆野苹果对NaCl胁迫的响应过程。     

山羊卵泡发育相关基因的筛选及分析

【背景】山羊第一卵泡波中的优势卵泡(dominant follicles, DF)和从属卵泡(subordinate follicles,SF)是整个卵泡发育过程中最为关键的两个阶段。随着卵泡的进一步发育,最终DF可能发育成为成熟卵泡,直到排卵;SF将走向闭锁,其中颗粒细胞的凋亡是导致卵泡发生闭锁的关键因素。然而目前对促进卵泡的优势化或导致其闭锁的分子机理尚不清楚。【目的】通过对山羊第一卵泡波中DF和SF颗粒细胞进行高通量测序,旨在筛选影响卵泡发育的关键基因,为深入探究卵泡发育的调控机制提供理论依据。【方法】选取10只1岁龄健康的贵州白山羊分别注射前列腺素F2α,使其同期发情,此后每天用B超检测并记录卵泡的生长情况,发情3 d后,统一屠宰并采集第一卵泡波中DF (直径4.5—6 mm)与SF (直径3—4.5 mm),分别分离其中的颗粒细胞,提取总RNA、构建文库后通过Illumina Hiseq 2500平台进行测序。利用FastQC对测序产出raw reads进行质量评估并经过过滤后,获得品质较高的clean reads;使用Trinity对得到的clean reads进行重新组装,从而获得unigenes;使用CLC Genomics Workbench将unigenes与山羊RefSeq数据库进行比对获得mRNA;使用DESeq2软件对获得的mRNA进行差异表达分析;分别采用goseq和kobas软件对得到的差异表达基因进行GO分析及KEGG信号通路分析;最终通过qRT-PCR对筛选出的可能影响卵泡发育的关键基因进行验证。【结果】分别对测序得到的raw reads进行过滤后,在DF颗粒细胞中获得43 217 934条clean reads,占raw reads的比例为95.19%;SF颗粒细胞中获得40 766 348条clean reads,占raw reads的比例为95.35%。将得到的unigenes与山羊的RefSeq数据库进行比对后,共得到33 896条带有注释的转录本,再通过设定FPKM1, q value0.05,共在两种卵泡颗粒细胞中获得13 644个基因。设定参数:FPKM≥1,SF-FPKM/DF-FPKM1,P0.05,获得695个差异表达mRNA,其中233个在SF颗粒细胞中表达显著上调,462个表达显著下调;对所获得695个差异表达mRNA进行GO功能富集分析,共分为三大类42组:其中生物学过程占47.6%,细胞组分占47.6%,分子功能占4.8%;KEGG信号通路分析,发现20条通路,其中与核糖体通路相关的基因富集最为显著。通过在Genecard中进行功能分析后,筛选6个可能与山羊卵泡发育密切相关的基因,其中PRLR、PTX3、RGN在SF颗粒细胞中表现为上调;DKK3、ALDH1A2、RARRES1则表现为下调。qRT-PCR显示PRLR、RGN、DKK3、ALDH1A2、RARRES1的表达趋势与高通量测序结果一致,且RGN在从属卵泡颗粒细胞中的表达量极显著地高于优势卵泡(P0.01);DKK3、ALDH1A2、RARRES1在优势卵泡颗粒细胞中的表达量极显著地高于从属卵泡(P0.01)。【结论】DKK3、ALDH1A2、RARRES1和RGN在优势卵泡和从属卵泡中表达量存在极显著差异,推测在山羊卵泡发育过程中可能促进卵泡的优势化或导致闭锁,对深入探究卵泡发育的调控机制具有重要意义。     

基于高通量测序分析肺炎支原体感染姜曲海猪肺组织miRNA表达谱

为探明姜曲海猪感染肺炎支原体后肺组织miRNA(microRNA)表达谱和分子机制,选取50日龄姜曲海猪为实验猪,随机分成感染组和对照组,人工感染肺炎支原体28 d后,采集肺部组织,进行高通量miRNA测序,采用生物信息学软件进行miRNA鉴定和靶基因预测。结果显示:感染组和对照组的肺组织分别筛选到14 265 786条和14 000 588条小RNA纯净序列(clean reads)。与对照组相比,感染组中筛选到73个显著差异表达的miRNAs,其中39个表达上调,34个表达下调,从中随机选取4个miRNAs进行定量PCR(quantitative real time PCR,qRT-PCR)验证,感染组和对照组的表达水平与测序结果基本一致。73个差异表达miRNAs预测到1 685个靶基因和4 220个靶位点,靶基因预测筛选到8个与免疫调控相关的miRNAs。靶基因GO(gene ontology)分析显示,miRNA广泛参与生物过程、细胞组成和分子功能的调控。靶基因KEGG(Kyoto encyclopedia of genes and genomes)分析显示,miRNA参与调控细胞凋亡、ECM受体相互作用、粘附斑通路等信号通路。本研究通过高通量测序获得姜曲海猪感染肺炎支原体后肺组织miRNA表达谱,通过生物信息学分析筛选到与免疫调控相关的miRNAs和信号通路,为进一步阐明姜曲海猪的肺炎支原体感染机制奠定了基础。     

育成期高能饲喂下开产与未开产蛋鸡肝脏miRNA表达谱差异分析

【目的】对育成期高能饲喂下开产与未开产蛋鸡的肝脏进行miRNA高通量测序分析,探明高能饮食状态下蛋鸡肝脏中影响其开产的miRNA,为提高产蛋性能打下基础。【方法】以代谢能水平为12.14 MJ/kg的高能饲粮饲喂育成期蛋鸡,通过Illumina NextSeq500高通量测序平台对开产与未开产蛋鸡肝脏进行small RNA测序,使用DESeq分析miRNA表达量,并以实时荧光定量PCR进行验证;采用miRanda对差异表达miRNA进行靶基因和靶位点预测,同时以超几何分布对差异表达的靶基因进行GO功能注释分析和KEGG信号通路富集分析。【结果】 6个样本（3个开产蛋鸡样本,3个未开产蛋鸡样本）能注释到的miRNA均超过300个,约占miRBase中已鉴定鸡miRNA的30.00%,且前体鉴定结果显示各样本注释到的miRNA有部分定位在同一前体上。与未开产组样本相比,开产组样本有9个miRNA表达上调、3个miRNA表达下调。其中,gga-miR-132c-5p、gga-miR-132b-5p、gga-miR-2184a-5p、gga-miR-132c-3p、gga-miR-132b-3p及gga-miR-2184a-3p属于同一miRNA基因簇,且gga-miR-1682、gga-miR-132b-3p和gga-miR-2184a-3p等3个上调miRNA在未开产蛋鸡样本中不表达。通过miRanda共预测得到129个差异表达潜在靶基因,以miR-203a预测得到的靶基因最多,为32个;其次是gga-miR-2184a-5p、gga-miR-148a-5p和gga-miR-12211-5p,分别预测到30、25和23个靶基因;而miR-132c-5p预测得到的靶基因最少,仅为9个。129个靶基因显著注释到8条GO功能条目上,生物学过程（Biological process）主要注释到脂质代谢过程、细胞脂质代谢过程、脂质生物合成过程、磷脂生物合成过程、磷脂代谢过程、甘油磷脂生物合成过程和解剖学结构发育,分子功能（Molecular funcion）仅注释到1条GO功能条目,即利钠肽受体活性,未发现涉及细胞组分（Cellular component）的GO功能条目;在注释到的8条GO功能条目中有6条与脂质代谢相关,涉及的靶基因有22个,占总潜在靶基因的17.10%。KEGG信号通路富集分析结果显示共显著富集到7条KEGG信号通路,其中脂肪酸降解、脂肪酸生物合成、酮体合成与降解及PPAR信号通路等4条信号通路与脂质代谢相关。【结论】育成期高能饲喂下开产与未开产蛋鸡中共存在12个差异表达miRNA,涉及129个差异表达潜在靶基因,且主要富集在肝脏脂质代谢相关过程和信号通路,说明miRNA是通过调控脂质代谢及其相关基因表达而影响蛋鸡开产。     

母猪健康和闭锁卵巢有腔卵泡的转录组测序比较分析

[目的]本文旨在通过对猪卵泡颗粒细胞进行转录组测序了解卵泡闭锁相关基因的表达谱并确定闭锁卵泡的潜在分子标记物。[方法]通过HE染色对卵泡及颗粒细胞进行形态学观察；抽取健康卵泡和晚期闭锁卵泡的颗粒细胞，进行转录组测序分析，最后采用实时荧光定量PCR检测差异表达基因的相对表达水平。[结果]转录组数据表明，与闭锁卵泡颗粒细胞相比，健康卵泡颗粒细胞中有771个差异基因，其中204个为上调基因，567个为下调基因；GO功能富集分析显示，差异基因显著富集到320个生物进程中，主要与细胞增殖和血管生成有关；KEGG通路分析表明，差异基因富集到48条信号通路中，主要参与氧化应激、激素合成和能量代谢等方面的信号通路；实时荧光定量PCR结果表明，CYP1B1、PDE10A、NUPR1、NOX4、LDLR、LHCGR和CYP19A1表达趋势与测序结果一致。[结论]氧化应激、血管生成、能量代谢以及激素合成等多方面因素综合调控猪卵巢中颗粒细胞的凋亡，继而影响卵泡闭锁现象的产生。     

山羊卵泡发育相关基因的筛选及分析

【背景】山羊第一卵泡波中的优势卵泡（dominant follicles, DF）和从属卵泡（subordinate follicles, SF）是整个卵泡发育过程中最为关键的两个阶段。随着卵泡的进一步发育,最终DF可能发育成为成熟卵泡,直到排卵;SF将走向闭锁,其中颗粒细胞的凋亡是导致卵泡发生闭锁的关键因素。然而目前对促进卵泡的优势化或导致其闭锁的分子机理尚不清楚。【目的】通过对山羊第一卵泡波中DF和SF颗粒细胞进行高通量测序,旨在筛选影响卵泡发育的关键基因,为深入探究卵泡发育的调控机制提供理论依据。【方法】选取10只1岁龄健康的贵州白山羊分别注射前列腺素F_2α,使其同期发情,此后每天用B超检测并记录卵泡的生长情况,发情3 d后,统一屠宰并采集第一卵泡波中DF (直径4.5—6 mm)与SF (直径3 —4.5 mm),分别分离其中的颗粒细胞,提取总RNA、构建文库后通过Illumina Hiseq 2500平台进行测序。利用FastQC对测序产出raw reads进行质量评估并经过过滤后,获得品质较高的clean reads;使用Trinity对得到的clean reads进行重新组装,从而获得unigenes;使用CLC Genomics Workbench将unigenes与山羊RefSeq数据库进行比对获得mRNA;使用DESeq2 软件对获得的mRNA进行差异表达分析;分别采用goseq和kobas软件对得到的差异表达基因进行GO分析及KEGG信号通路分析;最终通过qRT-PCR对筛选出的可能影响卵泡发育的关键基因进行验证。【结果】分别对测序得到的raw reads进行过滤后,在DF颗粒细胞中获得43 217 934条clean reads,占raw reads的比例为95.19%;SF颗粒细胞中获得40 766 348条clean reads,占raw reads的比例为95.35%。将得到的unigenes与山羊的RefSeq 数据库进行比对后,共得到33 896条带有注释的转录本,再通过设定FPKM>1, q value<0.05,共在两种卵泡颗粒细胞中获得13 644个基因。设定参数：FPKM≥1,SF-FPKM/DF-FPKM>1,P<0.05,获得695个差异表达mRNA,其中233个在SF颗粒细胞中表达显著上调,462个表达显著下调;对所获得695个差异表达mRNA进行GO功能富集分析,共分为三大类42组：其中生物学过程占47.6%,细胞组分占47.6%,分子功能占4.8%;KEGG信号通路分析,发现20条通路,其中与核糖体通路相关的基因富集最为显著。通过在Genecard中进行功能分析后,筛选6个可能与山羊卵泡发育密切相关的基因,其中PRLR、PTX3、RGN在SF颗粒细胞中表现为上调;DKK3、ALDH1A2、RARRES1则表现为下调。qRT-PCR显示PRLR、RGN、DKK3、ALDH1A2、RARRES1的表达趋势与高通量测序结果一致,且RGN在从属卵泡颗粒细胞中的表达量极显著地高于优势卵泡（P<0.01）;DKK3、ALDH1A2、RARRES1在优势卵泡颗粒细胞中的表达量极显著地高于从属卵泡（P<0.01）。【结论】DKK3、ALDH1A2、RARRES1和RGN在优势卵泡和从属卵泡中表达量存在极显著差异,推测在山羊卵泡发育过程中可能促进卵泡的优势化或导致闭锁,对深入探究卵泡发育的调控机制具有重要意义。     

基于高通量测序的山羊卵巢基质、卵泡转录组分析

【目的】比较山羊Capra hircus卵巢基质、大卵泡和小卵泡间的m RNA表达图谱,为探究卵泡发育机制提供一定的研究基础。【方法】采用高通量测序技术对发情期川中黑山羊的卵巢基质、大卵泡和小卵泡进行转录组测序,利用生物信息学方法检测mRNA表达谱,筛选差异基因,对其进行GO和KEGG通路分析,最后随机抽取5个差异基因进行荧光定量PCR验证。【结果】卵巢基质vs大卵泡、卵巢基质vs小卵泡和小卵泡vs大卵泡的差异表达基因分别为524、180和403个。筛选出INHA、TNFRSF19等15个与卵泡发育相关的基因,其主要富集在内质网蛋白质加工、类固醇生物合成、卵母细胞减数分裂等信号通路。通过q RT-PCR验证,定量结果与测序结果基本一致。【结论】川中黑山羊卵巢基质、大卵泡和小卵泡的基因表达模式不同,小卵泡与卵巢基质的基因表达模式更相近。     

持续感染口蹄疫病毒VP1蛋白对细胞转录组的影响

【目的】通过转录组测序（RNA-Seq）方法分析持续感染口蹄疫病毒结构蛋白VP1诱导的PK-15细胞基因表达谱的变化,探索VP1蛋白对宿主细胞信号通路的影响。为后续开展VP1影响FMDV持续感染相关机制研究提供基础数据和理论支持。【方法】PK-15细胞分别转染pCAGGS-HA（对照载体）、pCAGGS-HA-Pi-VP1质粒,转染36 h后,Western blot检测VP1蛋白表达,同时提取两组细胞总RNA进行转录组测序（每组3个重复,2组6个样本）。将测序所得的序列进行过滤比对,获得差异表达基因后,对差异表达基因进行了生物信息学分析。采用qRT-PCR验证关键差异表达基因的表达。【结果】RNA-Seq共鉴定出3 376个差异表达基因（DEGs）|log2(Fold Change)|>0和padj<0.05,其中上调基因1 744个,下调基因1 632个。许多参与免疫反应和炎症的效应基因差异表达。GO和KEGG富集分析显示,差异基因GO条目绝大部分与核糖体有关,而KEGG富集中涉及了与天然免疫相关的信号通路,如TNF信号通路、Toll受体信号通路。【结论】pCAGGS-...     

梅花鹿鹿茸快速生长期环状RNA鉴定及生物信息学

环状RNA(circRNA)是一类呈环状结构的内源性非编码RNA.通过对3只雄性梅花鹿鹿茸间充质组织的小鞍子时期和快速生长期进行高通量测序,根据测序结果鉴定其中的环状RNA,分析不同发育阶段环状RNA表达量差异,对差异表达的环状RNA来源基因进行GO富集分析和KEGG功能注释分类,最终构建了与鹿茸生长发育相关差异最显著的circRNA-miRNA-mRNA调控网络.结果表明:鉴定得到7790个环状RNA,其中,从鹿茸发育前期到鹿茸的快速生长期鉴定得到208个显著差异表达的环状RNA,包括39个显著上调差异表达的环状RNA以及169个显著下调差异表达的环状RNA.环状RNA显著富集到生物调控、细胞过程、代谢过程、细胞、催化活性等多种GO条目,显著富集到粘着斑、PI3K-Akt信号通路、MAPK信号通路、TNF信号通路、FoxO信号通路.最终筛选的最显著表达的10个环状RNA构建了与鹿茸发育相关的circRNA-miRNA-mRNA互作调控网络.     

禽致病性大肠杆菌phoP/Q基因缺失对雏鸡脾脏microRNAs功能的影响

从microRNA(miRNA)水平探究禽致病性大肠杆菌(APEC)phoP/Q基因缺失对雏鸡脾脏功能表达的影响,为APEC的防制提供参考资料。用禽致病性大肠杆菌和其phoP/Q基因的缺失株分别攻毒14日龄雏鸡,采集其脾脏组织为材料进行高通量测序获得miRNA表达谱,选择差异倍数大于8的miRNA进行Real-time PCR,对差异表达的miRNA进行靶基因预测以及GO、KEGG富集分析。测序结果获得10个差异表达的miRNA,其中1个miRNA表达量下调,9个miRNA表达量上调,Real-time PCR结果表明,miRNA变化趋势与测序结果一致,GO分析结果表明,差异表达基因主要富集在能量代谢、免疫系统、细胞生长与凋亡、糖合成与代谢等生物学过程。KEGG分析结果显示预测靶基因参与MAPK信号通路、糖代谢通路、Wnt信号通路等重要通路。本试验分析禽致病性大肠杆菌phoP/Q基因缺失后感染雏鸡脾脏miRNA的表达差异,从miRNA水平探讨了phoP/Q基因缺失对APEC致病性的影响,为APEC的防治提供参考资料。