山羊皮肤新microRNA的挖掘

已有证据表明,microRNA在动物皮毛发育和形态发生过程中发挥重要的作用.在山羊皮肤microRNA的研究中,数据库不完善仍是限制其功能研究的主要原因.本研究以山羊基因组为参考,对高通量测序中未注释的占总reads 5.40%的数据,使用Mireap软件预测的方法进行数据挖掘,并对新的microRNA的长度、表达量、基因组定位、靶基因预测及KEGG pathway进行分析.研究预测到28个新的microRNA,其长度分布与已知山羊、牛等动物的分布一致.对表达量统计表明,表达量在100以上的有4个,最高的是Novel-15,达到555;表达量在10~100的有19个;低于10的有5个.基因组定位研究表明,新microRNA在山羊染色体上分布均匀,在5、6、7、19和X染色体上分别有两个microRNA被定位.预测和分析得到6 519个靶基因.KEGG分析结果表明,有3 552个靶基因定位到51条通路中,占前3位的通路是Metabolic pathways、Pathways in cancer、MAPK signaling pathway.本研究发掘了更多山羊的microRNA,为今后山羊皮肤发育调控的研究奠定了基础.     

绿壳蛋鸡与白来航蛋鸡microRNA变异及其靶基因预测分析

本研究旨在分析绿壳蛋鸡与白来航蛋鸡microRNA变异及其靶基因,将已鉴定的microRNA与预测microRNA组成全基因组的microRNA,并通过与GGRS得到的两种鸡群体SNP位点进行映射,挖掘含SNP的microRNA。利用生物信息学的方法,对含SNP的microRNA进行靶基因预测,并对直接含有SNP的mature-microRNA进行聚焦。将这些靶基因进行富集,一共得到22个GO分类和10条KEGG信号通路与3个主要IPA调控网络,发现其在与生长相关的mTOR信号通路、Wnt信号通路、生长激素受体网络、胰岛素样生长因子Ⅰ受体网络得到了富集,与产蛋性状相关的卵母细胞减数分裂信号通路、黄体酮卵母细胞成熟信号通路得到了富集。此研究方法和结果可以给后期研究提供参考。     

热应激下齐兴肉兔睾丸组织的转录组分析

为研究公兔在急性热应激过程中相关基因的表达及响应的分子机制,本研究以齐兴肉兔公兔为实验对象,构建热应激公兔睾丸精子发生模型,提取热应激组和对照组睾丸组织总RNA,反转录建立cDNA文库,利用llumina Hiseq^TM测序平台进行转录组测序,对差异表达基因进行GO、KEGG富集分析及RT-PCR验证。结果表明:与对照组相比,热应激组已注释的有765个基因表达上调,37个基因表达下调。对差异显著基因进行KEGG通路富集分析,主要富集于胞吞作用、PI3K-Akt信号通路、核糖体、细胞黏附分子(CAMs)、MAPK信号通路等通路,最终筛选出3个与热应激反应、精子生成有关的基因进行RT-PCR检验,为精子发生的分子调控机制提供线索,并为耐热性家兔的育种提供理论基础。     

禽致病性大肠杆菌(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致病机制提供新思路。     

高能饲粮对育成期蛋鸡肝脏miRNA表达谱的影响

本研究利用高能和低能饲粮饲喂育成期蛋鸡,运用small RNA测序方法对肝脏中miR-NA进行检测,对其表达量进行分析,对差异表达miRNA进行靶基因预测,并对差异表达的预测靶基因进行基因本体论(GO)功能和京都基因与基因组百科全书(KEGG)信号通路富集分析.结果表明:与低能组相比,高能组有2个miRNA表达显著上调...     

小尾寒羊间情期和发情期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是通过代谢和免疫途径调控卵巢周期性活动。     

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

牦牛miR-101的靶基因预测及生物信息学分析

本研究利用相关生物信息学软件和数据库对牦牛miR-101进行了靶基因预测与生物信息学分析。通过高通量测序获取牦牛miR-101的序列,并分析其序列保守性;选取TargetScan、miRanda、PicTar预测结果的交集并结合miRTarbase中已证实的靶基因集合,取二者并集,采用BINGO和DAVID数据库获得靶基因集合的本体注释,进行GO功能富集及Pathway信号通路富集分析。结果显示,miR-101序列在各物种间高度保守,miR-101调控的靶基因GO功能富集主要包括多细胞生物体发育、发育过程、系统发育、解剖结构发育、器官发育等基本生物学过程和蛋白连接、转录因子活性等分子功能;靶基因存在于细胞各个组分中,包括细胞质、细胞核、细胞器中;信号转导通路显著富集于MAPK信号通路（MAPK signaling pathway）、黏着斑（focal adhesion）、ErbB信号通路（ErbB signaling pathway）、Wnt信号通路（Wnt signaling pathway）、TGF-beta信号通路（TGF-beta signaling pathway）和mTOR信号通路（mTOR signaling pathway）中（P<0.05）。通过对牦牛miR-101靶基因的生物信息学分析,为进一步研究miR-101在牦牛肌肉发育中的作用奠定基础。     

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

基于转录组测序数据的山羊下丘脑组织繁殖相关信号通路及候选基因的筛选

旨在通过对不同光照时长条件下河南槐山羊下丘脑组织进行转录组测序和生物信息学分析,筛选出与山羊繁殖相关的信号通路及候选基因。本研究利用Illumina二代高通量测序平台(NGST),采用PE150测序策略,对经自然光照(8 h光照:16 h黑暗)与人工光照(16 h光照:8 h黑暗)条件处理后的8只空怀母羊(每组各4只,平均年龄1周岁)下丘脑组织进行转录组测序,将组装得到的Unigene总数比对到参考基因组序列,进行差异表达基因分析、GO和KEGG富集分析以及候选基因筛选,并采用Real-time PCR方法对候选基因的相对表达量变化进行分析。结果显示,RNA-Seq共得到了约4.4亿条reads,平均每个样本的reads数约为5 249万条;自然光照与人工光照2组DESeq分析得到448个差异表达基因,并富集在KEGG数据库中的241个信号通路,其中包括5个与繁殖相关的信号通路;3个速激肽家族候选基因(TACR1、TACR2和TACR3)显著富集在Calcium signaling pathway(钙离子信号通路,chx04020)通路;Real-time PCR分析结果显示,转录组测序结果可靠。综上表明,Calcium signaling pathway(钙离子信号通路,chx04020)与TACR1、TACR2和TACR3基因,可能在山羊繁殖过程中发挥着重要作用。     

调控香猪繁殖候选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,可为从分子水平提高香猪产仔数提供理论基础。     

蒙古马高负荷运动训练前后转录组差异表达分析

旨在对蒙古马高负荷运动训练前后转录组差异表达进行分析。本试验对6匹进行了4个月高负荷运动训练的蒙古马进行研究,分别在训练前后两个时期采集肌肉样品,利用二代测序技术对肌肉样品建立转录组文库,对发生显著变化的差异基因进行GO功能富集和KEGG Pathway分析,寻找训练前后发生显著变化的代谢通路及相关候选基因。结果表明,在转录组筛选到的1 102个显著差异表达基因中,有299个在训练后发生下调,803个发生上调。这些差异基因被注释到3 398条GO term上。其中与运动学相关的生物学过程有:肌肉结构发育、心血管系统发育、循环系统发育、肌肉组织发育、肌肉器官发育、肌肉收缩和肌肉细胞分化等。对差异表达基因进行KEGG pathway分析,结果表明,差异基因被富集到如扩张型心肌病、肥厚型心肌病、心肌收缩、钙信号通路、肌动蛋白骨架调节等与运动相关的193条通路中。本研究结果为探究蒙古马强耐力的基本分子机理提供理论基础,同时也为理解蒙古马运动学特性提供新的思路。     

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

本研究旨在分析产肠毒素大肠杆菌（enterotoxigenic Escherichia coli，ETEC）感染猪小肠上皮细胞（IPEC-J2）后，在感染初期细胞内长链非编码RNA （long non coding RNA，lncRNA）的表达谱变化，探究lncRNA在ETEC感染初期所起到的调控作用。使用ETEC F41感染IPEC-J2，在感染前和感染初期（感染后4 h）收集细胞，通过Illumina Hiseq Xten平台进行高通量测序，共发现lncRNA 9 975条。感染初期与感染前相比，共发现100条差异表达lncRNA，其中40条表达上调，60条表达下调。通过miRanda-3.3a和psRobot_v1.2软件共同预测差异表达lncRNA的靶基因，并对差异表达lncRNA的靶基因进行GO功能和KEGG通路富集分析。结果显示，感染初期差异表达lncRNA的靶基因显著富集于细胞核、核仁、代谢过程调控及发育过程等GO功能条目中。KEGG分析表明，感染初期差异表达lncRNA的靶基因主要富集在PI3K-Akt信号通路、肌动蛋白细胞骨架调节、黏附斑及细胞周期等信号通路。利用实时荧光定量PCR随机验证了5条差异表达lncRNA，结果与测序分析中表达变化趋势一致。本研究对ETEC感染IPEC-J2细胞初期的lncRNA表达谱进行了差异分析，为深入探究lncRNA在ETEC感染初期中的作用机制提供了参考依据。     

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

兔感染日本血吸虫后血浆中宿主源循环microRNAs变化研究

循环性microRNAs不仅参与生命体的重要活动调控,还有望成为某些疾病诊断的新型标志分子。为深入探讨宿主与血吸虫之间的相互作用机制并发现新型诊断标志分子,我们对前期获得的感染和未感染新西兰大白兔血浆小RNA的Solexia测序结果进行生物信息学分析,对差异表达的microRNAs用实时定量PCR进行验证,同时对这些差异表达的microRNAs的靶基因进行了分析。在剔除虫源性序列后,我们发现了感染和未感染日本血吸虫的新西兰大白兔体内的循环性宿主源microRNAs 173条,选取了10个差异表达的microRNAs利用实时定量PCR进行验证分析,结果表明所选取的10个microRNAs表达与Solex测序结果一致。靶基因的预测与分析结果表明感染血吸虫兔血浆中差异表达microRNAs的靶基因主要参与细胞进程、代谢调控、生命调控、发育进程等生命进程调控。本研究结果为深入开展日本血吸虫与宿主的互作及其寄生的分子机理提供了基础。     

静原鸡miR-107的靶基因预测及生物信息学分析

为探究miR-107在静原鸡肌肉组织肌苷酸特异性沉积过程中的调控机制,本研究利用相关生物信息学软件和数据库对静原鸡miR-107进行了靶基因预测与生物信息学分析。通过转录组测序获取静原鸡miR-107序列,进行序列比对和保守性分析;通过TargetScan、miRDB和PicTar预测结果的交集并结合miRTarbase中已证实的靶基因集合,取二者并集,采用OmicShare工具对基因集合进行GO功能富集和信号通路富集分析,并讨论HMED数据库中miR-107在不同组织和疾病中的表达水平。结果显示,miR-107序列在各物种间高度保守,miR-107调控的靶基因GO功能富集主要包括细胞过程、生物调节、代谢过程、刺激反应等生物学过程和结合、催化活性、转录调节活性等分子功能;靶基因存在于细胞的多个组分中,包括细胞器、细胞膜、含蛋白质复合物。信号转导通路显著富集于癌症中的microRNAs途径（microRNAs in cancer）、人乳头瘤病毒感染（human papillomavirus infection）、癌症中途径（pathways in cancer）、PI3K-Akt信号通路（PI3K-Akt signaling pathway）中（P<0.05）。miRNA-靶基因互作分析发现,miR-107主要通过与靶基因mRNA的非翻译区靶向结合抑制靶基因的表达。根据HMED数据库中miR-107在不同组织和疾病中的表达数据,发现miR-107在黑素瘤（melanoma）中的表达水平最高,在扁桃体（tonsil）和肌肉（muscle）等组织中表达水平较低。通过对静原鸡miR-107靶基因的生物信息学分析,为miR-107功能及调控机制的深入研究提供参考依据,也为进一步研究miR-107在静原鸡肌肉发育中的作用奠定基础。     

Ssc-miR-133a-5p的靶基因预测及其相关信号通路的生物信息学分析

课题组前期通过高通量测序技术发现,microRNAssc-miR-133a-5p在肥胖猪脂肪组织中表达水平显著下调,推测其可能在脂肪沉积过程中发挥了重要作用。分析ssc-miR-133a-5p序列保守性,利用TargetScan,miRDB和miRWalk在线分析工具预测其候选靶基因,进一步对候选靶基因进行蛋白质互作分析以及KEGG分析,最后将预测的靶基因与课题组前期筛选出的与猪脂肪沉积能力相关的基因取交集。结果表明ssc-miR-133a-5p候选靶基因PPARGC1A与RORA等,参与AMPK,TNF,与胰岛素分泌等信号通路。结果提示ssc-miR-133a-5p的候选靶基因可能通过多种信号通路发挥重要作用。文章的结果可为microRNA参与猪的脂肪生成调节机制提供科学依据。     

初情期母牛下丘脑-垂体-卵巢转录组分析

通过对3头18月龄安格斯牛下丘脑-垂体-卵巢组织进行转录组测序和生物信息学比较分析,筛选和挖掘与肉牛繁殖活动相关候选基因和信号通路。通过转录组测序,进行表达基因分析、GO和KEGG富集分析以及信号通路筛选。结果显示,共获得61.92 Gb Clean Data,各样品Clean Data均达到6.88 Gb,Q30碱基百分比在92.42%及以上。有15154个基因在所有组织中都有表达,有4780个基因只在2种组织器官中共同表达,有5809个基因只在其中的一个组织器官中表达。下丘脑-垂体-卵巢各组织所表达的基因中,在GO数据库得到注释的分别为2168,2101,13777个,在KEGG数据库中得到注释分别有12477,12446,12176个。在所有KEGG通路中,Notch信号通路(notch signaling pathway)、神经营养因子信号通路(neurotrophin signaling pathway)和血管内皮生长因子信号通路(VEGF signaling pathway)的表达频率最高,说明这些信号通路可能也参与下丘脑-垂体-卵巢轴对母牛繁殖活动的调节。     

美国短毛黑水貂快速生长过程中转录组差异表达分析

旨在对美国短毛黑水貂快速生长过程中转录组差异表达进行分析。本试验将美国短毛黑水貂45、90日龄各3只健康公貂的胸肌组织作为试验材料,利用Illumina HiSeq~(TM)2500高通量测序平台建立转录组文库,对两个生长阶段差异显著性的基因进行GO功能富集和KEGG Pathway分析,寻找调控水貂快速生长期的相关候选基因和代谢通路,并利用qRT-PCR对转录组测序结果进行验证。结果表明,以P0.05,|log_2FC|1为条件筛选到279个显著差异基因,对这些差异基因进行GO和KEGG分析,筛选到与肌肉生长相关显著富集GO条目有66条、相关差异表达基因44个,其中上调20个,下调24个;与肌肉生长相关显著富集KEGG通路12条,如p53信号通路、PPAR信号通路、FOXO信号通路等。这些差异基因可能在水貂的快速生长过程中起到调控作用,可以作为后续研究水貂生长发育的候选基因。本研究结果为探究水貂生长发育调控机制及培育大体型水貂品种奠定基础。     

不同氧浓度下猪睾丸间质细胞转录组分析

【目的】研究不同氧浓度下猪睾丸间质细胞表达谱差异,筛选缺氧导致猪睾丸间质细胞损伤的关键基因。【方法】将猪睾丸间质细胞分为缺氧组(1%氧气,H24组)、正常组(15.75%氧气,N24组),每组3个重复,在相应条件下分别处理24 h,用CCK8法检测细胞存活率;利用转录组测序技术筛选出缺氧组和正常组睾丸间质细胞中的差异表达基因,利用DESeq 2.0软件进行差异表达基因分析,并对差异表达基因进行GO功能和KEGG通路富集分析,进一步筛选与氧含量相关的调控基因。【结果】缺氧刺激24 h后,与正常组相比,缺氧组猪睾丸间质细胞存活率极显著降低(P<0.01);转录组测序共获得1 654个差异表达基因,其中1 242个基因上调,412个基因下调。GO功能富集分析发现,共获得富集条目8条,其中生物过程1条,分子功能7条;KEGG通路富集分析发现,共获得10个显著富集的信号通路,包括HIF-1信号通路、JAK-STAT信号通路和糖尿病并发症中的AGE-RAGE信号通路等。通过GO功能和KEGG通路富集分析共获得6个与缺氧应激相关的基因,分别为信号传导及转录激活蛋白5(signal transd...