Identification of SNPs in interferon gamma, interleukin-22, and their receptors and associations with health and production-related traits in Canadian Holstein bulls

Genetic variants in a number of immunoregulatory genes have been previously associated with health and production traits in dairy cattle. Therefore, in the following study, the genes coding interferon gamma (IFNG), IFNG receptor 1 and 2 domains, interleukin-22 (IL22), and IL22 receptor alpha 1, were investigated for single nucleotide polymorphisms (SNPs) in Holstein bulls. These SNPs, along with SNPs previously identified in IL10, IL10 receptor, and transforming growth factor beta 1 (TGFB1) genes, were evaluated for statistical associations to estimated breeding values for milk somatic cell score (SCS), a trait highly correlated to mastitis incidence, and various production-related traits, including milk yield, protein yield, fat yield, and lactation persistency. While no significant associations were found between these SNPs and SCS, SNPs in IL10 receptor beta subunit showed a significant effect on protein yield and lactation persistency. While there is evidence that IL10 plays an important role during lactation, it is also likely that the effects of SNPs in IL10 receptor beta subunit on protein yield and lactation persistency are due to linkage disequilibrium with a neighboring QTL.     

Identification of single nucleotide polymorphisms in the bovine CCL2, IL8, CCR2 and IL8RA genes and their association with health and production in Canadian Holsteins

The aim of this study was to identify the presence of SNPs in the chemokine genes CCL2 and IL8 and the chemokine receptor genes IL8RA and CCR2, and assess their potential contribution to variation in estimated breeding values (EBVs) for somatic cell score (SCS) and four other traits in Canadian Holstein bulls. Pools of DNA for bulls with high (H) and low (L) EBVs for SCS were used for identification of 11 SNPs. Two unreported SNPs were found in the CCL2 gene and one SNP was found in the CCR2 gene. Previously reported SNPs (three in the IL8 gene and five in the IL8RA chemokine receptor) were also identified. Two SNPs in CCL2, three in IL8, one in IL8RA and one in CCR2 were genotyped in Canadian Holstein bulls (n = 338) using tetra primer ARMS-PCR. We investigated associations of these seven polymorphisms with three production traits (milk yield, fat yield and protein yield) and one conformation trait related to mastitis (udder depth). The allele substitution effect for the CCL2 rs41255713:T>C SNP was significant at an experimental-wise level for milk yield (247.5 +/- 79.9 kg) and protein yield (7.4 +/- 2.3 kg) EBVs (P T SNP on SCS was significant at the comparison-wise level (-0.04 +/- 0.02, P = 0.05), which might indicate a possible association in support of other published studies. Lastly, we assigned CCR2 to BTA22q24, where a previously QTL for SCS was identified.     

Associations of the bovine major histocompatibility complex DRB3 (BoLA-DRB3) with production traits in Canadian dairy cattle

Associations of two alleles of the bovine major histocompatibility complex DRB3 gene (BoLA-DRB3) with lowered somatic cell score (SCS) and occurrence of disease (BoLA-DRB3.2* 16 and *23, respectively) have previously been documented. The objective of this study was to evaluate potential relationships between BoLA-DRB3 alleles with production traits, namely 305-day milk, milk fat and milk protein yield, in a population of Canadian dairy cattle (Holstein, n = 835 and Jersey, n = 66) over the course of two lactations. No significant associations were detected between BoLA alleles and production traits in Jerseys. In Holsteins, alleles *16 and *23 also did not show associations with production traits but allele *8 was significantly associated with increased 305-day milk, fat and protein yields in the previous lactation (the lactation prior to immunization with a gram negative core antigen vaccine), and with increased protein production in the subsequent (with reference to the time of immunization) lactation. Allele *22 was associated with decreased milk and protein yield in both previous and subsequent lactations. Therefore, it can be concluded that increasing or decreasing the frequency of BoLA alleles *16 and *23 to reduce SCS or increase resistance to mastitis in this population would not have adverse effects on production in this population, and that certain BoLA alleles (*8 and *22) are associated with altered production traits in Canadian Holsteins.     

Single nucleotide polymorphisms in the growth hormone and insulin-like growth factor-1 genes are associated with milk production, body condition score and fertility traits in dairy cows

The somatotrophic axis (GH-IGF) is a key regulator of animal growth and development, affecting performance traits that include milk production, growth rate, body composition, and fertility. The aim of this study was to quantify the association of previously identified SNPs in bovine growth hormone (GH1) and insulin-like growth factor 1 (IGF-1) genes with direct performance trait measurements of lactation and fertility in Holstein-Friesian lactating dairy cows. Sixteen SNPs in both IGF-1 and GH1 were genotyped across 610 cows and association analyses were carried out with traits of economic importance including calving interval, pregnancy rate to first service and 305-day milk production, using animal linear mixed models accounting for additive genetic effects. Two IGF-1 SNPs, IGF1i1 and IGF1i2, were significantly associated with body condition score at calving, while a single IGF-1 SNP, IGF1i3, was significantly associated with milk production, including milk yield (means ± SEM; 751.3 ± 262.0 kg), fat yield (21.3 ± 10.2 kg) and protein yield (16.5 ± 8.0 kg) per lactation. Only one GH1 SNP, GH33, was significantly associated with milk protein yield in the second lactation (allele substitution effect of 9.8 ± 5.0 kg). Several GH1 SNPs were significantly associated with fertility, including GH32, GH35 and GH38 with calving to third parity (22.4 ± 11.3 days) (GH32 and GH38 only), pregnancy rate to first service (0.1%) and overall pregnancy rate (0.05%). The results of this study demonstrate the effects of variants of the somatotrophic axis on milk production and fertility traits in commercial dairy cattle.     

Associations between milk performance traits in Holstein cows and 16 candidate SNPs identified by arrayed primer extension (APEX) microarray

An oligonucleotide microarray-which allows for parallel genotyping of many SNPs in genes involved in cow milk protein biosynthesis-was used to identify which of the 16 candidate SNPs are associated with milk performance traits in Holstein cows. Four hundred cows were genotyped by the developed and validated microarray. Significant associations were found between four single SNPs, namely DGAT1 (acyloCoA:diacylglycerol acyltransferase), LTF (lactoferrin), CSN3 (kappa-casein), and GHR (growth hormone receptor) and with fat and protein yield and percentage. Many significant associations between combined genotypes (two SNPs) and milk performance traits were found. The associations between the combined genotypes DGAT1/LTF and DGAT1/LEPTIN analyzed traits are presented as examples. The microarray based on APEX (Arrayed Primer Extension) is a fast and reliable method for multiple SNP analysis of potential application in marker-assisted selection. After further development, the chip may prospectively be used for dairy cattle paternity analysis and evolutionary studies.     

Associations Between Milk Performance Traits in Holstein Cows and 16 Candidate SNPs Identified by Arrayed Primer Extension (APEX) Microarray

An oligonucleotide microarray—which allows for parallel genotyping of many SNPs in genes involved in cow milk protein biosynthesis—was used to identify which of the 16 candidate SNPs are associated with milk performance traits in Holstein cows. Four hundred cows were genotyped by the developed and validated microarray. Significant associations were found between four single SNPs, namely DGAT1 (acyloCoA:diacylglycerol acyltransferase), LTF (lactoferrin), CSN3 (kappa-casein), and GHR (growth hormone receptor) and with fat and protein yield and percentage. Many significant associations between combined genotypes (two SNPs) and milk performance traits were found. The associations between the combined genotypes DGAT1/LTF and DGAT1/LEPTIN analyzed traits are presented as examples.

The microarray based on APEX (Arrayed Primer Extension) is a fast and reliable method for multiple SNP analysis of potential application in marker-assisted selection. After further development, the chip may prospectively be used for dairy cattle paternity analysis and evolutionary studies.     

Bovine IFNGR2, IL12RB1, IL12RB2, and IL23R polymorphisms and MAP infection status

Mycobacterium avium ssp. paratuberculosis (MAP) infection causes a chronic granulomatous inflammatory condition of the bovine gut that is characterized by diarrhea, progressive weight loss, and emaciation, and ultimately leads to loss in productivity and profitability of dairy operations. The host cytokine machinery is known to play an important role in protecting against MAP infection. Therefore, the goal of the present study was to assess whether polymorphisms in candidate genes encoding important cytokines and cytokine receptors are associated with MAP infection status of dairy cattle. MAP infection status was evaluated based on serum and milk enzyme-linked immunosorbent assays (ELISAs) for MAP-specific antibodies. Twenty previously reported polymorphisms in genes encoding bovine interferon gamma (IFNG), IFNGR1, IFNGR2, IL22, IL22RA1, IL12RB1, IL12RB2, and IL23R were genotyped in a resource population of 446 dairy Holsteins with known MAP infection status, and logistic regression was used to assess the statistical association with a binomial MAP infection status phenotype. Four SNPs in IFNGR2, IL12RB1, IL12RB2, and IL23R were found to be associated with the MAP infection status of the resource population. These results underscore the importance of cytokines and their receptors in conferring protection against MAP infection and warrant further functional characterization of these associations.     

Genetic correlations between production and disease traits during first lactation in Holstein cows

The aim of this study was to estimate genetic correlations between milk yield, somatic cell score (SCS), mastitis, and claw and leg disorders (CLDs) during first lactation in Holstein cows by using a threshold–linear random regression test-day model. We used daily records of milk, fat and protein yields; somatic cell count (SCC); and mastitis and CLD incidences from 46 771 first-lactation Holstein cows in Hokkaido, Japan, that calved between 2000 and 2009. A threshold animal model for binary records (mastitis and CLDs) and linear animal model for yield traits were applied in our multiple trait analysis. For both liabilities and yield traits, additive genetic effects were used as random regression on cubic Legendre polynomials of days on milk. The highest positive genetic correlations between yields and disease incidences (0.36 for milk and mastitis, 0.56 for fat and mastitis, 0.24 for protein and mastitis, 0.32 for milk and CLD, 0.44 for fat and CLD and 0.31 for protein and CLD) were estimated at about the time of peak milk yield (36 to 65 days in milk). Selection focused on early lactation yield may therefore increase the risk of mastitis and CLDs. The positive genetic correlations of SCS with mastitis or CLD incidence imply that selection to reduce SCS in the early stages of lactation would decrease the incidence of both mastitis and CLD.     

Genome Wide Analysis of Fertility and Production Traits in Italian Holstein Cattle

A genome wide scan was performed on a total of 2093 Italian Holstein proven bulls genotyped with 50K single nucleotide polymorphisms (SNPs), with the objective of identifying loci associated with fertility related traits and to test their effects on milk production traits. The analysis was carried out using estimated breeding values for the aggregate fertility index and for each trait contributing to the index: angularity, calving interval, non-return rate at 56 days, days to first service, and 305 day first parity lactation. In addition, two production traits not included in the aggregate fertility index were analysed: fat yield and protein yield. Analyses were carried out using all SNPs treated separately, further the most significant marker on BTA14 associated to milk quality located in the DGAT1 region was treated as fixed effect. Genome wide association analysis identified 61 significant SNPs and 75 significant marker-trait associations. Eight additional SNP associations were detected when SNP located near DGAT1 was included as a fixed effect. As there were no obvious common SNPs between the traits analyzed independently in this study, a network analysis was carried out to identify unforeseen relationships that may link production and fertility traits.     

Pharmacogenomics of multiple sclerosis: association of immune response genes polymorphism with copaxone treatment efficacy

Complex association analysis of copaxone (glatiramer acetate) immunotherapy efficacy with allelic polymorphism in the number of immune response genes, which encode interferone beta (IFNB1), transforming growth factor beta1 (TGFB1), interferone gamma (IFNG), tumor necrosis factor (TNF), interferon alpha/beta receptor 1 (IFNAR1), CC chemokine receptor 5 (CCR5), interleukin 7 receptor alpha subunit (IL7RA), cytotoxic T-lymphocyte antigen 4 (CTLA4) and HLA class II histocompatibility antigen beta chain (DRB1) was performed with APSampler algorithm for 285 multiple sclerosis patients of Russian ethnicity. The results show evidence for the contribution of polymorphic variants in CCRS, DRB1, IFNG, TGFB1, IFNAR1, IL7RA and, probably, TNF and CTLA4 genes to copaxone treatment response. Single alleles of CCR5 and DRB1 genes are reliably associated with treatment efficacy. Carriage of allelic variants of other above mentioned genes contribute with reliable effect to copaxone treatment response as part of bi- and three-allelic combinations only. Present investigation may support basis toward the future possibility of prognostic test realization, which can provide a personal choice of immunomodulatory treatment for a patient with multiple sclerosis.     

Mapping quantitative trait loci (QTL) in sheep. IV. Analysis of lactation persistency and extended lactation traits in sheep

Background

In sheep dairy production, total lactation performance, and length of lactation of lactation are of economic significance. A more persistent lactation has been associated with improved udder health. An extended lactation is defined by a longer period of milkability. This study is the first investigation to examine the presence of quantitative trait loci (QTL) for extended lactation and lactation persistency in sheep.

Methods

An (Awassi × Merino) × Merino single-sire backcross family with 172 ewes was used to map QTL for lactation persistency and extended lactation traits on a framework map of 189 loci across all autosomes. The Wood model was fitted to data from multiple lactations to estimate parameters of ovine lactation curves, and these estimates were used to derive measures of lactation persistency and extended lactation traits of milk, protein, fat, lactose, useful yield, and somatic cell score. These derived traits were subjected to QTL analyses using maximum likelihood estimation and regression analysis.

Results

Overall, one highly significant (LOD > 3.0), four significant (2.0 < LOD < 3.0) and five suggestive (1.7 < LOD < 2.0) QTL were detected across all traits in common by both mapping methods. One additional suggestive QTL was identified using maximum likelihood estimation, and four suggestive (0.01 < P < 0.05) and two significant (P < 0.01) QTL using the regression approach only. All detected QTL had effect sizes in the range of 0.48 to 0.64 SD, corresponding to QTL heritabilities of 3.1 to 8.9%. The comparison of the detected QTL with results in cattle showed conserved linkage regions. Most of the QTL identified for lactation persistency and extended lactation did not coincide. This suggests that persistency and extended lactation for the same as well as different milk yield and component traits are not controlled by the same genes.

Conclusion

This study identified ten novel QTL for lactation persistency and extended lactation in sheep, but results suggest that lactation persistency and extended lactation do not have a major gene in common. These results provide a basis for further validation in extended families and other breeds as well as targeting regions for genome-wide association mapping using high-density SNP arrays.     

Associations between bovine IGFBP2 polymorphisms with fertility, milk production, and metabolic status in UK dairy cows

Insulin-like growth factor binding protein-2 (IGFBP2) is a key regulator of IGF activity that has been associated with insulin resistance and obesity. In cows, IGFBP2 mRNA expression is differentially regulated according to nutritional status in different tissues including the liver, reproductive tract, and mammary gland. This study investigated associations between single nucleotide polymorphisms (SNPs) in bovine IGFBP2 with fertility, milk production, and metabolic traits in Holstein-Friesian dairy cows. Fertility was assessed in heifers by measuring age at first service, age at first conception, and age at first calving. During the first and second lactation, the number of postpartum days for commencement of luteal activity (based on milk progesterone profiles), days to first service, days to conception, average milk production per day, 305-day milk yield, total milk yield, and total days in milk were recorded. Blood samples were taken at -1, +1, and +8 weeks relative to first and second calving for assessment of metabolic status (IGF1, insulin, beta-hydroxybutyrate, and glucose). Five novel SNPs were identified in IGFBP2, two of which had significant associations with fertility (age at conception in heifers and commencement of luteal activity) and 305-day milk yield in lactation 1. Trends of association were also observed with the peripartum metabolic status, in particular the glucose, insulin, and beta-hydroxybutyrate concentrations around second calving. These results indicate that IGFBP2 SNPs may influence tissue mobilization in dairy cows and may thus be of interest for marker assisted selection.     

Differential expression of genes in milk of dairy cattle during lactation

The milk fat globule (MFG) is one of the most representative of mammary gland tissues and can be utilized to study gene expression of lactating cows during lactation. In this study, RNA‐seq technology was employed to detect differential expression of genes in MFGs at day 10 and day 70 after calving between two groups of cows with extremely high (H group) and low (L group) 305‐day milk yield, milk fat yield and milk protein yield. In total, 1232, 81, 429 and 178 significantly differentially expressed genes (false discovery rate q < 0.05) were detected between H10 and L10, H70 and L70, H10 and H70, and L10 and L70 respectively. Gene Ontology enrichment and pathway analysis revealed that these differentially expressed genes were enriched in biological processes involved in mammary gland development, protein and lipid metabolism process, signal transduction, cellular process, differentiation and immune function. Among these differentially expressed genes, 178 (H10 vs. L10), 4 (H70 vs. L70), 68 (H10 vs. H70) and 22 (L10 vs. L70) were found to be located within previously reported QTL regions for milk production traits. Based on these results, some promising candidate genes for milk production traits in dairy cattle were suggested.     

Response profiles of dairy cows to a single 24-h milking interval in relation with milk proteolysis,udder expansion and immune traits

An extended milking interval of 24 h (24-h milking interval (24h-MI)) constitutes the acute phase of cow adaptation to once-daily milking (ODM). A recent trial including 724 24h-MI challenges demonstrated that milk yield responses to this acute phase of ODM are highly variable (from+22% to −52% of milk yield when switching to the 24h-MI, mean=−25.3%) and that factors such as stage of lactation parity and milk yield level influenced cows’ responses but did not account for all individual variability. Additional traits related to physiological, immune and behavioural adaptation were measured on a subset (96 observations) of this data set. This study aimed to determine (1) the relationship of these traits with cows’ milk yield responses, (2) their ability – combined with previously identified traits – to help predict milk yield responses to 24h-MI (adaptive profiles). The 24h-MI challenge consisted of three successive periods: one control week of twice-daily milking (cTDM), one single day of 24h-MI and then 13 days of TDM (pTDM). Milk yield responses to the 24h-MI (corrected for effects of stage of lactation, parity, milk yield level and milk yield) were related to physiological traits measured during cTDM (milk flow rate, presence or absence of interleukin-8) and to their changes during the 24h-MI (absolute increase in milk flow rate and relative udder distension). Analysis of associations between milk yield responses, stage of lactation, parity, milk yield level, proteolysis, udder expansion and immune traits found three adaptive cow profile clusters. Cows in cluster 1 had a less compliant udder than cows in cluster 2, and they lost more milk during the 24h-MI than cluster-2 and cluster-3 cows. After resuming twice daily-milking (TDM), cluster-2 cows fully recovered the milk they had lost during the 24h-MI. On the opposite, cluster-3 cows did not recover the milk they lost, likely due to udder inflammation during cTDM, as suggested by elevated concentrations of interleukin-8 in their milk. These results combining new traits with stage of lactation, parity and milk yield level constitute a first step towards predicting individual cow responses to a 24h-MI.     

Genetic parameters for milk urea concentration and milk traits in Polish Holstein-Friesian cows

Milk urea concentration (MU) used by dairy producers for management purposes can be affected by selection for milk traits. To assess this problem, genetic parameters for MU in Polish Holstein-Friesian cattle were estimated for the first three lactations. The genetic correlation of MU with milk production traits, lactose percentage, fat to protein ratio (FPR) and somatic cell score (SCS) were computed with two 5-trait random regression test-day models, separately for each lactation. Data used for estimation (159,044 daily observations) came from 50 randomly sampled herds. (Co)variance components were estimated with the Bayesian Gibbs sampling method. The coefficient of variation for MU in all three parities was high (40–41 %). Average daily heritabilities of MU were 0.22 for the first parity and 0.21 for the second and third lactations. Average genetic correlations for different days in milk in the first three lactations between MU and other traits varied. They were small and negative for protein percentage (from ?0.24 to ?0.11) and for SCS (from ?0.14 to ?0.09). The weakest genetic correlation between MU and fat percentage, and between MU and lactose percentage were observed (from ?0.10 to 0.10). Negative average genetic correlation with the fat to protein ratio was observed only in the first lactation (?0.14). Genetic correlations with yield traits were positive and ranged from low to moderate for protein (from 0.09 to 0.33), fat (from 0.16 to 0.35) and milk yield (from 0.20 to 0.42). These results suggest that the selection on yield traits and SCS tends to increase MU slightly.     

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.