A post-GWAS confirming the SCD gene associated with milk medium- and long-chain unsaturated fatty acids in Chinese Holstein population

Phenotypic variation in milk fatty acid composition and its association with stearoyl-CoA desaturase 1 (SCD1) gene polymorphisms in Gir cows

Individual variation in milk fatty acid (FA) composition has been partially attributed to stearoyl-CoA desaturase 1 (SCD1) gene polymorphisms in taurine breeds, but much less is known for Zebu breeds. This study investigated the phenotypic variation in milk FA composition, and the influence of SCD1 variants on this trait and on milk fat desaturase indices (DI) in Gir cows. The functional impact of SCD1 variants was predicted using bioinformatics tools. Milk and blood samples were collected from 312 cows distributed in 10 herds from five states of Brazil. SCD1 variants were identified through target sequencing, and milk FA composition was determined by gas chromatography. Phenotypic variation in milk FA composition fell within the range reported for taurine breeds, with SCD18 index showing the lowest variation among the DI. Fourteen SCD1 variants were identified, six of which not previously described. Regarding the A293V polymorphism, all cows were homozygous for the C allele (coding for alanine), whereas all genotypes were detected for the second SNP affecting the 293 codon (G > A), with compelling evidence for functional effects. Significant associations (based on raw p-values) were found between this SNP and C12:0, cis-9, trans-11 CLA and short-chain FA, and between another SNP (rs523411937) and C15:0 and odd-chain linear FA. A new SNP on Chr26:21277069 was associated with trans-11 C18:1, cis-9, trans-11 CLA, C18:3 n-3 and n-3 FA. These findings indicate that SCD1 polymorphisms also contributes to the phenotypic variation in milk FA composition of Gir cows, with potential use in their breeding programmes.

Diversity of Mitochondrial DNA Haplogroups and Their Association with Bovine Antral Follicle Count

Maternal origins based on the bovine mitochondrial D-loop region are proven to have two main origins: Bos taurus and Bos indicus. To examine the association between the maternal origins of bovine and reproductive traits, the complete mitochondrial D-loop region sequences from 501 Chinese Holstein cows and 94 individuals of other breeds were analyzed. Based on the results obtained from the haplotype analysis, 260 SNPs (single nucleotide polymorphism), 32 indels (insertion/deletion), and 219 haplotypes were identified. Moreover, the nucleotide diversity (π) and haplotype diversity (Hd) were 0.024 ± 0.001 and 0.9794 ± 0.003, respectively, indicating the abundance of genetic resources in Chinese Holstein cows. The results of the median-joining network analysis showed two haplogroups (HG, including HG1 and HG2) that diverged in genetic distance. Furthermore, the two haplogroups were significantly (p < 0.05) correlated with the antral follicle (diameter ≥ 8 mm) count, and HG1 individuals had more antral follicles than HG2 individuals, suggesting that these different genetic variants between HG1 and HG2 correlate with reproductive traits. The construction of a neighbor-joining phylogenetic tree and principal component analysis also revealed two main clades (HG1 and HG2) with different maternal origins: Bos indicus and Bos taurus, respectively. Therefore, HG1 originating from the maternal ancestors of Bos indicus may have a greater reproductive performance, and potential genetic variants discovered may promote the breeding process in the cattle industry.

Identification of genetic effects and potential causal polymorphisms of CPM gene impacting milk fatty acid traits in Chinese Holstein

Our previous GWAS revealed 83 significant SNPs and 20 promising candidate genes associated with milk fatty acid traits in dairy cattle. Out of them, the carboxypeptidase M (CPM) gene contains a genome-wide significant SNP, Hapmap49848-BTA-106779, which is strongly associated with myristic acid (C14:0; P = 0.0064). Herein, we aimed to confirm the genetic effects of CPM on milk fatty acids in Chinese Holstein. Seven SNPs were detected by re-sequencing the sequences of entire exons and 3000 bp of up-/downstream flanking regions of the CPM gene, of which three were in 5' flanking region, one in the 3' UTR and three were in the 3' flanking region. Using the Haploview 4.1, we estimated the LD among the identified SNPs and found two haplotype blocks. With the animal model, we performed the SNP- and haplotype-based association analyses, and observed that these SNPs and haplotype blocks mainly had strong genetic associations with medium-chain saturated fatty acids (caproic acid, C6:0; caprylic acid, C8:0; capric acid, C10:0; and lauric acid, C12:0) (P < 0.0001-0.0257). In addition, using the Genomatix software, we predicted that three SNPs in the 5' flanking region of CPM (g.45079507A>G, g.45080228C>A and g.45080335C>G) changed the transcription factor binding sites for PREF (progesterone receptor biding site), ZBRK1 (transcription factor with eight central zinc fingers and an N-terminal KRAB domain), SOX9 (sex-determining region Y-box 9, dimeric binding sites), SOX6 (sex-determining region Y-box 6) and FOXP1-ES (alternative splicing variant of FOXP1, activated in ESCs). Further, the dual-luciferase reporter assay showed these three SNPs altered the transcriptional activity of CPM gene (P ≤ 0.0006). In summary, using the post-GWAS strategy, we first confirmed the significant genetic effects of CPM with milk fatty acids in dairy cattle, and identified three potential causal mutations.

A targeted genotyping approach to enhance the identification of variants for lactation persistency in dairy cows

Lactation persistency (LP), defined as the ability of a cow to maintain milk production at a high level after milk peak, is an important phenotype for the dairy industry. In this study, we used a targeted genotyping approach to scan for potentially functional single nucleotide polymorphisms (SNPs) within 57 potential candidate genes derived from our previous genome wide association study on LP and from the literature. A total of 175,490 SNPs were annotated within 10-kb flanking regions of the selected candidate genes. After applying several filtering steps, a total of 105 SNPs were retained for genotyping using target genotyping arrays. SNP association analyses were performed in 1,231 Holstein cows with 69 polymorphic SNPs using the univariate liner mixed model with polygenic effects using DMU package. Six SNPs including rs43770847, rs208794152, and rs208332214 in ADRM1; rs209443540 in C5orf34; rs378943586 in DDX11; and rs385640152 in GHR were suggestively significantly associated with LP based on additive effects and associations with 4 of them (rs43770847, rs208794152, rs208332214, and rs209443540) were based on dominance effects at P < 0.05. However, none of the associations remained significant at false discovery rate adjusted P (FDR) < 0.05. The additive variances explained by each suggestively significantly associated SNP ranged from 0.15% (rs43770847 in ADRM1) to 5.69% (rs209443540 in C5orf34), suggesting that these SNPs might be used in genetic selection for enhanced LP. The percentage of phenotypic variance explained by dominance effect ranged from 0.24% to 1.35% which suggests that genetic selection for enhanced LP might be more efficient by inclusion of dominance effects. Overall, this study identified several potentially functional variants that might be useful for selection programs for higher LP. Finally, a combination of identification of potentially functional variants followed by targeted genotyping and association analysis is a cost-effective approach for increasing the power of genetic association studies.

Identification of genetic associations of ECHS1 gene with milk fatty acid traits in dairy cattle

Our previous genome-wide association study identified 83 genome-wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl-CoA hydratase, short chain 1 (ECHS1) and enoyl-CoA hydratase and 3-hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re-sequencing the full-length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3' flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single-locus association analysis showed that it was significantly associated with three milk fatty acids (P = 0.0002-0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (P = <0.0001-0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (P = <0.0001-0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long-chain unsaturated and medium-chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in-depth validation.

A post-GWAS confirming effects of PRKG1 gene on milk fatty acids in a Chinese Holstein dairy population

BACKGROUND

We previously conducted a genome-wide association study (GWAS) strategy for milk fatty acids in Chinese Holstein, and identified 83 genome-wide significant single nucleotide polymorphisms (SNPs) and 314 suggestive significant SNPs. Among them, two SNPs, BTB-01077939 and BTA-11275-no-rs associated with C10:0, C12:0, and C14 index (P = 0.000014 ~ 0.000024), were within and close to (0.85 Mb) protein kinase, cGMP-dependent, type І (PRKG1) gene on BTA26, respectively. PRKG1 gene plays a key role in lipolysis to release fatty acids and glycerol through the hydrolysis of triacyglycerol in adipocytes. We herein considered it as a promising candidate for milk fatty acids. The purpose of this study was to investigate whether PRKG1 had effects on milk fatty acids.

RESULTS

By direct sequencing the PCR products of pooled DNA, we identified a total of six SNPs, including one in 5' flanking region, four in 3' untranslated region (UTR), and one in 3' flanking region. The single-locus association analysis was carried out, and showed that the six SNPs mainly had significant associations with C6:0, C8:0 and C17:1 (P < 0.0001 ~ 0.0035). In addition, we observed a haplotype block formed by g.6903810G > A and g.6904047G > T with Haploview 4.1, and it was strongly associated with C8:0, C10:0, C16:1, C17:1, C20:0 and C16 index (P = < 0.0001 ~ 0.0123). The SNP, g.8344262A > T, was predicted to alter the binding site (BS) of transcription factor (TF) GAGA box with Genomatix software, and the subsequent luciferase assay verified that it really changed the transcriptional activity of PRKG1 gene (P = 0.0009).

CONCLUSION

In conclusion, to our best of knowledge, we are the first who identified the significant effects of PRKG1 on milk fatty acids in dairy cattle.

Determination of Genetic Effects of LIPK and LIPJ Genes on Milk Fatty Acids in Dairy Cattle

In our previous genome-wide association study (GWAS) on milk fatty acids (FAs) in Chinese Holstein, we discovered 83 genome-wide significant single nucleotide polymorphisms (SNPs) associated with milk FAs. Two of them were close to lipase family member K (LIPK) and lipase family member J (LIPJ), respectively. Hence, this study is a follow-up to verify whether the LIPK and LIPJ have significant genetic effects on milk FAs in dairy cattle. By re-sequencing the entire exons, and 3 kb of 5′ and 3′ flanking regions, two and seven SNPs were identified in LIPK and LIPJ, respectively, including a novel SNP, ss158213049726. With the Haploview 4.1 software, we found that five of the SNPs in LIPJ formed a haplotype block (D′ = 0.96 ~ 1.00). Single-locus association analyses revealed that each SNP in LIPK and LIPJ was significantly associated with at least one milk FA (p = < 1.00 × 10⁻⁴ ~ 4.88 × 10⁻²), and the haplotype-based association analyses showed significant genetic effects on nine milk FAs (p = < 1.00 × 10⁻⁴ ~ 3.98 × 10⁻²). Out of these SNPs, the missense mutation in LIPK gene, rs42774527, could change the protein secondary structure and function predicted by SOPMA, SIFT, and PROVEAN softwares. With the Genomatix software, we predicted that two SNPs, rs110322221 in LIPK and rs211373799 in LIPJ, altered the transcription factors binding sites (TFBSs), indicating their potential regulation on promoter activity of the genes. Furthermore, we found that both LIPK and LIPJ had relatively high expressions in the mammary gland. In conclusion, our research is the first to demonstrate that LIPK and LIPJ genes have significant associations with milk FAs, and the identified SNPs might be served as genetic markers to optimize breeding programs for milk FAs in dairy cattle. This research deserves in-depth verification.