Exploration of Key Functional Genes Affecting Milk Production Traits in Dairy Cattle Based on RNA-seq

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

Abstract

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

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.

References

[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.