Search

Background: Bovine lactoferrin (Lf) is an iron absorbing whey protein with antibacterial, antiviral, and antifungal activity. Lactoferrin is economically valuable and has an extremely variable concentration in milk, partly driven by environmental influences such as milking frequency, involution, or mastitis. A significant genetic influence has also been previously observed to regulate lactoferrin content in milk. Here, we conducted genetic mapping of lactoferrin protein concentration in conjunction with RNA-seq, ChIP-seq, and ATAC-seq data to pinpoint candidate causative variants that regulate lactoferrin concentrations in milk., Results: We identified a highly-significant lactoferrin protein quantitative trait locus (pQTL), as well as a cis lactotransferrin (LTF) expression QTL (cis-eQTL) mapping to the LTF locus. Using ChIP-seq and ATAC-seq datasets representing lactating mammary tissue samples, we also report a number of regions where the openness of chromatin is under genetic influence. Several of these also show highly significant QTL with genetic signatures similar to those highlighted through pQTL and eQTL analysis. By performing correlation analysis between these QTL, we revealed an ATAC-seq peak in the putative promotor region of LTF, that highlights a set of 115 high-frequency variants that are potentially responsible for these effects. One of the 115 variants (rs110000337), which maps within the ATAC-seq peak, was predicted to alter binding sites of transcription factors known to be involved in lactation-related pathways., Conclusions: Here, we report a regulatory haplotype of 115 variants with conspicuously large impacts on milk lactoferrin concentration. These findings could enable the selection of animals for high-producing specialist herds., (© 2024. The Author(s).)

Fourier-transform mid-infrared (FT-MIR) spectroscopy is a high-throughput and inexpensive methodology used to evaluate concentrations of fat and protein in dairy cattle milk samples. The objective of this study was to compare the genetic characteristics of FT-MIR predicted fatty acids and individual milk proteins with those that had been measured directly using gas and liquid chromatography methods. The data used in this study was based on 2,005 milk samples collected from 706 Holstein-Friesian × Jersey animals that were managed in a seasonal, pasture-based dairy system, with milk samples collected across 2 consecutive seasons. Concentrations of fatty acids and protein fractions in milk samples were directly determined by gas chromatography and high-performance liquid chromatography, respectively. Models to predict each directly measured trait based on FT-MIR spectra were developed using partial least squares regression, with spectra from a random selection of half the cows used to train the models, and predictions for the remaining cows used as validation. Variance parameters for each trait and genetic correlations for each pair of measured/predicted traits were estimated from pedigree-based bivariate models using REML procedures. A genome-wide association study was undertaken using imputed whole-genome sequence, and quantitative trait loci (QTL) from directly measured traits were compared with QTL from the corresponding FT-MIR predicted traits. Cross-validation prediction accuracies based on partial least squares for individual and grouped fatty acids ranged from 0.18 to 0.65. Trait prediction accuracies in cross-validation for protein fractions were 0.53, 0.19, and 0.48 for α-casein, β-casein, and κ-casein, 0.31 for α-lactalbumin, 0.68 for β-lactoglobulin, and 0.36 for lactoferrin. Heritability estimates for directly measured traits ranged from 0.07 to 0.55 for fatty acids; and from 0.14 to 0.63 for individual milk proteins. For FT-MIR predicted traits, heritability estimates were mostly higher than for the corresponding measured traits, ranging from 0.14 to 0.46 for fatty acids, and from 0.30 to 0.70 for individual proteins. Genetic correlations between directly measured and FT-MIR predicted protein fractions were consistently above 0.75, with the exceptions of C18:0 and C18:3 cis-3, which had genetic correlations of 0.72 and 0.74, respectively. The GWAS identified trait QTL for fatty acids with likely candidates in the DGAT1, CCDC57, SCD, and GPAT4 genes. Notably, QTL for SCD were largely absent in the FT-MIR predicted traits, and QTL for GPAT4 were absent in directly measured traits. Similarly, for directly measured individual proteins, we identified QTL with likely candidates in the CSN1S1, CSN3, PAEP, and LTF genes, but the QTL for CSN3 and LTF were absent in the FT-MIR predicted traits. Our study indicates that genetic correlations between directly measured and FT-MIR predicted fatty acid and protein fractions are typically high, but that phenotypic variation in these traits may be underpinned by differing genetic architecture., (The Authors. Published by Elsevier Inc. and Fass Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).)

Accurate and timely pregnancy diagnosis is an important component of effective herd management in dairy cattle. Predicting pregnancy from Fourier-transform mid-infrared (FT-MIR) spectroscopy data is of particular interest because the data are often already available from routine milk testing. The purpose of this study was to evaluate how well pregnancy status could be predicted in a large data set of 1,161,436 FT-MIR milk spectra records from 863,982 mixed-breed pasture-based New Zealand dairy cattle managed within seasonal calving systems. Three strategies were assessed for defining the nonpregnant cows when partitioning the records according to pregnancy status in the training population. Two of these used records for cows with a subsequent calving only, whereas the third also included records for cows without a subsequent calving. For each partitioning strategy, partial least squares discriminant analysis models were developed, whereby spectra from all the cows in 80% of herds were used to train the models, and predictions on cows in the remaining herds were used for validation. A separate data set was also used as a secondary validation, whereby pregnancy diagnosis had been assigned according to the presence of pregnancy-associated glycoproteins (PAG) in the milk samples. We examined different ways of accounting for stage of lactation in the prediction models, either by including it as an effect in the prediction model, or by pre-adjusting spectra before fitting the model. For a subset of strategies, we also assessed prediction accuracies from deep learning approaches, utilizing either the raw spectra or images of spectra. Across all strategies, prediction accuracies were highest for models using the unadjusted spectra as model predictors. Strategies for cows with a subsequent calving performed well in herd-independent validation with sensitivities above 0.79, specificities above 0.91 and area under the receiver operating characteristic curve (AUC) values over 0.91. However, for these strategies, the specificity to predict nonpregnant cows in the external PAG data set was poor (0.002-0.04). The best performing models were those that included records for cows without a subsequent calving, and used unadjusted spectra and days in milk as predictors, with consistent results observed across the training, herd-independent validation and PAG data sets. For the partial least squares discriminant analysis model, sensitivity was 0.71, specificity was 0.54 and AUC values were 0.68 in the PAG data set; and for an image-based deep learning model, the sensitivity was 0.74, specificity was 0.52 and the AUC value was 0.69. Our results demonstrate that in pasture-based seasonal calving herds, confounding between pregnancy status and spectral changes associated with stage of lactation can inflate prediction accuracies. When the effect of this confounding was reduced, prediction accuracies were not sufficiently high enough to use as a sole indicator of pregnancy status., (The Authors. Published by Elsevier Inc. and Fass Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).)

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., (© 2022. The Author(s).)

Background: Deleterious recessive conditions have been primarily studied in the context of Mendelian diseases. Recently, several deleterious recessive mutations with large effects were discovered via non-additive genome-wide association studies (GWAS) of quantitative growth and developmental traits in cattle, which showed that quantitative traits can be used as proxies of genetic disorders when such traits are indicative of whole-animal health status. We reasoned that lactation traits in cattle might also reflect genetic disorders, given the increased energy demands of lactation and the substantial stresses imposed on the animal. In this study, we screened more than 124,000 cows for recessive effects based on lactation traits., Results: We discovered five novel quantitative trait loci (QTL) that are associated with large recessive impacts on three milk yield traits, with these loci presenting missense variants in the DOCK8, IL4R, KIAA0556, and SLC25A4 genes or premature stop variants in the ITGAL, LRCH4, and RBM34 genes, as candidate causal mutations. For two milk composition traits, we identified several previously reported additive QTL that display small dominance effects. By contrasting results from milk yield and milk composition phenotypes, we note differing genetic architectures. Compared to milk composition phenotypes, milk yield phenotypes had lower heritabilities and were associated with fewer additive QTL but had a higher non-additive genetic variance and were associated with a higher proportion of loci exhibiting dominance., Conclusions: We identified large-effect recessive QTL which are segregating at surprisingly high frequencies in cattle. We speculate that the differences in genetic architecture between milk yield and milk composition phenotypes derive from underlying dissimilarities in the cellular and molecular representation of these traits, with yield phenotypes acting as a better proxy of underlying biological disorders through presentation of a larger number of major recessive impacts., (© 2022. The Author(s).)

Background: Fourier-transform mid-infrared (FT-MIR) spectroscopy provides a high-throughput and inexpensive method for predicting milk composition and other novel traits from milk samples. While there have been many genome-wide association studies (GWAS) conducted on FT-MIR predicted traits, there have been few GWAS for individual FT-MIR wavenumbers. Using imputed whole-genome sequence for 38,085 mixed-breed New Zealand dairy cattle, we conducted GWAS on 895 individual FT-MIR wavenumber phenotypes, and assessed the value of these direct phenotypes for identifying candidate causal genes and variants, and improving our understanding of the physico-chemical properties of milk., Results: Separate GWAS conducted for each of 895 individual FT-MIR wavenumber phenotypes, identified 450 1-Mbp genomic regions with significant FT-MIR wavenumber QTL, compared to 246 1-Mbp genomic regions with QTL identified for FT-MIR predicted milk composition traits. Use of mammary RNA-seq data and gene annotation information identified 38 co-localized and co-segregating expression QTL (eQTL), and 31 protein-sequence mutations for FT-MIR wavenumber phenotypes, the latter including a null mutation in the ABO gene that has a potential role in changing milk oligosaccharide profiles. For the candidate causative genes implicated in these analyses, we examined the strength of association between relevant loci and each wavenumber across the mid-infrared spectrum. This revealed shared association patterns for groups of genomically-distant loci, highlighting clusters of loci linked through their biological roles in lactation and their presumed impacts on the chemical composition of milk., Conclusions: This study demonstrates the utility of FT-MIR wavenumber phenotypes for improving our understanding of milk composition, presenting a larger number of QTL and putative causative genes and variants than found from FT-MIR predicted composition traits. Examining patterns of significance across the mid-infrared spectrum for loci of interest further highlighted commonalities of association, which likely reflects the physico-chemical properties of milk constituents., (© 2021. The Author(s).)

Background: The DGAT1 gene encodes an enzyme responsible for catalysing the terminal reaction in mammary triglyceride synthesis, and underpins a well-known pleiotropic quantitative trait locus (QTL) with a large influence on milk composition phenotypes. Since first described over 15 years ago, a protein-coding variant K232A has been assumed as the causative variant underlying these effects, following in-vitro studies that demonstrated differing levels of triglyceride synthesis between the two protein isoforms., Results: We used a large RNAseq dataset to re-examine the underlying mechanisms of this large milk production QTL, and hereby report novel expression-based functions of the chr14 g.1802265AA > GC variant that encodes the DGAT1 K232A substitution. Using expression QTL (eQTL) mapping, we demonstrate a highly-significant mammary eQTL for DGAT1, where the K232A mutation appears as one of the top associated variants for this effect. By conducting in vitro expression and splicing experiments in bovine mammary cell culture, we further show modulation of splicing efficiency by this mutation, likely through disruption of an exon splice enhancer as a consequence of the allele encoding the 232A variant., Conclusions: The relative contributions of the enzymatic and transcription-based mechanisms now attributed to K232A remain unclear; however, these results suggest that transcriptional impacts contribute to the diversity of lactation effects observed at the DGAT1 locus.

The objective of this study was to estimate genetic correlations among milk fatty acid (FA) concentrations in New Zealand dairy cattle. Concentrations of each of the most common FA, expressed as a percentage of the total FA, were determined by gas chromatography on a specific cohort of animals. Using this data set, prediction equations were derived using mid-infrared (MIR) spectroscopy data collected from the same samples. These prediction equations were applied to a large data set of MIR measurements in 34,141 milk samples from 3,445 Holstein-Friesian, 2,935 Jersey, and 3,609 crossbred Holstein-Friesian × Jersey cows, sampled an average of 3.42 times during the 2007-2008 season. Data were analyzed using univariate and bivariate repeatability animal models. Heritability of predicted FA concentration in milk fat ranged from 0.21 to 0.42, indicating that genetic selection could be used to change the FA composition of milk. The de novo synthesized FA (C6:0, C8:0, C10:0, C12:0, and C14:0) showed strong positive genetic correlations with each other, ranging from 0.24 to 0.99. Saturated FA were negatively correlated with unsaturated (-0.93) and polyunsaturated (-0.84) FA. The saturated FA were positively correlated with milk fat yield and fat percentage, whereas the unsaturated FA were negatively associated with fat yield and fat percentage. Our results indicate that bovine milk FA composition can be changed through genetic selection using MIR as a phenotypic proxy., (Copyright © 2020 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

The use of Fourier-transform mid-infrared (FTIR) spectroscopy is of interest to the dairy industry worldwide for predicting milk composition and other novel traits that are difficult or expensive to measure directly. Although there are many valuable applications for FTIR spectra, noise from differences in spectral responses between instruments is problematic because it reduces prediction accuracy if ignored. The purpose of this study was to develop strategies to reduce the impact of noise and to compare methods for standardizing FTIR spectra in order to reduce between-instrument variability in multiple-instrument networks. Noise levels in bands of the infrared spectrum caused by the water content of milk were characterized, and a method for identifying and removing outliers was developed. Two standardization methods were assessed and compared: piecewise direct standardization (PDS), which related spectra on a primary instrument to spectra on 5 other (secondary) instruments using identical milk-based reference samples (n = 918) analyzed across the 6 instruments; and retroactive percentile standardization (RPS), whereby percentiles of observed spectra from routine milk test samples (n = 2,044,094) were used to map and exploit primary- and secondary-instrument relationships. Different applications of each method were studied to determine the optimal way to implement each method across time. Industry-standard predictions of milk components from 2,044,094 spectra records were regressed against predictions from spectra before and after standardization using PDS or RPS. The PDS approach resulted in an overall decrease in root mean square error between industry-standard predictions and predictions from spectra from 0.190 to 0.071 g/100 mL for fat, from 0.129 to 0.055 g/100 mL for protein, and from 0.143 to 0.088 g/100 mL for lactose. Reductions in prediction error for RPS were similar but less consistent than those for PDS across time, but similar reductions were achieved when PDS coefficients were updated monthly and separate primary instruments were assigned for the North and South Islands of New Zealand. We demonstrated that the PDS approach is the most consistent method to reduce prediction errors across time. We also showed that the RPS approach is sensitive to shifts in milk composition but can be used to reduce prediction errors, provided that secondary-instrument spectra are standardized to a primary instrument with samples of broadly equivalent milk composition. Appropriate implementation of either of these approaches will improve the quality of predictions based on FTIR spectra for various downstream applications., (The Authors. Published by FASS Inc. and Elsevier Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).)

In cattle, the X chromosome accounts for approximately 3 and 6% of the genome in bulls and cows, respectively. In spite of the large size of this chromosome, very few studies report analysis of the X chromosome in genome-wide association studies and genomic selection. This lack of genetic interrogation is likely due to the complexities of undertaking these studies given the hemizygous state of some, but not all, of the X chromosome in males. The first step in facilitating analysis of this gene-rich chromosome is to accurately identify coordinates for the pseudoautosomal boundary (PAB) to split the chromosome into a region that may be treated as autosomal sequence (pseudoautosomal region) and a region that requires more complex statistical models. With the recent release of ARS-UCD1.2, a more complete and accurate assembly of the cattle genome than was previously available, it is timely to fine map the PAB for the first time. Here we report the use of SNP chip genotypes, short-read sequences, and long-read sequences to fine map the PAB (X chromosome:133,300,518) and simultaneously determine the neighboring regions of reduced homology and true pseudoautosomal region. These results greatly facilitate the inclusion of the X chromosome in genome-wide association studies, genomic selection, and other genetic analysis undertaken on this reference genome., (The Authors. Published by FASS Inc. and Elsevier Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).)

Background: Over many years, artificial selection has substantially improved milk production by cows. However, the genes that underlie milk production quantitative trait loci (QTL) remain relatively poorly characterised. Here, we investigate a previously reported QTL located at the CSF2RB locus on chromosome 5, for several milk production phenotypes, to better understand its underlying genetic and molecular causes., Results: Using a population of 29,350 taurine dairy cows, we conducted association analyses for milk yield and composition traits, and identified highly significant QTL for milk yield, milk fat concentration, and milk protein concentration. Strikingly, protein concentration and milk yield appear to show co-located yet genetically distinct QTL. To attempt to understand the molecular mechanisms that might be mediating these effects, gene expression data were used to investigate eQTL for 11 genes in the broader interval. This analysis highlighted genetic impacts on CSF2RB and NCF4 expression that share similar association signatures to those observed for lactation QTL, strongly implicating one or both of these genes as responsible for these effects. Using the same gene expression dataset representing 357 lactating cows, we also identified 38 novel RNA editing sites in the 3' UTR of CSF2RB transcripts. The extent to which two of these sites were edited also appears to be genetically co-regulated with lactation QTL, highlighting a further layer of regulatory complexity that involves the CSF2RB gene., Conclusions: This locus presents a diversity of molecular and lactation QTL, likely representing multiple overlapping effects that, at a minimum, highlight the CSF2RB gene as having a causal role in these processes.

Inflammation of the mammary gland following bacterial infection, commonly known as mastitis, affects all mammalian species. Although the aetiology and epidemiology of mastitis in the dairy cow are well described, the genetic factors mediating resistance to mammary gland infection are not well known, due in part to the difficulty in obtaining robust phenotypic information from sufficiently large numbers of individuals. To address this problem, an experimental mammary gland infection experiment was undertaken, using a Friesian-Jersey cross breed F2 herd. A total of 604 animals received an intramammary infusion of Streptococcus uberis in one gland, and the clinical response over 13 milkings was used for linkage mapping and genome-wide association analysis. A quantitative trait locus (QTL) was detected on bovine chromosome 11 for clinical mastitis status using micro-satellite and Affymetrix 10 K SNP markers, and then exome and genome sequence data used from the six F1 sires of the experimental animals to examine this region in more detail. A total of 485 sequence variants were typed in the QTL interval, and association mapping using these and an additional 37 986 genome-wide markers from the Illumina SNP50 bovine SNP panel revealed association with markers encompassing the interleukin-1 gene cluster locus. This study highlights a region on bovine chromosome 11, consistent with earlier studies, as conferring resistance to experimentally induced mammary gland infection, and newly prioritises the IL1 gene cluster for further analysis in genetic resistance to mastitis.

The genetic merit of a herd is a key determinant in productivity for dairy farmers. However, making breeding decisions to maximize the rate of genetic gain can be complex because there is no certainty about which cows will become pregnant with a heifer calf. In this study, breeding worth (BrW) was used as a measure of genetic merit, and several mating strategies were evaluated. These strategies included randomly mating whole herds to the entire bull team, excluding low-ranked cows from producing replacement heifers, and nominating high-ranked cows to the most highly ranked bulls. Simulations were undertaken using 4 bull teams generated from bulls currently marketed in New Zealand and a selection of New Zealand dairy herds. Average replacement heifer BrW was calculated for 1,000 iterations of each combination of mating strategy, herd, and bull team (scenario). Variation in resulting average replacement heifer BrW within scenarios was due to random sampling of which cows became pregnant with a heifer calf. Relative to mating the whole herd to an entire bull team, excluding the lowest ranked cows from producing replacements resulted in the greatest increase in average replacement heifer BrW across all herds and bull teams, with a gain of approximately 0.4 BrW point for each 1% of cows excluded. Nominating top-ranking cows to the highest ranking bulls in the team had little effect (0.06-0.13 BrW increase for each 1% of top cows nominated) in improving BrW of replacement heifers. The number of top bulls nominated had a variable effect depending on the BrW spread of the entire bull team. Although excluding cows with the lowest BrW from producing replacement heifers is most effective for improving BrW, it is important to ensure that the number of heifers born is sufficient to replace cows leaving the herd. It is likely that optimal strategies for improving BrW will vary from farm to farm depending not only on the BrW structure of the herd, the bull team available, and the reproduction success on farm but also on farm management practices. This simulation study provides expected outcomes from a variety of mating strategies to allow informed decision making on farm., (The Authors. Published by FASS Inc. and Elsevier Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).)

Background: Lactose provides an easily-digested energy source for neonates, and is the primary carbohydrate in milk in most species. Bovine lactose is also a key component of many human food products. However, compared to analyses of other milk components, the genetic control of lactose has been little studied. Here we present the first GWAS focussed on analysis of milk lactose traits., Results: Using a discovery population of 12,000 taurine dairy cattle, we detail 27 QTL for lactose concentration and yield, and subsequently validate the effects of 26 of these loci in a distinct population of 18,000 cows. We next present data implicating causative genes and variants for these QTL. Fine mapping of these regions using imputed, whole genome sequence-resolution genotypes reveals protein-coding candidate causative variants affecting the ABCG2, DGAT1, STAT5B, KCNH4, NPFFR2 and RNF214 genes. Eleven of the remaining QTL appear to be driven by regulatory effects, suggested by the presence of co-locating, co-segregating eQTL discovered using mammary RNA sequence data from a population of 357 lactating cows. Pathway analysis of genes representing all lactose-associated loci shows significant enrichment of genes located in the endoplasmic reticulum, with functions related to ion channel activity mediated through the LRRC8C, P2RX4, KCNJ2 and ANKH genes. A number of the validated QTL are also found to be associated with additional milk volume, fat and protein phenotypes., Conclusions: Overall, these findings highlight novel candidate genes and variants involved in milk lactose regulation, whose impacts on membrane transport mechanisms reinforce the key osmo-regulatory roles of lactose in milk.

X chromosome inactivation (XCI) is a process by which 1 of the 2 copies of the X chromosomes present in female mammals is inactivated. The transcriptional silencing of one X chromosome achieves dosage compensation between XX females and XY males and ensures equal expression of X-linked genes in both sexes. Although all mammals use this form of dosage compensation, the complex mechanisms that regulate XCI vary between species, tissues, and development. These mechanisms include not only varying levels of inactivation, but also the nature of inactivation, which can range from being random in nature to driven by parent of origin. To date, no data describing XCI in calves or adult cattle have been reported and we are reliant on data from mice to infer potential mechanisms and timings for this process. In the context of dairy cattle breeding and genomic prediction, the implications of X chromosome inheritance and XCI in the mammary gland are particularly important where a relatively small number of bulls pass their single X chromosome on to all of their daughters. We describe here the use of RNA-seq, whole genome sequencing and Illumina BovineHD BeadChip (Illumina, San Diego, CA) genotypes to assess XCI in lactating mammary glands of dairy cattle. At a population level, maternally and paternally inherited copies of the X chromosome are expressed equally in the lactating mammary gland consistent with random inactivation of the X chromosome. However, average expression of the paternal chromosome ranged from 10 to 90% depending on the individual animal. These results suggest that either the mammary gland arises from 1 or 2 stem cells, or a nongenetic mechanism that skews XCI exists. Although a considerable amount of future work is required to fully understand XCI in cattle, the data reported here represent an initial step in ensuring that X chromosome variation is captured and used in an appropriate manner for future genomic selection., (The Authors. Published by the Federation of Animal Science Societies and Elsevier Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).)

Single nucleotide polymorphisms have been the DNA variant of choice for genomic prediction, largely because of the ease of single nucleotide polymorphism genotype collection. In contrast, structural variants (SV), which include copy number variants (CNV), translocations, insertions, and inversions, have eluded easy detection and characterization, particularly in nonhuman species. However, evidence increasingly shows that SV not only contribute a substantial proportion of genetic variation but also have significant influence on phenotypes. Here we present the discovery of CNV in a prominent New Zealand dairy bull using long-read PacBio (Pacific Biosciences, Menlo Park, CA) sequencing technology and the Sniffles SV discovery tool (version 0.0.1; https://github.com/fritzsedlazeck/Sniffles). The CNV identified from long reads were compared with CNV discovered in the same bull from Illumina sequencing using CNVnator (read depth-based tool; Illumina Inc., San Diego, CA) as a means of validation. Subsequently, further validation was undertaken using whole-genome Illumina sequencing of 556 cattle representing the wider New Zealand dairy cattle population. Very limited overlap was observed in CNV discovered from the 2 sequencing platforms, in part because of the differences in size of CNV detected. Only a few CNV were therefore able to be validated using this approach. However, the ability to use CNVnator to genotype the 557 cattle for copy number across all regions identified as putative CNV allowed a genome-wide assessment of transmission level of copy number based on pedigree. The more highly transmissible a putative CNV region was observed to be, the more likely the distribution of copy number was multimodal across the 557 sequenced animals. Furthermore, visual assessment of highly transmissible CNV regions provided evidence supporting the presence of CNV across the sequenced animals. This transmission-based approach was able to confirm a subset of CNV that segregates in the New Zealand dairy cattle population. Genome-wide identification and validation of CNV is an important step toward their inclusion in genomic selection strategies., (The Authors. Published by the Federation of Animal Science Societies and Elsevier Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).)

Increasing environmental temperatures are a threat to the sustainability of livestock production and, because of the high metabolic demands of lactation, to dairy production in particular. Summer heat waves in temperate climates reduce feed intake, milk production, and cow comfort. In extreme heat events, there is an increase in cow mortality. In tropical climates, dairy cattle are mostly (zebu) type or zebu crossbred with temperate dairy breeds. Crossbreeding is undertaken to combine the heat tolerance and tick resistance of zebu with the productivity of temperate dairy breeds. In the absence of improved heat tolerance, milk production and fertility of temperate cattle is severely impaired. We have recently identified a key role for the prolactin pathway in regulating heat tolerance. A de novo mutation in prolactin that impairs prolactin activity was discovered in hairy and heat intolerant, New Zealand dairy cattle. The phenotypes produced were remarkably similar to those seen in fescue toxicosis, a syndrome seen in grazing cattle in the U.S. where ingestion of ergovaline, a fungal toxin from infected pasture, inhibits prolactin secretion. Recognition of the role of prolactin in hairy cattle led us to identify a deletion in exon 10 of the long-form of the prolactin receptor in Senepol cattle that causes truncation of the protein and determines the slick coat and heat tolerance traits found in this , beef breed. The short form of the prolactin receptor is predicted to be unaffected by the deletion. Knowledge of this dominant mutation has provided the impetus to begin a crossbreeding program to investigate performance and heat tolerance of temperate dairy cattle carrying the slick, prolactin receptor variant. The perceived opportunity is to introgress this variant into temperate dairy cattle to enable performance and welfare improvement in hot climates. Heat tolerance of cattle with slick coats appears to be mostly associated with coat type although sweating ability may also be enhanced. Further investigation is required of performance traits in cows homozygous for the slick variant because the published data are almost exclusively from heterozygous animals. Combination of the slick mutation with other favorable genes for heat tolerance, especially those for coat color, will be particularly enabled by gene editing technologies, offering opportunities for further improvement in bovine thermotolerance.

A major pleiotropic quantitative trait locus (QTL) located at ~25 Mbp on bovine chromosome 14 affects a myriad of growth and developmental traits in Bos taurus and indicus breeds. These QTL have been attributed to two functional variants in the bidirectional promoter of PLAG1 and CHCHD7. Although PLAG1 is a good candidate for mediating these effects, its role remains uncertain given that these variants are also associated with expression of five additional genes at the broader locus. In the current study, we conducted expression QTL (eQTL) mapping of this region using a large, high depth mammary RNAseq dataset representing 375 lactating cows. Here we show that of the seven previously implicated genes, only PLAG1 and LYN are differentially expressed by QTL genotype, and only PLAG1 bears the same association signature of the growth and body weight QTLs. For the first time, we also report significant association of PLAG1 genotype with milk production traits, including milk fat, volume, and protein yield. Collectively, these data strongly suggest PLAG1 as the causative gene underlying this diverse range of traits, and demonstrate new effects for the locus on lactation phenotypes.

The mammary gland is a prolific lipogenic organ, synthesising copious amounts of triglycerides for secretion into milk. The fat content of milk varies widely both between and within species, and recent independent genome-wide association studies have highlighted a milk fat percentage quantitative trait locus (QTL) of large effect on bovine chromosome 5. Although both EPS8 and MGST1 have been proposed to underlie these signals, the causative status of these genes has not been functionally confirmed. To investigate this QTL in detail, we report genome sequence-based imputation and association mapping in a population of 64,244 taurine cattle. This analysis reveals a cluster of 17 non-coding variants spanning MGST1 that are highly associated with milk fat percentage, and a range of other milk composition traits. Further, we exploit a high-depth mammary RNA sequence dataset to conduct expression QTL (eQTL) mapping in 375 lactating cows, revealing a strong MGST1 eQTL underpinning these effects. These data demonstrate the utility of DNA and RNA sequence-based association mapping, and implicate MGST1, a gene with no obvious mechanistic relationship to milk composition regulation, as causally involved in these processes.

The main objectives of this study were to establish the relative value of milk yields under twice-daily milking (TDM) as a predictor of yield and yield loss under once-daily milking (ODM), and to understand the role of residual milk and udder storage capacity-related traits in regulating yield and yield loss during ODM. A Holstein-Friesian × Jersey crossbred herd was established over 2 seasons (years), as 2 individual cohorts on the same farm, managed on a pasture-based system over 4 lactations. Short-term (1-wk) ODM studies, with a starting total of 690 cows, were undertaken in mid- and late-lactation in lactation 2 and in mid-lactation in lactation 3 for each cohort. A 10-wk study of ODM performance began in mid-lactation in lactation 3, whereas lactation 4 was a full-lactation assessment of ODM. In the short-term studies, milk yield under ODM was well predicted (R(2)=0.7 to 0.8 in 5 of 6 studies) by the daily yield under TDM in the week before ODM. Yield loss (kg/d) increased with increasing milk yield and with increasing somatic cell count (SCC), although predictions were relatively poor (R(2)=0.09 to 0.30). Yield loss (%) decreased with increasing TDM yield in 3 of the 6 studies and was positively correlated with SCC during ODM. Nevertheless, ODM yield loss, in absolute or percentage terms, was a poorly repeatable trait in grazing cows. Part of the variation in yield loss percentage (30%) was positively associated with residual milk (%), measured pretrial, during measurement of functional udder capacity in lactation 3. Total production (kg of milk) over the full-lactation ODM study in lactation 4 was correlated with total production in the 10-wk trial in lactation 3 (r=0.72 and 0.63 for cohorts 1 and 2, respectively). Identifying the highest- and lowest-producing 10% of animals during the full lactation of ODM indicated that poor production was associated with high yields of residual milk (measured in lactation 3) and, conversely, high production was associated with low yields of residual milk, relative to the other 80% of animals. These same "high" and "low" production groups from lactation 4 had similar differences in performance in the earlier short-term studies and a larger or smaller percentage yield loss associated with the residual milk measurement. Breeding strategies for ODM may benefit, therefore, from greater emphasis on selecting for a low residual milk fraction to optimize milking performance. Nevertheless, the level of milk production under TDM is a strong phenotypic predictor of milk production under ODM., (Copyright © 2015 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

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.

Lactation, hair development and homeothermy are characteristic evolutionary features that define mammals from other vertebrate species. Here we describe the discovery of two autosomal dominant mutations with antagonistic, pleiotropic effects on all three of these biological processes, mediated through the prolactin signalling pathway. Most conspicuously, mutations in prolactin (PRL) and its receptor (PRLR) have an impact on thermoregulation and hair morphology phenotypes, giving prominence to this pathway outside of its classical roles in lactation.

The objective of this study was to estimate heritability and crossbreeding parameters (breed and heterosis effects) of various fatty acid (FA) concentrations in milk fat of New Zealand dairy cattle. For this purpose, calibration equations to predict concentration of each of the most common FAs were derived with partial least squares (PLS) using mid-infrared (MIR) spectral data from milk samples (n=850) collected in the 2003-04 season from 348 second-parity crossbred cows during peak, mid and late lactation. The milk samples produced both, MIR spectral data and concentration of the most common FAs determined using gas chromatography (GC). The concordance correlation coefficients (CCC) between the concentration of a FA determined by GC and the PLS equation ranged from 0.63 to 0.94, suggesting that some prediction equations can be considered to have substantial predictive ability. The PLS calibration equations were then used to predict the concentration of each of the fatty acids in 26,769 milk samples from 7385 cows that were herd-tested during the 2007-08 season. Data were analysed using a single-trait repeatability animal model. Shorter chain FA (16:0 and below) were significantly higher (P<0.05) in Jersey cows, while longer chain, including unsaturated longer chain FA were higher in Holstein-Friesian cows. The estimates of heritabilities ranged from 0.17 to 0.41 suggesting that selective breeding could be used to ensure milk fat composition stays aligned to consumer, market and manufacturing needs.

A genomic prediction for residual feed intake (RFI) developed in growing dairy heifers (RFIgro) was used to predict and test breeding values for RFI in lactating cows (RFIlac) from an independent, industry population. A selection of 3,359 cows, in their third or fourth lactation during the study, of above average genetic merit for milk production, and identified as at least 15/16ths Holstein-Friesian breed, were selected for genotyping from commercial dairy herds. Genotyping was carried out using the bovine SNP50 BeadChip (Illumina Inc., San Diego, CA) on DNA extracted from ear-punch tissue. After quality control criteria were applied, genotypes were imputed to the 624,930 single nucleotide polymorphisms used in the growth study. Using these data, genomically estimated breeding values (GEBV) for RFIgro were calculated in the selected cow population based on a genomic prediction for RFIgro estimated in an independent group of growing heifers. Cows were ranked by GEBV and the top and bottom 310 identified for possible purchase. Purchased cows (n=214) were relocated to research facilities and intake and body weight (BW) measurements were undertaken in 99 "high" and 98 "low" RFIgro animals in 4 consecutive groups [beginning at d 61 ± 1.0 standard error (SE), 91 ± 0.5 SE, 145 ± 1.3 SE, and 191 ± 1.5 SE d in milk, respectively] to measure RFI during lactation (RFIlac). Each group of ~50 cows (~25 high and ~25 low RFIgro) was in a feed intake facility for 35 d, fed pasture-alfalfa cubes ad libitum, milked twice daily, and weighed every 2 to 3 d. Milk composition was determined 3 times weekly. Body weight change and BW at trial mid-point were estimated by regression of pre- and posttrial BW measurements. Residual feed intake in lactating cows was estimated from a linear model including BW, BW change, and milk component yield (as MJ/d); RFIlac differed consistently between the high and low selection classes, with the overall means for RFIlac being +0.32 and -0.31 kg of dry matter (DM) per day for the high and low classes, respectively. Further, we found evidence of sire differences for RFIlac, with one sire, in particular, being highly represented in the low RFIgro class, having a mean RFIlac of -0.83 kg of DM per day in 47 daughters. In conclusion, genomic prediction of RFIgro based on RFI measured during growth will discriminate for RFIlac in an independent group of lactating cows., (Copyright © 2014 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Combining data from research herds may be advantageous, especially for difficult or expensive-to-measure traits (such as dry matter intake). Cows in research herds are often genotyped using low-density single nucleotide polymorphism (SNP) panels. However, the precision of quantitative trait loci detection in genome-wide association studies and the accuracy of genomic selection may increase when the low-density genotypes are imputed to higher density. Genotype data were available from 10 research herds: 5 from Europe [Denmark, Germany, Ireland, the Netherlands, and the United Kingdom (UK)], 2 from Australasia (Australia and New Zealand), and 3 from North America (Canada and the United States). Heifers from the Australian and New Zealand research herds were already genotyped at high density (approximately 700,000 SNP). The remaining genotypes were imputed from around 50,000 SNP to 700,000 using 2 reference populations. Although it was not possible to use a combined reference population, which would probably result in the highest accuracies of imputation, differences arising from using 2 high-density reference populations on imputing 50,000-marker genotypes of 583 animals (from the UK) were quantified. The European genotypes (n=4,097) were imputed as 1 data set, using a reference population of 3,150 that included genotypes from 835 Australian and 1,053 New Zealand females, with the remainder being males. Imputation was undertaken using population-wide linkage disequilibrium with no family information exploited. The UK animals were also included in the North American data set (n=1,579) that was imputed to high density using a reference population of 2,018 bulls. After editing, 591,213 genotypes on 5,999 animals from 10 research herds remained. The correlation between imputed allele frequencies of the 2 imputed data sets was high (>0.98) and even stronger (>0.99) for the UK animals that were part of each imputation data set. For the UK genotypes, 2.2% were imputed differently in the 2 high-density reference data sets used. Only 0.025% of these were homozygous switches. The number of discordant SNP was lower for animals that had sires that were genotyped. Discordant imputed SNP genotypes were most common when a large difference existed in allele frequency between the 2 imputed genotype data sets. For SNP that had ≥ 20% discordant genotypes, the difference between imputed data sets of allele frequencies of the UK (imputed) genotypes was 0.07, whereas the difference in allele frequencies of the (reference) high-density genotypes was 0.30. In fact, regions existed across the genome where the frequency of discordant SNP was higher. For example, on chromosome 10 (centered on 520,948 bp), 52 SNP (out of a total of 103 SNP) had ≥ 20% discordant SNP. Four hundred and eight SNP had more than 20% discordant genotypes and were removed from the final set of imputed genotypes. We concluded that both discordance of imputed SNP genotypes and differences in allele frequencies, after imputation using different reference data sets, may be used to identify and remove poorly imputed SNP., (Copyright © 2014 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Residual feed intake (RFI), as a measure of feed conversion during growth, was estimated for around 2,000 growing Holstein-Friesian heifer calves aged 6 to 9 mo in New Zealand and Australia, and individuals from the most and least efficient deciles (low and high RFI phenotypes) were retained. These animals (78 New Zealand cows, 105 Australian cows) were reevaluated during their first lactation to determine if divergence for RFI observed during growth was maintained during lactation. Mean daily body weight (BW) gain during assessment as calves had been 0.86 and 1.15 kg for the respective countries, and the divergence in RFI between most and least efficient deciles for growth was 21% (1.39 and 1.42 kg of dry matter, for New Zealand and Australia, respectively). At the commencement of evaluation during lactation, the cows were aged 26 to 29 mo. All were fed alfalfa and grass cubes; it was the sole diet in New Zealand, whereas 6 kg of crushed wheat/d was also fed in Australia. Measurements of RFI during lactation occurred for 34 to 37 d with measurements of milk production (daily), milk composition (2 to 3 times per week), BW and BW change (1 to 3 times per week), as well as body condition score (BCS). Daily milk production averaged 13.8 kg for New Zealand cows and 20.0 kg in Australia. No statistically significant differences were observed between calf RFI decile groups for dry matter intake, milk production, BW change, or BCS; however a significant difference was noted between groups for lactating RFI. Residual feed intake was about 3% lower for lactating cows identified as most efficient as growing calves, and no negative effects on production were observed. These results support the hypothesis that calves divergent for RFI during growth are also divergent for RFI when lactating. The causes for this reduced divergence need to be investigated to ensure that genetic selection programs based on low RFI (better efficiency) are robust., (Copyright © 2014 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Milk is composed of a complex mixture of lipids, proteins, carbohydrates and various vitamins and minerals as a source of nutrition for young mammals. The composition of milk varies between individuals, with lipid composition in particular being highly heritable. Recent reports have highlighted a region of bovine chromosome 27 harbouring variants affecting milk fat percentage and fatty acid content. We aimed to further investigate this locus in two independent cattle populations, consisting of a Holstein-Friesian x Jersey crossbreed pedigree of 711 F2 cows, and a collection of 32,530 mixed ancestry Bos taurus cows. Bayesian genome-wide association mapping using markers imputed from the Illumina BovineHD chip revealed a large quantitative trait locus (QTL) for milk fat percentage on chromosome 27, present in both populations. We also investigated a range of other milk composition phenotypes, and report additional associations at this locus for fat yield, protein percentage and yield, lactose percentage and yield, milk volume, and the proportions of numerous milk fatty acids. We then used mammary RNA sequence data from 212 lactating cows to assess the transcript abundance of genes located in the milk fat percentage QTL interval. This analysis revealed a strong eQTL for AGPAT6, demonstrating that high milk fat percentage genotype is also additively associated with increased expression of the AGPAT6 gene. Finally, we used whole genome sequence data from six F1 sires to target a panel of novel AGPAT6 locus variants for genotyping in the F2 crossbreed population. Association analysis of 58 of these variants revealed highly significant association for polymorphisms mapping to the 5'UTR exons and intron 1 of AGPAT6. Taken together, these data suggest that variants affecting the expression of AGPAT6 are causally involved in differential milk fat synthesis, with pleiotropic consequences for a diverse range of other milk components.

Aim: To estimate genetic and crossbreeding parameters for the incidence of recorded clinical lameness in New Zealand dairy cattle., Methods: Herd records from 76,357 cows, collected during the 2005/06 to 2008/09 milking seasons from 155 herds in the Livestock Improvement Corporation young sire progeny test scheme, were used to estimate genetic parameters and breed effects for incidence of recorded clinical lameness in HolsteinFriesian, Jersey and crossbred dairy cattle. Recorded clinical lameness was coded "1" for cows that presented at least one event of clinical lameness at any day during the season and "0" for unaffected cows. Genetic parameters were estimated using an animal model across breeds considering all and then only first lactation records. Heritability and repeatability of recorded clinical lameness were calculated from the variance component estimates both with and without logit transformation., Results: The mean incidence of recorded clinical lameness per herd was 6.3 (min 2, max 34)%. The incidence of recorded clinical lameness in Holstein Friesian cows (mean 6.8, SE 0.24%) was higher than the incidence of recorded clinical lameness in crossbred (mean 6.1, SE 0.19%) and Jersey cows (mean 6.0, SE 0.28%) (p=0.0002). There was no difference in incidence between crossbred and Jersey cows (p=0.96). Estimates of the heritability of recorded clinical lameness as an untransformed trait were 0.053 (SE 0.014) for first lactation records and 0.016 (SE 0.003) for all lactation records. As a transformed (logit) trait heritabilities were 0.067 (SE 0.024) and 0.044 (SE 0.016) for first and all lactation records, respectively. The repeatability estimates of recorded clinical lameness were 0.071 (SE 0.005) and 0.107 (SE 0.011) for untransformed and logit transformed lactation records, respectively. Sire estimated breeding values for recorded clinical lameness showed the lowest values in Jersey sires, and ranged between -5 and 8%., Conclusions: Despite the low heritability of recorded clinical lameness, this study provided evidence that there is significant exploitable animal genetic variation. Selection of specific sires across and within breeds could be an option for increasing genetic resistance to lameness in New Zealand dairy cattle.

Feed makes up a large proportion of variable costs in dairying. For this reason, selection for traits associated with feed conversion efficiency should lead to greater profitability of dairying. Residual feed intake (RFI) is the difference between actual and predicted feed intakes and is a useful selection criterion for greater feed efficiency. However, measuring individual feed intakes on a large scale is prohibitively expensive. A panel of DNA markers explaining genetic variation in this trait would enable cost-effective genomic selection for this trait. With the aim of enabling genomic selection for RFI, we used data from almost 2,000 heifers measured for growth rate and feed intake in Australia (AU) and New Zealand (NZ) genotyped for 625,000 single nucleotide polymorphism (SNP) markers. Substantial variation in RFI and 250-d body weight (BW250) was demonstrated. Heritabilities of RFI and BW250 estimated using genomic relationships among the heifers were 0.22 and 0.28 in AU heifers and 0.38 and 0.44 in NZ heifers, respectively. Genomic breeding values for RFI and BW250 were derived using genomic BLUP and 2 bayesian methods (BayesA, BayesMulti). The accuracies of genomic breeding values for RFI were evaluated using cross-validation. When 624,930 SNP were used to derive the prediction equation, the accuracies averaged 0.37 and 0.31 for RFI in AU and NZ validation data sets, respectively, and 0.40 and 0.25 for BW250 in AU and NZ, respectively. The greatest advantage of using the full 624,930 SNP over a reduced panel of 36,673 SNP (the widely used BovineSNP50 array) was when the reference population included only animals from either the AU or the NZ experiment. Finally, the bayesian methods were also used for quantitative trait loci detection. On chromosome 14 at around 25 Mb, several SNP closest to PLAG1 (a gene believed to affect stature in humans and cattle) had an effect on BW250 in both AU and NZ populations. In addition, 8 SNP with large effects on RFI were located on chromosome 14 at around 35.7 Mb. These SNP may be associated with the gene NCOA2, which has a role in controlling energy metabolism., (Copyright © 2012 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Selection for divergence between individuals for efficiency of feed utilization (residual feed intake, RFI) has widespread application in the beef industry and is usually undertaken when animals are fed diets based on silages with grain. The objective of this research was to develop a feeding system (using Gallagher, Hamilton, New Zealand, electronics) to measure RFI for growth in Holstein-Friesian heifers (aged 5-9 mo), and identify divergent individuals to be tested for RFI during lactation. A dry forage diet (alfalfa cubes) was fed because intakes could be measured accurately, and the New Zealand dairy industry (4.4 million milking cows in lactation) relies heavily on forage feeding. The evaluation was undertaken over 3 yr with 1,052 animals fed in a facility for 7 wk, and weighed 3 times weekly. The mean age at the start of measurements was 215 d, body weight (BW) 189 kg, and mean daily dry matter intakes averaged 6.7 kg. Body weight gain (all animals) averaged 0.88 kg/d. The RFI was determined as the residuals from the regression of mean intake on mean BW(0.75) and daily BW gain of individuals. Actual and fitted intakes were strongly related (R(2) = 0.82). In terms of gross efficiency (feed intake/BW gain), RFI+year explained 43% of the variation, BW gain+year explained 66%, and RFI+BW gain+year explained 79% of the variation (all P<0.001). Daily BW gains (kg) of the most and least efficient 10% averaged (± standard deviation) 0.88 ± 0.15 and 0.88 ± 0.12 (P = 0.568), respectively, and the divergence between mean intakes was 1.46 kg of dry matter/d. The most and least efficient animals will be tested for RFI during lactation and genetic markers will be identified for the trait., (Copyright © 2012 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

We report mapping of a quantitative trait locus (QTL) with a major effect on bovine stature to a ∼780-kb interval using a Hidden Markov Model-based approach that simultaneously exploits linkage and linkage disequilibrium. We re-sequenced the interval in six sires with known QTL genotype and identified 13 clustered candidate quantitative trait nucleotides (QTNs) out of >9,572 discovered variants. We eliminated five candidate QTNs by studying the phenotypic effect of a recombinant haplotype identified in a breed diversity panel. We show that the QTL influences fetal expression of seven of the nine genes mapping to the ∼780-kb interval. We further show that two of the eight candidate QTNs, mapping to the PLAG1-CHCHD7 intergenic region, influence bidirectional promoter strength and affect binding of nuclear factors. By performing expression QTL analyses, we identified a splice site variant in CHCHD7 and exploited this naturally occurring null allele to exclude CHCHD7 as single causative gene.

Regulation of milk synthesis and secretion is controlled mostly through local (intramammary) mechanisms. To gain insight into the molecular pathways comprising this response, an analysis of mammary gene expression was conducted in 12 lactating cows shifted from twice daily to once daily milking. Tissues were sampled by biopsy from adjacent mammary quarters of these animals during the two milking frequencies, allowing changes in gene expression to be assessed within each animal. Using bovine-specific, oligonucleotide arrays representing 21,495 unique transcripts, a range of differentially expressed genes were found as a result of less frequent milk removal, constituting transcripts and pathways related to apoptotic signaling (NF-kappaB, JUN, ATF3, IGFBP5, TNFSF12A) mechanical stress and epithelial tight junction synthesis (CYR61, CTGF, THBS1, CLDN4, CLDN8), and downregulated milk synthesis (LALBA, B4GALT1, UGP2, CSN2, GPAM, LPL). Quantitative real-time PCR was used to assess the expression of 13 genes in the study, and all 13 of these were correlated (P < 0.05) with values derived from array analysis. It can be concluded that the physiological changes that occur in the bovine mammary gland as a result of reduced milk removal frequency likely comprise the earliest stages of the involution response and that mechano-signal transduction cascades associated with udder distension may play a role in triggering these events.

Aim: 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., Methods: 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)., Results: 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., Conclusions: 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.

beta-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 beta-carotene and subsequently identified a premature stop codon in bovine beta-carotene oxygenase 2 (BCO2), which also affects serum beta-carotene content. The BCO2 enzyme is thereby identified as a key regulator of beta-carotene metabolism.

The imprints of domestication and breed development on the genomes of livestock likely differ from those of companion animals. A deep draft sequence assembly of shotgun reads from a single Hereford female and comparative sequences sampled from six additional breeds were used to develop probes to interrogate 37,470 single-nucleotide polymorphisms (SNPs) in 497 cattle from 19 geographically and biologically diverse breeds. These data show that cattle have undergone a rapid recent decrease in effective population size from a very large ancestral population, possibly due to bottlenecks associated with domestication, selection, and breed formation. Domestication and artificial selection appear to have left detectable signatures of selection within the cattle genome, yet the current levels of diversity within breeds are at least as great as exists within humans.

New platforms utilizing single nucleotide polymorphisms (SNP) offer operational advantages over the conventional microsatellite-based ones, making them a promising alternative for parentage exclusion. Through simulation and empirical data, a 40-SNP panel (where the minor allele frequency was 0.35 on average) was shown to be a comparable or better diagnostic tool than the current 14-microsatellite panel that is used to parentage test New Zealand dairy animals. The 40 SNP alone did not have sufficient power of exclusion to match more than 75% of the progeny to the correct sire and dam. Utilizing mating records and grouping progeny and dams by birth and calving dates, respectively, decreased the number of sire-dam combinations that each progeny was tested against and dramatically increased the utility of the SNP. These results highlight the importance of combining genotypes with on-farm data to maximize the ability to assign parentage in the New Zealand dairy herd.

When a genetic marker and a quantitative trait locus (QTL) are in linkage disequilibrium (LD) in one population, they may not be in LD in another population or their LD phase may be reversed. The objectives of this study were to compare the extent of LD and the persistence of LD phase across multiple cattle populations. LD measures r and r(2) were calculated for syntenic marker pairs using genomewide single-nucleotide polymorphisms (SNP) that were genotyped in Dutch and Australian Holstein-Friesian (HF) bulls, Australian Angus cattle, and New Zealand Friesian and Jersey cows. Average r(2) was approximately 0.35, 0.25, 0.22, 0.14, and 0.06 at marker distances 10, 20, 40, 100, and 1000 kb, respectively, which indicates that genomic selection within cattle breeds with r(2) >or= 0.20 between adjacent markers would require approximately 50,000 SNPs. The correlation of r values between populations for the same marker pairs was close to 1 for pairs of very close markers (<10 kb) and decreased with increasing marker distance and the extent of divergence between the populations. To find markers that are in LD with QTL across diverged breeds, such as HF, Jersey, and Angus, would require approximately 300,000 markers.

The aim of this paper was to characterize the diversity among haplotypes based on 22 single nucleotide polymorphisms (SNPs) and one deletion within four casein genes in two geographically distant goat populations, the Sicilian Girgentana breed and the Norwegian goat breed. Forty Girgentana goats were genotyped for the aforementioned polymorphisms and the resulting data set was compared with 436 goats from the Norwegian population previously genotyped for these markers. Several casein gene polymorphisms were not in Hardy-Weinberg equilibrium either in Girgentana, or in the Norwegian breed. The SNP haplotype frequencies for the four casein genes were calculated and despite the large geographical distance and phenotypic divergence between these breeds, a proportion of casein loci haplotypes were found to be identical between both Norwegian and Girgentana goats. However, for the CSN2 gene there were no common haplotypes between the two populations. The level of linkage disequilibrium between the casein genes was less in the Girgentana population than in the Norwegian population.

In seasonal dairy production systems a high level of reproductive performance is required to maintain a 365-day calving interval. Compared with many countries New Zealand dairy cattle are fertile, with an average calving interval of 368 days and first service conception rates of around 55%. However, fertility has been demonstrated to have declined over the past 15-years. An extended interval from calving to first luteal activity (CLA) has previously been identified as a major cause of infertility in New Zealand dairy cows. This paper investigates the phenotypic relationship between CLA and reproductive performance in a population of second cross FriesianxJersey dairy cattle. The data consists of 1780 lactations, from 801 individuals that were managed in a seasonal, pasture-based dairy production system. Overall, reproductive performance in the herd was high, with first service conception rates of over 60% and more than 80% of the herd pregnant within 42-days of the start of the breeding period. Animals were divided into six groups of at least 270 animals per group based on CLA (50 days). First service conception rates and the interval from the start of mating to conception were not statistically different between the CLA groups (p>0.05), despite the number of oestrous cycles prior to first mating being inversely proportional to CLA. However, only 47% of animals with a CLA of >70 days (n=88) were pregnant within 21-day of the start of mating, compared with 65% of animals with CLA of <70 days. Animals that had not ovulated by 1-week prior to the start of mating also had reduced reproductive performance (45% versus 67% pregnant within 21-days of the start of mating, p<0.001). The poor reproductive performance of animals with CLA>70 days is likely to be due to a lower proportion of animals having ovulated by the start of mating (0.28 versus 0.95, p<0.001). Contrary to previous findings in non-seasonal production systems the length of CLA was not related to fertility, except when CLA was greater than 70 days.

A mechanobullous skin disorder was identified in the progeny of a 3-y-old Friesian-Jersey crossbred bull. The condition presented as loss of skin and mucosa from contact areas and inflammation. Examination of skin samples under light microscopy revealed separation of the epidermis from the dermis. Electron microscopic analysis refined the site of cleavage to above the basement membrane involving lysis of basal keratinocytes. These observations were consistent with the simplex form of epidermolysis bullosa (EB) in humans. Candidate genes based on human gene mutations were assessed, resulting in keratin 5 being identified as the most likely candidate gene. The sequence of bovine keratin 5 was established and sequencing led to identification of a G to A substitution in all affected animals. This mutation leads to an amino acid change of glutamic acid to lysine in the final E (478) of the KLLEGE motif of the protein. The sire carried a de novo mutation and was mosaic, explaining his asymptomatic status and the less than expected frequency of affected offspring. Remarkably, the same mutation has been previously described in EB simplex in humans.

The discovery of markers linked to genes that are responsible for traits of interest to the dairy industry might prove useful because they could aid in selection and breeding decisions. We have developed a selective DNA pooling methodology to allow us to efficiently screen the bovine genome in order to find genes responsible for production traits. Using markers on chromosome 14 as a test case, we identified a gene (DGAT1) previously known to affect three traits (fat yield, protein yield and total milk yield). Furthermore, we predicted similar effects to those previously shown for DGAT1 in a New Zealand Holstein-Friesian herd. Additionally, we showed a low error rate (1.6%) for the pooling procedure. Hence we are confident that we can apply this procedure to an entire genome scan in the search for quantitative trait loci (QTL).

Recently, DGAT1 was identified as the gene that underlies the QTL for bovine milk production on chromosome 14. This study investigated the effect of the reported polymorphism in three dairy breeds in New Zealand. Statistically significant results were identified for milk fat, milk protein, and volume for Jersey and Holstein-Friesian breeds, and only milk volume for Ayrshires. The average allele substitution effects were 2 to 3 kg of protein and 120 to 130 l milk for both the Jersey and Holstein-Friesian breeds. For milk fat, the average allele substitution effect was 6 kg for Holstein-Friesians and 3 kg for Jerseys. In all breeds, where the polymorphism increased milk fat yield, it decreased milk protein yield and milk volume.

The aim of this study was to evaluate the profitability of dairy herds under three mating systems involving the Holstein-Friesian, Jersey, and Ayrshire breeds. Mating systems were straight breeding and rotational cross-breeding using two or three breeds. A deterministic model was developed to simulate the nutritional, biological, and economic performance of dairy herds under New Zealand conditions. Expected performances per cow were obtained using estimates of breed group and heterosis effects, age effects, and age distribution in the herd. Requirements for dry matter in feed were estimated per cow for maintenance, lactation, pregnancy, and growth of the replacements. Stocking rate was calculated by assuming 12,000 kg of dry matter utilized annually per hectare. Productivity per hectare was calculated as performance per cow multiplied by stocking rate. Profitability was the difference between income (sale of milk and salvage value of animals) and costs (related to the number of cows in the herd and the land area farmed). Under current market values for milk and meat, all of the rotational crossbred herds showed superior profitability to the straightbred herds (Holstein-Friesian x Jersey, NZ$505/ha; Holstein-Friesian x Jersey x Ayrshire NZ$493/ha; Jersey x Ayrshire, NZ$466/ha; Holstein-Friesian x Ayrshire, NZ$430/ha; Jersey, NZ$430/ha; Holstein-Friesian, NZ$398/ha; and Ayrshire, NZ$338/ha). Changes in the value for fat relative to protein affected profitability more significantly in herds using the Jersey breed, and changes in the value for meat affected profitabiity more significantly in herds using the Holstein-Friesian and Ayrshire breeds. Results suggested that, under New Zealand conditions, the use of rotational crossbreeding systems could increase profitability of dairy herds under the conceivable market conditions.

The effects of selection and crossbreeding on the New Zealand dairy industry net income were evaluated with a deterministic model over a 25 yr planning horizon. Several mating strategies involving Holstein-Friesian, Jersey and Ayrshire cattle were evaluated. These strategies were straight breeding, upgrading to Holstein-Friesian, upgrading to Jersey, upgrading to Ayrshire, use of the best bulls irrespective of breed and two- and three-breed rotational crossbreeding. Industry productions of milk, fat, protein, and lactose were calculated assuming that 12,000 kg of dry matter per hectare was utilized from 1,224,911 hectares of pasture. Profitability was the difference between income (international sale of whole milk powder, casein, butter, and beef from salvage animals) and costs (farm expenses, milk collection, manufacture, and marketing). Casein and whole milk powder were valued at NZ$8.345 and NZ$3.306/ per kilogram, respectively, over 25 yr. Butter was valued at NZ$2.995/kilogram for base year production levels and NZ$0.45/kilogram for marginal increases in production. Upgrading to Holstein-Friesian resulted in the highest industry net income (NZ$1119 million) followed by straight breeding (NZ$1086 million) and two-breed rotational Holstein-Friesian x Jersey (NZ$1076 million). However, if the marginal value of extra butter production was assumed equal to the average base value, then upgrading to Jersey resulted in the highest industry net income (NZ$1185 million) followed by two-breed rotational Holstein-Friesian x Jersey (NZ$1177 million) and use of the best bulls (NZ$1173 million). Future costs and prices of dairy products have major impact on mating strategies.

Selective genotyping (i.e., genotyping of individuals that are extreme for a quantitative trait) has been undertaken in many studies but has usually been limited to one trait of interest. This paper outlines the implications of selective genotyping when applied to a daughter design. Formulas are presented that enable unbiased estimation of QTL effects for the selectively genotyped and correlated traits. Formulas were verified using simulation. Further, formulas are presented to enable power calculations to be undertaken for the selectively genotyped trait and correlated traits. These algorithms are demonstrated with a numerical example.

Seven New Zealand Holstein-Friesian families and two Jersey families, a total of 274 sires, were analyzed in a granddaughter design for marker-quantitative trait loci associations. For 17 nonproduction traits (management, size, and conformation traits), an across-family analysis was undertaken using multimarker regression procedures. Threshold levels were set empirically by permuting the data. A quantitative trait locus for stature was identified on chromosome 14, which was significant at the 15% experimentwise level (suggestive linkage). The quantitative trait locus was identified to be segregating in three of the nine families.