水牛MFSD2A基因多态性及其与产奶性状的关联分析

doi:10.16431/j.cnki.1671-7236.2022.10.018

摘要/Abstract

摘要： 【目的】研究水牛主要促进因子超家族成员2a(major facilitator superfamily domain containing 2a,MFSD2A)基因多态性,并筛选与水牛产奶性状相关联的SNP。【方法】采集383头健康水牛血液DNA,利用PCR扩增和Sequenom MassARRAY飞行时间质谱技术对SNP位点进行筛选及基因型检测,并进行遗传多样性、连锁不平衡(LD)和单倍型分析,以及对不同基因型和单倍型与水牛产奶性状进行关联分析。【结果】在水牛MFSD2A基因中检出7个SNPs位点,多态信息含量(PIC)均在0.25~0.40之间,属于中度多态,除g.8290 T＞C位点外,其余位点均处于Hardy-Weinberg平衡状态(P＞0.05);关联分析结果表明,这7个SNPs位点与305 d产奶量均显著关联(P＜0.05),且突变杂合子基因型的305 d产奶量是3种基因型中最低的,g.568 C＞G和g.701 A＞G位点与乳脂率显著关联(P＜0.05),乳脂率从大到小依次为突变纯合子＞野生纯合子＞突变杂合子,g.568 C＞G、g.3326 G＞A、g.8290 T＞C和g.8523 C＞G位点与乳蛋白率显著关联(P＜0.05),7个SNPs位点与乳尿素氮关联均不显著(P＞0.05);连锁不平衡和单倍型分析结果明,g.568 C＞G、g.701 A＞G和g.3203 T＞G位点可组成一个单倍型模块(Block 1),g.8523 C＞G和g.9138 G＞T位点可组成另一个单倍型模块(Block 2)。其中g.568 C＞G和g.701 A＞G、g.568 C＞G和g.3326 G＞A、g.8523 C＞G和g.9138 G＞T位点的r²分别为0.87、0.37和0.48,处于强连锁不平衡状态。Block 1中的H1(CAG)单倍型乳蛋白率显著高于H2(CAT)和H3(GCT)(P＜0.05),Block 2中的D3(GG)单倍型305 d产奶量极显著高于D1(CG)和D2(CT)(P＜0.01)。【结论】筛选出的MFSD2A基因7个SNPs位点与305 d产奶量均显著关联,g.568 C＞G和g.701 A＞G位点与乳脂率显著关联,g.568 C＞G、g.701 A＞G、g.3326 G＞A、g.8290 T＞C和g.8523 C＞G位点与乳蛋白率显著关联,这7个SNPs位点与乳尿素氮均无显著关联。纯合子的基因型在高305 d产奶量和乳脂率的选择上更有优势,H1(CAG)和D3(GG)单倍型是提高水牛305 d产奶量和乳蛋白率的有利单倍型。

关键词: 水牛; MFSD2A基因; 连锁不平衡; 产奶性状; 关联分析

Abstract: 【Objective】 In order to study the polymorphism of major facilitator superfamily domain containing 2a (MFSD2A) gene of buffalo and screen the SNP associated with milk production traits of buffalo.【Method】 The blood DNA was collected from 383 healthy buffaloes,and the SNP were screened and genotyped by PCR amplification and Se-quenom MassARRAY SNP time-of-flight mass spectrometry.The genetic diversity,linkage disequilibrium (LD) and haplotype analysis were analyzed,and the association between different genotypes and haplotypes and milk production traits was analyzed.【Result】 Seven SNPs were detected in MFSD2A gene of buffalo,and the polymorphism information content (PIC) ranged from 0.25 to 0.40,which belonged to moderate polymorphism.Except g.8290 T＞C,all the other sites were in Hardy-Weinberg equilibrium(P＞0.05).The results of association analysis showed that the seven SNPs were significantly associated with 305 d milk yield (P＜0.05),and the 305 d milk yield of mutant heterozygote was the lowest among the three genotypes.g.568 C＞G and g.701 A＞G were significantly associated with milk fat percentage (P＜0.05).The order of milk fat percentage genotypes was mutant homozygote＞wild homozygote＞mutant heterozygote.g.568 C＞G,g.3326 G＞A,g.8290 T＞C and g.8523 C＞G were significantly associated with milk protein percentage (P＜0.05),but these seven SNPs were not significantly associated with milk urea nitrogen (P＞0.05).The results of LD and haplotype analysis showed that g.568 C＞G,g.701 A＞G and g.3203 T＞G could form a haplotype module (Block 1),while g.8523 C＞G and g.9138 G＞T could form another haplotype module (Block 2).The r² of g.568 C＞G and g.701 A＞G,g.568 C＞G and g.3326 G＞A,g.8523 C＞G and g.9138 G＞T were 0.87,0.37 and 0.48,respectively,which were in a state of strong linkage disequilibrium.The milk protein percentage of H1 (CAG) haplotype in Block 1 was significantly higher than that of H2 (CAT) and H3 (GCT) haplotype (P＜0.05),and the 305 d milk yield of D3 (GG) haplotype in Block 2 was significantly higher than that of D1 (CG) and D2 (CT) haplotype (P＜0.01).【Conclusion】 The seven SNPs of MFSD2A gene were significantly associated with 305 d milk yield,g.568 C＞G and g.701 A＞G were significantly associated with milk fat percentage,g.568 C＞G,g.701 A＞G,g.3326 G＞A,g.8290 T＞C and g.8523 C＞G were significantly associated with milk protein percentage were significantly associated,and none of these seven SNPs were significantly associated with milk urea nitrogen.Homozygous genotypes had more advantages in the selection of high 305 d milk yield and milk fat percentage.H1 (CAG) and D3 (GG) haplotypes were favorable haplotypes to increase 305 d milk yield and milk protein rate in buffalo.

Key words: buffalo; MFSD2A gene; linkage disequilibrium; milk production traits; association analysis

中图分类号:

S823.8⁺3

梁莎莎, 庞春英, 邓廷贤, 陆杏蓉, 段安琴, 马小娅, 方艳艳, 梁贤威. 水牛MFSD2A基因多态性及其与产奶性状的关联分析[J]. 中国畜牧兽医, 2022, 49(10): 3856-3865.

LIANG Shasha, PANG Chunying, DENG Tingxian, LU Xingrong, DUAN Anqin, MA Xiaoya, FANG Yanyan, LIANG Xianwei. Polymorphism of MFSD2A Gene and Its Association with Milk Production Traits in Buffalo[J]. China Animal Husbandry and Veterinary Medicine, 2022, 49(10): 3856-3865.

参考文献

[1] ZICARELLI L.Buffalo milk:Its properties,dairy yield and mozzarella production[J].Veterinary Research Communications,2004,28(Suppl 1):127-135.
[2] ETHAYATHULLA A S,YOUSEF M S,AMIN A,et al.Structure-based mechanism for Na⁺/melibiose symport by MelB[J].Nature Communications,2014,5:3009.
[3] 王文广,匡德宣,李娜,等.树鼩Mfsd2a基因的克隆分析和不同组织表达量的检测[J].实验动物与比较医学,2019,39(3):178-186. WANG W G,KUANG D X,LI N,et al.Cloning and analysis of tree shrew Mfsd2a gene and detection of its expression in different tissues[J].Laboratory Animal and Comparative Medicine,2019,39(3):178-186.(in Chinese)
[4] SANCHEZ-CAMPILLO M,RUIZ-PALACIOS M,RUIZ-ALCARAZ A J,et al.Child head circumference and placental MFSD2a expression are associated to the level of MFSD2a in maternal blood during pregnancy[J].Front Endocrinol (Lausanne),2020,11:38.
[5] NGUYEN L N,MA D,SHUI G,et al.Mfsd2a is a transporter for the essential omega-3 fatty acid docosahexaenoic acid[J].Nature,2014,509(7501):503-506.
[6] YAN N.Structural biology of the major facilitator superfamily transporters[J].Annual Review of Biophysics,2015,44:257-283.
[7] BERGER J H,CHARRON M J,SILVER D L.Major facilitator superfamily domain-containing protein 2a (MFSD2A) has roles in body growth,motor function,and lipid metabolism[J].PLoS One,2012,7(11):e50629.
[8] 章纯熙.中国水牛科学[M].南宁:广西科学技术出版社,2000. ZHANG C X.Science and Technology in Chinese Buffaloes[M].Nanning:Guangxi Science & Technology Publishing House,2000.(in Chinese)
[9] 杨龙.利用cDNA/EST序列大规模开发内含子多态性标记的研究[D].杭州:浙江大学,2008. YANG L.Studies on large-scale development of intron polymorphism markers using cDNA/EST sequences[D].Hangzhou:Zhejiang University,2008.(in Chinese)
[10] JEONG Y M,MUN J H,LEE I,et al.Distinct roles of the first introns on the expression of Arabidopsis profilin gene family members[J].Plant Physiology,2006,140(1):196-209.
[11] 唐馨,苟萍.内含子的功能[J].生命的化学,2019,39(4):772-777. TANG X,GOU P.The introns function[J].Chemistry of Life,2019,39(4):772-777.(in Chinese)
[12] SLATKIN M.Linkage disequilibrium--understanding the evolutionary past and mapping the medical future[J].Nature Reviews Genetics,2008,9(6):477-485.
[13] DE LA VEGA F M,DAILEY D,ZIEGLE J,et al.New generation pharmacogenomic tools:A SNP linkage disequilibrium map,validated SNP assay resource,and high-throughput instrumentation system for large-scale genetic studies[J].Biotechniques,2002,Suppl:48-50.
[14] SVED J A,HILL W G.One hundred years of linkage disequilibrium[J].Genetics,2018,209(3):629-636.
[15] LU D,SARGOLZAEI M,KELLY M,et al.Linkage disequilibrium in Angus,Charolais,and crossbred beef cattle[J].Frontiers in Genetics,2012,3:152.
[16] ARDLIE K G,KRUGLYAK L,SEIELSTAD M.Patterns of linkage disequilibrium in the human genome[J].Nature Reviews Genetics,2002,3(4):299-309.
[17] MOKRY F B,BUZANSKAS M E,DE ALVARENGA M M,et al.Linkage disequilibrium and haplotype block structure in a composite beef cattle breed[J].BMC Genomics,2014,15(Suppl 7):S6.
[18] LI X,LI J.Haplotype inference[J].Methods in Molecular Biology,2012,850:411-421.
[19] 许金根,赵磊,顾有方,等.猪CD4基因多态性及生物信息学分析[J].基因组学与应用生物学,2018,37(11):4736-4741. XU J G,ZHAO L,GU Y F,et al.Polymorphisms and bioinformatics analysis of porcine CD4 gene[J].Genomics and Applied Biology,2018,37(11):4736-4741.(in Chinese)
[20] PAO S S,PAULSEN I T,SAIER M J.Major facilitator superfamily[J].Microbiology and Molecular Biology Reviews,1998,62(1):1-34.
[21] REDDY V S,SHLYKOV M A,CASTILLO R,et al.The major facilitator superfamily (MFS) revisited[J].The FEBS Journal,2012,279(11):2022-2035.
[22] HUANG B,LI X.The role of Mfsd2a in nervous system diseases[J].Frontiers Neuroscience,2021,15:730534.
[23] QUEK D Q,NGUYEN L N,FAN H,et al.Structural insights into the transport mechanism of the human sodium-dependent lysophosphatidylcholine transporter MFSD2A[J]. The Journal of Biological Chemistry,2016,291(18):9383-9394.