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

【目的】研究水牛主要促进因子超家族成员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~2分别为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产奶量和乳蛋白率的有利单倍型。