中国畜牧兽医 ›› 2021, Vol. 48 ›› Issue (8): 2713-2726.doi: 10.16431/j.cnki.1671-7236.2021.08.004

• 生物技术 • 上一篇    下一篇

沐川乌骨黑鸡肌内脂肪沉积相关microRNAs的筛查与鉴定

喻世刚1, 王钢1, 廖娟1, 谢伟2   

  1. 1. 乐山师范学院, 川南地方品种鸡产业化四川省高等学校工程研究中心, 乐山 614000;
    2. 乐山市动物疫病预防控制中心, 乐山 614000
  • 收稿日期:2020-12-01 出版日期:2021-08-20 发布日期:2021-08-19
  • 通讯作者: 王钢 E-mail:Lswanggang@163.com
  • 作者简介:喻世刚(1985-),男,重庆人,副教授,研究方向:家禽遗传育种与繁殖,E-mail:shigang_yu@163.com
  • 基金资助:
    四川省科技计划资助项目(2018JY0452、2019YJ0304);乐山市科技项目(20NZD005);乐山师范学院科技项目(DGZZ202002)

Screening and Identification of the microRNAs Associated with Intramuscular Fat Deposition in Muchuan Black-bone Chicken

YU Shigang1, WANG Gang1, LIAO Juan1, XIE Wei2   

  1. 1. Engineering Research Center of Sichuan Province Higher School of Local Chicken Breeds Industrialization in Southern Sichuan, Leshan Normal University, Leshan 614000, China;
    2. Leshan Animal Disease Prevention and Control Center, Leshan 614000, China
  • Received:2020-12-01 Online:2021-08-20 Published:2021-08-19

摘要: 试验旨在挖掘高肌内脂肪(IMF)含量和低IMF含量鸡胸肌中差异表达的microRNA(miRNA)(DEM),以期从miRNA角度探究鸡IMF沉积的分子调控机制。以120只沐川乌骨黑鸡为研究对象,采用高通量测序分别完成3个IMF含量极高(H组)和极低(L组)鸡胸肌组织的miRNA测序,测序结果经生物信息学分析,鉴定H和L组中差异表达的miRNAs,随机选择6个miRNAs进行实时荧光定量PCR验证。IMF含量测定结果表明,每只鸡平均100 g胸肌中IMF含量为4.08 g,H和L组IMF含量存在极显著差异(P<0.01);6个测序文库结果表明,每个测序文库获得Clean reads均超过10 000 000条,Clean reads占Raw reads的百分比>93%,21~23 nt的small RNA(sRNA)数量占总sRNA的百分比>50%,其中22 nt长度的sRNA所占比例最高,平均GC碱基含量为48.95%,获得的测序数据可靠;6个测序文库中sRNA注释为rRNA、tRNA、snRNA、snoRNA以及重复序列的比例较低,分别为2.74%、1.85%、6.21%、1.58%、2.82%和2.81%;6个测序文库共鉴定已知miRNAs成熟体578个,miRNAs前体500个,鉴定差异表达的miRNAs 23个,其中H组中上调表达的miRNAs有16个,下调表达的miRNAs有7个;miRNA靶基因预测分析结果表明,23个差异表达miRNAs共预测靶基因628个;GO分析结果表明,共鉴定了30个显著富集GO条目,包括6个细胞成分、13个生物学过程和11个分子功能;KEGG富集分析表明,差异表达miRNA的靶基因显著富集于胰岛素信号通路、半乳糖代谢、脂肪细胞因子信号通路、鞘糖脂生物合成、糖酵解及淀粉和蔗糖代谢;实时荧光定量PCR结果表明,gga-miR-124a-3p、gga-let-7a-5p在H组中表达量极显著或显著高于L组(P<0.01;P<0.05);gga-miR-19a-3p在H组中表达量极显著低于L组(P<0.01),gga-miR-6553-3p和gga-miR-128-3p在H组中的表达量显著低于L组(P<0.05),其表达趋势与测序结果一致。以上结果表明,差异表达的miRNA gga-miR-124a-3p、gga-let-7a-5p、gga-miR-6553-3p、gga-miR-128-3p和gga-miR-19a-3p参与脂肪生成,是鸡IMF含量调控重要的候选miRNAs,为阐释miRNA调控鸡IMF沉积分子机制提供理论基础。

关键词: 鸡; 肌内脂肪; miRNA; 高通量测序

Abstract: This study was aimed to identify differentially expressed microRNA(miRNA) (DEM)in the breast muscle between high intramuscular fat (IMF) and low IMF content chickens, and explore the miRNA regulatory mechanisms involved in IMF deposition in chicken.120 Muchuan Black-bone chickens were employed in this study.High-throughput RNA sequencing (RNA-Seq) was performed to identify miRNA in the breast muscle between three high IMF (H group) and low IMF (L group) content chickens.The sequencing data was analyzed to identify DEMs.Six miRNAs were randomly selected for quantitative analysis using Real-time quantitative PCR.The results of IMF contents determination showed that the average IMF level of 120 chickens was 4.08 g per 100 g muscle, and there was extremely significant difference in IMF content between H and L groups (P<0.01).Results of 6 RNA-Seq libraries showed that more than 10 000 000 Clean reads were obtained for each library, and the Clean reads ratio was greater than 93%, more than 50% of small RNA (sRNA) ranged from 21 to 23 nt in length, and the 22 nt sRNA was the most abundant.The average GC content was 48.95%.These results indicated that the sequencing data were reliably produced in our study.The proportion of sRNA annotated as rRNA, tRNA, snRNA, snoRNA and repetitive sequence in the 6 RNA-Seq libraries was 2.74%, 1.85%, 6.21%, 1.58%, 2.82% and 2.81%, respectively.A total of 578 known miRNAs and 500 pre-miRNAs were identified from the 6 chicken sequence libraries, with 23 DEMs were identified between two groups, of which 16 were upregulated and 7 were downregulated.The prediction analysis of miRNA target genes showed that 23 miRNAs which were differentially expressed between two chicken groups targeted a total of 628 genes.GO analysis results showed that a total of 30 significantly enriched GO terms were identified, including 6 cellular components, 13 biological process and 11 molecular function.KEGG pathway analysis showed that the target genes were mainly involved with insulin signaling pathway, galactose metabolism, adipocytokine signaling pathway, glycosphingolipid biosynthesis, glycolysis and gtarch and sucrose metabolism.Real-time quantitative PCR results indicated that the expression of gga-miR-124a-3p and gga-let-7a-5p in H group was significantly higher than that in L group (P<0.01;P<0.05), the expression of gga-miR-19a-3p in H group was significantly lower than that in L group (P<0.01), the expression of gga-miR-6553-3p and gga-miR-128-3p in H group were significantly lower than that in L group (P<0.05), and 5 DEMs expression pattern were consistent with sequencing data.The results showed that gga-miR-124a-3p, gga-let-7a-5p, gga-miR-6553-3p, gga-miR-128-3p and gga-miR-19a-3p were associated with adipogenesis in animal, which might be important for IMF deposition in chicken and would provide a foundation for clarifying the molecular regulatory mechanisms of miRNA involved in IMF deposition in chicken.

Key words: chicken; intramuscular fat; miRNA; high throughput sequencing

中图分类号: