Holstein Friesian dairy cattle edited for diluted coat color as a potential adaptation to climate change.

BACKGROUND: High-producing Holstein Friesian dairy cattle have a characteristic black and white coat, often with large proportions of black. Compared to a light coat color, black absorbs more solar radiation which is a contributing factor to heat stress in cattle. To better adapt dairy cattle to rapidly warming climates, we aimed to lighten their coat color by genome editing. RESULTS: Using gRNA/Cas9-mediated editing, we introduced a three bp deletion in the pre-melanosomal protein 17 gene (PMEL) proposed as causative variant for the semi-dominant color dilution phenotype observed in Galloway and Highland cattle. Calves generated from cells with homozygous edits revealed a strong color dilution effect. Instead of the characteristic black and white markings of control calves generated from unedited cells, the edited calves displayed a novel grey and white coat pattern. CONCLUSION: This, for the first time, verified the causative nature of the PMEL mutation for diluting the black coat color in cattle. Although only one of the calves was healthy at birth and later succumbed to a naval infection, the study showed the feasibility of generating such edited animals with the possibility to dissect the effects of the introgressed edit and other interfering allelic variants that might exist in individual cattle and accurately determine the impact of only the three bp change.

The evolving role of Fourier-transform mid-infrared spectroscopy in genetic improvement of dairy cattle.

Over the last 100 years, significant advances have been made in the characterisation of milk composition for dairy cattle improvement programs. Technological progress has enabled a shift from labour intensive, on-farm collection and processing of samples that assess yield and fat levels in milk, to large-scale processing of samples through centralised laboratories, with the scope extended to include quantification of other traits. Fourier-transform mid-infrared (FT-MIR) spectroscopy has had a significant role in the transformation of milk composition phenotyping, with spectral-based predictions of major milk components already being widely used in milk payment and animal evaluation systems globally. Increasingly, there is interest in analysing the individual FT-MIR wavenumbers, and in utilising the FT-MIR data to predict other novel traits of importance to breeding programs. This includes traits related to the nutritional value of milk, the processability of milk into products such as cheese, and traits relevant to animal health and the environment. The ability to successfully incorporate these traits into breeding programs is dependent on the heritability of the FT-MIR predicted traits, and the genetic correlations between the FT-MIR predicted and actual trait values. Linking FT-MIR predicted traits to the underlying mutations responsible for their variation can be difficult because the phenotypic expression of these traits are a function of a diverse range of molecular and biological mechanisms that can obscure their genetic basis. The individual FT-MIR wavenumbers give insights into the chemical composition of milk and provide an additional layer of granularity that may assist with establishing causal links between the genome and observed phenotypes. Additionally, there are other molecular phenotypes such as those related to the metabolome, chromatin accessibility, and RNA editing that could improve our understanding of the underlying biological systems controlling traits of interest. Here we review topics of importance to phenotyping and genetic applications of FT-MIR spectra datasets, and discuss opportunities for consolidating FT-MIR datasets with other genomic and molecular data sources to improve future dairy cattle breeding programs.

Strategies for noise reduction and standardization of milk mid-infrared spectra from dairy cattle.

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.

Detection and assessment of copy number variation using PacBio long-read and Illumina sequencing in New Zealand dairy cattle.

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.

Bovine mammary gland X chromosome inactivation.

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.

BREEDING AND GENETICS SYMPOSIUM:Breeding heat tolerant dairy cattle: the case for introgression of the "slick" prolactin receptor variant into dairy breeds.

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.

Leveraging genetically simple traits to identify small-effect variants for complex phenotypes.

BACKGROUND: Polymorphisms underlying complex traits often explain a small part (less than 1 %) of the phenotypic variance (σ2P). This makes identification of mutations underling complex traits difficult and usually only a subset of large-effect loci are identified. One approach to identify more loci is to increase sample size of experiments but here we propose an alternative. The aim of this paper is to use secondary phenotypes for genetically simple traits during the QTL discovery phase for complex traits. We demonstrate this approach in a dairy cattle data set where the complex traits were milk production phenotypes (fat, milk and protein yield; fat and protein percentage in milk) measured on thousands of individuals while secondary (potentially genetically simpler) traits are detailed milk composition traits (measurements of individual protein abundance, mineral and sugar concentrations; and gene expression). RESULTS: Quantitative trait loci (QTL) were identified using 11,527 Holstein cattle with milk production records and up to 444 cows with milk composition traits. There were eight regions that contained QTL for both milk production and a composition trait, including four novel regions. One region on BTAU1 affected both milk yield and phosphorous concentration in milk. The QTL interval included the gene SLC37A1, a phosphorous antiporter. The most significant imputed sequence variants in this region explained 0.001 σ2P for milk yield, and 0.11 σ2P for phosphorus concentration. Since the polymorphisms were non-coding, association mapping for SLC37A1 gene expression was performed using high depth mammary RNAseq data from a separate group of 371 lactating cows. This confirmed a strong eQTL for SLC37A1, with peak association at the same imputed sequence variants that were most significant for phosphorus concentration. Fitting any of these variants as covariables in the association analysis removed the QTL signal for milk production traits. Plausible causative mutations in the casein complex region were also identified using a similar strategy. CONCLUSIONS: Milk production traits in dairy cows are typical complex traits where polymorphisms explain only a small portion of the phenotypic variance. However, here we show that these mutations can have larger effects on secondary traits, such as concentrations of minerals, proteins and sugars in the milk, and expression levels of genes in mammary tissue. These larger effects were used to successfully map variants for milk production traits. Genetically simple traits also provide a direct biological link between possible causal mutations and the effect of these mutations on milk production.

Gene expression in liver and adipose tissue is altered during and after temporary changes to postpartum milking frequency.

Short-term changes to milking frequency can alter the metabolic status of dairy cows depending on the duration, magnitude, and stage of lactation at which the milking frequency changes occur. Additionally, effects of altered milking frequency that are subsequent to cows returning to a normal twice-daily (2×) milking regimen are not well established. This study tested the hypothesis that plasma concentrations of key hormones and metabolites and transcription of genes involved in the somatotropic axis and lipid metabolism would be altered in liver and subcutaneous adipose tissue from cows milked with different frequencies. Multiparous Holstein-Friesian dairy cows were allocated to 2× milking for the whole lactation, or once-(1×) or 3 times-(3×) daily milking for 3 or 6 wk, immediately postpartum, and then 2× milking for the remainder of the lactation. Liver and subcutaneous fat were biopsied at wk 1 (liver only), 3, 6, and 9 postpartum, and transcription of genes involved in the somatotropic axis and lipid metabolism were measured. At wk 3, cows milked 3× had lower hepatic expression of growth hormone receptor (GHR1A) compared with cows milked 2× or 1×, and lower IGF1 expression compared with cows milked 1×, indicating greater uncoupling of the somatotropic axis. At wk 6, reduced transcription of total GHR and GHR1B occurred in the adipose tissue of cows milked 3×. Cows milked 1× had greater transcription in adipose tissue of lipogenesis genes at wk 3 and 6, and lipolysis genes at wk 6, compared with cows milked 2×, indicating a period of increased fatty acid storage, followed by increased fatty acid reesterification. At wk 9, cows previously milked 3× for 6 wk maintained lower transcription of genes involved in lipogenesis, lipolysis, and ketolysis in adipose tissue compared with cows milked 2×, indicating that the effects of 3× milking persist for at least 3 wk after switching to 2× milking. Results indicate that alterations to milking frequency affect the transcription of genes involved in lipid mobilization and storage, enabling the animal to manage the energy demands associated with the change in milk production. Some of these gene transcription changes were maintained in cows previously milked 3×, indicating that the adipose tissue gene expression changes were still required even after 3 wk of the less-demanding 2× milking regimen.

The expression of genes involved in hepatic metabolism is altered by temporary changes to milking frequency.

Changes to milking frequency (MF) affect the metabolic and energetic status of dairy cows. However, the duration of altered MF necessary to modify hepatic transcription during early lactation is less clear. Additionally, long-term responses to short-term alterations in MF have not been established. Holstein-Friesian dairy cows (n = 120) were allocated to 3 or 6 wk of either once-daily (1 ×) or thrice-daily (3 ×) milking, immediately postpartum. Following treatment, cows were switched to twice-daily (2 ×) milking. These 4 treatment groups were compared with cows milked 2 × (n = 30) for the whole lactation. Liver tissue was collected by biopsy at 1, 3, 6, and 9 wk postpartum from 12 cows per treatment, RNA was extracted, and transcript abundance of genes involved in hepatic metabolism was quantified. Milking frequency altered the expression of most of the genes measured; however, we observed no effects caused by the length of time on the alternative milking frequency and no interactions between MF and length. During the MF treatment, mRNA expression of some, but not all, genes involved in gluconeogenesis (G6PC, PCK1), fatty acid ß-oxidation (CPT1A, CPT2), ketogenesis (HMGCS2), lipid transport (APOA1), and lipolysis (PNPLA2) were lower for cows milked 1 × and plasma glucose and insulin concentrations were greater. Cows milked 3 × had reduced mRNA expression for some of the genes involved in fatty acid synthesis (ACACA) and lipid transport (APOB) and had greater plasma NEFA concentrations at wk 1. At 9 wk postpartum, expression data indicated that cows previously milked 3 × had a greater capacity for gluconeogenesis (PCK1), ketogenesis (HMGCS2), and urea cycling (ASL, CPS1) and lower glucose concentrations than cows previously milked 1 ×, because some of the genes involved in these processes were still altered. Milking cows 1 × relative to 2 ×, however, did not result in significant carryover effects on the expression of the genes measured in this study, indicating that metabolic changes are not sustained beyond the period of reduced MF. Changes to MF altered the hepatic response during early lactation; however, this was not dependent on the duration of MF change. Although we observed only minimal carryover effects on hepatic metabolism from short periods of reduced MF postpartum, there may be long-term effects on urea cycling (ASL, CPS1) and ketogenesis (HMGCS2) when 3 × milking occurs immediately postpartum.

Expression analysis of key somatotropic axis and liporegulatory genes in ghrelin- and obestatin-infused dairy cows.

Ghrelin, an orexigenic hormone, is the endogenous ligand for the growth hormone secretagogue receptor (GHSR). Obestatin is produced from the same precursor peptide as ghrelin, and although obestatin was initially thought to promote actions opposite to those of ghrelin, many studies have failed to confirm this hypothesis. In the current study, multiparous cows were continuously infused with ghrelin (n = 10) or obestatin (n = 10) for 8 wk and compared to an untreated group (n = 10) to examine the effects of these hormones on somatotropic and liporegulatory gene expression. The expression of key genes was measured by quantitative real-time polymerase chain reaction. Growth hormone secretagogue receptor mRNA expression was altered in ghrelin- and obestatin-infused cows in a similar manner, as expression was increased at 4 wk, however it had decreased by 8 wk. Obestatin-infused cows presented with a significant decrease in the expression of ATP-binding cassette A1 (ABCA1) in adipose tissue, suggesting changes in cholesterol transport. Liver insulin-like growth factor (IGF) binding protein-3 mRNA displayed a week-by-treatment interaction, as expression was increased in control and obestatin-infused cows; however, expression decreased in ghrelin-infused cows. Adipose expression of hormone sensitive lipase (LIPE) mRNA was not altered by treatment or time, suggesting hormone infusion is not initiating lipolysis. The expression of lipogenic genes in adipose tissue increased with time in all groups, consistent with the general lactational profile of lipogenesis in dairy cows. These data indicate that continuous infusion of ghrelin or obestatin does not alter the expression of key somatotropic or liporegulatory genes in the lactating dairy cow, although obestatin infusion may alter cholesterol transport.

Effects of reduced frequency of milk removal on gene expression in the bovine mammary gland.

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.