Pregnancy status predicted using milk mid-infrared spectra from dairy cattle.

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.

Genetic correlations of milk fatty acid contents predicted from milk mid-infrared spectra in New Zealand dairy cattle.

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.

Short communication: Identification of the pseudoautosomal region in the Hereford bovine reference genome assembly ARS-UCD1.2.

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.

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.

Mating strategies to maximize genetic merit in dairy cattle herds.

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.

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.

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.

Lactation traits associated with short- and long-term once-daily milking performance in New Zealand crossbred dairy cattle.

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.

Residual feed intake of lactating Holstein-Friesian cows predicted from high-density genotypes and phenotyping of growing heifers.

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.

Holstein-Friesian calves selected for divergence in residual feed intake during growth exhibited significant but reduced residual feed intake divergence in their first lactation.

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.

Imputation of genotypes from low density (50,000 markers) to high density (700,000 markers) of cows from research herds in Europe, North America, and Australasia using 2 reference populations.

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.

Measuring residual feed intake in dairy heifers fed an alfalfa (Medicago sativa) cube diet.

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.

Accuracy of genomic predictions of residual feed intake and 250-day body weight in growing heifers using 625,000 single nucleotide polymorphism markers.

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.

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.

Mutation in bovine beta-carotene oxygenase 2 affects milk color.

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.

Linkage disequilibrium and persistence of phase in Holstein-Friesian, Jersey and Angus cattle.

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 number of single nucleotide polymorphisms and on-farm data required for whole-herd parentage testing in dairy cattle herds.

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.

Phenotypic relationships between the interval from calving to first luteal activity and fertility in a seasonal dairy production system.

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.

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.

Quantitative trait loci analysis for five milk production traits on chromosome six in the Dutch Holstein-Friesian population.

Twenty Dutch Holstein-Friesian families, with a total of 715 sires, were evaluated in a granddaughter experiment design for marker-QTL associations. Five traits-milk, fat and protein yield and fat and protein percent-were analyzed. Across-family analysis was undertaken using multimarker regression principles. One and two QTL models were fitted. Critical values for the test statistic were calculated empirically by permuting the data. Individual trait distributions of permuted test statistics differed and, thus distributions, had to be calculated for each trait. Experimentwise critical values, which account for evaluating marker-QTL associations on all 29 autosomal-bovine chromosomes and for five traits, were calculated. A QTL for protein percent was identified in one, and two QTL models and was significant at the 1 and 2% level, respectively. Extending the multimarker regression approach to an analysis including two QTL was limited by families not being informative at all markers, which resulted in singularity. Below average heterozygosity for the first and last marker lowered information content for the first and last marker bracket. Highly informative markers at the ends of the mapped chromosome would overcome the decrease in information content in the first and last marker bracket and singularity for the two QTL model.