基于转录组测序的奶牛产奶性状重要功能基因挖掘

doi:10.13560/j.cnki.biotech.bull.1985.2020-0255

摘要/Abstract

摘要： 挖掘奶牛产奶性状功能基因和重要遗传变异历来是我国奶牛分子育种的关键。本研究分别采集高产奶牛和低产奶牛泌乳期的血液组织,利用转录组测序进行高低产组差异表达基因分析、遗传变异的检测和功能注释。基因表达分析结果显示,在两组间共检测到431个差异表达基因,这些基因富集在代谢过程、激素调节、细胞迁移、免疫等生物学过程。研究同时检测到373 049个SNPs和42 273个Indels,其中75 363个是新发现的遗传变异。将遗传变异在高低产组中的等位基因频率分布情况与基因表达分析结果进行整合分析,最终筛选出12个影响奶牛产奶性状的重要候选基因：ASS1、CKB、GGT1、MGAM、SDSL、HP、LTF、MMP9、UPP1、DEFB4A、PGLYRP1、MS4A8和14个位于这些基因内的重要遗传变异。

关键词: 产奶性状, 转录组, 基因表达, 功能基因, SNP

Abstract: Exploring genetic variants affecting milk production traits is very crucial for molecular breeding in dairy cattle in China. In current study,RNA-seq technology was used to detect differential expression and genetic variation of genes in the blood samples between 2 groups of cows with high and low milk production. As a result,431 differentially expressed genes were detected between the 2 groups. GO and pathway analysis showed that these differentially expressed genes were enriched in biological processes such as metabolic processes,hormone regulation,cell migration and immune function. In addition,a total of 373 049 SNPs and 42 273 Indels were detected,of which 75 363 were novel genetic mutations. By integrating the allele frequency distribution of genetic variations with the expression level of genes in the 2 groups,12 candidate functional genes affecting milk production traits were suggested,namely ASS1,CKB,GGT1,MGAM,SDSL,HP,LTF,MMP9,UPP1,DEFB4A,PGLYRP1 and MS4A8,and in them 14 important genetic variations occurred.

Key words: milk production traits, transcriptome, gene expression, functional genes, single nucleotide polymorphism

陈一丹, 张昱, 杨洁, 张勤, 姜力. 基于转录组测序的奶牛产奶性状重要功能基因挖掘[J]. 生物技术通报, 2020, 36(9): 244-252.

CHEN Yi-dan, ZHANG Yu, YANG Jie, ZHANG Qin, JIANG Li. Exploration of Key Functional Genes Affecting Milk Production Traits in Dairy Cattle Based on RNA-seq[J]. Biotechnology Bulletin, 2020, 36(9): 244-252.

参考文献

[1] 郑靖宇. 牛乳营养价值浅析[J]. 农民致富之友, 2011, 20:85.
Zheng JY.Analysis on the nutritional value of milk[J]. Nong Min Zhi Fu Zhi You, 2011, 20:85.
[2] Zhan B, Fadista J, Thomsen B, et al.Global assessment of genomic variation in cattle by genome resequencing and high-throughput genotyping[J]. BMC Genomics, 2011, 12(1):557.
[3] Suarez-Vega A, Gutierrez-gil B, Klopp C, et al. Variant discovery in the sheep milk transcriptome using RNA sequencing[J]. BMC Genomics, 2017, 18(1):1-13.
[4] Pham-quoc C, Kieu-do B, Thinh N. A high-performance FPGA-based BWA-MEM DNA sequence alignment[J]. Concurrency Computation, 2019:1-12.
[5] Li H, Handsaker B, Wysoker A, et al.The sequence alignment/map format and SAMtools[J]. Bioinformatics, 2009, 25(16):2078-2079.
[6] Danecek P, Mccarthy A.BCFtools/csq:Haplotype-aware variant consequences[J]. Bioinformatics, 2017, 33(13):2037-2039.
[7] Cingolani P, Platts A, Wang L, et al.A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff[J]. Fly, 2012, 6:80-92.
[8] Liao Y, Smyth K, Shi W.FeatureCounts:An efficient general purpose program for assigning sequence reads to genomic features[J]. Bioinformatics, 2014, 30(7):923-930.
[9] Love I, Huber W, Anders S.Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome Biology, 2014, 15(12):1-21.
[10] Jiang L, Liu J, Sun D, et al.Genome wide association studies for milk production traits in Chinese Holstein population[J]. PLoS One, 2010, 5(10):e13661.
[11] Yang J, Jiang J, Liu X, et al.Differential expression of genes in milk of dairy cattle during lactation[J]. Animal Genetics, 2016, 47(2):174-180.
[12] Cui X, Hou Y, Yang S, et al.Transcriptional profiling of mammary gland in Holstein cows with extremely different milk protein and fat percentage using RNA sequencing[J]. BMC Genomics, 2014, 15(1):226.
[13] Zheng X, Ning C, Dong Y, et al.Quantitative proteome analysis of bovine mammary gland reveals protein dynamic changes involved in peak and late lactation stages[J]. Biochemical and Biophysical Research Communications, 2017, 494(1-2):292-297.
[14] Liang R, Han B, Li Q, et al.Using RNA sequencing to identify putative competing endogenous RNAs(ceRNAs)potentially regulating fat metabolism in bovine liver[J]. Scientific Reports, 2017, 7(1):6396.
[15] Wang L, Li X, Wang J, et al.Polymorphism in PGLYRP-1 gene by PCR-RFLP and its association with somatic cell score in Chinese Holstein[J]. Research in Veterinary Science, 2013, 95(2):508-514.
[16] Engel K, Hohne W, Haberle J.Mutations and polymorphisms in the human argininosuccinate synthetase(ASS1)gene[J]. Human Mutation, 2009, 30(3):300-307.
[17] Bertrand F, Hewan-lowe K, Catrou P, et al. Abstract 3640:Potential cooperation of monoclonal antibody based targeting of gamma-glutamyl transferase(GGT1)with chemotherapy in kidney tumor xenografts[J]. Cancer Research, 2011, 71(S8):3640.
[18] Rose R, Chaudet M, Jones K.Structural studies of the intestinal α-glucosidases, maltase-glucoamylase and sucrase-isomaltase[J]. Journal of Pediatric Gastroenterology and Nutrition, 2018, 66:S11-S13.
[19] Farnaud S, Evans R.Lactoferrin - A multifunctional protein with antimicrobial properties[J]. Molecular Immunology, 2003, 40:395-405.
[20] Legrand D, Mazurier J.A critical review of the roles of host lactoferrin in immunity[J]. Biometals, 2010, 23:365-376.
[21] Bullen J, Rogers J, Leigh L.Iron-binding proteins in milk and resistance to Escherichia coli infection in infants[J]. Br Med J, 1972, 1(5792):69-75.
[22] Kawakami H, Hiratsuka M, Dosako S.Effects of iron-saturated lactoferrin on iron absorption.[J]. Agricultural and Biological Chemistry, 1988, 52(4):903-908.
[23] Li Y, Li P, Yang G, et al.Investigation and comparison of the anti-tumor activities of lactoferrin, α-lactalbumin, and β-lactoglobulin in A549, HT29, HepG2, and MDA231-LM2 tumor models[J]. Journal of Dairy Science, 2019, 102(11):9586-9597.
[24] O’halloran F, Berry P, Bahar B, et al.Polymorphisms in the bovine lactoferrin promoter are associated with reproductive performance and somatic cell count[J]. Journal of Dairy Science, 2010, 93(3):1253-1259.
[25] Mao Y, Zhu X, Xing S, et al.Polymorphisms in the promoter region of the bovine lactoferrin gene influence milk somatic cell score and milk production traits in Chinese Holstein cows[J]. Research in Veterinary Science, 2015, 103:107-112.
[26] Viale E, Tiezzi F, Maretto F, et al.Association of candidate gene polymorphisms with milk technological traits, yield, composition, and somatic cell score in Italian Holstein-Friesian sires[J]. Journal of Dairy Science, 2017, 100(9):7271-7281.
[27] Raschia A, Nani P, Maizon O, et al.Single nucleotide polymorphisms in candidate genes associated with milk yield in argentinean holstein and holstein x jersey cows[J]. Journal of Animal Science and Technology, 2018, 60(1):1-10.
[28] 涂茂林, 杨戬, 刘猛等. 乳铁蛋白检测方法研究进展[J]. 乳业科学与技术, 2015(1):32-37.
Tu MY, Yang J, Liu M, et al.Advances in detection methods for lactoferrin[J]. Journal of Dairy Science and Technology, 2015(1):32-37.
[29] 肖奕昕, 李大鹏, 赵善仓等. 乳铁蛋白生物学功能及其检测方法研究进展[J]. 中国乳品工业, 2018, 46(3):29-34.
Xiao YX, Li DP, Zhao SC, et al.Advances in biological functions and determination methods of lactoferrin[J]. China Dairy Industry, 2018, 46(3):29-34.
[30] Iglesias-figueroa B, Espinoza-scanchez E, Siqueiros T, et al. Lactoferrin as a nutraceutical protein from milk, an overview[J]. International Dairy Journal, 2018, 89:37-41.
[31] Bagnicka E, Strzalkowska N, Flisikowsi K, et al.The polymorphism in the β4-defensin gene and its association with production and somatic cell count in Holstein-Friesian cows[J]. Journal of Animal Breeding and Genetics, 2007, 124(3):150-156.
[32] Cochran D, Cole B, Null J, et al.Discovery of single nucleotide polymorphisms in candidate genes associated with fertility and production traits in Holstein cattle[J]. BMC Genetics, 2013, 14(1):1.
[33] Zoncu R, Efeyan A, Sabatini M.MTOR:From growth signal integration to cancer, diabetes and ageing[J]. Nature Reviews Molecular Cell Biology, 2011, 12(1):21-35.
[34] Yamada T, Komoto J, Kasuya T, et al.A catalytic mechanism that explains a low catalytic activity of serine dehydratase like-1 from human cancer cells:Crystal structure and site-directed mutagenesis studies[J]. Biochimica et Biophysica Acta, 2008, 1780(5):809-818.