As a member of the fatty acid desaturase family, fatty acid desaturase 2 (FADS2) gene is a rate‐limiting enzyme in the synthesis of unsaturated fatty acids and within/near to the reported QTL regions for milk‐production traits. We previously found that FADS2 is differentially expressed during different lactations of Chinese Holstein cows, and participates in lipid metabolic processes by influencing the insulin, PI3K‐Akt, MAPK, AMPK, mTOR and PPAR signaling pathways. Therefore, we considered this gene as a candidate gene for milk‐production traits. In this study, we identified 12 SNPs in FADS2 by re‐sequencing, including two SNPs in the 5′ flanking region, one in the seventh exon, five in introns, two in the 3′ untranslated region and two in the 3′ flanking region. The 29:g.40378819C>T is a missense mutation that causes alanine (GCG) to be replaced with valine (GTG). Through single marker association analysis, we found that all of the 12 SNPs were significantly associated with 305 day milk yield, fat yield, fat percentage, protein yield or protein percentage (p < 0.0493). The results of the subsequent haplotype association analysis also confirmed the associations between the gene and milk‐production traits. In summary, this study suggests that there is a significant genetic association between FADS2 and milk‐production traits, and that the SNPs with significant genetic effects can provide important molecular information for the development of a genomic selection chip in dairy cattle.

Background Our preliminary work confirmed that, SLC22A7 (solute carrier family 22 member 7), NGFR (nerve growth factor receptor), ARNTL (aryl hydrocarbon receptor nuclear translocator like) and PPP2R2B (protein phosphatase 2 regulatory subunit Bβ) genes were differentially expressed in dairy cows during different stages of lactation, and involved in the lipid metabolism through insulin, PI3K-Akt, MAPK, AMPK, mTOR, and PPAR signaling pathways, so we considered these four genes as the candidates affecting milk production traits. In this study, we detected polymorphisms of the four genes and verified their genetic effects on milk yield and composition traits in a Chinese Holstein cow population. Results By resequencing the whole coding region and part of the flanking region of SLC22A7, NGFR, ARNTL and PPP2R2B, we totally found 20 SNPs, of which five were located in SLC22A7, eight in NGFR, three in ARNTL, and four in PPP2R2B. Using Haploview4.2, we found three haplotype blocks including five SNPs in SLC22A7, eight in NGFR and three in ARNTL. Single-SNP association analysis showed that 19 out of 20 SNPs were significantly associated with at least one of milk yield, fat yield, fat percentage, protein yield or protein percentage in the first and second lactations (P < 0.05). Haplotype-based association analysis showed that the three haplotypes were significantly associated with at least one of milk yield, fat yield, fat percentage, protein yield or protein percentage (P < 0.05). Further, we used SOPMA software to predict a SNP, 19:g.37095131C > T in NGFR, changed the structure of NGFR protein. In addition, we used Jaspar software to found that four SNPs, 19:g.37113872C > G,19:g.37113157C > T, and 19:g.37112276C > T in NGFR and 15:g.39320936A > G in ARNTL, could change the transcription factor binding sites and might affect the expression of the corresponding genes. These five SNPs might be the potential functional mutations for milk production traits in dairy cattle. Conclusions In summary, we proved that SLC22A7, NGFR, ARNTL and PPP2R2B have significant genetic effects on milk production traits. The valuable SNPs can be used as candidate genetic markers for genomic selection of dairy cattle, and the effects of these SNPs on other traits need to be further verified.

SEC13 homolog, nuclear pore and COPII coat complex component (SEC13) is the core component of the cytoplasmic COPII complex, which mediates material transport from the endoplasmic reticulum to the Golgi complex. Our preliminary work found that SEC13 gene was differentially expressed in dairy cows during different stages of lactation, and involved in metabolic pathways of milk synthesis such as citric acid cycle, fatty acid, starch and sucrose metabolisms, so we considered that the SEC13 might be a candidate gene affecting milk production traits. In this study, we detected the polymorphisms of SEC13 gene and verified their genetic effects on milk yield and composition traits in a Chinese Holstein cow population. By sequencing the whole coding and partial flanking regions of SEC13, we found four single nucleotide polymorphisms (SNPs). Subsequent association analysis showed that these four SNPs were significantly associated with milk yield, fat yield, protein yield or protein percentage in the first and second lactations (p ≤.0351). We also found that two SNPs in SEC13 formed one haplotype block by Haploview4.2, and the block was significantly associated with milk yield, fat yield, fat percentage, protein yield or protein percentage (p ≤ .0373). In addition, we predicted the effect of SNP on 5′region on transcription factor binding sites (TFBSs), and found that the allele A of 22:g.54362761A&gt;G could bind transcription factors (TFs) GATA5, GATA3, HOXD9, HOXA10, CDX1 and Hoxd13; and further dual-luciferase reporter assay verified that the allele A of this SNP inhibited the fluorescence activity. We speculate that the A allele of 22:g.54362761A&gt;G might inhibit the transcriptional activity of SEC13 gene by binding the TFs, which may be a cause mutation affecting the formation of milk production traits in dairy cows. In summary, we proved that SEC13 has a significant genetic effect on milk production traits and the identified significant SNPs could be used as candidate genetic markers for GS SNP chips development; on the other hand, we verified the transcriptional regulation of 22:g.54362761A&gt;G on SEC13 gene, providing research direction for further function validation tests.