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.

Phenotypic population screen identifies a new mutation in bovine DGAT1 responsible for unsaturated milk fat.

Selective breeding has strongly reduced the genetic diversity in livestock species, and contemporary breeding practices exclude potentially beneficial rare genetic variation from the future gene pool. Here we test whether important traits arising by new mutations can be identified and rescued in highly selected populations. We screened milks from 2.5 million cows to identify an exceptional individual which produced milk with reduced saturated fat content, and improved unsaturated and omega-3 fatty acid concentrations. The milk traits were transmitted dominantly to her offspring, and genetic mapping and genome sequencing revealed a new mutation in a previously unknown splice enhancer of the DGAT1 gene. Homozygous carriers show features of human diarrheal disorders, and may be useful for the development of therapeutic strategies. Our study demonstrates that high-throughput phenotypic screening can uncover rich genetic diversity even in inbred populations, and introduces a novel strategy to develop novel milks with improved nutritional properties.

A triad of highly divergent polymeric immunoglobulin receptor (PIGR) haplotypes with major effect on IgA concentration in bovine milk.

The aim of this study was to determine a genetic basis for IgA concentration in milk of Bos taurus. We used a Holstein-Friesian x Jersey F2 crossbred pedigree to undertake a genome-wide search for QTL influencing IgA concentration and yield in colostrum and milk. We identified a single genome-wide significant QTL on chromosome 16, maximising at 4.8 Mbp. The polymeric immunoglobulin receptor gene (PIGR) was within the confidence interval of the QTL. In addition, mRNA expression analysis revealed a liver PIGR expression QTL mapping to the same locus as the IgA quantitative trait locus. Sequencing and subsequent genotyping of the PIGR gene revealed three divergent haplotypes that explained the variance of both the IgA QTL and the PIGR expression QTL. Genetic selection based on these markers will facilitate the production of bovine herds producing milk with higher concentrations of IgA.