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.

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. 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. 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.

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.

The number of secretory cells in the mammary gland is often cited as a major determinant of milk production. However, literature data for proxy measures of secretory cell number do not fully support such a claim. In particular, measurements of total mammary DNA in livestock explain only <25% of variation in milk yield, probably because of tissue heterogeneity for DNA concentration. Relative to BW, measurements of udder size in dairy cattle, as total DNA or udder weight, are approximately double those seen in most other livestock classes. Therefore, selection for dairy production, not surprisingly, has resulted in cows with greater secretory capacity. There is limited evidence that genetic selection is still increasing udder size in some cattle populations, but more recent data are needed. It is contended that the most important period of mammary growth for determination of milk yield is that occurring in pregnancy and early lactation. Mammary development is largely complete, at term, in sheep, goats, and cattle, but in pigs, the udder continues to grow during the first 3 wk of lactation, depending, in part, on litter size. Increased litter size in sheep and goats will enhance the extent of mammary development at the end of gestation (and milk yield) by 20 to 25%. However, twinning in dairy cattle does not affect milk production and, by inference only, is not likely to affect numbers of secretory cells at term. Milking frequency and suckling intensity in very early lactation will increase milk yield in cows and increase mammogenesis and milk yield in sheep, indicating that even at term, the ruminant gland retains some capacity for further development, if demand requires it. There is limited understanding of the hormonal signals in pregnancy that regulate mammary development relative to the number of young carried. Furthermore, the genetic differences between dairy and beef cattle that lead to substantially greater udder size in the dairy breeds have not been identified. During lactation, the drivers for secretory cell loss in relation to milking frequency and nutritional status are still not known. Measurement of mammary development and using this phenotype in genomewide association studies to identify key genetic variants for mammogenesis will provide knowledge that is fundamental to understanding the quantitative regulation of milk production. Copyright © 2017. American Society of Animal Science .

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 Bos indicus (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 Bos taurus, 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. Copyright © 2017. American Society of Animal Science .

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.