牛多个组织microRNASolexa测序与生物信息学分析

本研究旨在发掘牛更多的microRNA(miRNA)信息及新的miRNA序列,为进一步研究miRNA的生物学功能奠定理论基础.运用高通量测序技术对来自一头西门塔尔公牛和一头荷斯坦母牛的多个组织的小RNA进行混池测序,对测序结果进行了生物信息学分析.随机选取了一条成熟miRNA:bta-miR-2346-5p和两条新预测的miRNAs:novel_miR-48和novel_miR-86,采用茎环RT-PCR进行验证.共鉴定了604条miRNAs,其中429条为牛的已知miRNAs,175条为新预测的miRNAs.通过茎环RT-PCR发现选取的这3条miRNAs在背最长肌、心肌、肝脏、脾脏、肺脏、肾脏、大肠、小肠中均有表达,证明了高通量测序结果的准确性.本试验得到的数据在一定程度上丰富了牛miRNA的数量和种类,为后续牛或其他物种的miRNA相关生物学研究提供了更加详实的信息.     

牛胎儿骨骼肌来源成肌细胞分化相关microRNA的筛查与鉴定

旨在发掘牛胎儿骨骼肌来源的成肌细胞分化过程中差异表达的miRNA,为进一步研究牛胎儿期成肌细胞分化的分子机理提供理论基础。本研究体外培养3头4月龄的牛胎儿骨骼肌来源的成肌细胞,运用高通量测序技术对成肌细胞分化起始和终末两个不同时期进行miRNA测序,并对测序结果进行生物信息学分析,鉴定差异表达的miRNA。随机选取差异表达miRNA,进行荧光定量PCR验证。结果表明,本研究从6个测序文库中共鉴定出已知miRNA 619个,其中差异表达miRNA 150个。差异表达miRNA中bta-miR-199a-5p等多个miRNA已经被证实与肌肉发育相关。通过靶基因预测富集到MAPK等多个与肌肉分化相关的通路。随机选取其中的8个miRNA进行荧光定量PCR验证,定量结果与测序结果完全一致。本研究鉴定出了150个与牛胎儿骨骼肌来源的成肌细胞分化相关的差异表达miRNAs。本研究结果将为研究miRNA调控牛胎儿骨骼肌发育提供科学依据。     

应用生物信息学预测牛新microRNA及验证

本研究旨在通过比较基因组学和结构序列分析方法预测牛基因组中的新miRNAs并进行试验验证。根据miRNA分子序列具有一定保守性,将人、小鼠、绵羊、猪和狗5种哺乳动物已知的miRNA分子与NCBI中牛的全基因组序列(UMD3.1)对比,获得牛miRNAs;随机选取部分预测的miRNAs进行实时荧光定量PCR(qRTPCR)验证。结果,基于物种间序列同源比对共计预测到44条新的牛miRNAs,选取其中4条miRNAs,通过qRTPCR方法,发现他们在泌乳期中国荷斯坦牛乳腺、心、子宫和肝组织中均有表达。结果表明,基于比较基因组学和生物信息学预测新的miRNA分子可行,为奶牛基因表达调控及其性状形成机制研究提供了前期基础。     

不同季节夏南牛外周血新microRNA的挖掘

研究表明microRNA在生物体生长、发育和疾病发生等过程中起着重要作用。本研究运用高通量测序技术对热应激组与非热应激组的夏南牛血液小RNA进行混池测序,并对测序结果进行了生物信息学分析。共鉴定了1 227条microRNAs,其中已知microRNAs为509条,新预测的microRNAs为718条。对新的microRNA的长度、表达量、靶基因预测及KEGG pathway进行分析。预测和分析得到16 077个靶基因,通过KEGG定位到219条通路中,占前4位的通路是Pathways in cancer、MAPK signaling pathway、Cytokine-cytokine receptor interaction、Focal adhesion。本研究旨在发掘更多牛的microRNA信息及新的microRNA序列,为进一步研究microRNA的生物学功能奠定理论基础。     

Solexa测序技术分析鸡与鹌鹑属间杂交雌性和雄性胚胎的差异microRNAs

分别从孵化3d的已鉴定性别的鸡与鹌鹑属间杂交早期雌性和雄性胚胎组织中提取小RNA,并构建各自的cDNA文库,通过Illumina/Solexa二代测序平台对两样本Small RNA进行高通量深度测序,结合生物信息学的方法对其属间杂交早期的雌性和雄性胚胎中已知的和预测的miRNAs的类型、长度、丰度以及参与的KEGG通路等进行了初步的统计和分析。结果显示,鸡与鹌鹑属间杂交早期雌性胚胎组织的比对序列有16 058 009条,雄性组织有17 943 294条;雌性胚胎预测得到132种miRNA,雄性胚胎预测得到460种miRNA;雌性和雄性样本间差异显著的miRNA有117个(P<0.01);KEGG通路分析表明,这些差异miRNAs的靶基因主要参与了Wnt信号通路、MAPK信号通路、TGF-beta信号通路等与胚胎生长发育相关的途径。结果表明,就已知的miRNAs而言,在鸡与鹌鹑属间杂交雌性和雄性胚胎间存在着一定的差异,这些差异的miRNAs将为后续的功能和调控机理实验研究提供依据。     

Solexa测序技术分析鸡与鹌鹑属间杂交雌性和雄性胚胎的差异microRNAs

分别从孵化3d的已鉴定性别的鸡与鹌鹑属间杂交早期雌性和雄性胚胎组织中提取小RNA,并构建各自的cDNA文库,通过Illumina／Solexa二代测序平台对两样本SmallRNA进行高通量深度测序,结合生物信息学的方法对其属间杂交早期的雌性和雄性胚胎中已知的和预测的miRNAs的类型、长度、丰度以及参与的KEGG通路等进行了初步的统计和分析。结果显示,鸡与鹌鹑属间杂交早期雌性胚胎组织的比对序列有16058009条,雄性组织有17943294条;雌性胚胎预测得到132种miRNA,雄性胚胎预测得到460种miRNA;雌性和雄性样本间差异显著的miRNA有117个（P〈0．01）;KEGG通路分析表明,这些差异miRNAs的靶基因主要参与了wnt信号通路、MAPK信号通路、TGF-beta信号通路等与胚胎生长发育相关的途径。结果表明,就已知的miRNAs而言,在鸡与鹌鹑属间杂交雌性和雄性胚胎问存在着一定的差异,这些差异的miRNAs将为后续的功能和调控机理实验研究提供依据。     

鸡性成熟启动前后下丘脑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功能的保守性,该研究结果将为深入开展动物(包括人类)性成熟的分子调控机制研究提供新的理论依据。     

牛基因组中新的microRNA预测及分析

MicroRNA(miRNA)是一类调控真核基因表达的非编码小分子RNA,目前已证实miRNA在生物体生长、发育和疾病发生等过程中起着重要的作用。用同源搜索的计算分析法预测新的miRNAs是当前发现和鉴定miRNA的有效方法之一,本研究根据NCBI数据库中的牛EST、GSS信息以及猪、家犬、人、大猩猩和小家鼠5种哺乳动物的已经注册miRNA分子信息,预测牛新的候选miRNA,通过碱基错配、二级结构、A+U含量、自由能大小等分析候选序列特征,得到了17条牛新的候选miRNA。这对进一步寻找牛的miRNA和遗传育种研究提供了新的思路,具有一定的理论和实际意义。     

牦牛卵巢小RNA高通量测序及生物信息学分析

本研究旨在利用高通量测序技术描绘牦牛卵巢sRNA图谱,为探析牦牛繁殖性能提供基础。以牦牛卵巢组织作为研究对象,构建牦牛sRNA文库,进行高通量测序和生物信息学分析,最后利用RT-qPCR技术验证miRNA在牦牛卵巢组织中的表达。经测序后得到包含12 126 208条clean reads,其中特异序列为212 757条;比对分析显示,只有7 383 536(60.89%)条总序列及80 981(38.06%)条特异序列能与牦牛基因组匹配。与黄牛相比,牦牛sRNA中有56个表达量上调,其中miR-135a表达上调最显著;有33个表达下调,其中miR-2316表达下调最显著。RT-qPCR验证6个miRNA在牦牛卵巢组织中的表达情况,定量结果和测序结果基本一致。与牦牛EST序列信息比对后共预测出5个新miRNA。本研究成功绘制牦牛卵巢sRNA图谱,同时证实RNA-Seq高通量测序技术在完善sRNA信息及挖掘新miRNA方面的优势,对研究sRNA在牦牛遗传育种以及以牦牛为重要高原动物模型来进行高原适应性和抗逆性等相关领域的研究具有重要意义。     

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

乏情和发情初产母猪卵巢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数据资源,为进一步深入研究初产母猪的繁殖性能提供了理论依据。     

乏情和发情初产母猪卵巢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数据资源,为进一步深入研究初产母猪的繁殖性能提供了理论依据。     

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

外源皮质酮对雏鸡法氏囊microRNA表达的影响

MicroRNA(miRNA)是在真核生物中发现的一类内源性的具有调控功能的非编码小RNA,在鸡的许多生物学功能中发挥重要作用,但在鸡的应激性免疫抑制的分子调控机制中尚不完全清楚。本试验通过在日粮中添加外源皮质酮激素,构建雏鸡应激性免疫抑制模型,高通量Solexa测序技术建立雏鸡法氏囊差异miRNA表达谱,共鉴定出77个差异表达的miRNAs,34个上调,43个下调,生物学分析预测差异表达miRNA的靶基因有282个,对这些预测的靶基因进行GO和KEGG通路富集分析,结果显示这些通路主要富集到心肌收缩、精氨酸和脯氨酸代谢、RNA聚合酶、肾素-血管紧张素系统。试验结果为揭示家禽应激性免疫抑制转录后分子调控机制提供了依据。     

成年滩羊和小尾寒羊皮肤毛囊差异表达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可能通过调控其靶基因最终参与了绵羊皮肤毛囊的发育。     

猪血清和血浆外泌体miRNA的鉴定与比较

为了解血液中外泌体miRNA的表达情况,利用高通量测序技术检测了来源于相同个体的猪血清和血浆外泌体miRNA表达谱。结果显示,血清外泌体中提取的RNA浓度显著高于血浆外泌体(P=0.02),从2个文库中共检测到187个成熟体miRNAs,对应165个前体。鉴定出了在2个文库中丰富表达的17个miRNAs,暗示这些miRNAs在血清和血浆中发挥着重要的作用。筛选出99个差异表达miRNAs,其中54个在血清外泌体中高表达,45个在血浆外泌体中高表达。通过靶基因预测和功能注解分析,丰富表达的miRNAs广泛地参与到机体的生理学过程中。本研究对猪血清和血浆外泌体miRNA进行了鉴定及比较,为深入研究外泌体miRNA在生物体内的功能提供理论依据。     

罗伊氏乳杆菌对仔猪空肠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表达来调控有关肠道健康的信号通路。     

基于实时荧光定量反转录PCR技术的动物miRNA表达分析策略

在人体内有约60%的蛋白质编码基因受到到miRNA调控。近年来,miRNA参与调控的与动物生长相关的细胞周期、发育、分化和新陈代谢等多种生物学过程也越来越受到关注。目前针对miRNA表达检测的常用方法有很多,其中基于反转录PCR(RT-PCR)的放大检测技术是应用范围最广,实施手段最灵活的一种研究策略。主要包括miRNA的全定量PCR法(miR-Q)、基于茎环引物的实时荧光定量PCR法、Poly(A) 加尾性通用反转录法、miRNA表达谱的高通量RT-PCR分析和miRNA扩增谱系(mRAP)法。作者对这几种基于实时荧光定量RT-PCR技术的动物miRNA表达分析策略分别进行了介绍,对其原理、方法和应用范畴进行了概括性总结。     

西藏牛和三江牛大脑组织中低氧适应相关circRNAs的比较分析

旨在从circRNA层面加深哺乳动物对高原低氧环境适应性的认识,探索西藏牛和三江牛大脑组织中差异表达的circRNA及相关调控网络,为研究circRNA在西藏牛低氧适应性方面的分子调控机制奠定理论基础。本研究分别采集3头4.5岁健康、雌性西藏牛和三江牛的大脑组织作为试验样本,构建cDNA文库进行高通量测序,利用生物信息学方法对宿主基因进行GO和KEGG分析,预测差异表达circRNA下游靶向基因,构建circRNA-miRNA和circRNA-miRNA-mRNA可视化调控网络,辅以RNaseR酶抗性检测和实时荧光定量PCR（qRT-PCR）验证测序数据的可靠性。结果,在西藏牛和三江牛的6个样本中,共筛选获得858个显著差异的circRNAs,其中上调表达394个,下调表达464个。其宿主基因共参与49个功能亚类,注释到31条显著富集的信号通路（P<0.05）,主要涉及MAPK信号通路、谷氨酸能突触、Rap1信号通路、磷脂酶D信号通路及cGMP-PKG信号通路等生物学过程。靶基因预测结果显示,350个上调和364个下调差异表达circRNAs分别靶向结合492和508个miRNAs。其中,miRNA靶点数最多的是novel_circ_017825和novel_circ_046762,说明这两个circRNAs可能在西藏牛脑功能调控方面发挥重要作用。circRNA-bta-miR-2284z-mRNA调控网络展示了bta-miR-2284z与circRNA、mRNA间的靶向关系,推测上述circRNAs可能通过充当bta-miR-2284z海绵的方式,间接影响西藏牛脑组织中相关靶向mRNA的表达水平。随机挑选6个circRNA进行RNaseR酶抗性试验和qRT-PCR验证,表达趋势与测序结果一致,说明所获得的circRNA为环状转录本。本研究利用高通量测序技术获得了西藏牛和三江牛大脑组织中circRNA的表达谱及信号通路,并初步构建了circRNA-miRNA-mRNA网络互作模型,为后续进一步探索circRNA参与西藏牛低氧适应性的生物学过程和分子功能,研究circRNA在大型哺乳动物低氧适应性方面的调控机制提供了可靠的数据支持。     

西藏牛和三江牛大脑组织中低氧适应相关circRNAs的比较分析

旨在从circRNA层面加深哺乳动物对高原低氧环境适应性的认识,探索西藏牛和三江牛大脑组织中差异表达的circRNA及相关调控网络,为研究circRNA在西藏牛低氧适应性方面的分子调控机制奠定理论基础。本研究分别采集3头4.5岁健康、雌性西藏牛和三江牛的大脑组织作为试验样本,构建cDNA文库进行高通量测序,利用生物信息学方法对宿主基因进行GO和KEGG分析,预测差异表达circRNA下游靶向基因,构建circRNA-miRNA和circRNA-miRNA-mRNA可视化调控网络,辅以RNaseR酶抗性检测和实时荧光定量PCR(qRT-PCR)验证测序数据的可靠性。结果,在西藏牛和三江牛的6个样本中,共筛选获得858个显著差异的circRNAs,其中上调表达394个,下调表达464个。其宿主基因共参与49个功能亚类,注释到31条显著富集的信号通路(P0.05),主要涉及MAPK信号通路、谷氨酸能突触、Rap1信号通路、磷脂酶D信号通路及cGMP-PKG信号通路等生物学过程。靶基因预测结果显示,350个上调和364个下调差异表达circRNAs分别靶向结合492和508个miRNAs。其中,miRNA靶点数最多的是novel_circ_017825和novel_circ_046762,说明这两个circRNAs可能在西藏牛脑功能调控方面发挥重要作用。circRNA-bta-miR-2284z-mRNA调控网络展示了bta-miR-2284z与circRNA、mRNA间的靶向关系,推测上述circRNAs可能通过充当bta-miR-2284z海绵的方式,间接影响西藏牛脑组织中相关靶向mRNA的表达水平。随机挑选6个circRNA进行RNaseR酶抗性试验和qRT-PCR验证,表达趋势与测序结果一致,说明所获得的circRNA为环状转录本。本研究利用高通量测序技术获得了西藏牛和三江牛大脑组织中circRNA的表达谱及信号通路,并初步构建了circRNA-miRNA-mRNA网络互作模型,为后续进一步探索circRNA参与西藏牛低氧适应性的生物学过程和分子功能,研究circRNA在大型哺乳动物低氧适应性方面的调控机制提供了可靠的数据支持。