Effects of SNPs in AANAT and ASMT Genes on Milk and Peripheral Blood Melatonin Concentrations in Holstein Cows (Bos taurus)

Aralkylamine N-acetyltransferase (AANAT) and acetylserotonin O-methyltransferase (ASMT), the two rate-limiting enzymes for melatonin synthesis, regulate melatonin production in mammals. Through analysis of the milk melatonin level and dairy herd improvement (DHI) index, it was found that the melatonin concentration in milk was significantly negatively correlated with the 305 day milk yield (305M) and peak milk yield (PeakM) (p < 0.05), while it was significantly positively correlated with the serum melatonin concentration (p < 0.05). The full-length of AANAT and ASMT were sequenced and genotyped in 122 cows. Three SNPs in AANAT and four SNPs in ASMT were significantly related to MT levels in the milk and serum (p < 0.05). The SNPs in AANAT were temporarily denoted as N-SNP1 (g.55290169 T>C), N-SNP2 (g.55289357 T>C), and N-SNP3 (g.55289409 C>T). The SNPs in ASMT were temporarily denoted as M-SNP1 (g.158407305 G>A), M-SNP2 (g.158407477 A>G), M-SNP3 (g.158407874 G>A), and M-SNP4 (g.158415342 T>C). The M-SNP1, M-SNP2, and M-SNP3 conformed to the Hardy−Weinberg equilibrium (p > 0.05), while other SNPs deviated from the Hardy−Weinberg equilibrium (p < 0.05). The potential association of MT production and each SNP was statistically analyzed using the method of linkage disequilibrium (LD). The results showed that N-SNP2 and N-SNP3 had some degree of LD (D′ = 0.27), but M-SNP1 and M-SNP2 had a strong LD (D′ = 0.98). Thus, the DHI index could serve as a prediction of the milk MT level. The SNPs in AANAT and ASMT could be used as potential molecular markers for screening cows to produce high melatonin milk.