福清山羊和戴云山羊夏季发情期卵巢组织比较转录组分析

试验旨在从分子遗传学角度探讨戴云山羊夏季维持较高繁殖力的遗传机制,从中发掘有用的遗传变异信息。利用转录组测序方法对福清山羊和戴云山羊夏季发情期卵巢组织进行研究,筛选品种间的差异表达基因（differentially expressed genes,DEGs）,利用基因本体（gene ontology,GO）、蛋白相邻类的聚簇（cluster of orthologous groups of proteins,COG）和京都基因与基因组百科全书（kyoto encyclopedia of genes and genomes,KEGG）数据库对筛选的差异表基因进行功能注释,并进行差异表达基因聚类分析,通过实时荧光定量PCR验证测序结果。结果显示,福清山羊和戴云山羊发情期卵巢的DEGs共357个,其中上调基因221个,下调基因136个。357个DEGs中的287个基因能够被GO数据库注释,78个DEGs能够被COG数据库注释,219个DEGs能够被KEGG数据库注释。KEGG数据库注释结果显示,这些差异基因显著富集的信号通路为蛋白消化吸收（protein digestion and absorption）、吞噬体（phagosome）、肺结核（tuberculosis）、胞外基质受体互作（ECM-receptor interaction）、金黄色葡萄球菌感染（Staphylococcus aureus infection）、同种异体移植物排斥（allograft rejection）等通路。经实时荧光定量PCR验证,所选基因（转录本）表达变化模式与转录组测序结果一致,表明测序结果可靠。本试验利用高通量测序获得大量山羊卵巢转录组信息,有助于从分子水平对山羊卵巢进行深入研究。     

基于转录组测序筛选高盐培养下双峰驼肾皮质细胞差异表达基因

为探讨双峰驼的高盐适应机制并筛选出受高盐调控的基因，本试验采用RNA-Seq技术对双峰驼肾皮质细胞进行转录组测序分析。试验分为2个组，等渗培养基处理的肾皮质细胞为对照组（Control），高渗培养基处理的肾皮质细胞为高渗组（HS），进行转录组测序，从而筛选出一系列差异表达基因（DEGs），并用GO富集和KEGG富集分析对DEGs进行基因功能和途径的注释。结果显示，在高盐环境下双峰驼肾皮质细胞中共有4 854个显著DEGs;GO富集分析显示，DEGs显著富集在G蛋白偶联受体活性、受体结合、核酸结合转录因子活性、质膜、离子通道的活动和免疫反应等功能中；KEGG通路富集显示，DEGs显著富集在信号传导、辅助因子和维生素代谢、信号分子相互作用、运输和分解代谢、免疫系统等通路中；采用实时荧光定量PCR技术随机检测了参与高盐代谢的3个差异表达基因，其表达趋势与RNA-Seq一致，可以验证RNA-Seq技术的准确性。因此，双峰驼肾皮质细胞的高盐耐受性可能与这些途径和通路中重要基因的调控有关。本试验结果将为进一步揭示双峰驼耐盐性的分子调控机制提供重要的参考依据。     

基于高通量测序的水貂胚胎滞育期和激活期卵巢转录组分析

试验旨在获取胚胎滞育期与激活期水貂卵巢组织的转录组信息，挖掘滞育的胚胎激活前后水貂卵巢差异表达基因（DEGs）及其功能信息，为探讨卵巢信号调控水貂胚胎滞育的分子机制提供参考。随机采集8只健康雌性水貂（胚胎滞育期和激活期各4只）的卵巢组织为样本，利用Illumina HiSeq平台RNA-Seq技术对其进行转录组测序，筛选在水貂胚胎滞育期和激活期卵巢中的DEGs并进行生物信息学分析。结果显示，测序获得661 195 568个raw reads，经过滤后获得650 834 900个clean reads，组装后得到389 895个unigenes，通过与Nr、GO、KOG和KEGG数据库比对，对unigenes进行注释，其中Nr数据库注释了156 419个unigenes；GO数据库注释了122 657个unigenes；KOG数据库注释了58 320个unigenes；KEGG数据库注释了72 653个unigenes。滞育期和激活期卵巢中有1 797个DEGs，与胚胎滞育期水貂卵巢相比，激活期卵巢有1 298个DEGs显著上调，499个DEGs显著下调。GO功能分析发现，DEGs显著富集的生物学过程主要有跨膜信号受体活性、信号受体活性、G蛋白偶联受体活性、酶联受体蛋白信号通路、细胞表面受体信号通路、细胞周期阻滞、芳香酯酶活性。KEGG通路分析发现，水貂滞育期和激活期卵巢中的DEGs显著富集于神经活动配体-受体相互作用信号通路。本研究利用高通量测序技术获得水貂胚胎滞育期与胚胎激活期卵巢的转录组信息，为深入探究水貂卵巢调节胚胎滞育的分子机制奠定了基础。     

基于RNA-Seq技术的塔里木马鹿毛色相关差异表达基因筛选及初步分析

旨在基于RNA-Seq技术对塔里木马鹿毛色相关基因进行筛选及分析。采用Illumina Hi Seq TM2000测序平台对塔里木马鹿和天山马鹿的皮肤组织进行转录组测序,所得序列经质控、组装后比对到NR、Swiss-Prot、COG、KOG、KEGG、GO和Pfam数据库中注释,并对差异表达基因进行筛选、功能注释和富集分析。结果表明,测序获得25 038个有注释信息的Unigenes,比对分析显示,塔里木马鹿与天山马鹿有922个差异表达基因,其中上调表达基因495个,下调表达基因427个;GO功能富集分析结果显示,568个差异表达基因富集到61个GO条目上,分别参与了生物学过程、细胞组分及分子功能;KEGG代谢通路富集分析发现,在差异表达基因中富集最显著的代谢通路是ECM-受体相互作用。利用实时荧光定量PCR（qRT-PCR）方法分析与塔里木马鹿毛色相关的7个候选基因的转录水平变化来验证转录组测序结果的准确性和可靠性,这些基因的表达趋势与转录组测序结果相一致。ECM-受体相互作用、蛋白质消化与吸收、PI3K-Akt信号通路及与黑色素合成相关的酪氨酸等通路可能与塔里木马鹿的毛色有关;候选基因MITF、Ggt1、VDR、PTPRF、CⅡTA、ARPC5L、POMC等可能在塔里木马鹿毛色形成过程中发挥重要作用。本研究结果为今后塔里木马鹿毛色相关基因的分子调控机制方面及挖掘潜在的新基因提供了丰富的试验数据。     

基于转录组测序挖掘广西麻鸡生长发育相关基因

【目的】通过对1和60日龄广西麻鸡腿肌进行转录组测序(RNA-Seq),筛选出广西麻鸡生长发育关键基因。【方法】选择1和60日龄健康母鸡各3只,分别采集腿肌组织,利用Illumina Novaseq^TM 6000平台进行mRNA转录组测序,利用DESeq2软件进行差异表达基因的分析,并对差异表达基因进行GO功能和KEGG通路富集分析,进一步筛选出生长发育相关基因,并用实时荧光定量PCR验证转录组测序结果。【结果】转录组测序结果显示,6个样本分别获得34 296 198～41 647 642条clean reads。与1日龄腿肌相比,60日龄腿肌共获得2 304个差异表达基因,其中998个基因上调,1 306个基因下调。GO功能富集分析发现,共获得富集条目572条,其中生物过程381条,细胞组分70条,分子功能121条。KEGG通路富集分析发现,共获得34个显著富集的信号通路,其中心肌收缩、MAPK信号通路、紧密连接、肌动蛋白细胞骨架信号通路、心肌细胞中的肾上腺素能信号传导等与生长发育相关。通过GO功能和KEGG通路富集分析获得5个生长发育相关基因,分别为肌球蛋白重...     

基于单细胞转录组测序技术筛选牛体内囊胚性别特异性基因

【目的】探究牛体内雌性和雄性囊胚在mRNA水平上的特异性差异,挖掘雌性和雄性囊胚发育差异相关候选基因。【方法】参考GenBank中牛牙釉质基因（AMEL）序列（AMELX基因,登录号：NM_001014984.1;AMELY基因,登录号：NM_174240.2）设计引物,以胚胎内细胞DNA为模板,采用巢式PCR扩增技术鉴定单个牛体内囊胚的性别。选用确定性别的单个雌性和雄性囊胚为试验材料,采用Smart-Seq2扩增技术构建单个胚胎测序文库,应用Illumina HiSeq Xten高通量测序平台对单个胚胎进行单细胞转录组测序（scRNA-Seq）,并进行了基因差异表达、GO功能和KEGG通路分析。【结果】通过巢式PCR扩增胚胎内细胞中的AMEL基因,确定了胚胎性别;基因表达分析筛选出两组之间差异表达基因（DEGs）6 160个,其中雌性特异性基因675个,雄性特异性基因305个;GO功能注释发现雌性特异性基因显著注释在核苷酸结合、信号转导调控、多细胞生物发育、细胞骨架等条目,雄性特异性基因显著注释在氧化磷酸化、线粒体、线粒体内膜、核糖体等条目;KEGG通路富集分析发现7条与雌性和雄性囊胚发育差异相关的通路,分别为代谢途径、糖酵解、磷酸戊糖途径、细胞衰老、氧化磷酸化、调节干细胞多能性的信号通路和Wnt信号通路;并利用GO和KEGG结果筛选出5个在牛体内囊胚期胚胎发育过程中具有直接或间接作用的基因：FBP1、GADD45G、FHL2、FOSB和WNT2B。【结论】牛体内雌性和雄性囊胚之间存在广泛的转录差异,性别特异性基因FBP1、GADD45G、FHL2、FOSB和WNT2B可能是直接或间接影响胚胎发育的关键基因。     

鞭毛蛋白诱导表达差异基因的分析

试验旨在分析鞭毛蛋白刺激机体基因表达谱的特征，揭示鞭毛蛋白佐剂的作用机制。从GEO数据库筛选目标微阵列GSE46421和GSE72016及其差异表达基因（DEGs），经DAVID分析DEGs参与的GO功能注释和KEGG信号通路、String构建DEGs的蛋白质-蛋白质相互作用（PPI）网络、Cytoscape可视化PPI网络并筛选高连通度的PPI子模块和hub基因。结果显示，共筛选出182个DEGs，包含上调基因161个、下调基因121个。DEGs参与GO生物过程主要有信号转导、免疫应答、炎症反应、细胞凋亡等，参与分子功能包含激酶激活、受体结合、细胞因子激活、膜信号转导。DEGs参与KEGG通路涵盖细胞因子-细胞因子受体相互作用、TNF信号通路、NF-κB信号通路、Toll样受体信号通路、NOD样受体信号通路等免疫相关途径。分析DEGs的PPI网络发现，2个重要的PPI子模块和10个关键基因（IL-10、IL-6、CCL2、CXCL2、NFKBIA、MyD88等）。研究表明，鞭毛蛋白通过识别TLR5和NOD样受体，触发MyD88依赖性和MyD88非依赖性途径发出信号，激活转录因子NF-...     

不同鸡种胚胎腿肌差异表达mRNA和lncRNA的筛选及其竞争性调控网络的构建

旨在筛选藏鸡和大恒肉鸡胚胎期骨骼肌差异表达的mRNA和lncRNA,并构建lncRNA-miRNA-mRNA竞争性调控网络,为进一步探讨两个鸡品种骨骼肌生长发育速度的差异机制奠定基础。随机选取已孵化18 d的藏鸡和大恒肉鸡胚胎各3个,采集胚胎腿肌组织进行转录组测序。筛选两个品种中差异表达的mRNAs和lncRNAs,对差异表达基因进行GO和KEGG富集分析,并对lncRNA的靶标miRNA以及靶向mRNA的miRNA进行预测,筛选与肌肉发育相关的mRNA和lncRNA构建lncRNA-miRNA-mRNA竞争性调控网络,最后利用实时荧光定量PCR (qRT-PCR)对测序结果进行验证。结果显示,共有106个mRNAs在两个品种中差异表达,其中上调表达mRNAs 48个,下调表达mRNAs 58个。差异表达lncRNAs共有28个,其中10个上调,18个下调。差异表达mRNA的GO富集结果显示,骨骼肌细胞分化条目被显著富集,KEGG富集分析中脂质代谢通路被显著富集。lncRNA靶基因的KEGG富集结果显示,在10个显著富集通路中,有基因显著富集于类固醇生物合成和脂肪酸生物合成等与脂质代谢...     

外源注射miR-499-5p对肉鸡腓肠肌转录组的影响

旨在探讨外源注射miR-499-5p对肉鸡腓肠肌肌纤维直径及基因表达的影响，完善miR-499-5p对肉鸡骨骼肌发育的分子调控机制。本研究选用14日龄健康大恒肉鸡公鸡16只，随机分为2组，每组8个重复，均正常饲喂日粮，分别在腓肠肌部位肌肉注射agomiR-NC（对照组，5 nmol）和agomiR-miR-499-5p（注射组，5 nmol），为保证试验效果，每3 d注射1次，共注射6次。注射完成后将16只肉鸡屠宰，采集腓肠肌组织进行肌纤维直径测定和转录组学（RNA-seq）分析。结果显示，注射miR-499-5p可导致肉鸡腓肠肌肌纤维直径显著减小（P<0.05）。转录组测序结果显示，注射组和对照组分别获得308 263 164和316 696 152条clean reads，与鸡参考基因组的比对率分别为90.29%和89.57%。与对照组相比，注射组共有差异表达基因18个，其中上调基因12个，下调基因6个（P<0.05）。其中差异表达最大的两个基因分别是原癌基因FOS（fos proto oncogene，FOS）和早期生长应答蛋白基因（EGR1）。GO富集分析发现，差异表达基因显著富集于跨膜受体蛋白丝氨酸/苏氨酸激酶信号通路、转录调控、脂质合成和骨骼肌再生等过程。KEGG富集分析结果显示，差异表达基因显著富集于丝裂原活化蛋白激酶（MAPK）信号通路和Wnt信号通路。利用生物信息学软件预测到miR-499-5p的靶基因共78个，GO富集分析发现，靶基因也显著富集于丝氨酸/苏氨酸激酶活性的负调控和转录因子活性等过程。本研究表明，注射miR-499-5p可改变肉鸡腓肠肌肌纤维直径从而影响肌纤维特性，差异表达基因中的FOS、EGR1、CYR61和WNT7A基因可能在miR-499-5p的作用下通过上述通路参与了脂肪沉积和肌肉发育的调控。     

共轭亚油酸影响边鸡胸肌脂类代谢的相关差异基因的筛选

旨在通过转录组测序技术筛选出共轭亚油酸(conjugated linoleic acids,CLA)影响边鸡胸肌脂类代谢的相关基因,为日粮中CLA调控边鸡脂类代谢的作用机制提供理论依据。本试验选取90日龄健康边鸡180只,随机分为5组,每组3个重复,每组日粮中分别添加不同比例的CLA 0%(对照组)、0.5%、1.0%、1.5%和2.0%,预饲期1周,正饲期6周。饲养试验结束后采集胸肌组织进行转录组测序,使用DESeq2软件对测序数据进行差异表达分析,对差异表达基因进行GO功能和KEGG通路富集分析,筛选出与胸肌脂类代谢相关的差异表达基因,利用qPCR对差异表达基因进行验证。通过转录组测序数据分析,共筛选出1 229个差异表达基因,其中594个基因上调,635个基因下调。随机选取9个差异表达基因进行qPCR验证,其相对表达量变化趋势与测序结果一致。GO功能分析发现,差异表达基因主要集中在生物学过程中的细胞过程、单一生物过程、生物调节和代谢过程。KEGG通路富集分析发现,各试验组中差异表达基因显著富集在不同的信号通路中,主要富集在ECM-受体相互作用、黏着斑、吞噬体、细胞凋亡、肌动蛋白细胞骨架调节和细胞因子-细胞因子受体相互作用等通路中。通过GO和KEGG功能注释分析,共筛选出MCAT、APOA1、PTGDS、ALDH3A2、PLTP、FABP3、FOXO1、UCP3和FABP4等18个主要的参与脂类代谢相关的差异表达基因,可能在调控边鸡胸肌脂类代谢中发挥重要的作用。本研究通过转录组测序筛选出CLA影响边鸡胸肌脂类代谢过程中的主要调控基因,发现了日粮中CLA影响差异表达基因作用的主要生物学过程和信号通路,为今后研究CLA影响边鸡胸肌脂类代谢的分子作用机制奠定了基础。     

藏鸡不同发育阶段腿部肌肉组织转录组及microRNA联合分析

旨在从转录组和miRNA角度探讨藏鸡肌肉发育的调控机制,了解藏鸡肌肉发育的特殊性。本研究对120和150日龄藏鸡的肌肉组织进行转录组和small RNA测序,筛选出两个日龄阶段差异表达的基因和miRNA,并利用qRT-PCR技术对测序结果进行验证。结果,共筛选出1 691个差异表达基因,其中上调基因330个,下调基因1 361个。差异表达miRNA共有22个,其中9个上调,13个下调。随机选择的5个基因和miRNAs的qPCR验证结果表明表达趋势与测序结果一致。GO富集分析显示,在富集前10的条目中,与免疫相关的条目占较大比例。KEGG分析结果显示,在19个显著富集通路中,与免疫相关的通路占较大比例,且83个miRNA的靶基因出现在显著富集通路上。转录组和small RNA测序数据联合分析表明,343个miRNA-mRNA对为负相关调控模式。本研究从多组学层面为进一步理解藏鸡的肌肉发育提供了理论基础。     

解淀粉芽孢杆菌DGL1促燕麦生长分子机制及代谢通路探究

分离自青海大格勒干旱沙地的解淀粉芽孢杆菌DGL1(Bacillus.amyloliquefaciens)被证实能够促进青海本地燕麦(Avena sativa)品种‘青燕1号’的生长。为了揭示DGL1促进燕麦生长的分子机制,本研究通过Illumina高通量转录组测序技术分析燕麦根部与菌株DGL1菌悬液分别互作2h,4h,8h,12h的处理组与对照组(CK)间的转录组差异表达基因,并通过GO富集、KEGG Pathway富集分析燕麦地上部分差异表达基因的相关代谢通路及关键基因。结果表明：燕麦根部与DGL1菌悬液互作2h,4h,8h和12h,燕麦叶部分别有1894,6130,8033和12215个差异表达基因,显著富集在氨基酸代谢、植物光合作用、次级产物代谢通路和与燕麦生长发育调控等相关代谢途径;通过KEGG富集分析发现IAA合成的色氨酸代谢途径中生长素早期响应蛋白编码基因AUXI、茉莉酸信号途径中核心受体编码基因COI1,铵转运蛋白编码基因AMT等基因显著上调,推测菌株DGL1对燕麦的促生机制是通过促进燕麦光合作用、激素代谢、次生代谢物合成、氨基酸代谢等多个途径相互协调的结果。     

基于转录组测序技术研究饲粮赖氨酸水平调控贵州黄鸡采食量的差异表达基因和信号通路

本试验旨在利用转录组测序(RNA?Seq)技术研究饲粮赖氨酸(Lys)水平调控贵州黄鸡采食量的分子机制,挖掘与采食量相关的候选基因和信号通路,进一步提高家禽采食量以及促进养殖效益.选取健康的1日龄贵州黄鸡雏鸡400羽,随机分为4组,每组5个重复,每个重复20羽.各组饲粮Lys水平分别为0.7％、0.9％、1.1％和1....     

大型迪庆藏猪不同生长阶段背脂与腹脂脂质代谢差异基因及调控网络分析

旨在筛选大型迪庆藏猪不同生长阶段、不同部位脂质代谢差异的关键基因。本研究选择胎次相同、出生日期相近、体重10 kg左右的大型迪庆藏猪36头,随机分为3组,相同条件育肥,分别在平均体重达40、80和120 kg时屠宰,测定胴体性能,每组采集3头猪的背脂和腹脂进行高通量转录组测序,测序数据经拼接、比对,筛选与脂质代谢相关的显著差异基因并进行GO、KEGG分析、基因互作网络分析。结果表明,40、80和120 kg大型迪庆藏猪腹脂vs.背脂分别筛到486、765和339个差异表达显著基因,随机挑选的EGR2、SOD3等5个差异表达显著基因的qPCR结果与转录组测序结果一致,差异表达显著基因主要富集在肌肉收缩、细胞黏附、间充质细胞增殖正调节等GO条目,心肌收缩、PI3K-Akt信号通路、Hippo信号通路等KEGG通路;40 kg组EGR2、RARRES2、TMOD4和SFRP2基因位于网络核心,EGR2、RARRES2、SFRP2在腹脂上调,TMOD4下调;80 kg组THBS1、PPARA、NRIP1和LPL基因位于网络核心,4个核心基因均在腹脂上调;120 kg组HTRA1、TSHR、LR...     

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

通过对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)的表达频率最高,说明这些信号通路可能也参与下丘脑-垂体-卵巢轴对母牛繁殖活动的调节。     

C型产气荚膜梭菌外毒素致小鼠肠道损伤的转录组分析

旨在探索产气荚膜梭菌外毒素造成机体炎性损伤及免疫调控紊乱的毒性机制,为产气荚膜梭菌病的致病机理研究提供理论基础。将C型产气荚膜梭菌强毒株C59-2培养上清腹腔注射BALB/c小鼠,采集小肠样本进行转录组测序,筛选差异表达基因,并对其进行GO功能注释和KEGG通路富集分析。结果显示,共获得40.99 Gb有效碱基,筛选后共得到795个差异表达基因,其中,229个基因表达上调,566个基因表达下调,对随机选取的10个基因进行q-PCR验证,其相对表达量与转录表达谱一致。GO功能注释主要涉及G蛋白偶联核苷酸受体活性和G蛋白偶联嘌呤核苷酸受体活性等。KEGG通路富集分析发现,主要富集在TNF信号通路、IL-17信号通路、p53信号通路、FOXO信号通路、Toll样受体信号通路、NF-κB信号通路等。产气荚膜梭菌外泌毒素侵入机体,会激活TNF等炎性信号通路,进而造成肠道发生炎性损伤甚至坏死。