本研究旨在检测水牛STAT1基因的单核苷酸多态性(SNP),探究中国水牛多态性位点的群体遗传特征,寻找与产奶性状相关的分子标记。采用PCR和测序法筛选水牛STAT1基因序列的SNPs,以357头水牛为研究对象,采用基质辅助激光解吸电离飞行时间质谱法(MALDI-TOF-MS)对群体进行了基因型检测,运用SAS(9.3)软件一般线性模型(GLM)对STAT1基因型与产奶性状的关联程度进行统计分析。结果表明,在水牛STAT1基因中检测出3个突变位点,分别位于5′UTR(g.68G>A)、外显子10(g.986G>T)和3′UTR(g.2881C>A)。关联分析结果显示,水牛STAT1基因g.986G>T位点与305天产奶量和乳脂率显著关联(P<0.05)。Bonferroni t检验多重比较结果显示,基因型TT和GG个体的305天产奶量显著高于GT基因型(P<0.05),基因型TT个体的乳脂率显著高于GG和GT基因型(P<0.05),且305天的产奶量随着胎次的增加呈现上升趋势,直到第3或4胎次时达到最大值,表明胎次是影响水牛产奶量的重要环境因素之一。本研究表明,水牛STAT1基因的g.986G>T位点可能是影响其产奶性状的候选分子遗传标记,为今后建立中国水牛标记辅助选择育种研究奠定基础。