Search

Additional file 1: Figure S1. Correlation between the concentrations of BLG A and B protein isoforms obtained by the mass spectrometry (MS) screening method with those obtained by a high-performance liquid chromatography (HPLC) method [18]. Cows (n = 40) were a mix of heterozygous AB and homozygous AA and BB individuals.

Additional file 3 Supplementary 3 Ensembl Variant Effect Predictor (VEP) calculated variant consequences on the 340 variants identified from whole genome sequencing of 302 goats

Additional file 1 Supplementary 1 Full protein sequence alignment between species across several taxa indicating the degree of conservation around the SENP3 E89del (A), PSMB6 V222I (B) and SHBG S267I/S277I/S281I/S166I (C) variant sites (outlined in red). Multiple sequence alignment was performed with ClustalW. NCBI or UniProt accession numbers provided in brackets. Sequence logos were created using WebLogo web based application [34]

Background Economically important milk production traits including milk volume, milk fat and protein yield vary considerably across dairy goats in New Zealand. A significant portion of the variation is attributable to genetic variation. Discovery of genetic markers linked to milk production traits can be utilised to drive selection of high-performance animals. A previously reported genome wide association study across dairy goats in New Zealand identified a quantitative trait locus (QTL) located on chromosome 19. The most significantly associated single nucleotide polymorphism (SNP) marker for this locus is located at position 26,610,610 (SNP marker rs268292132). This locus is associated with multiple milk production traits including fat, protein and volume. The predicted effect of selection for the beneficial haplotype would result in an average production increase of 2.2 kg fat, 1.9 kg protein and 73.6 kg milk yield. An outstanding question was whether selection for the beneficial allele would co-select for any negative pleiotropic effects. An adverse relationship between milk production and udder health traits has been reported at this locus. Therefore, a genome wide association study was undertaken looking for loci associated with udder traits. Results The QTL and production associated marker rs268292132 was identified in this study to also be associated with several goat udder traits including udder depth (UD), fore udder attachment (FUA) and rear udder attachment (RUA). Our study replicates the negative relationship between production and udder traits with the high production allele at position 19:26,610,610 (SNP marker rs268292132) associated with an adverse change in UD, FUA and RUA. Conclusions Our study has confirmed the negative relationship between udder traits and production traits in the NZ goat population. We have found that the frequency of the high production allele is relatively high in the NZ goat population, indicating that its effect on udder conformation is not significantly detrimental on animal health. It will however be important to monitor udder conformation as the chromosome 19 locus is progressively implemented for marker assisted selection. It will also be of interest to determine if the gene underlying the production QTL has a direct effect on mammary gland morphology or whether the changes observed are a consequence of the increased milk volume.

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.

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.

To identify quantitative trait loci (QTL) affecting the concentration of beta-lactoglobulin in milk, and to evaluate the effect of beta-lactoglobulin genetic variants on the concentration of fat, protein and casein in bovine milk.A herd of 850 F2 Holstein-Friesian x Jersey crossbred cows was produced through mating six Holstein-Friesian x Jersey F1 bulls of high genetic merit with F1 cows from the national herd. A total of 1,610 herd-test records from 556 second-parity crossbreds were analysed. The concentration of fat, protein and casein in milk was measured at peak, mid- and late lactation, during the production seasons of 2003-2004 and 2004-2005. Liveweight was measured daily. DNA from the F2 animals, their F1 dams and sires, and selected grandsires was genotyped across the genome, initially with 285 microsatellite markers, and subsequently with 6,634 single nucleotide polymorphisms (SNP).A highly significant QTL for the concentration of beta-lactoglobulin in milk was identified, which coincided with the position of the beta-lactoglobulin gene on bovine Chromosome 11. No other consistently significant QTL for the concentration of beta-lactoglobulin in milk were detected. Cows with the BB beta-lactoglobulin genotype produced milk with a 30% lower concentration of beta-lactoglobulin than cows with the AA genotype. The beta-lactoglobulin polymorphism also explained variation in the proportion of casein in total protein. In addition, the percentage of fat was higher for BB than AA animals, whereas the percentage of total protein, mean daily milk yield and liveweight did not differ between AA and BB animals.A significant QTL determining the concentration of beta-lactoglobulin in milk was identified. Selection of animals for the beta-lactoglobulin B-allele may enable the production of milk naturally enriched for casein, thus allowing a potential increase in the yield of cheese. There may be additional future value in production of bovine milk more like human milk, where decreasing the concentration of beta-lactoglobulin is desirable.

β-Carotene biochemistry is a fundamental process in mammalian biology. Aberrations either through malnutrition or potentially through genetic variation may lead to vitamin A deficiency, which is a substantial public health burden. In addition, understanding the genetic regulation of this process may enable bovine improvement. While many bovine QTL have been reported, few of the causative genes and mutations have been identified. We discovered a QTL for milk β-carotene and subsequently identified a premature stop codon in bovine β-carotene oxygenase 2 (BCO2), which also affects serum β-carotene content. The BCO2 enzyme is thereby identified as a key regulator of β-carotene metabolism.