Casein haplotypes and their association with milk production traits in the Finnish Ayrshire cattle

Polymorphism of casein genes was studied in half-sib families of artificial insemination bulls of the Finnish Ayrshire dairy breed. Ten grandsires and 300 of their sons were genotyped for the following polymorphisms: α_S1-casein (B, C), β-casein (A¹, A²), the microsatellite within the K-casein gene (ms⁵, ms⁴) and K-casein (A, B, E). Nine different combinations of these alleles, casein haplotypes, were found. Associations between casein haplotypes and milk production traits (milk and protein yield, fat and protein percentage and milking speed) were studied with ordinary least-squares analysis to find a direct effect of the haplotypes or an association within individual grandsire families using the granddaughter design. Estimated breeding values of sons were obtained from cow evaluation by animal model. No direct effect of the casein haplotypes on the traits was found. Within grandsire families, in one out of four families the chromosomal segment characterized by haplotype 3 (B-A²-ms⁴-A) was associated with an increase in milk yield ( P <0.01) and a decrease in fat percentage ( P < 0.01) when contrasted with haplotype 8 (B-A¹-ms⁴-E). The results provide evidence that in the Finnish Ayrshire breed at least one quantitative trait locus affecting the genetic variation in yields traits is segregating linked to either haplotype 3 (B-A²-ms⁴-A) or 8 (B-A¹-ms⁴-E).     

Genetic variation in BoLA microsatellite loci in Portuguese cattle breeds

Major histocompatibility complex (MHC) typing based on microsatellites can be a valuable approach to understanding the selective processes occurring at linked or physically close MHC genes and can provide important information on variability and relationships of populations. Using microsatellites within or in close proximity with bovine lymphocyte antigen (BoLA) genes, we investigated the polymorphisms in the bovine MHC, known as the BoLA, in eight Portuguese cattle breeds. Additional data from non-BoLA microsatellite loci were also used to compare the variability between these regions. Diversity was higher in BoLA than in non-BoLA microsatellites, as could be observed by the number of alleles, allelic richness and observed heterozygosity. Brava de Lide, a breed selected for aggressiveness and nobility, presented the lowest values of observed heterozygosity and allelic richness in both markers. Results from neutrality tests showed few statistically significant differences between the observed Hardy–Weinberg homozygosity ( F ) and the expected homozygosity ( F _E), indicating the apparent neutrality of the BoLA microsatellites within the analysed breeds. Nevertheless, we detected a trend of lower values of observed homozygosity compared with the expected one. We also detected some differences in the levels of allelic variability among the four BoLA microsatellites. Our data showed a higher number of alleles at the BoLA-DRB3 locus than at the BoLA-DRBP1 locus. These differences could be related to their physical position in the chromosome and may reflect functional requirements for diversity.     

Genetic diversity and differentiation in Portuguese cattle breeds using microsatellites

Genotype data from 30 microsatellites were used to assess genetic diversity and relationships among 10 native Portuguese cattle breeds, American Charolais and the Brazilian Caracú. Hardy–Weinberg equilibrium was observed for all loci/population combinations except for five loci in Brava de Lide and one locus in Alentejana that exhibited heterozygote deficiency. Estimates of average observed and expected heterozygosities, total number of alleles (TNA) per breed and mean number of alleles (MNA) per locus/population were obtained. A total of 390 alleles were detected. TNA among Iberian cattle ranged from 170 to 237 and MNA ranged from 5.67 to 8.07. The highest observed heterozygosities were found in the Caracú, Maronesa, Garvonesa and Arouquesa and the lowest in Brava de Lide and Mirandesa. Estimation of population subdivision using Wright's F_ST index showed that the average proportion of genetic variation explained by breed differences was 9%. Neighbour‐joining phylogenetic trees based on D_A distances showed that the genetic relationships of present‐day Portuguese native breeds are consistent with historical origins in the Brown Concave (Arouquesa, Mirandesa, Marinhoa) and Red Convex (Mertolenga, Alentejana, Garvonesa, Minhota) evolutionary groups. The Iberian Black Orthoide group, represented by Brava de Lide and Maronesa, and the Barrosã breed appeared to be more closely related to the Brown Concave group but may represent a separate lineage. The Caracú breed was not found to be closely associated with any of the native Portuguese breeds.     

Comparison of milk protein allele frequencies in Nordic cattle breeds

Allele frequencies at four milk protein loci were studied in five modern and 17 old Nordic cattle breeds in order to reveal variants that are characteristic for these populations. The B allele of CSN3, which has been associated with improved manufacturing properties of milk, showed significantly lower frequencies in modern production breeds than in old breeds of interest for conservation purposes. Characteristic frequencies of CSN1S1 (C), CSN2 (A2) and CSN3 (B) were found in Icelandic cattle, Swedish Mountain cattle, Northern Finncattle and Western Fjord cattle, which indicate a common origin of these populations. Further comparisons of allele frequencies in old Nordic breeds suggest sorting of these breeds into two groups with a northern and southern geographic location.     

Contributions of Portuguese cattle breeds to genetic diversity using marker-estimated kinships

The quantitative assessment of genetic diversity within and between populations is important for decision-making in genetic conservation plans. In our study, we applied the livestock core set method to define the contribution of 15 cattle breeds, 11 of which are Portuguese indigenous cattle breeds, to genetic diversity. In livestock core set theory genetic diversity is defined as the maximum genetic variance that can be obtained in a random-mating population that is bred from the populations present in that core set. Two methods to estimate marker-estimated kinships to obtain the contributions to the core set were used in this study: the weighted log-linear model (WLM) and the weighted log-linear mixed model (WLMM). The breeds that contributed most to diversity in the core set were Holstein-Friesian followed by the Portuguese Mertolenga and Cachena for both WLM and WLMM methods. The ranking of relative contributions of cattle breeds was maintained when we considered only the Portuguese cattle breeds. Furthermore, we were able to identify the marginal contributions and respective losses of diversity for each of the 11 Portuguese cattle breeds when we considered a subset of populations that are not threatened of being lost (the Safe set composed of the four exotic breeds present in this study). When WLM was used losses in genetic diversity ranged from 2.68 to 0.65% while the loss in founder genome equivalents ranged from 37.37 to 8.43% for Mertolenga and Brava de Lide breeds respectively. When WLMM was used losses in genetic diversity and founder genome equivalents were less extreme than for the WLM method, ranging from 1.27 to 0.69 and 26.8 to 12.99 respectively.     

Assignment of the casein kinase II gene family to cattle chromosomes

Theileriosis, or East Coast fever, a parasitic disease in cattle, is associated with overexpression of casein kinase II. Casein kinase II is composed of two catalytic subunits (α or α') and two regulatory β subunits. The genes encoding these subunits of casein kinase II were mapped to bovine chromosomes by polymerase chain reaction analysis of a well-characterized bovine × rodent somatic hybrid cell panel. The α-subunit (CSNK2A1) was mapped to bovine chromosome 13, the α'-subunit (CSNK2A2) to chromosome 5 and the β-subunit (CSNK2B) to chromosome 23. Both CSNK2A1 and CSNK2B mapped to known regions of conserved synteny between human and cattle, while CSNK2A2 defined a new homology segment between the human and bovine genomes.     

Goat casein genotypes are associated with milk production traits in the Sarda breed

The aim of the current work was to analyze, in the Sarda breed goat, genetic polymorphism within the casein genes and to assess their influence on milk traits. Genetic variants at the CSN1S1, CSN2, CSN1S2 and CSN3 gene loci were investigated using PCR‐based methods, cloning and sequencing. Strong alleles prevailed at the CSN1S1 gene locus and defective alleles also were revealed. Null alleles were evidenced at each calcium‐sensitive gene locus. At the CSN3 gene locus, we observed a prevalence of the CSN3 A and B alleles; the occurrence of rare alleles such as CSN3 B'', C, C', D, E and M; and the CSN3 S allele (GenBank KF644565 ) described here for the first time in Capra hircus. Statistical analysis showed that all genes, except CSN3, significantly influenced milk traits. The CSN1S1 BB and AB genotypes were associated with the highest percentages of protein (4.41 and 4.40 respectively) and fat (5.26 and 5.34 respectively) (P < 0.001). A relevant finding was that CSN2 and CSN1S2 genotypes affected milk protein content and yield. The polymorphism of the CSN2 gene affected milk protein percentage with the highest values recorded in the CSN2 AA goats (4.35, at P < 0.001). The CSN1S2 AC goats provided the highest fat (51.02 g/day) and protein (41.42 g/day) (P < 0.01) production. This information can be incorporated into selection schemes for the Sarda breed goat.     

Validation of fatty acid predictions in milk using mid-infrared spectrometry across cattle breeds

The aim of this study was to investigate the accuracy to predict detailed fatty acid (FA) composition of bovine milk by mid-infrared spectrometry, for a cattle population that partly differed in terms of country, breed and methodology used to measure actual FA composition compared with the calibration data set. Calibration equations for predicting FA composition using mid-infrared spectrometry were developed in the European project RobustMilk and based on 1236 milk samples from multiple cattle breeds from Ireland, Scotland and the Walloon Region of Belgium. The validation data set contained 190 milk samples from cows in the Netherlands across four breeds: Dutch Friesian, Meuse-Rhine-Yssel, Groningen White Headed (GWH) and Jersey (JER). The FA measurements were performed using gas–liquid partition chromatography (GC) as the gold standard. Some FAs and groups of FAs were not considered because of differences in definition, as the capillary column of the GC was not the same as used to develop the calibration equations. Differences in performance of the calibration equations between breeds were mainly found by evaluating the standard error of validation and the average prediction error. In general, for the GWH breed the smallest differences were found between predicted and reference GC values and least variation in prediction errors, whereas for JER the largest differences were found between predicted and reference GC values and most variation in prediction errors. For the individual FAs 4:0, 6:0, 8:0, 10:0, 12:0, 14:0 and 16:0 and the groups’ saturated FAs, short-chain FAs and medium-chain FAs, predictions assessed for all breeds together were highly accurate (validation R² > 0.80) with limited bias. For the individual FAs cis-14:1, cis-16:1 and 18:0, the calibration equations were moderately accurate (R² in the range of 0.60 to 0.80) and for the individual FA 17:0 predictions were less accurate (R² < 0.60) with considerable bias. FA concentrations in the validation data set of our study were generally higher than those in the calibration data. This difference in the range of FA concentrations, mainly due to breed differences in our study, can cause lower accuracy. In conclusion, the RobustMilk calibration equations can be used to predict most FAs in milk from the four breeds in the Netherlands with only a minor loss of accuracy.     

Evidence of positive selection towards Zebuine haplotypes in the BoLA region of Brangus cattle

The Brangus breed was developed to combine the superior characteristics of both of its founder breeds, Angus and Brahman. It combines the high adaptability to tropical and subtropical environments, disease resistance, and overall hardiness of Zebu cattle with the reproductive potential and carcass quality of Angus. It is known that the major histocompatibility complex (MHC, also known as bovine leucocyte antigen: BoLA), located on chromosome 23, encodes several genes involved in the adaptive immune response and may be responsible for adaptation to harsh environments. The objective of this work was to evaluate whether the local breed ancestry percentages in the BoLA locus of a Brangus population diverged from the estimated genome-wide proportions and to identify signatures of positive selection in this genomic region. For this, 167 animals (100 Brangus, 45 Angus and 22 Brahman) were genotyped using a high-density single nucleotide polymorphism array. The local ancestry analysis showed that more than half of the haplotypes (55.0%) shared a Brahman origin. This value was significantly different from the global genome-wide proportion estimated by cluster analysis (34.7% Brahman), and the proportion expected by pedigree (37.5% Brahman). The analysis of selection signatures by genetic differentiation (F_st) and extended haplotype homozygosity-based methods (iHS and Rsb) revealed 10 and seven candidate regions, respectively. The analysis of the genes located within these candidate regions showed mainly genes involved in immune response-related pathway, while other genes and pathways were also observed (cell surface signalling pathways, membrane proteins and ion-binding proteins). Our results suggest that the BoLA region of Brangus cattle may have been enriched with Brahman haplotypes as a consequence of selection processes to promote adaptation to subtropical environments.     

Genetic heterogeneity at the bovine KIT gene in cattle breeds carrying different putative alleles at the spotting locus

According to classical genetic studies, piebaldism in cattle is largely influenced by the allelic series at the spotting locus (S), which includes the S^H (Hereford pattern), S⁺ (non‐spotted) and s (spotted) alleles. The S locus was mapped on bovine chromosome 6 in the region containing the KIT gene. We investigated the KIT gene, analysing its variability and haplotype distribution in cattle of three breeds (Angus, Hereford and Holstein) with different putative alleles (S⁺, S^H and s respectively) at the S locus. Resequencing of a whole of 0.485 Mb revealed 111 polymorphisms. The global nucleotide diversity was 0.087%. Tajima’s D‐values were negative for all breeds, indicating putative directional selection. Of the 28 inferred haplotypes, only five were observed in the Hereford breed, in which one was the most frequent. Coalescent simulation showed that it is highly unlikely (P < 10E‐6) to obtain this low number of haplotypes conditionally on the observed number of segregating SNPs. Therefore, the neutral model could be rejected for the Hereford breed, suggesting that a selection sweep occurred at the KIT locus. Twelve haplotypes were inferred in Holstein and Angus. For these two breeds, the neutral model could not be rejected. High heterogeneity of the KIT gene was confirmed from a phylogenetic analysis. Our results suggest a role of the KIT gene in determining the S^H allele(s) in the Hereford, but no evidence of selective sweep was obtained in Holstein, suggesting that complex mechanisms (or other genes) might be the cause of the spotted phenotype in this breed.     

Detection of a new k-casein variant in milk of Pinzgauer cattle

A new k casein variant (k-CN G) with a frequency of 0.003 was demonstrated in Pinzgauer cattle from Austria and Bavaria, Germany by iso-electric focusing in polyacrylamide gels and by alkaline polyacrylamide gel electrophoresis. k-CN G was not present in milk samples of Limpurger, another endangered breed.     

A meta-analysis on association between CSN3 gene variants and milk yield and composition in cattle

A meta-analysis on the effects of A and B alleles, the most frequent alleles of CSN3 gene, on milk yield and composition traits was conducted by pooling a large dataset consisting of 30 471 genotyped cattle. Four genetic models, comprising dominant (AA + AB vs. BB), recessive (AA vs. AB + BB), additive (AA vs. BB) and co-dominant (AA + BB vs. AB), were employed to analyze data. Standardized mean difference (SMD) was used to measure the size of the effects of A and B alleles of CSN3 on studied traits. Effect sizes of 0.2, 0.5 and 0.8 represent small, medium and large effects, respectively. The results indicate that B allele, in the form of BB genotype, has a significant, but medium effect on lactation yield under dominant (SMD = 0.259, P-value = 0.006) and additive (SMD = 0.279, P-value = 0.035) models. Moreover, a small decrease in the fat percentage occurred in cows having A allele under dominant (SMD = −0.077, P-value = 0.006) and additive (SMD = −0.106, P-value = 0.035) models. Furthermore, CSN3 variants significantly but slightly affect protein percentage under dominant (SMD = −0.146, P-value = 0.000), recessive (SMD = −0.077, P-value = 0.000) and additive (SMD = −0.219, P-value = 0.000) models, showing the negative effect of A allele on this trait. Meta-analysis results reveal that daily milk yield is slightly affected by CSN3 variants under recessive (SMD = 0.056, P-value = 0.033) and additive (SMD = 0.061, P-value = 0.013) genetic models. There is no effect of CSN3 variants on either protein or fat yield. In addition, the effects of CSN3 variants on milk-related traits were not observed under the co-dominant model. Sensitivity and publication bias analyses were carried out to confirm the stability of meta-analyses results.     

Frequencies of genes for coat colour and horns in Nordic cattle breeds

Gene frequencies of coat colour and horn types were assessed in 22 Nordic cattle breeds in a project aimed at establishing genetic profiles of the breeds under study. The coat colour loci yielding information on genetic variation were: extension, agouti, spotting, brindle, dun dilution and colour sided. The polled locus was assessed for two alleles. A profound variation between breeds was observed in the frequencies of both colour and horn alleles, with the older breeds generally showing greater variation in observed colour, horn types and segregating alleles than the modern breeds. The correspondence between the present genetic distance matrix and previous molecular marker distance matrices was low (r = 0.08 – 0.12). The branching pattern of a neighbour-joining tree disagreed to some extent with the molecular data structure. The current data indicates that 70% of the total genetic variation could be explained by differences between the breeds, suggesting a much greater breed differentiation than typically found at protein and microsatellite loci. The marked differentiation of the cattle breeds and observed disagreements with the results from the previous molecular data in the topology of the phylogenetic trees are most likely a result of selection on phenotypic characters analysed in this study.     

Bioavailability of zinc and its binding to casein in milks and formulas

Differences in zinc bioavailability among milk and formulas may be attributed to binding of zinc to various ligands. We determined the distribution of zinc and protein at different pHs and zinc and calcium concentrations. We used radiolabelled cow's milk, human milk, whey-predominant (WPF) and casein-predominant (CPF) infant formula. Lowering the pH changed zinc and protein distribution: zinc shifted from pellet (casein) to whey in cow's milk, from fat to whey in human milk and from fat and pellet to whey in formulas. Protein shifted from whey to pellet in human milk and from whey and pellet to fat in formulas. Increasing zinc and calcium concentrations shifted protein and zinc from pellet to whey for cow's milk and from whey and pellet to fat for the formulas. Protein distribution was not affected by calcium or zinc addition in human milk or CPF, while zinc shifted from whey to fat in human milk and from fat and pellet to whey in CPF. Zinc and calcium binding to isolated bovine or human casein increased with pH. At 500 mg/L of zinc, bovine casein bound 32.0 +/- 1.8 and human casein 10.0 +/- 0.9 mg zinc/g protein. At 500 mg/L of calcium, calcium was preferentially bound over zinc. Adding calcium and zinc resulted in 32.0 +/- 1.8 mg zinc/g bound to bovine casein and 17.0 +/- 0.8 mg zinc/g to human casein, while calcium binding was low. Suckling rat pups dosed with 65Zn labelled infant diets were killed and individual tissues were gamma counted. Lower zinc bioavailability was found for bovine milk at pH = 4.0 (%65Zn in liver = 18.7+1.4) when compared to WPF (22.8 +/- 1.6) or human milk (26.9 +/- 0.8). Lowering the pH further decreased zinc bioavailability from human milk, but not from cow's milk or WPF. Knowledge of the compounds binding minerals and trace elements in infant formulas is essential for optimizing zinc bioavailability.     

Extractive bioconversions in aqueous two-phase systems: enzymatic hydrolysis of casein proteins

Partition coefficients in poly(ethylene glycol)/dextran aqueous two-phase systems are reported for mixed-casein and its components, alpha, beta and kappa casein. Rates of casein proteolysis by alpha-chymotrypsin and by trypsin are reported in single-phase and aqueous two-phase reactor systems. The advantages resulting from selective partitioning of substrates, enzymes, and products are examined in terms of relative volumetric reaction rates.     

Linkage between amylase-1 locus and a major gene for milk fat content in cattle

Summary
Linkage between the amylase-1 ( Am-1 ) locus and a quantitative trait locus influencing fat content in milk was studied in offspring from heterozygous sires of the Swedish Red and White dairy breed. The effect on bull breeding values for fat content was estimated as interactions between sire and paternal Am-1 allele using a model eliminating the direct effects of sire and Am-1 allele. There were strong indications of linkage, confirming results of previous studies. The interaction was caused by strong associations in 7 out of 14 sire families. A test for within-family variance heterogeneity performed on the whole population of breeding bulls also supported the presence of a major gene for fat content in milk. The results indicate that there is genetic linkage between the Am-1 locus and a locus with large effect on milk fat content.     

Research on association between variants and haplotypes of TRPV1 and TRPA1 genes with growth traits in three cattle breeds

The objective of this study was to examine the association of transient receptor potential vanilloid 1 (TRPV1) and transient receptor potential ankyrin 1 (TRPA1) genes polymorphisms with growth traits in three Chinese cattle breeds (Jiaxian red cattle, Qinchuan cattle and Luxi cattle). Through experiments we identified three single nucleotide polymorphisms (SNPs) in these three cattle breeds TRPV1 and TRPA1 genes using PCR-SSCP, (forced) PCR-RFLP methods. Three of these polymorphisms are all synonymous mutation which includes (NW_003104493.1: 30327 C?>?T), (NW_003104493.1: 33394 A?>?G) and (NW_003104493.1: 38471?G?>?A) are in exons. The other three polymorphisms are located at 3'UTR. Furthermore, we evaluated the haplotype frequency and the statistical analyses indicated that these SNPs of TRPV1 and TRPA1 genes were associated with bovine body height, body length, waist angle width, hucklebone width, cross ministry height, chest qingwidth (p?<?0.05) and recommendation height, cannon circumference (p?<?0.01) of Qingchuan cattle; body length, waist angle width (p?<?0.05) of Jiaxian red cattle; body weight, Body length, cannon circumference, chest circumference (p?<?0.05) and body height (p?<?0.01) of Luxi cattle. Our result confirms the polymorphisms in the TRPV1 and TRPA1 genes are associated with growth traits that may be used for marker-assisted selection (MAS) in three beef breeds programs.     

Assessing priorities for conservation in Tuscan cattle breeds using microsatellites

Preservation of rare genetic stocks requires assessment of within-population genetic diversity and between-population differentiation to make inferences on their degree of uniqueness. A total of 194 Tuscan cattle (44 Calvana, 35 Chianina, 25 Garfagnina, 31 Maremmana, 31 Mucca Pisana and 28 Pontremolese) individuals were genotyped for 34 microsatellite markers. Moreover, 56 samples belonging to Argentinean Creole and Asturiana de la Montaña cattle breeds were used as an outgroup. Genetic diversity was quantified in terms of molecular coancestry and allelic richness. STRUCTURE analyses showed that the Tuscan breeds have well-differentiated genetic backgrounds, except for the Calvana and Chianina breeds, which share the same genetic ancestry. The between-breed Nei's minimum distance (D_m) matrices showed that the pair Calvana–Chianina was less differentiated (0.049 ± 0.006). The endangered Tuscan breeds (Calvana, Garfagnina, Mucca Pisana and Pontremolese) made null or negative contributions to diversity, except for the Mucca Pisana contribution to allelic richness (C_T = 1.8%). The Calvana breed made null or negative within-breed contributions (f¯_W = 0.0%; C_W = −0.4%). The Garfagnina and Pontremolese breeds made positive contributions to between-breed diversity but negative and high within-breed contributions, thus suggesting population bottleneck with allelic losses and increase of homozygosity in the population. Exclusion of the four endangered Tuscan cattle breeds did not result in losses in genetic diversity (f¯_T = −0.7%; C_T = −1.2%), whereas exclusion of the non-endangered breeds (Chianina and Maremmana) did (f¯_T = 2.1%; C_T = 3.9%); the simple exclusion of the Calvana breed from the former group led to losses in genetic diversity (f¯_T = 0.47%; C_T = 2.34%), indicating a diverse significance for this breed. We showed how quantifying both within-population diversity and between-population differentiation in terms of allelic frequencies and allelic richness provides different and complementary information on the genetic backgrounds assessed and may help to implement priorities and strategies for conservation in livestock.     

Estimating the frequency of Asian cytochrome B haplotypes in standard European and local Spanish pig breeds

Mitochondrial DNA has been widely used to perform phylogenetic studies in different animal species. In pigs, genetic variability at the cytochrome B gene and the D-loop region has been used as a tool to dissect the genetic relationships between different breeds and populations. In this work, we analysed four SNP at the cytochrome B gene to infer the Asian (A1 and A2 haplotypes) or European (E1 and E2 haplotypes) origins of several European standard and local pig breeds. We found a mixture of Asian and European haplotypes in the Canarian Black pig (E1, A1 and A2), German Piétrain (E1, A1 and A2), Belgian Piétrain (E1, A1), Large White (E1 and A1) and Landrace (E1 and A1) breeds. In contrast, the Iberian (Guadyerbas, Ervideira, Caldeira, Campanario, Puebla and Torbiscal strains) and the Majorcan Black pig breeds only displayed the E1 haplotype. Our results show that the introgression of Chinese pig breeds affected most of the major European standard breeds, which harbour Asian haplotypes at diverse frequencies (15–56%). In contrast, isolated local Spanish breeds, such as the Iberian and Majorcan Black pig, only display European cytochrome B haplotypes, a feature that evidences that they were not crossed with other Chinese or European commercial populations. These findings illustrate how geographical confinement spared several local Spanish breeds from the extensive introgression event that took place during the 18th and 19th centuries in Europe.