Identification of Genetic Associations and Functional Polymorphisms of SAA1 Gene Affecting Milk Production Traits in Dairy Cattle

Our initial RNA sequencing (RNA-seq) revealed that the Serum amyloid A1 (SAA1) gene was differentially expressed in the mammary glands of lactating Holstein cows with extremely high versus low phenotypic values of milk protein and fat percentage. To further validate the genetic effect and potential molecular mechanisms of SAA1 gene involved in regulating milk production traits in dairy cattle, we herein performed a study through genotype-phenotype associations. Six identified SNPs were significantly associated with one or more milk production traits (0.00002< P < 0.0025), providing additional evidence for the potential role of SAA1 variants in milk production traits in dairy cows. Subsequently, both luciferase assay and electrophoretic mobility shift assay (EMSA) clearly demonstrated that the allele A of g.-963C>A increased the promoter activity by binding the PARP factor while allele C did not. Bioinformatics analysis indicated that the secondary structure of SAA protein changed by the substitution A/G in the locus c. +2510A>G. Our findings were the first to reveal the significant associations of the SAA1 gene with milk production traits, providing basis for further biological function validation, and two identified SNPs, g.-963C>A and c. +2510A>G, may be considered as genetic markers for breeding in dairy cattle.