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

本研究旨在发掘牛更多的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相关生物学研究提供了更加详实的信息。     

牛多个组织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预测及分析

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

利用生物信息学方法预测家蚕核型多角体病毒基因组编码的miRNA

microRNAs(miRNAs)是长度为18～25 nt的非蛋白质编码小RNA分子,可以通过抑制有关基因mRNA的翻译从而调控机体的生长发育、疾病发生等过程。在很多病毒中发现有miRNA的存在,这些miRNA可能在病毒基因的表达以及与宿主的相互作用中发挥作用。利用vMir软件预测家蚕核型多角体病毒(BmNPV)基因组可能编码的miRNA前体序列,对BmNPV基因组扫描分析后,得到7条能形成发夹结构的可能是miRNA的候选前体序列。以BmNPV添食感染5龄起蚕,72 h后取幼虫血淋巴提取总RNA,对这7条序列进行反转录PCR验证,其中扩增出的5条序列在BmNPV基因组中有相同的序列,推测这5条序列为BmNPV基因组编码的miRNA前体序列。进一步分析5条miRNA前体序列的发夹结构,并用反转录PCR的方法验证其成熟体序列,结果在5条miRNA前体序列共扩增出6条在BmNPV基因组中有相同序列的成熟体序列,它们分别位于不同的ORF之间,是BmNPV基因组编码的miRNA。利用在线靶基因预测程序RNAhybrid和RNA22预测到6条miRNA在BmNPV中的13个可能的靶基因,其中有6个靶基因是病毒的结构蛋白基因,2个为晚期表达因子基因,5个为未知功能基因。推测这些BmNPV基因组编码的miRNA可能与病毒感染宿主有关。     

鸡性成熟启动前后下丘脑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的筛查与鉴定

旨在发掘牛胎儿骨骼肌来源的成肌细胞分化过程中差异表达的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转录组

本研究采集的Texel绵羊胎儿的背最长肌样品包括5个发育阶段:70、85、100、120和135 d。通过采用Solexa测序技术对混池样品进行了microRNA(miRNA)深度测序,获得了16532850条原始序列读数。使用ACGT101-miR v4.2分析软件对miRNA测序数据进行了深度发掘;通过与最新的哺乳动物成熟miRNA序列(miRNA数据库:miRBase v17.0)、miRNA前体序列、绵羊基因组序列(绵羊基因组数据库,2010年2月)的比对分析,对绵羊microRNA转录组数据库进行了大幅更新,pre-miRNA(miRNA前体)序列增至1529条,编码的miRNA成熟体序列增至1999条。通过实时荧光定量PCR技术对深度测序获得的5个miRNA在混池样品中进行了验证。绵羊miRNA的研究起步较晚,而miRNA的发现与鉴定将促进基因调控机制及miRNA功能的研究。     

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

基于UniGene寻找山羊新的miRNA及其试验验证

根据miRNA进化的保守性,以已知动物的miRNA为搜索序列与山羊的UniGene进行BLAST比对,预测山羊新的miRNA,并通过RT-PCR和克隆测序进行验证。对新发现的山羊miRNA进行实时荧光定量,检测其在12个月皮肤及肌肉中的表达量;用基于UniGene的方法获得了5条山羊miRNA候选分子,通过RT-PCR和克隆测序验证,发现这5条miRNA分子在山羊的皮肤和肌肉中均有表达,其中由chi-miR-374b、chi-miR-421和chi-miR-421*构成了山羊miRNA基因簇,chi-miR-2284n是牛和山羊特异的。实时荧光定量结果显示,chi-miR-1839、chi-miR-374b和chi-miR-2284n在2月份皮肤中的表达量明显高于其他月份,chi-miR-421和chi-miR-421*在10月份皮肤中的表达量最高;发现了5条山羊新的miRNA(chi-miR-2284n、chi-miR-421*、chi-miR-421、chi-miR-1839和chi-miR-374b,登录号:JQ002550-JQ002554),补充了山羊miRNA数据库的不足。     

牛肺泡巨噬细胞微小RNA的克隆与鉴定

MicroRNA(miRNA)是一类长19～26个核苷酸的内源性非编码单链小分子RNA,其异常表达将会导致生理的异常和疾病的产生.为了研究miRNAs在牛肺泡巨噬细胞抗菌机制中的作用,构建了牛肺泡巨噬细胞miRNA cDNA文库,从438个克隆中筛选出22个miRNAs,其中有19个miRNAs在牛的miRBase数据库中没有发现,但是有8个miRNAs与人和(或)鼠中的序列完全同源,因此,将它们命名为Bta-miR-141、Bta-miR-187、Bta-miR-191、Bta-miR-448、Bta-miR-589、Bta-miR-873、Bta-miR-463和Bta-miR-562;另外有11个miRNAs在所有miR-Base数据库中都没有发现,有待进一步研究.这些数据为研究miRNAs调控牛肺泡巨噬细胞抗菌机制奠定了基础.     

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

饲粮中添加共轭亚油酸对猪肌肉组织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在牛早期妊娠中的作用及其调控机制。【方法】 以3~6岁、体重480~600 kg的夏南牛作为研究对象, 选取10头供体牛进行同期发情、超数排卵和人工授精, 23头受体牛只做同期发情处理。在人工授精后第7天, 冲洗供体牛子宫以获取囊胚, 选取3头囊胚数相近的供体牛及3头与供体牛体重和年龄均相近的受体牛, 颈静脉采血, 进行血浆外泌体的分离与鉴定; 然后提取血浆外泌体miRNA, 并检测其表达量; 采用R语言中的DESeq差异算法计算P值, 并筛选出P<0.05的miRNA, 对差异表达的miRNA进行靶基因预测、GO功能富集分析和KEGG信号通路分析。【结果】 6个样本的囊泡粒径均在135 nm左右, 符合外泌体的特征。与供体牛相比, 受体牛中有9个miRNAs表达显著上调(P<0.05), 13个miRNAs显著下调(P<0.05);22个差异表达的miRNAs中, 有15个miRNAs预测出无重复的靶基因2 990个。GO功能富集分析和KEGG信号通路分析的结果表明, 这些靶基因主要富集在与生物黏附(biological adhesion)、定位(localization)、细胞连接(cell junction)功能有关的通路上, 显著富集的信号通路与黏着斑(focal adhesion)、黏着连接(adherens junction)有关, 提示血浆外泌体miRNA可能参与调控胚胎着床。【结论】 研究结果可为筛选和探究影响胚胎着床的血浆外泌体miRNA提供参考, 并为进一步阐明血浆外泌体miRNA在母牛早期妊娠调控中的作用提供依据。     

家蚕环形RNA的预测与鉴定

环形RNA(circRNA)是具有闭合环状结构的单链RNA,其功能涉及对亲本基因转录的调控、miRNA的吸附和作为某些疾病的生物标志物。以家蚕幼虫、蛹和蛾3个发育时期蚕体的混合样本提取总RNA,利用高通量测序共获得112 271 320个测序read。采用软件Map Splice和Find_circ从获得的测序read中预测家蚕的circRNAs,其中用Map Splice预测到的circRNA共1 659条,用Find_circ预测到的circRNA共9 777条,2个软件共同预测到的circRNA有1 598条。利用circRNA预测程序CircRNApredictor.pl从家蚕转录组read(SRR315361)和家蚕基因组序列中预测到12 955条circRNA,进一步将高通量测序read匹配预测的circRNA序列,结果显示序列首尾重叠,为环状分子。通过软件miRanda预测了部分家蚕circRNA的靶miRNA。设计特异性的正、反向PCR引物对,选取其中11条SRPBM(spliced reads per billion mapping)值最高的circRNA进行PCR,产物测序鉴定为家蚕潜在的circRNA分子,其中BmcircR_3621、BmcircR_8955、BmcircR_8926、BmcircR_6123、BmcircR_8349、BmcircR_1198、BmcircR_9535、BmcircR_6215、BmcircR_7937在家蚕酪氨酸蛋白激酶Abl(tyrosine-protein kinase Abl)和异常细胞迁移蛋白10(abnormal cell migration protein 10)等的基因序列中以环的形式存在,并且主要成环位点与生物信息学预测结果一致,而BmcircR_2196和BmcircR_6226的成环位点与预测存在差异。研究结果可供家蚕乃至昆虫新型RNA分子的鉴定和研究参考。     

中国山羊痘弱毒疫苗AV41株基因组分子标记的比较分析

为鉴别我国山羊痘病毒(GTPV)弱毒疫苗AV41株与牛结节性皮肤病病毒(LSDV)野毒株,运用生物信息学方法将AV41株基因组全序列与GenBank中已发表的LSDV、GTPV和绵羊痘病毒(SPPV)基因组进行比较,筛选AV41株基因组分子标记.采用我国不同厂家来源的AV41株疫苗测序验证分子标记相关序列,并与我国LS...     

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

基于全基因组信息鉴定中国荷斯坦牛产奶性状基因及功能注释

前期研究基于中国荷斯坦牛女儿设计资源群体,采用Illumina公司Bovine50K微珠芯片对产奶性状进行了全基因组关联分析(GWAS),利用传递不平衡(L1-TDT)和回归分析(MMRA)2种统计分析方法共同检测到35个显著SNPs位点。本研究旨在基于该GWAS结果进一步对产奶性状基因进行鉴定及功能注释。基于牛基因组序列草图(Btau4.2),采用生物信息学和比较基因组学方法进行显著SNPs位置候选基因筛查和功能预测。分析发现,12个SNPs位点位于基因内部,23个位于基因侧翼,最终鉴定到28个位置候选基因,并确定了其物理位置、基因类型及潜在功能。基因功能可归纳为6种类型:调节机体营养成分代谢和平衡、细胞骨架或基质成分、调节细胞增殖和周期及凋亡、参与细胞信号转导和盐离子通道构成、具有激酶活性、参与mRNA转录调控或翻译调控。该研究为进一步鉴定中国荷斯坦牛产奶性状主效基因及功能验证打下了基础。     

罗伊氏乳杆菌对仔猪空肠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组图谱研究

家猪对农业有重要的意义,也是研究人类疾病良好的模式动物,而猪microRNA(miRNA)组图谱研究滞后,有必要建立一个汇总micRNA调节功能和相互作用的文库以便对其信息进行解读。本研究收集了猪从胚胎到成年发育的全部代表性时期共10个阶段的样品,建立了miRNA序列库,并对其miRNA进行了研究。测序结果与哺乳动物miRNA、前体发夹序列(pre-miRNA)、首次公布的猪基因组序列,以及表达序列标签(EST)进行比对分析。该研究结果使猪miRNA组序列数增至867条(623条与基因组序列相匹配),可编码1004个miRNA,其中777个是特异性序列。该研究还运用实时定量PCR技术从47个特异组织样品中选择30种miRNA进行试验,定量结果与测序结果完全一致。本研究成果提供了很多关于miRNA的末端序列变异、miRNA前体结构、染色体定位、特定发育阶段的miRNA表达和保守序列的信息;为猪miRNA组图谱、miRNA种类及异构体和miRNA特征的研究提供了很多宝贵的数据资源;也将大大促进猪生物学和以猪作为模式动物进行的人类生物学和生物医学的发展。     

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将为后续的功能和调控机理实验研究提供依据。