Screening for phenotypic outliers identifies an unusually low concentration of a ß-lactoglobulin B protein isoform in bovine milk caused by a synonymous SNP.

BACKGROUND: Milk samples from 10,641 dairy cattle were screened by a mass spectrometry method for extreme concentrations of the A or B isoforms of the whey protein, ß-lactoglobulin (BLG), to identify causative genetic variation driving changes in BLG concentration. RESULTS: A cohort of cows, from a single sire family, was identified that produced milk containing a low concentration of the BLG B protein isoform. A genome-wide association study (GWAS) of BLG B protein isoform concentration in milk from AB heterozygous cows, detected a group of highly significant single nucleotide polymorphisms (SNPs) within or close to the BLG gene. Among these was a synonymous G/A variation at position + 78 bp in exon 1 of the BLG gene (chr11:103256256G > A). The effect of the A allele of this SNP (which we named B') on BLG expression was evaluated in a luciferase reporter assay in transfected CHO-K1 and MCF-7 cells. In both cell types, the presence of the B' allele in a plasmid containing the bovine BLG gene from -922 to + 898 bp (relative to the transcription initiation site) resulted in a 60% relative reduction in mRNA expression, compared to the plasmid containing the wild-type B sequence allele. Examination of a mammary RNAseq dataset (n = 391) identified 14 heterozygous carriers of the B' allele which were homozygous for the BLG B protein isoform (BB'). The level of expression of the BLG B' allele was 41.9 ± 1.0% of that of the wild-type BLG B allele. Milk samples from three cows, homozygous for the A allele at chr11:103,256,256 (B'B'), were analysed (HPLC) and showed BLG concentrations of 1.04, 1.26 and 1.83 g/L relative to a mean of 4.84 g/L in milk from 16 herd contemporaries of mixed (A and B) BLG genotypes. The mechanism by which B' downregulates milk BLG concentration remains to be determined. CONCLUSIONS: High-throughput screening and identification of outliers, enabled the discovery of a synonymous G > A mutation in exon 1 of the B allele of the BLG gene (B'), which reduced the milk concentration of ß-lactoglobulin B protein isoform, by more than 50%. Milk from cows carrying the B' allele is expected to have improved processing characteristics, particularly for cheese-making.

Metabolic proteomics of the liver and mammary gland during lactation.

The liver and the mammary gland have complementary metabolic roles during lactation. Glucose synthesized by the liver is released into the circulation and is taken up by the mammary gland where major metabolic products of glucose include milk sugar (lactose) and the glycerol backbone of milk fat (triglycerides). Hepatic synthesis of glucose is often accompanied by ß-oxidation in that organ to provide energy for glucose synthesis, while mammary gland synthesizes rather than oxidizes fat during lactation. We have therefore compared enzyme abundances between the liver and mammary gland of lactating Friesian cows where metabolic output is well established. Quantitative differences in protein amount were assessed using two-dimensional differential in-gel electrophoresis. As predicted, the abundances of enzymes catalysing gluconeogenesis and ß-oxidation were greatest in the liver, and enzyme abundances in mammary tissue were consistent with fat synthesis rather than ß-oxidation.

Profiling the metabolic proteome of bovine mammary tissue.

2-DE and MALDI mass fingerprinting were used to analyse mammary tissue from lactating Friesian cows. The goal was detection of enzymes in metabolic pathways for synthesis of milk molecules including fatty acids and lactose. Of 418 protein spots analysed by PMF, 328 were matched to database sequences, resulting in 215 unique proteins. We detected 11 out of the 15 enzymes in the direct pathways for conversion of glucose to fatty acids, two of the pentose phosphate pathway enzymes and two of the enzymes for lactose synthesis from glucose. We did not detect enzymes that catalyse the first three reactions of glycolysis. Our results are typical of enzyme detection using 2-DE of mammalian tissues. We therefore advocate caution when relating enzyme abundances measured by 2-DE to metabolic output as not all relevant proteins are detected. 2-D DIGE was used to measure interindividual variation in enzyme abundance from eight animals. We extracted relative protein abundances from 2-D DIGE data and used a logratio transformation that is appropriate for compositional data of the kind represented in many proteomics experiments. Coefficients of variation for abundances of detected enzymes were 3-8%. We recommend use of this transformation for DIGE and other compositional data.