牦牛皱胃全长转录组测序分析

doi:10.16431/j.cnki.1671-7236.2022.05.002

摘要/Abstract

摘要： 【目的】获得牦牛(Bos grunniens)皱胃全长转录组数据库，深入挖掘牦牛皱胃功能基因。【方法】采用PacBio Sequel高通量测序系统，使用单分子实时(single molecular real time，SMRT)测序技术对成年牦牛皱胃全长转录组进行测序，对原始数据进行质控和去冗余分析，再比对参考基因组获取过滤后的非重复序列(Unigenes)，使用多种生物信息软件对牦牛皱胃全长转录本数据进行功能注释、转录因子注释、编码区预测、简单重复序列(simple sequence repeats，SSR)分析及可变剪接分析。【结果】通过测序共获得14467420条子序列，平均子序列长度为3 344.23 bp，质控得到循环一致性序列(CCS)有296840条，全长非嵌合(FLNC)序列有277 402条，过滤和去冗余后对比参考基因组最终获得8 556条Unigenes。通过与NR、Swiss-Prot、KEGG、KOG、eggNOG、GO和Pfam数据库比对，对Unigenes进行注释，其中NR数据库注释了8 544条Unigenes;Swiss-Prot数据库注释了8 475条Unigenes;KEGG数据库注释了1721条Unigenes;KOG数据库注释了6 572条Unigenes;eggNOG数据库注释了8491条Unigenes;GO数据库注释了7 725条Unigenes;Pfam数据库注释了8 162条Unigenes。此外，经鉴定或预测，还获得了943个转录因子、8544个编码区片段、3596个SSR位点和1 825个可变剪接事件。【结论】本研究获得了较为可靠的牦牛皱胃全长转录组数据，可为进一步研究牦牛皱胃生物学特性、相关代谢途径、信号通路及其分子机制提供数据支持。

关键词: 牦牛; 皱胃; 全长转录组; 第3代测序; 功能注释

Abstract: 【Objective】 In order to obtain the full-length transcriptome database of abomasum in yak (Bos grunniens),the functional genes of abomasum of yak were deeply excavated.【Method】 The full-length transcriptome of abomasum of adult yak was sequenced by single molecular real time (SMRT) sequencing technology using high-throughput sequencing PacBio Sequel high-throughput sequencing system.The original data were quality controlled and de redundant analysis,and then compared with the reference genome to obtain filtered non repetitive sequences (Unigenes).A variety of bioinformatics software were used to compare the full-length transcriptome data of abomasum in yak for functional annotation,transcription factor annotation,coding region prediction,simple sequence repeats (SSR) analysis and alternative splices analysis.【Result】 A total of 14 467 420 subsequences were obtained by sequencing,with an average length of 3 344.23 bp.296 840 CCS sequences and 277 402 FLNC sequences were obtained by quality control.After filtering and de redundancy,the reference genome was compared,and 8 556 Unigenes were finally obtained.By comparing with NR,Swiss-Prot,KEGG,KOG,eggNOG,GO and Pfam databases,Unigenes were annotated,of which 8 544 Unigenes were annotated in NR database.Swiss-Prot database annotated 8 475 Unigenes.The KEGG database annotated 1 721 unigenes.The KOG database annotated 6 572 Unigenes.The eggNOG database annotated 8 491 Unigenes.GO database annotated 7 725 Unigenes;The Pfam database annotated 8 162 Unigenes.In addition,943 transcription factors,8 544 CDS fragments,3 596 SSR sites and 1 825 alternative splicing events were identified or predicted.【Conclusion】 In this study,reliable full-length transcriptome data of yak abomasum were obtained,it could provide data support for the further study of the biological characteristics,related metabolic pathways,signal pathways and molecular mechanisms of abomasum in yak.

Key words: yak; abomasum; full-length transcriptome; the third generation sequencing; functional notes

中图分类号:

S823

谢书琼, 马士龙, 唐娇, 刘益丽, 江明锋. 牦牛皱胃全长转录组测序分析[J]. 中国畜牧兽医, 2022, 49(5): 1610-1620.

XIE Shuqiong, MA Shilong, TANG Jiao, LIU Yili, JIANG Mingfeng. Sequencing Analysis of Full-length Transcriptome of Abomasum in Yak[J]. China Animal Husbandry and Veterinary Medicine, 2022, 49(5): 1610-1620.

参考文献

[1] 孟祥琼,胡学敏,王悦如,等.不同发育阶段牦牛肺组织内肥大细胞的分布[J].畜牧兽医科技信息,2020,9:43-44. MENG X Q,HU X M,WANG Y R,et al.The distribution of mast cells in lung tissue of yak at different developmental stages[J].Livestock and Veterinary Science and Technology Information,2020,9:43-44.(in Chinese)
[2] 胡泉军.牦牛基因组数据库建设[D].兰州:兰州大学,2014. HU Q J.Construction of yak genome database[D].Lanzhou:Lanzhou University,2014.(in Chinese)
[3] 张剑搏,丁学智,AHMAD A A,等.高原土著动物适应性进化的研究进展[J].畜牧兽医学报,2019,50(9):1723-1736. ZHANG J B,DING X Z,AHMAD A A,et al.Advances in research on adaptive evolution of native animals of Tibetan plateau[J].Acta Veterinaria et Zootechnica Sinica,2019,50(9):1723-1736.(in Chinese)
[4] QI X,ZHANG Q,HE Y,et al.The transcriptomic landscape of yaks reveals molecular pathways for high altitude adaptation[J].Genome Biology and Evolution,2019,11(1):72-85.
[5] WANG H,ZHONG J,ZHANG C,et al.The whole-transcriptome landscape of muscle and adipose tissues reveals the ceRNA regulation network related to intramuscular fat deposition in yak[J].BMC Genomics,2020,21:1-15.
[6] WU S,MIPAM T,XU C,et al.Testis transcriptome profiling identified genes involved in spermatogenic arrest of cattle yak[J].PloS One,2020,15(2):e0229503.
[7] XU S R,WEI P,YANG Q L,et al.Transcriptome analysis revealed key signaling networks regulating ovarian activities in the domestic yak[J].Theriogenology,2020,147:50-56.
[8] 鲁延刚,高世杰,方梅,等.牦牛皱胃组织结构及黏膜免疫相关细胞研究[J].畜牧兽医学报,2012,43(3):469-475. LU Y G,GAO S J,FANG M,et al.Study on the tissue structure and mucosal lmmunity-associated cells in Abomasum of the yak[J].Acta Veterinaria et Zootechnica Sinica,2012,43(3):469-475.(in Chinese)
[9] 杨晓晴.牦牛皱胃分泌细胞的分布及GSDMA/HSP70在胃的表达与变化规律[D].兰州:甘肃农业大学,2018. YANG X Q.The distribution of gastric decretory cells and the expression and variation of GSDMA/HSP70 in yak stomach[D].Lanzhou:Gansu Agricultural University,2018.(in Chinese)
[10] 鲁延刚.牦牛皱胃组织结构及黏膜免疫相关细胞的研究[D].兰州:甘肃农业大学,2012. LU Y G.Study on the tissue structure and mucosal immunity-associated cells in abomasum of the yak[D].Lanzhou:Gansu Agricultural University,2012.(in Chinese)
[11] DOBSON D E,PRAGER E M,WILSON A.Stomach lysozymes of ruminants.Ⅰ.Distribution and catalytic properties[J].Journal of Biological Chemistry,1984,259(18):11607-11616.
[12] JOLLÈS J,JOLLES P,BOWMAN B H,et al.Episodic evolution in the stomach lysozymes of ruminants[J].Journal of Molecular Evolution,1989,28(6):528-535.
[13] SALK J J,SCHMITT M W,LOEB L A.Enhancing the accuracy of next-generation sequencing for detecting rare and subclonal mutations[J].Nature Reviews Genetics,2018,19(5):269-285.
[14] MCCARTHY A.Third generation DNA sequencing:Pacific biosciences' single molecule real time technology[J].Chemistry & Biology,2010,17(7):675-676.
[15] CHIN C-S,PELUSO P,SEDLAZECK F J,et al.Phased diploid genome assembly with single-molecule real-time sequencing[J].Nature Methods,2016,13(12):1050-1054.
[16] GORDON S P,TSENG E,SALAMOV A,et al.Widespread polycistronic transcripts in fungi revealed by single-molecule mRNA sequencing[J].PloS One,2015,10(7):e0132628.
[17] ABDEL-GHANY S E,HAMILTON M,JACOBI J L,et al.A survey of the sorghum transcriptome using single-molecule long reads[J].Nature Communications,2016,7(1):1-11.
[18] TARDAGUILA M,DE LA FUENTE L,MARTI C,et al.SQANTI:Extensive characterization of long-read transcript sequences for quality control in full-length transcriptome identification and quantification[J].Genome Research,2018,28(3):396-411.
[19] LI H.Minimap2:Pairwise alignment for nucleotide sequences[J].Bioinformatics,2018,34(18):3094-3100.
[20] ISELI C,JONGENEEL C V,BUCHER P.ESTScan:A program for detecting,evaluating,and reconstructing potential coding regions in EST sequences[C].ISMB,1999,99:138-148.
[21] BEIER S,THIEL T,MVNCH T,et al.MISA-web:A web server for microsatellite prediction[J].Bioinformatics,2017,33(16):2583-2585.
[22] KORESSAAR T,REMM M.Enhancements and modifications of primer design program Primer3[J].Bioinformatics,2007,23(10):1289-1291.
[23] FOISSAC S,SAMMETH M.ASTALAVISTA:Dynamic and flexible analysis of alternative splicing events in custom gene datasets[J].Nucleic Acids Research,2007,35(suppl_2):W297-W299.
[24] RHOADS A,AU K F.PacBio sequencing and its applications[J].Genomics,Proteomics & Bioinformatics,2015,13(5):278-289.
[25] 王明成.牦牛基因组组装与注释的改进[D].兰州:兰州大学,2018. WANG M C.Yak genome assembly and annotation improvements[D].Lanzhou:Lanzhou University,2018.(in Chinese)
[26] 张春兰.小尾寒羊和杜泊羊臂二头肌转录组及肌球蛋白轻链基因家族结构特征分析[D].泰安:山东农业大学,2014. ZHANG C L.Transcriptome analysis of Small-tailed Han sheep and Dorper's biceps brachii and structure characteristics of myosin light chain gene families[D].Tai'an:Shandong Agricultural University,2014.(in Chinese)
[27] MILLEN,DANILO D,ARRIGONI M D B,et al.Rumenology[M].Cham,Switzerland:Springer International Publishing,2016.
[28] MYRZABEKOVA M,LABEIT S,NIYAZOVA R Y.Features of miRNAs binding sites within the C2H2 ZNF family:A Bos taurus,Eqqus caballus,and Ovies aries comparative approach[J].International Journal of Biology and Chemistry,2020,13(1):33-46.
[29] 李枫梅.肉用牛bHLH转录因子家族成员的鉴定与分析[D].阜阳:阜阳师范学院,2017. LI F M.Identification and functional analysis of the basic helix-loop-helix (bHLH) transcription factors in the meat cattle[D].Fuyang:Fuyang Normal University,2017.(in Chinese)
[30] ZAKANY J,DUBOULE D.Hox genes and the making of sphincters[J].Nature,1999,401(6755):761-762.
[31] 鲍晶晶,浦亚斌,马月辉,等.绵羊不同发育阶段背最长肌组织中可变剪接的鉴定与分析[J].生物技术通报,2019,35(7):32-38. BAO J J,PU Y B,MA Y H,et al.ldentification and analysis of alternative splicing in longissimus dorsi of sheep at different development stages[J].Biotechnology Bulletin,2019,35(7):32-38.(in Chinese)