A QTL with major effect on milk yield and composition maps to bovine chromosome 14

Dynamic Profile of the Yak Mammary Transcriptome during the Lactation Cycle

The objective of this study was to assess the transcriptome of the mammary tissue of four yaks during the whole lactation cycle. For this purpose, biopsies of the mammary gland were performed at -30, -15, 1, 15, 30, 60, 120, 180, and 240 days relative to parturition (d). The transcriptome analysis was performed using a commercial bovine microarray platform and the results were analyzed using several bioinformatic tools. The statistical analysis using an overall false discovery rate ≤ 0.05 for the effect of whole lactation and p < 0.05 for each comparison identified >6000 differentially expressed genes (DEGs) throughout lactation, with a large number of DEGs observed at the onset (1 d vs. -15 d) and at the end of lactation (240 d vs. 180 d). Bioinformatics analysis revealed a major role of genes associated with BTA3, BTA4, BTA6, BTA9, BTA14, and BTA28 in lactation. Functional analysis of DEG underlined an overall induction of lipid metabolism, suggesting an increase in triglycerides synthesis, likely regulated by PPAR signaling. The same analysis revealed an induction of amino acid metabolism and secretion of protein, with a concomitant decrease in proteasome, indicating a major role of amino acid handling and reduced protein degradation in the synthesis and secretion of milk proteins. Glycan biosynthesis was induced for both N-glycan and O-glycan, suggesting increased glycan content in the milk. The cell cycle and immune response, especially antigen processing and presentation, were strongly inhibited during lactation, suggesting that morphological changes are minimized during lactation, while the mammary gland prevents immune hyper-response. Transcripts associated with response to radiation and low oxygen were enriched in the down-regulated DEG affected by the stage of lactation. Except for this last finding, the functions affected by the transcriptomic adaptation to lactation in mammary tissue of yak are very similar to those observed in dairy cows.

Polymorphism detection of <i>DGAT1</i> and <i>Lep</i> genes in Anatolian water buffalo (<i>Bubalus bubalis</i>) populations in Turkey

Acyl-CoA: diacylglycerol–acyltransferase 1 (DGAT1) enzyme plays a key role in controlling the synthesis rate triglyceride from diacylglycerol. Leptin (LP, OB, obese) is an important hormone that synthesizes mostly from adipose tissue and regulates glucose metabolism and homeostasis. DGAT1 and Lep genes are closely related to reproduction, growth, milk yield and composition in water buffalo breeds. This study aimed to identify genetic variation in the DGAT1 and Lep gene regions in 150 water buffalo individuals from five different provinces of Turkey using DNA sequencing. A total of 38 nucleotide variations and indels have identified 761 bp long partial intron 2 and exon 3 and 5′ UTR regions of the Lep gene in Anatolian water buffalo populations; 422 bp long partial exon 7–9 and exon 8 regions of DGAT1 gene were amplified and two mutations were defined in the point of 155 and 275 nucleotide that is three genotypes for S allele and Y allele of DGAT1 gene in intron 7 in Anatolian buffalo populations, respectively. These SNPs may have an effect on reproduction, growth, milk yield and composition in water buffalo populations and may prove to be useful for water buffalo breeding.

Genome-Wide Identification of Candidate Genes for Milk Production Traits in Korean Holstein Cattle

Simple Summary

Milk production traits that are economically important in the dairy industry have been considered the main selection criteria for breeding. The present genome-wide association study was performed to identify chromosomal loci and candidate genes with potential effects on milk production phenotypes in a Korean Holstein population. A total of eight significant quantitative trait locus regions were identified for milk yield (Bos taurus autosome (BTA) 7 and 14), adjusted 305-d fat yield (BTA 3, 5, and 14), adjusted 305-d protein yield (BTA 8), and somatic cell score (BTA 8 and 23) of milk production traits. Furthermore, we discovered three main candidate genes (diacylglycerol O-acyltransferase 1 (DGAT1), phosphodiesterase 4B (PDE4B), and anoctamin 2 (ANO2)) through bioinformatics analysis. These genes may help to understand better the underlying genetic and molecular mechanisms for milk production phenotypes in the Korean Holstein population.

Abstract

We performed a genome-wide association study and fine mapping using two methods (single marker regression: frequentist approach and Bayesian C (BayesC): fitting selected single nucleotide polymorphisms (SNPs) in a Bayesian framework) through three high-density SNP chip platforms to analyze milk production phenotypes in Korean Holstein cattle (n = 2780). We identified four significant SNPs for each phenotype in the single marker regression model: AX-311625843 and AX-115099068 on Bos taurus autosome (BTA) 14 for milk yield (MY) and adjusted 305-d fat yield (FY), respectively, AX-428357234 on BTA 18 for adjusted 305-d protein yield (PY), and AX-185120896 on BTA 5 for somatic cell score (SCS). Using the BayesC model, we discovered significant 1-Mb window regions that harbored over 0.5% of the additive genetic variance effects for four milk production phenotypes. The concordant significant SNPs and 1-Mb window regions were characterized into quantitative trait loci (QTL). Among the QTL regions, we focused on a well-known gene (diacylglycerol O-acyltransferase 1 (DGAT1)) and newly identified genes (phosphodiesterase 4B (PDE4B), and anoctamin 2 (ANO2)) for MY and FY, and observed that DGAT1 is involved in glycerolipid metabolism, fat digestion and absorption, metabolic pathways, and retinol metabolism, and PDE4B is involved in cAMP signaling. Our findings suggest that the candidate genes in QTL are strongly related to physiological mechanisms related to the fat production and consequent total MY in Korean Holstein cattle.

Single nucleotide polymorphisms in the growth hormone receptor gene and Alu1 polymorphisms in the diacylglycerol acyltransferase 1 gene as related to meat production in sheep

Aim:

This study aimed to investigate the polymorphisms in genes related to meat production, including growth hormone receptor (GHR) and diacylglycerol acyltransferase 1 (DGAT1) genes, in different breeds of sheep, including Barki, Najdi, and Harri.

Materials and Methods:

Blood samples were collected from 75 randomly selected healthy Barki, Najdi, and Harri breeds of sheep, with 25 samples per breed. GHR and DGAT1 genes were identified using a single nucleotide polymorphism assay followed by digestion with the restriction enzyme Alu1.

Results:

The analysis of the GHR gene sequence showed nucleotide substitutions at nt 69 in exon 10 (c.69 G > A); this mutation is considered a transition mutation. The sequences of detected SNPs in the GHR gene in the different sheep breeds were submitted to the GenBank database with accession numbers MG906773 to MG906781. The substitutions at exon 10 (c.69 G > A) results in an alteration to the amino acid (p. Lysine > Arginine). At c.69, the A allele frequency was 0.61, 0.59, and 0.54, while the G allele frequency was 0.39, 0.41, and 0.46, for Barki, Najdi, and Harri breeds, respectively. The genotype AG at nt 69 locus had the highest frequency in the Najdi and Harri sheep. The frequency of AG was 0.62, 0.61, and 0.64, while the frequency of AA was 0.30, 0.28, and 0.22, for Barki, Najdi, and Harri sheep, respectively. After digestion with the restriction enzyme AluI, the DGAT1 locus had two genotypes, CC and CT. The highest frequency, 0.88, was found for allele C, which was detected in Barki breed. The lowest frequency, 0.75, for the same allele was found for Harri.

Conclusions:

The detected CT genotype may explain the moderate intramuscular fat content and muscle marbling in the Barki sheep breed.

Exploring novel single nucleotide polymorphisms and haplotypes of the diacylglycerol O-acyltransferase 1 (DGAT1) gene and their effects on protein structure in Iranian buffalo

BACKGROUND

Diacylglycerol O-acyltransferase 1 (DGAT1) plays a key role in the synthesis of triglycerides. Recent studies have shown that a transition mutation resulting in substitutions of guanine by adenine in the DGAT1 gene in cattle has considerable effects on milk yield and composition. Currently, there is no systematic research reporting on the utilization of this gene segment in Iranian buffalo (Bubalus bubalis).

OBJECTIVE

In this study, the genetic differentiation of three indigenous Iranian buffalo populations was investigated in the region spanning exon 3 to exon 17 of the DGAT1 gene.

METHODS

A total of 200 buffaloes were genotyped, all the samples were sequenced directly in both directions with forward and reverse sequencing primers.

RESULTS

Sequence analysis showed novel SNPs compared to the reference GenBank sequence (DQ886485) at nucleotide positions g.6097A>G, g.7036C>T, g.7338G>A, g.7710C>T, g.8087C>T, g.8259G>A, g.8275G>A, g.8367C>T, and g.8426C>T. No polymorphisms were found within exon 8. Therefore, the K232A position was thought to be a conserved and fixed region for high milk fat content (K allele) in Bos indicus and all buffalo breeds. Comparison with Indian buffalo revealed three exonic SNPs, one of which was nonsynonymous. A unique 22 bp insertion was observed in intron 10 of DGAT1. Linkage disequilibrium analysis allowed the identification of nine haplotypes among the sampled animals. To our knowledge, this is the first report of sequencing analysis of the DGAT1 gene in Iranian buffalo.

CONCLUSION

Our results suggest that genetic diversity exists and could be useful in examining the association between the DGAT1 gene and milk production traits in buffalo.

Genome-wide association study of conformation and milk yield in mixed-breed dairy goats

Identification of genetic markers that affect economically important traits is of high value from a biological point of view, enabling the targeting of candidate genes and providing practical benefits for the industry such as wide-scale genomic selection. This study is one of the first to investigate the genetic background of economically important traits in dairy goats using the caprine 50K single nucleotide polymorphism (SNP) chip. The aim of the project was to perform a genome-wide association study for milk yield and conformation of udder, teat, and feet and legs. A total of 137,235 milk yield records on 4,563 goats each scored for 10 conformation traits were available. Out of these, 2,381 goats were genotyped with the Illumina Caprine 50K BeadChip (Illumina Inc., San Diego, CA). A range of pseudo-phenotypes were used including deregressed breeding values and pseudo-estimated breeding values. Genome-wide association studies were performed using the multi-locus mixed model (MLMM) algorithm implemented in SNP & Variation Suite v7.7.8 (Golden Helix Inc., Bozeman, MT). A genome-wise significant [-log₁₀(P-value) > 5.95] SNP for milk yield was identified on chromosome 19, with additional chromosome-wise significant (-log₁₀(P-value) > 4.46] SNP on chromosomes 4, 8, 14, and 29. Three genome-wise significant SNP for conformation of udder attachment, udder depth, and front legs were identified on chromosome 19, and chromosome-wise SNP were found on chromosomes 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 21, 23, and 27. The proportion of variance explained by the significant SNP was between 0.4 and 7.0% for milk yield and between 0.1 and 13.8% for conformation traits. This study is the first attempt to identify SNP associated with milk yield and conformation in dairy goats. Two genome-wise significant SNP for milk yield and 3 SNP for conformation of udder attachment, udder depth, and front legs were found. Our results suggest that conformation traits have a polygenic background because, for most of them, we did not identify any quantitative trait loci with major effect.

Regulation of DNA methylation on EEF1D and RPL8 expression in cattle

Dynamic changes to the epigenome play a critical role in a variety of biology processes and complex traits. Many important candidate genes have been identified through our previous genome wide association study (GWAS) on milk production traits in dairy cattle. However, the underlying mechanism of candidate genes have not yet been clearly understood. In this study, we analyzed the methylation variation of the candidate genes, EEF1D and RPL8, which were identified to be strongly associated with milk production traits in dairy cattle in our previous studies, and its effect on protein and mRNA expression. We compared DNA methylation profiles and gene expression levels of EEF1D and RPL8 in five different tissues (heart, liver, mammary gland, ovary and muscle) of three cows. Both genes showed the highest expression level in mammary gland. For RPL8, there was no difference in the DNA methylation pattern in the five tissues, suggesting no effect of DNA methylation on gene expression. For EEF1D, the DNA methylation levels of its first CpG island differed in the five tissues and were negatively correlated with the gene expression levels. To further investigate the function of DNA methylation on the expression of EEF1D, we collected blood samples of three cows at early stage of lactation and in dry period and analyzed its expression and the methylation status of the first CpG island in blood. As a result, the mRNA expression of EEF1D in the dry period was higher than that at the early stage of lactation, while the DNA methylation level in the dry period was lower than that at the early stage of lactation. Our result suggests that the DNA methylation of EEF1D plays an important role in the spatial and temporal regulation of its expression and possibly have an effect on the milk production traits.

A genome scan for milk production traits in dairy goats reveals two new mutations in Dgat1 reducing milk fat content

The quantity of milk and milk fat and proteins are particularly important traits in dairy livestock. However, little is known about the regions of the genome that influence these traits in goats. We conducted a genome wide association study in French goats and identified 109 regions associated with dairy traits. For a major region on chromosome 14 closely associated with fat content, the Diacylglycerol O-Acyltransferase 1 (DGAT1) gene turned out to be a functional and positional candidate gene. The caprine reference sequence of this gene was completed and 29 polymorphisms were found in the gene sequence, including two novel exonic mutations: R251L and R396W, leading to substitutions in the protein sequence. The R251L mutation was found in the Saanen breed at a frequency of 3.5% and the R396W mutation both in the Saanen and Alpine breeds at a frequencies of 13% and 7% respectively. The R396W mutation explained 46% of the genetic variance of the trait, and the R251L mutation 6%. Both mutations were associated with a notable decrease in milk fat content. Their causality was then demonstrated by a functional test. These results provide new knowledge on the genetic basis of milk synthesis and will help improve the management of the French dairy goat breeding program.

A post-GWAS replication study confirming the PTK2 gene associated with milk production traits in Chinese Holstein

Our initial genome-wide association study (GWAS) demonstrated that two SNPs (ARS-BFGL-NGS-33248, UA-IFASA-9288) within the protein tyrosine kinase 2 (PTK2) gene were significantly associated with milk production traits in Chinese Holstein dairy cattle. To further validate if the statistical evidence provided in GWAS were true-positive findings, a replication study was performed herein through genotype-phenotype associations. The two tested SNPs were found to show significant associations with milk production traits, which confirmed the associations observed in the original study. Specifically, SNPs lying in the PTK2 gene were also detected by sequencing 14 unrelated sires in Chinese Holsteins and a total of thirty-three novel SNPs were identified. Thirteen out of these identified SNPs were genotyped and tested for association with milk production traits in an independent resource population. After Bonferroni correction for multiple testing, twelve SNPs were statistically significant for more than two milk production traits. Analyses of pairwise D' measures of linkage disequilibrium (LD) between all SNPs were also explored. Two haplotype blocks were inferred and the association study at haplotype level revealed similar effects on milk production traits. In addition, the RNA expression analyses revealed that a non-synonymous coding SNP (g.4061098T>G) was involved in the regulation of gene expression. Thus the findings presented here provide strong evidence for associations of PTK2 variants with dairy production traits and may be applied in Chinese Holstein breeding program.

A novel dynamic impact approach (DIA) for functional analysis of time-course omics studies: validation using the bovine mammary transcriptome

The overrepresented approach (ORA) is the most widely-accepted method for functional analysis of microarray datasets. The ORA is computationally-efficient and robust; however, it suffers from the inability of comparing results from multiple gene lists particularly with time-course experiments or those involving multiple treatments. To overcome such limitation a novel method termed Dynamic Impact Approach (DIA) is proposed. The DIA provides an estimate of the biological impact of the experimental conditions and the direction of the impact. The impact is obtained by combining the proportion of differentially expressed genes (DEG) with the log2 mean fold change and mean –log P-value of genes associated with the biological term. The direction of the impact is calculated as the difference of the impact of up-regulated DEG and down-regulated DEG associated with the biological term. The DIA was validated using microarray data from a time-course experiment of bovine mammary gland across the lactation cycle. Several annotation databases were analyzed with DIA and compared to the same analysis performed by the ORA. The DIA highlighted that during lactation both BTA6 and BTA14 were the most impacted chromosomes; among Uniprot tissues those related with lactating mammary gland were the most positively-impacted; within KEGG pathways ‘Galactose metabolism’ and several metabolism categories related to lipid synthesis were among the most impacted and induced; within Gene Ontology “lactose biosynthesis” among Biological processes and “Lactose synthase activity” and “Stearoyl-CoA 9-desaturase activity” among Molecular processes were the most impacted and induced. With the exception of the terms ‘Milk’, ‘Milk protein’ and ‘Mammary gland’ among Uniprot tissues and SP_PIR_Keyword, the use of ORA failed to capture as significantly-enriched (i.e., biologically relevant) any term known to be associated with lactating mammary gland. Results indicate the DIA is a biologically-sound approach for analysis of time-course experiments. This tool represents an alternative to ORA for functional analysis.

RNA interference-mediated knockdown of DGAT1 decreases triglyceride content of bovine mammary epithelial cell line

Diacylglyceroltransferase-1 (DGAT1) expresses in nearly all tissues, including the mammary gland. Mice lacking DGAT1 exhibit decreased triglyceride content in mammary tissue, and are resistant to diet-induced obesity and diabetes mellitus. Thus, DGAT1 has received considerable attention. In the present study, the function of DGAT1 was examined by liposome mediated RNA interference (RNAi) to knockdown the expression of endogenous DGAT1 expression in bovine mammary epithelial cells (BMEC) and the changes of the biological functions of cells were analyzed. The mRNA and protein levels, intracellular triglyceride (TG) content, and total protein of BMECs were analyzed by real-time PCR, Western blot, TG kit, and ultraviolet spectrophotometer, respectively, before and after RNAi treatment. The results indicated that knockdown of DGAT1 expression significantly reduced TG content in BMECs. This study further confirmed the importance of DGAT1 in triglyceride synthesis in bovine mammary tissue.

Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature

We report mapping of a quantitative trait locus (QTL) with a major effect on bovine stature to a ∼780-kb interval using a Hidden Markov Model-based approach that simultaneously exploits linkage and linkage disequilibrium. We re-sequenced the interval in six sires with known QTL genotype and identified 13 clustered candidate quantitative trait nucleotides (QTNs) out of >9,572 discovered variants. We eliminated five candidate QTNs by studying the phenotypic effect of a recombinant haplotype identified in a breed diversity panel. We show that the QTL influences fetal expression of seven of the nine genes mapping to the ∼780-kb interval. We further show that two of the eight candidate QTNs, mapping to the PLAG1-CHCHD7 intergenic region, influence bidirectional promoter strength and affect binding of nuclear factors. By performing expression QTL analyses, we identified a splice site variant in CHCHD7 and exploited this naturally occurring null allele to exclude CHCHD7 as single causative gene.

Genome wide association studies for milk production traits in Chinese Holstein population

Genome-wide association studies (GWAS) based on high throughput SNP genotyping technologies open a broad avenue for exploring genes associated with milk production traits in dairy cattle. Motivated by pinpointing novel quantitative trait nucleotide (QTN) across Bos Taurus genome, the present study is to perform GWAS to identify genes affecting milk production traits using current state-of-the-art SNP genotyping technology, i.e., the Illumina BovineSNP50 BeadChip. In the analyses, the five most commonly evaluated milk production traits are involved, including milk yield (MY), milk fat yield (FY), milk protein yield (PY), milk fat percentage (FP) and milk protein percentage (PP). Estimated breeding values (EBVs) of 2,093 daughters from 14 paternal half-sib families are considered as phenotypes within the framework of a daughter design. Association tests between each trait and the 54K SNPs are achieved via two different analysis approaches, a paternal transmission disequilibrium test (TDT)-based approach (L1-TDT) and a mixed model based regression analysis (MMRA). In total, 105 SNPs were detected to be significantly associated genome-wise with one or multiple milk production traits. Of the 105 SNPs, 38 were commonly detected by both methods, while four and 63 were solely detected by L1-TDT and MMRA, respectively. The majority (86 out of 105) of the significant SNPs is located within the reported QTL regions and some are within or close to the reported candidate genes. In particular, two SNPs, ARS-BFGL-NGS-4939 and BFGL-NGS-118998, are located close to the DGAT1 gene (160bp apart) and within the GHR gene, respectively. Our findings herein not only provide confirmatory evidences for previously findings, but also explore a suite of novel SNPs associated with milk production traits, and thus form a solid basis for eventually unraveling the causal mutations for milk production traits in dairy cattle.

Assessment of DGAT1 and LEP gene polymorphisms in three Nelore (Bos indicus) lines selected for growth and their relationship with growth and carcass traits

The aim of this study was to analyze LEP and DGAT1 gene polymorphisms in 3 Nelore lines selected for growth and to evaluate their effects on growth and carcass traits. Traits analyzed were birth, weaning, and yearling weight, rump height, LM area, backfat thickness, and rump fat thickness obtained by ultrasound. Two SNP in the LEP gene [LEP 1620(A/G) and LEP 305(T/C)] and the K232A mutation in the DGAT1 gene were analyzed. The sample consisted of 357 Nelore heifers from 2 lines selected for yearling weight and a control line, established in 1980, at the Estação Experimental de Zootecnia de Sertãozinho (Sertãozinho, Brazil). Three genotypes were obtained for each marker. Differences in allele frequencies among the 3 lines were only observed for the DGAT1 K232A polymorphism, with the frequency of the A allele being greater in the control line than in the selected lines. The DGAT1 K232A mutation was associated only with rump height, whereas LEP 1620(A/G) was associated with weaning weight and LEP 305(T/C) with birth weight and backfat thickness. However, more studies, with larger data sets, are necessary before these makers can be used for marker-assisted selection.

Genome assembly anchored QTL map of bovine chromosome 14

Bovine chromosome 14 (BTA14) has been widely explored for quantitative trait loci (QTL) and genes related to economically important traits in both dairy and beef cattle. We reviewed more than 40 investigations and anchored 126 QTL to the current genome assembly (Btau 4_0). Using this anchored QTL map, we observed that, in dairy cattle, the region spanning 0 – 10 Mb on BTA14 has the highest density QTL map with a total of 56 QTL, mainly for milk production traits. It is very likely that both somatic cell score (SCS) and clinical mastitis share some common QTL in two regions: 61.48 Mb - 73.84 Mb and 7.86 Mb – 39.55 Mb, respectively. As well, both ovulation rate and twinning rate might share a common QTL region from 34.16 Mb to 65.38 Mb. However, there are no common QTL locations in three pregnancy related phenotypes: non-return rate, pregnancy rate and daughter pregnancy rate. In beef cattle, the majority of QTL are located in a broad region of 15 Mb – 45 Mb on the chromosome. Functional genes, such as CRH, CYP11B1, DGAT1, FABP4 and TG, as potential candidates for some of these QTL, were also reviewed. Therefore, our review provides a standardized QTL map anchored within the current genome assembly, which would enhance the process of selecting positional and physiological candidate genes for many important traits in cattle.

QTL detection for milk production traits in goats using a longitudinal model

Eight paternal half-sib families were used to identify chromosomal regions associated with variation in the lactation curves of dairy goats. DNA samples from 162 animals were amplified by PCR for 37 microsatellite markers, from Capra hircus autosomes CHI3, CHI6, CHI14 and CHI20. Milk samples were collected during 6 years, and there were 897 records for milk yield (MY) and 814 for fat (FP) and protein percentage (PP). The analysis was conducted in two stages. First, a random regression model with several fixed effects was fitted to describe the lactation function, using a scale (alpha) plus four shape parameters: beta and gamma, both associated with a decrease in the slope of the curve, and delta and phi that are related to the increase in slope. Predictions of alpha, beta, gamma, delta and phi were regressed using an interval mapping model, and F-tests were used to test for quantitative trait loci (QTL) effects. Significant (p < 0.05) QTLs were found for: (i) MY: CHI6 at 70-80 cM for all parameters; CHI14 at 14 cM for delta and phi; (ii) FP: CHI14, at 63 cM was associated with beta; CHI20, at 72 cM, showed association with alpha; (iii) PP: chromosomal regions associated with beta were found at 59 cM in CHI3 and at 55 cM in CHI20 with alpha and gamma. Analyses using more families and more animals will be useful to confirm or to reject these findings.

Effect of DGAT1 and TG gene polymorphisms on intramuscular fat and on milk production traits in different cattle breeds in Hungary

The objective of this study was to estimate the effect of the thyroglobulin (TG) locus on beef quality traits in some beef cattle breeds and to investigate the effect of the DGAT1 locus on milk production traits in the Hungarian Holstein Friesian population. TG and DGAT1 genotypes were determined by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay. At the TG locus TT bulls showed the highest fat percentage values in the longissimus dorsi muscle (m. longissimus dorsi); the difference between CC and TT genotypes was significant. DGAT1 GC/GC cows had the highest milk, fat and protein yield values. Due to the relatively small number of GC/GC cows the difference proved to be significant only between AA/AA and AA/GC genotypes.

Alleles of the bovine DGAT1 variable number of tandem repeat associated with a milk fat QTL at chromosome 14 can stimulate gene expression

A quantitative trait locus (QTL) affecting milk fat percentage has been mapped to the centromeric end of the bovine chromosome 14 (BTA14). This genomic area includes the DGAT1 gene, which encodes acyl-CoA:diacylglycerol acyltransferase 1, the key enzyme of triglyceride biosynthesis. Genetic and biochemical studies led to the identification of the nonconservative DGAT1-K232A polymorphism as a causal mutation for the QTL. In addition to this, another polymorphism in the 5'-regulatory region of this gene, the DGAT1 variable number of tandem repeat (VNTR), also showed a strong association with milk fat percentage. This promoter VNTR polymorphism affects the number of potential Sp1 binding sites and therefore might have an impact on DGAT1 expression and also milk fat content. Hence, the DGAT1 VNTR polymorphism might be another causal mutation for the BTA14 QTL. However, evidence for Sp1 binding to this polymorphic site and for the capability of DGAT1 VNTR alleles to stimulate gene expression was lacking. In the current work Sp1-VNTR interactions were analyzed by EMSA. In addition, effects of DGAT1 VNTR alleles on gene expression were measured with reporter gene analyses. Conclusions from the results are that 1) the DGAT1 VNTR sequence is indeed a target for Sp1 binding; 2) DGAT1 VNTR alleles can stimulate gene expression in vitro and probably in vivo as well; and 3) although the stimulating effects of the different DGAT1 VNTR alleles did not show significant differences in vitro, their effects on transcription might be different in the chromatin context existing in vivo.

From QTL to QTN identification in livestock--winning by points rather than knock-out: a review

Many quantitative trait loci (QTL) affecting economic traits in livestock have now been identified. However, the confidence interval (CI) of individual QTL as determined by linkage analysis often spans tens of map units, containing hundreds of genes. Linkage disequilibrium (LD) mapping can reduce the CI to individual map units, but this reduced interval will still contain tens of genes. Methods suitable for model animals to find and validate specific quantitative trait nucleotides (QTN) underlying the QTL cannot be easily applied to livestock species because of their long generation intervals, the cost of maintaining each animal and the difficulty of producing transgenics or 'knock-outs'. Considering these limitations, we review successful approaches for identifying QTN in livestock and outline a schematic strategy for QTN determination and verification. In addition to linkage and LD mapping, the methods include positional cloning, selection of candidate genes, DNA sequencing and statistical analyses. Concordance determination and functional assays are the critical tests for validation of a QTN; we provide a generalized formula for the probability of concordance by chance. Three genes that meet the burden of proof for QTN identification--DGAT1 in cattle, IGF2 in swine and GDF8 in sheep--are discussed in detail. The genetic and economic ramifications of identified QTN and the horizon for selection and introgression are also considered.

Statistical modeling of candidate gene effects on milk production traits in dairy cattle

A major objective of dairy cattle genomic research is to identify genes underlying the variability of milk production traits that could be useful in breeding programs. The candidate gene approach provides tools for searching for causative polymorphisms affecting quantitative traits. Genes with a possible effect on milk traits in cattle can be involved in different physiological pathways, such as triglyceride synthesis [acyl-CoA:diacylglycerol acyltransferase 1 gene (DGAT1)], fat secretion from the mammary epithelial tissue (butyrophilin), or entire-body energy homeostasis regulation (leptin and leptin receptor). In this study, based on data from 252 Black and White bulls from the active Polish dairy population, effects and potential interactions of 9 single nucleotide polymorphisms in the butyrophilin, DGAT1, leptin, and leptin receptor genes were investigated. Additionally, the effect of the number of additive, dominance, and epistatic genetic effects fitted into the model on the estimates of model parameters and model selection was illustrated. Phenotypic records were daughter yield deviations for milk, fat, and protein yields, obtained from a routine national genetic evaluation. Out of all the analyzed polymorphisms, DGAT1 K232A had a much larger effect on milk traits than the other single nucleotide polymorphisms considered. Estimates of the additive genetic effect of K232A expressed as half of the difference between Lys- and Ala-encoding variants were -107.4 kg of milk, 5.4 kg of fat, and -1.6 kg of protein at first parity, as well as -120 kg of milk and 6.8 kg of fat at second parity. In terms of model selection, it was demonstrated that the modified version of Bayesian information criterion selects models with the parameterization reflecting the genetic background of the analyzed trait, while the Bayesian information criterion chooses models that are too highly parameterized.

Molecular Cloning and Single Nucleotide Polymorphism Detection of Buffalo DGAT1 Gene

We report molecular cloning and single nucleotide polymorphism detection of the buffalo DGAT1 gene. Diacylglycerol acyltransferase (DGAT1) is considered the key enzyme in controlling the rate of synthesis of triglycerides. The DGAT1 gene was recently identified as a strong functional candidate gene affecting milk yield and composition in cattle. A full-length buffalo DGAT1 genomic DNA was amplified by iterative PCR based on homolog cloning. The buffalo DGAT1 gene comprises 17 exons and spans approximately 8.3 kb. The genomic structures of DGAT1 are highly conserved among mammal species. The deduced protein of buffalo DGAT1 contains 489 amino acids, showing high-sequence similarity with mammal homologs. Through PCR-SSCP analysis and sequencing, seven polymorphic positions were detected in the complete genomic region of buffalo DGAT1, and their frequencies were observed from a collection of 117 buffalo. The SNP (C/T) detected at position 11785 in exon 17 creates a substitution change for the amino acid sequence, resulting in an Ala residue (GCG) transition to a Val residue (GTG) in position 484 of buffalo DGAT1 protein. Information provided in this study will be useful in further studies to determine the role DGAT1 plays in the regulation of milk fat synthesis and quality improvement for milk in buffalo.

Characterization of the DGAT1 K232A and Variable Number of Tandem Repeat Polymorphisms in French Dairy Cattle

A quantitative trait locus (QTL) underlying different milk production traits has been identified with a high significance threshold value in the genomic region containing the acylCoA:diacylglycerol acyltransferase (DGAT1) gene, in the 3 main French dairy cattle breeds: French Holstein, Normande, and Montbéliarde. Previous studies have confirmed that the K232A polymorphism in DGAT1 is responsible for a major QTL underlying several milk production traits in Holstein dairy cattle and several other bovine breeds. In this study, we estimate the frequency of the 2 alternative alleles, K and A, of the K232A polymorphism in French Holstein, Normande, and Montbéliarde breeds. Although the K allele segregates in French Holstein and Normande breeds with a similar effect on production traits, the existence of additional mutations contributing to the observed QTL effect is strongly suggested in both breeds by the existence of sires heterozygous at the QTL but homozygous at the K232A polymorphism. One allele at a variable number of tandem repeats (VNTR) locus in the 5' noncoding region of DGAT1 has been recently proposed as a putative causative variant. In our study, this marker was found to present a high mutation rate of 0.8% per gamete and per generation, making the allele diversity observed compatible with that expected under neutrality. Moreover, among the sires homozygous at the K232A polymorphism, no allele at the VNTR can fully explain their QTL status. Finally, no allele at the VNTR was found to be significantly associated with the fat percentage variation in the 3 breeds simultaneously after correction for the effect of the K232A polymorphism. Therefore, our results suggest the existence of at least one other causative polymorphism not yet described. Because the A allele is nearly fixed in the Montbéliarde breed, this breed represents an interesting model to identify and confirm other mutations that have a strong effect on milk production traits.

Identification of bovine QTL for growth and carcass traits in Japanese Black cattle by replication and identical-by-descent mapping

To map quantitative trait loci (QTL) for growth and carcass traits in a purebred Japanese Black cattle population, we conducted multiple QTL analyses using 15 paternal half-sib families comprising 7860 offspring. We identified 40 QTL with significant linkages at false discovery rates of less than 0.1, which included 12 for intramuscular fat deposition called marbling and 12 for cold carcass weight or body weight. The QTL each explained 2%-13% of the phenotypic variance. These QTL included many replications and shared hypothetical identical-by-descent (IBD) alleles. The QTL for CW on BTA14 was replicated in five families with significant linkages and in two families with a 1% chromosome-wise significance level. The seven sires shared a 1.1-Mb superior Q haplotype as a hypothetical IBD allele that corresponds to the critical region previously refined by linkage disequilibrium mapping. The QTL for marbling on BTA4 was replicated in two families with significant linkages. The QTL for marbling on BTA6, 7, 9, 10, 20, and 21 and the QTL for body weight on BTA6 were replicated with 1% and/or 5% chromosome-wise significance levels. There were shared IBD Q or q haplotypes in the marbling QTL on BTA4, 6, and 10. The allele substitution effect of these haplotypes ranged from 0.7 to 1.2, and an additive effect between the marbling QTL on BTA6 and 10 was observed in the family examined. The abundant and replicated QTL information will enhance the opportunities for positional cloning of causative genes for the quantitative traits and efficient breeding using marker-assisted selection.

Analysis of a polymorphism in the DGAT1 gene in 14 cattle breeds through PCR-SSCP methods

The diacylglycerol O-acyltransferase (DGAT1) is a microsomal enzyme that catalyzes the final step of triglyceride synthesis. Recent work have evidenced a significant association between lysine at amino acid position 232 with elevated milk fat content, while an alanine at this position is associated with lowered milk fat content. The aim of the present work was to develop a simple and inexpensive PCR-SSCP assay in order to discriminate the CG/AA alleles in exon 8 of the DGAT1 gene. In addition, this method was used to analyze the polymorphism of the DGAT1 through PCR-SSCP methods in 14 populations of cattle from Argentine, Bolivia and Uruguay. The PCR primers were designed from GenBank reported sequences. In this study, we found three PCR-SSCP variants, which were denominated from "A" to "C". However, DNA sequencing analysis showed that "A" variant corresponded with the A allele, while both "B" and "C" observed pattern have the motif AA at positions 10,433-10,434 (K allele), being two alternative conformations of the same DNA sequence. Both variants were detected within each breed with the exception of Hereford, and the heterozygosity varied between 0.000 and 0.524. The gene frequency analysis evidenced significant differences among the studied breeds (F(ST) = 0.325, p = 0.000). European Bos taurus breeds, with the exception of Jersey breed, showed the lowest frequency of the K allele, while highest K allele frequencies were harboured by Bos indicus type cattle. In addition, unselected South American Creole cattle breeds and the synthetic Brangus breed had intermediate allele frequencies.

Multitrait Quantitative Trait Loci Mapping for Milk Production Traits in Danish Holstein Cattle

The aims of this study were (1) to confirm previously identified quantitative trait loci (QTL) on bovine chromosomes 6, 11, 14, and 23 in the Danish Holstein cattle population, (2) to assess the pleiotropic nature of each QTL on milk production traits by building multitrait and multi-QTL models, and (3) to include pedigree information on nongenotyped individuals to improve the estimation of genetic parameters underlying the random QTL model. Nineteen grandsire families were analyzed by single-trait (ST) and multitrait (MT) QTL mapping methods. The variance component-based QTL mapping model was implemented via restricted maximum likelihood (REML) to estimate QTL position and parameters. Segregation of the previously identified QTL was confirmed on bovine chromosomes 6, 11, and 14, but not on 23. A highly significant (1% chromosome-wise level) QTL was found on chromosome 6, between 37 and 73 cM. This QTL had a strong effect on protein percentage (PP) and fat percentage (FP) according to ST analyses, and effects on PP, FP, milk yield (MY), fat yield (FY), and protein yield (PY) in MT analyses. A QTL affecting PP was detected on chromosome 11 (at 70 cM) using ST analysis. The MT analysis revealed a second QTL (at 67 cM) approaching significance with an effect on MY. The ST analysis identified a QTL for MY and FP on chromosome 14, between 10 and 24 cM. The extended pedigree (nongenotyped animals) was included to estimate genetic parameters underlying the random QTL model; that is, additive polygenic and QTL variances. In general, the estimates of the QTL variance components were smaller but more precise when the extended pedigree was considered in the analysis.

Linkage disequilibrium on the bovine X chromosome: characterization and use in quantitative trait locus mapping

We herein demonstrate that in the Holstein-Friesian dairy cattle population, microsatellites are as polymorphic on the X chromosome as on the autosomes but that the level of linkage disequilibrium between these markers is higher on the X chromosome than on the autosomes. The latter observation is not compatible with the small male-to-female ratio that prevails in this population and results in a higher gonosomal than autosomal effective population size. It suggests that the X chromosome undergoes distinct selective or mutational forces. We describe and characterize a novel Markovian approach to exploit this linkage disequilibrium to compute the probability that two chromosomes are identical-by-descent conditional on flanking marker data. We use the ensuing probabilities in a restricted maximum-likelihood approach to search for quantitative trait loci (QTL) affecting 48 traits of importance to the dairy industry and provide evidence for the presence of QTL affecting 5 of these traits on the bovine X chromosome.

Whole-genome scan to detect QTL for milk production, conformation, fertility and functional traits in two US Holstein families

A genome scan was conducted in two US Holstein half-sib families to identify quantitative trait loci (QTL) affecting milk production and conformation traits using the granddaughter design. The sires of the two studied families were related as sire and son and had 96 and 212 sons respectively. A total of 221 microsatellite loci were scored in both families. Statistical analysis was performed using two different analytical methods; half-sib least squares regression and Bayesian Monte Carlo Markov Chain. Traits analysed included five traditional milk production traits, somatic cell count, daughter pregnancy rate, male fertility and 20 conformation traits. A total of 47 tests achieved at least genome-wise significance. However, results from the two methods of analysis were only concordant for QTL location and level of significance in eight instances.

Linkage disequilibrium in the domesticated pig

This study investigated the extent of linkage disequilibrium (LD) in two genomic regions (on chromosomes 4 and 7) in five populations of domesticated pigs. LD was measured with D' and tested for significance with the Fisher exact test. Effects of genetic (linkage) distance, chromosome, population, and their interactions on D' were tested both through a linear model analysis of covariance and by a theoretical nonlinear model. The overall result was that (1) the distance explained most of the variability of D', (2) the effect of chromosome was significant, and (3) the effect of population was significant. The significance of the chromosome effect may have resulted from selection and the significance of the population effect illustrates the effects of population structures and effective population sizes on LD. These results suggest that mapping methods based on LD may be valuable even with only moderately dense marker spacing in pigs.

QTL fine-mapping with recombinant-inbred heterogeneous stocks and in vitro heterogeneous stocks

We compared strategies to fine-map Quantitative Trait Loci (QTL) in mice with heterogeneous stocks (HS). We showed that a panel of about 100 Recombinant Inbred Lines (RIL) derived from an HS, and which we called an RIHS, was ideally suited to fine-map QTL to very high resolution, without the cost of additional genotyping. We also investigated a strategy based on in vitro fertilization of large numbers of F(1) offspring of HS males crossed with an inbred line (IVHS). This method required some additional genotyping but avoided the breeding delays and costs associated with the construction of an RI panel. We showed that QTL detection was higher by using RIHS than with IVHS and that it was independent of the number of RI lines, provided the total number of animals phenotyped was constant. However, fine-mapping accuracy was slightly better with IVHS. We also investigated the effects of varying the number of HS generations and using multiallelic microsatellites instead of SNPs. We found that quite modest generation times of 10-20 generations were optimal. Microsatellites were superior to SNPs only when the generation time was 30 or more and when the markers were widely spaced.

Genetic and functional confirmation of the causality of the DGAT1 K232A quantitative trait nucleotide in affecting milk yield and composition

We recently used a positional cloning approach to identify a nonconservative lysine to alanine substitution (K232A) in the bovine DGAT1 gene that was proposed to be the causative quantitative trait nucleotide underlying a quantitative trait locus (QTL) affecting milk fat composition, previously mapped to the centromeric end of bovine chromosome 14. We herein generate genetic and functional data that confirm the causality of the DGAT1 K232A mutation. We have constructed a high-density single-nucleotide polymorphism map of the 3.8-centimorgan BULGE30-BULGE9 interval containing the QTL and show that the association with milk fat percentage maximizes at the DGAT1 gene. We provide evidence that the K allele has undergone a selective sweep. By using a baculovirus expression system, we have expressed both DGAT1 alleles in Sf9 cells and show that the K allele, causing an increase in milk fat percentage in the live animal, is characterized by a higher Vmax in producing triglycerides than the A allele.

The DGAT1 K232A mutation is not solely responsible for the milk production quantitative trait locus on the bovine chromosome 14

The gene, acyl-CoA:diacylglycerol acyltransferase1 (DGAT1), was recently identified as the one underlying the quantitative trait locus (QTL) for milk production traits in the centromeric region of the bovine chromosome 14. Until now, 2 alleles, the lysine variant (increasing fat yield, fat and protein percentage) and the alanine variant (increasing protein and milk yield), were postulated at DGAT1. This study investigated whether the diallelic DGAT1 polymorphism is responsible for all the genetic variation at the centromeric region of this chromosome for milk, fat, and protein yield and fat and protein percentage. A statistical model was applied to a granddaughter design to analyze 16 German Holstein families. The model included the diallelic DGAT1 effect and the QTL transition probability estimated for each chromosomal position by a multiple marker approach. Because the regression coefficient of this probability was corrected for the diallelic DGAT1 polymorphism, it represented a putative conditional QTL effect. The effect of the DGAT1 gene was always highly significant. The conditional QTL effect was significant genomewise for fat percentage at the proximal end of the chromosome and for protein percentage at a more distal chromosomal region. Additional chromosomewise significance was found for fat and protein yield. Our results suggest an additional source of genetic variance on this chromosome for these traits; either one or more additional alleles segregating at DGAT1 that were not previously detected, a second quantitative trait locus affecting these traits, or both.

Assessment of the gene content of the chromosomal regions flanking bovine DGAT1

As a first step towards verifying the candidate status of DGAT1 as the causal gene for milk fat percentage in cattle, we constructed a bovine BAC contig spanning 576 kb of the chromosomal region containing DGAT1. High content of NotI sites (21 within the contig) indicated that the region is gene-rich. Twenty-three genes neighboring DGAT1 were mapped, including two bovine cDNA sequences that have no orthologous sequences within the NCBI sequence databases. On average, 2015 bp for each of the 23 neighboring genes were sequenced and entered into EMBL. Likewise, 10 new STS markers were established by BAC-end sequencing. Within the genes and STS markers, 55 polymorphisms were discovered. These will form the basis of future linkage disequilibrium studies to test whether any genes neighboring DGAT1 are associated with variation in milk fat percentage, thereby testing the candidate status of DGAT1.

Fine mapping of quantitative trait loci and assessment of positional candidate genes for backfat on bovine chromosome 14 in a commercial line of Bos taurus

Backfat thickness is one of the major quantitative traits that affect carcass quality in beef cattle. In this study, we have fine mapped a QTL for backfat EBV on bovine chromosome 14, using an identical-by-descent haplotype-sharing analysis, in a commercial line of Bos taurus. We also examined the association between gene-specific single nucleotide polymorphism (SNP) markers of the genes diacylgcerol acyltransferase 1 (DGAT1) and thyroglobulin (TG) and the backfat EBV. The results indicate that the QTL region for backfat identified on chromosome 14 is in agreement with previous studies. However, neither of the two polymorphisms of candidate genes tested, DGAT1 nor TG, showed a significant (P > 0.10) association with the backfat EBV in the cattle populations examined. However, a strong association (P = 0.0058) was detected between a microsatellite marker (CSSM66) lying approximately mid-way between the two candidate genes and the backfat EBV. These results suggest that other SNP of DGAT1, TG, or other gene(s) in the chromosomal region should be examined to test whether they have a significant effect on lipid metabolism.

Population-wide analysis of a QTL affecting milk-fat production in the Israeli Holstein population

A population-wide linkage disequilibrium on bovine chromosome 14 between microsatellite ILSTS039 and DGAT1, a putative quantitative trait locus affecting milk production traits, was found in the Israeli Holstein population. A total of 394 bulls were genotyped for both DGAT1 and ILSTS039, and 1747 cows were genotyped for ILSTS039. The ILSTS039 allele termed "225," and the DGAT1 K allele (substitution of a lysine residue with alanine), were associated with decreased milk production, and increased fat production and fat and protein percent. The number of 225 ILSTS039 and K DGAT1 alleles per individual were the same for 80% of the bulls genotyped. From the effects associated with cows homozygous for the 225 allele, the effect of the quantitative trait locus appears to be approximately codominant. The substitution effect was 0.16% fat. Genotype probabilities for the quantitative gene were determined for the entire Israeli Holstein milk-recorded population, including 507,725 cows and 1442 bulls, using segregation analysis. Overall frequency of the allele that increased fat percent was 8.9% in cows and 15.5% in bulls. The frequency of this allele decreased from 1981 until 1990, from 15 to 5%, and since has increased to 10%. The effects estimated on the population-wide analyses of both cows and bulls were similar to the effect associated with DGAT1 in the daughters of genotyped bulls. Modified animal model evaluations were computed for the entire population with the effect of this gene included in the model. The correlations between the modified and standard animal model evaluations for all traits were > 0.99.

Quantitative trait loci affecting milk production traits in Finnish Ayrshire dairy cattle

A whole genome scan of Finnish Ayrshire was conducted to map quantitative trait loci (QTL) affecting milk production. The analysis included 12 half-sib families containing a total of 494 bulls in a granddaughter design. The families were genotyped with 150 markers to construct a 2764 cM (Haldane) male linkage map. In this study interval mapping with multiple-marker regression approach was extended to analyse multiple chromosomes simultaneously. The method uses identified QTL on other chromosomes as cofactors to increase mapping power. The existence of multiple QTL on the same linkage group was also analyzed by fitting a two-QTL model to the analysis. Empirical values for chromosome-wise significance thresholds were determined using a permutation test. Two genome-wise significant QTL were identified when chromosomes were analyzed individually, one affecting fat percentage on chromosome (BTA) 14 and another affecting fat yield on BTA12. The cofactor analysis revealed in total 31 genome-wise significant QTL. The result of two-QTL analysis suggests the existence of two QTL for fat percentage on BTA3. In general, most of the identified QTL confirm results from previous studies of Holstein-Friesian cattle. A new QTL for all yield components was identified on BTA12 in Finnish Ayrshire.

Finding genes that underlie complex traits

Phenotypic variation among organisms is central to evolutionary adaptations underlying natural and artificial selection, and also determines individual susceptibility to common diseases. These types of complex traits pose special challenges for genetic analysis because of gene-gene and gene-environment interactions, genetic heterogeneity, low penetrance, and limited statistical power. Emerging genome resources and technologies are enabling systematic identification of genes underlying these complex traits. We propose standards for proof of gene discovery in complex traits and evaluate the nature of the genes identified to date. These proof-of-concept studies demonstrate the insights that can be expected from the accelerating pace of gene discovery in this field.