A post-GWAS confirming the genetic effects and functional polymorphisms of AGPAT3 gene on milk fatty acids in dairy cattle

Background: People are paying more attention to the healthy and balanced diet with the improvement of their living standards. Milk fatty acids(FAs) have been reported that they were related to some atherosclerosis and coronary heart diseases in human. In our previous genome-wide association study(GWAS) on milk FAs in dairy cattle, 83 genomewide significant single nucleotide polymorphisms(SNPs) were detected. Among them, two SNPs, ARS-BFGL-NGS-109493 and BTA-56389-no-rs associated with C18 index(P = 0.0459), were located in the upstream of 1-acylglycerol-3-phosphate O-acyltransferase 3(AGPAT3) gene. AGPAT3 is involved in glycerol-lipid, glycerol-phospholipid metabolism and phospholipase D signaling pathways. Hence, it was inferred as a candidate gene for milk FAs. The aim of this study was to further confirm the genetic effects of the AGPAT3 gene on milk FA traits in dairy cattle.Results: Through re-sequencing the complete coding region, and 3000 bp of 5′ and 3′ regulatory regions of the AGPA T3 gene, a total of 17 SNPs were identified, including four in 5′ regulatory region, one in 5′ untranslated region(UTR),three in introns, one in 3′ UTR, and eight in 3′ regulatory region. By the linkage disequilibrium(LD) analysis with Haploview4.1 software, two haplotype blocks were observed that were formed by four and 12 identified SNPs,respectively. Using SAS9.2, we performed single locus-based and haplotype-based association analysis on 24 milk FAs in 1065 Chinese Holstein cows, and discovered that all the SNPs and the haplotype blocks were significantly associated with C6:0, C8:0 and C10:0(P < 0.0001–0.0384). Further, with Genomatix, we predicted that four SNPs in 5′ regulatory region(g.146702957 G > A, g.146704373 A > G, g.146704618 A > G and g.146704699 G > A) changed the transcription factor binding sites(TFBSs) for transcription factors SMARCA3, REX1, VMYB, BRACH, NKX26, ZBED4, SP1, USF1, ARNT and FOXA1. Out of them, two SNPs were validated to impact transcriptional activity by performing luciferase assay that the alleles A of both SNPs, g.146704373 A > G and g.146704618 A > G, increased the transcriptional activities of AGPAT3 promoter compared with alleles G(P = 0.0004).Conclusions: In conclusion, our findings first demonstrated the significant genetic associations of the AGPAT3 gene with milk FAs in dairy cattle, and two potential causal mutations were detected.

Identification of single nucleotide polymorphisms of PIK3R1 and DUSP1 genes and their genetic associations with milk production traits in dairy cows

Background:Previously,phosphoinositide-3-kinase regulatory subunit 1(PIK3R1)and dual specificity phosphatase 1(DUSP1)were identified as promising candidate genes for milk production traits due to their being differentially expressed between the dry period and the peak of lactation in livers of dairy cows.Hence,in this study,the single nucleotide polymorphisms(SNPs)of PIK3R1 and DUSP1 genes were identified and their genetic associations with milk yield,fat yield,fat percentage,protein yield,and protein percentage,were investigated using 1067 Chinese Holstein cows from 40 sire families.Results:By re-sequencing the entire coding region and 2000 bp of the 5′and 3′flanking regions of the two genes,one SNP in the 5′untranslated region(UTR),three in the 3′UTR,and two in the 3′flanking region of PIK3R1 were identified,and one in the 5′flanking region,one in the 3′UTR,and two in the 3′flanking region of DUSP1 were found.Subsequent single-locus association analyses showed that five SNPs in PIK3 R1,rs42590258,rs210389799,rs208819656,rs41255622,rs133655926,and rs211408208,and four SNPs in DUSP1,rs207593520,rs208460068,rs209154772,and rs210000760,were significantly associated with milk,fat and protein yields in the first or second lactation(P values≤0.0001 and 0.0461).In addition,by the Haploview 4.2 software,the six and four SNPs in PIK3R1 and DUSP1 respectively formed one haplotype block,and the haplotype-based association analyses showed significant associations between their haplotype combinations and the milk traits in both two lactations(P values≤0.0001 and 0.0364).One SNP,rs207593520(T/G),was predicted to alter the transcription factor binding sites(TFBSs)in the 5′flanking region of DUSP1.Further,the dual-luciferase assay showed that the transcription activity of allele T in rs207593520 was significantly higher than that of allele G,suggesting the activation of transcriptional activity of DUSP1 gene by allele T of rs207593520.Thus,the rs207593520 SNP was highlighted as a potential causal mutation that should be further verified.Conclusions:We demonstrated novel and significant genetic effects of the PIK3R1 and DUSP1 genes on milk production traits in dairy cows,and our findings provide information for use in dairy cattle breeding.

Weighted single-step GWAS and RNA sequencing reveals key candidate genes associated with physiological indicators of heat stress in Holstein cattle

Background: The study of molecular processes regulating heat stress response in dairy cattle is paramount for developing mitigation strategies to improve heat tolerance and animal welfare. Therefore, we aimed to identify quantitative trait loci(QTL) regions associated with three physiological indicators of heat stress response in Holstein cattle, including rectal temperature(RT), respiration rate score(RS), and drooling score(DS). We estimated genetic parameters for all three traits. Subsequently, a weighted single-step genome-wide association study(WssGWAS) was performed based on 3200 genotypes, 151,486 phenotypic records, and 38,101 animals in the pedigree file. The candidate genes located within the identified QTL regions were further investigated through RNA sequencing(RNA-seq) analyses of blood samples for four cows collected in April(non-heat stress group) and four cows collected in July(heat stress group).Results: The heritability estimates for RT, RS, and DS were 0.06, 0.04, and 0.03, respectively. Fourteen, 19, and 20 genomic regions explained 2.94%, 3.74%, and 4.01% of the total additive genetic variance of RT, RS, and DS, respectively. Most of these genomic regions are located in the Bos taurus autosome(BTA) BTA3, BTA6, BTA8, BTA12, BTA14, BTA21, and BTA24. No genomic regions overlapped between the three indicators of heat stress, indicating the polygenic nature of heat tolerance and the complementary mechanisms involved in heat stress response. For the RNA-seq analyses, 2627 genes were significantly upregulated and 369 downregulated in the heat stress group in comparison to the control group. When integrating the WssGWAS, RNA-seq results, and existing literature, the key candidate genes associated with physiological indicators of heat stress in Holstein cattle are: PMAIP1, SBK1, TMEM33, GATB, CHORDC1, RTN4IP1, and BTBD7.Conclusions: Physiological indicators of heat stress are heritable and can be improved through direct selection. Fifty-three QTL regions associated with heat stress indicators confirm the polygenic nature and complex genetic determinism of heat tolerance in dairy cattle. The identified candidate genes will contribute for optimizing genomic evaluation models by assigning higher weights to genetic markers located in these regions as well as to the design of SNP panels containing polymorphisms located within these candidate genes.

Effects of feeding untreated, pasteurized and acidified waste milk and bunk tank milk on the performance, serum metabolic profiles, immunity, and intestinal development in Holstein calves

Background: The present experiment was performed to assess the effects of different sources of milk on the growth performance, serum metabolism, immunity, and intestinal development of calves. Eighty-four Holstein male neonatal calves were assigned to one of the following four treatment groups: those that received bunk tank milk(BTM), untreated waste milk(UWM), pasteurized waste milk(PWM), and acidified waste milk(AWM) for 21 d.Results: Calves in the BTM and AWM groups consumed more starter(P < 0.05) than those in the UWM group.Average daily gain in the UWM group was the highest(P < 0.05). Calves exhibited the highest(P < 0.05) serum total protein, albumin, total cholesterol, high density lipoprotein, triglycerides, growth hormone, immunoglobulin(Ig) A and Ig M concentrations in the UWM group, highest malondialdehyde and tumor necrosis factor-α in the PWM group(P < 0.05), and highest glutathione peroxidase and Ig G in the BTM group(P < 0.05). The jejunum and ileum of the calves in all treatments presented a slight inflammatory response. The jejunal inflammation scores were higher(P < 0.05) in the UWM and AWM groups than the BTM group; the ileal inflammation scores increased more(P < 0.05) in the AWM group than the BTM group. Jejunal immunohistochemical scores(IHS) were higher(P < 0.05) in the PWM and AWM groups than the BTM group. Compared to the other three groups, calves feeding on BTM had lower(P < 0.05) ileal IHS. Jejunal interleukin(IL)-1β, IL-8, and IL-10 mRNA expression in the UWM group was the highest(P < 0.05). Calves fed AWM increased(P < 0.05) mRNA expression of IL-8 and toll like receptor 4(TLR-4) in the jejunum and IL-8, IL-1β, IL-6, IL-8, and IL-10 in the mesenteric lymph nodes.Conclusions: Overall, bunk tank milk is the best choice for calf raising compared to waste milk. The efficiency of feeding pasteurized and acidified waste milk are comparable, and the acidification of waste milk is an acceptable labor-saving and diarrhea-preventing feed for young calves.