Genome-wide association study of age at puberty and its (co)variances with fertility and stature in growing and lactating Holstein-Friesian dairy cattle.

Reproductive performance is a key determinant of cow longevity in a pasture-based, seasonal dairy system. Unfortunately, direct fertility phenotypes such as intercalving interval or pregnancy rate tend to have low heritabilities and occur relatively late in an animal's life. In contrast, age at puberty (AGEP) is a moderately heritable, early-in-life trait that may be estimated using an animal's age at first measured elevation in blood plasma progesterone (AGEP4) concentrations. Understanding the genetic architecture of AGEP4 in addition to genetic relationships between AGEP4 and fertility traits in lactating cows is important, as is its relationship with body size in the growing animal. Thus, the objectives of this research were 3-fold. First, to estimate the genetic and phenotypic (co)variances between AGEP4 and subsequent fertility during first and second lactations. Second, to quantify the associations between AGEP4 and height, length, and BW measured when animals were approximately 11 mo old (standard deviation = 0.5). Third, to identify genomic regions that are likely to be associated with variation in AGEP4. We measured AGEP4, height, length, and BW in approximately 5,000 Holstein-Friesian or Holstein-Friesian × Jersey crossbred yearling heifers across 54 pasture-based herds managed in seasonal calving farm systems. We also obtained calving rate (CR42, success or failure to calve within the first 42 d of the seasonal calving period), breeding rate (PB21, success or failure to be presented for breeding within the first 21 d of the seasonal breeding period) and pregnancy rate (PR42, success or failure to become pregnant within the first 42 d of the seasonal breeding period) phenotypes from their first and second lactations. The animals were genotyped using the Weatherby's Versa 50K SNP array (Illumina, San Diego, CA). The estimated heritabilities of AGEP4, height, length, and BW were 0.34 (90% credibility interval [CRI]: 0.30, 0.37), 0.28 (90% CRI: 0.25, 0.31), 0.21 (90% CRI: 0.18, 0.23), and 0.33 (90% CRI: 0.30, 0.36), respectively. In contrast, the heritabilities of CR42, PB21 and PR42 were all <0.05 in both first and second lactations. The genetic correlations between AGEP4 and these fertility traits were generally moderate, ranging from 0.11 to 0.60, whereas genetic correlations between AGEP4 and yearling body-conformation traits ranged from 0.02 to 0.28. Our GWAS highlighted a genomic window on chromosome 5 that was strongly associated with variation in AGEP4. We also identified 4 regions, located on chromosomes 14, 6, 1, and 11 (in order of decreasing importance), that exhibited suggestive associations with AGEP4. Our results show that AGEP4 is a reasonable predictor of estimated breeding values for fertility traits in lactating cows. Although the GWAS provided insights into genetic mechanisms underpinning AGEP4, further work is required to test genomic predictions of fertility that use this information.

Genome-wide association study of anogenital distance and its (co)variances with fertility in growing and lactating Holstein-Friesian dairy cattle.

Anogenital distance (AGD) is a moderately heritable trait that can be measured at a young age that may provide an opportunity to indirectly select for improved fertility in dairy cattle. In this study, we characterized AGD and its genetic and phenotypic relationships with a range of body stature and fertility traits. We measured AGD, shoulder height, body length, and body weight in a population of 5,010 Holstein-Friesian and Holstein-Friesian × Jersey crossbred heifers at approximately 11 mo of age (AGD1). These animals were born in 2018 across 54 seasonal calving, pasture-based dairy herds. A second measure of AGD was collected in a subset of herds (n = 17; 1,956 animals) when the animals averaged 29 mo of age (AGD2). Fertility measures included age at puberty (AGEP), then time of calving, breeding, and pregnancy during the first and second lactations. We constructed binary traits reflecting the animal's ability to calve during the first 42 d of their herd's seasonal calving period (CR42), be presented for breeding during the first 21 d of the seasonal breeding period (PB21) and become pregnant during the first 42 d of the seasonal breeding period (PR42). The posterior mean of sampled heritabilities for AGD1 was 0.23, with 90% of samples falling within a credibility interval (90% CRI) of 0.20 to 0.26, whereas the heritability of AGD2 was 0.29 (90% CRI 0.24 to 0.34). The relationship between AGD1 and AGD2 was highly positive, with a genetic correlation of 0.89 (90% CRI 0.82 to 0.94). Using a GWAS analysis of 2,460 genomic windows based on 50k genotype data, we detected a region on chromosome 20 that was highly associated with variation in AGD1, and a second region on chromosome 13 that was moderately associated with variation in AGD1. We did not detect any genomic regions associated with AGD2 which was measured in fewer animals. The genetic correlation between AGD1 and AGEP was 0.10 (90% CRI 0.00 to 0.19), whereas the genetic correlation between AGD2 and AGEP was 0.30 (90% CRI 0.15 to 0.44). The timing of calving, breeding, and pregnancy (CR42, PB21, and PR42) during first or second lactations exhibited moderate genetic relationships with AGD1 (0.19 to 0.52) and AGD2 (0.46 to 0.63). Genetic correlations between AGD and body stature traits were weak (≤0.16). We conclude that AGD is a moderately heritable trait, which may have value as an early-in-life genetic predictor for reproductive success during lactation.

Animal- and herd-level factors associated with onset of puberty in grazing dairy heifers.

AIMS: To explore animal- and herd-level risk factors influencing age at puberty in predominantly Holstein-Friesian dairy heifers managed in seasonal, pasture-based systems. METHODS: Heifers born in spring 2018 (n = 5,010) from 54 commercial dairy herds in New Zealand were visited on three occasions when the mean heifer age, within herd, was 10 (visit 1; V1), 11 (V2) and 12 (V3) months old. Blood samples were collected on each visit and liveweight, stature and anogenital distance (AGD) were measured at V2. Heifers were defined as having reached puberty at the first visit where blood progesterone was elevated (≥ 1 ng/mL). Animal-level response variables included pubertal status by V1, V2 and V3, and age at puberty (or age at V3 plus 31 days for those that had not attained puberty by V3). To explore herd-level management factors, farmers answered a questionnaire relating to animal location, land type, health, feeding, and management between weaning and mating. A partial least squares regression was undertaken to identify herd-level factors associated with the greatest influence on puberty rate within herd. RESULTS: The mean age at puberty was 352 (SD 34.9) days. Heavier animals at a greater proportion of expected mature liveweight based on their breeding value for liveweight, or animals with a higher breed proportion of Jersey and lower breed proportion of Holstein, were associated with earlier puberty. Herd puberty rates varied widely among enrolled herds, and averaged 20%, 39% and 56% by V1, V2 and V3, respectively. Liveweight, followed by breed and land type, had the greatest influence on the herd puberty rate. Heifer herds with a greater mean liveweight (absolute and proportion of expected mature weight) or greater Jersey proportion had more animals that reached puberty at any visit, whereas herds located on steep land or with greater Holstein breed proportions had lower puberty rates. Management-related factors such as vaccinations, provision of feed supplements, and weighing frequency were also herd-level risk factors of puberty but had less influence. CONCLUSIONS AND CLINICAL RELEVANCE: This study highlights the importance of having well-grown heifers for increasing the chances of earlier puberty onset and the effect of breed and youngstock management to achieve growth targets. These outcomes have important implications for the optimal management of heifers to achieve puberty before their maiden breeding and for the timing of measurements to potentially incorporate a puberty trait in genetic evaluations.

Effect of nonsteroidal anti-inflammatory drugs on the inflammatory response of bovine endometrial epithelial cells in vitro.

Chronic postpartum uterine infection detrimentally affects subsequent fertility. Nonsteroidal anti-inflammatory drugs (NSAID) are used to alleviate pain and treat inflammatory conditions in transition dairy cows with varying success. To screen the efficacy of NSAID in the absence of animal experiments, we have established an in vitro model to study uterine inflammation. Inflammation was induced in cultured bovine endometrial epithelial cells by challenging cells with an inflammation cocktail: lipopolysaccharide and proinflammatory cytokines, interleukin-1ß (IL1ß) and tumor necrosis factor α (TNFα). Release of the inflammation markers, serum amyloid A (SAA) and α-1-acid glycoprotein (αAGP), was measured by ELISA. Concentration of these markers was used to indicate the effectiveness in dampening inflammation of 5 NSAID: meloxicam, flunixin meglumine, aspirin, ketoprofen, and tolfenamic acid. Three NSAID, meloxicam, flunixin meglumine, and tolfenamic acid, were successful at dampening the release of SAA and αAGP into cell-culture supernatant, and the corresponding treated cells were selected for down-stream mRNA expression analysis. Expression of 192 genes involved in regulation of inflammatory pathways were investigated using Nanostring. Of the genes investigated, 81 were above the mRNA expression-analysis threshold criteria and were included in expression analysis. All SAA genes investigated (SAA2, SAA3, M-SAA3.2) were upregulated in response to the inflammation cocktail, relative to mRNA expression in control cells; however, AGP mRNA expression was below the expression analysis threshold and was, therefore, excluded from analysis. Treatment with NSAID downregulated genes involved in regulating chemokine signaling (e.g., CXCL2, CXCR4, CXCL5, and CXCL16) and genes that regulate the eicosanoid pathway (e.g., LTA4H, PTGS2, PLA2G4A, and PTGDS). Of the 5 NSAID investigated, meloxicam, flunixin meglumine, and tolfenamic acid are recommended for further investigation into treatment of postpartum uterine inflammation. The results from this study confirm the immunomodulatory properties of the endometrial epithelium in response to inflammatory stimuli and suggest that NSAID may be beneficial in alleviating uterine inflammation.

Variance parameter estimation for age at puberty phenotypes under 2 levels of phenotype censorship.

Age at puberty (AGEP) is a moderately heritable trait in cattle that may be predictive of an animal's genetic merit for reproductive success later in life. In addition, under some mating strategies (for example, where mating begins before all animals have attained puberty) animals that attain puberty at a relatively young age will also likely conceive earlier than their herd mates, and thus begin their productive life earlier. Unfortunately, AGEP is challenging to measure because animals must be observed over a period of several months. Our objectives for this study were twofold. The first objective was to produce variance components for AGEP. The second objective was to investigate the implications of a simplified phenotyping strategy for AGEP, when the interval between repeated blood plasma progesterone measures was extended from weekly to monthly, increasing the extent of left, interval, and right censoring. We measured AGEP in a closely monitored population of around 500 Holstein-Friesian heifers, born in 2015 and managed under a seasonal, pasture-based dairy system. Animals were blood tested weekly from approximately 240 to 440 d of age and were deemed to have reached puberty when blood plasma progesterone elevation (>1 ng/mL) was detected in 2 of 3 consecutive blood tests (AGEP_Weekly). To simulate a simplified phenotyping strategy based on monthly herd visits (AGEP_Monthly), we selectively disregarded data from all but 3 blood test events, when animals were around 300, 330, and 360 d of age (standard deviation = 14.5 d). The posterior mean of estimated heritabilities for AGEP_Weekly was 0.54, with a 90% credibility interval (90% CRI) of 0.41 to 0.66, whereas it was 0.44 (90% CRI 0.32 to 0.57) for AGEP_Monthly. The correlation between EBVs for AGEP_Weekly and AGEP_Monthly was 0.87 (90% CRI, 0.84 to 0.89). We conclude that in this population, AGEP is a moderately heritable trait. Further, increasing phenotype censorship from weekly to monthly observations would not have altered the main conclusions of this analysis. Our results support the strategic use of censoring to reduce costs and animal ethics considerations associated with collection of puberty phenotypes.

The first week following insemination is the period of major pregnancy failure in pasture-grazed dairy cows.

A 60% pregnancy success for inseminations is targeted to optimize production efficiency for dairy cows within a seasonal, pasture-grazed system. Routine measures of pregnancy success are widely available but are limited, in practice, to a gestation stage beyond the first 28 d. Although some historical data exist on embryonic mortality before this stage, productivity of dairy systems and genetics of the cows have advanced significantly in recent decades. Accordingly, the aim was to construct an updated estimate of pregnancy success at key developmental stages during the first 70 d after insemination. Blood samples were collected for progesterone concentrations on d 0 and 7. A temporal series of 4 groups spanning fertilization through d 70 were conducted on 4 seasonal, pasture-grazed dairy farms (n = 1,467 cows) during the first 21 d of the seasonal breeding period. Morphological examination was undertaken on embryos collected on d 7 (group E7) and 15 (group E15), and pregnancy was diagnosed via ultrasonography on approximately d 28 and 35 (group E35) as well as d 70 (group E70). Fertilization, embryo, and fetal evaluation for viability established a pregnancy success pattern. Additionally, cow and on-farm risk factor variables associated with pregnancy success were evaluated. We estimated pregnancy success rates of 70.9%, 59.1%, 63.8%, 62.3%, and 56.7% at d 7, 15, 28, 35, and 70, respectively. Fertilization failure (15.8%) and embryonic arrest before the morula stage (10.3%) were the major developmental events contributing to first-week pregnancy failures. Embryo elongation failure of 7% contributed to pregnancy failure during the second week. The risk factors for pregnancy success that were related to the cows included interval between calving and insemination, and d-7 plasma progesterone concentrations, whereas insemination sire was associated with pregnancy outcome. Most pregnancy failure occurs during the first week among seasonal-calving pasture-grazed dairy cows.

Adaptive immune response ranking is associated with reproductive phenotypes in grazing dairy cows divergent in genetic merit for fertility traits.

A strong adaptive immune response has been reported to have positive effects on fertility; therefore, we investigated antibody- and cell-mediated adaptive immune responses (AMIR and CMIR, respectively) and their associations with reproductive phenotypes using a population of animals that differed in their estimated genetic merit for fertility traits (fertility breeding value; FertBV). Holstein-Friesian heifers (n = 528) grazed on pasture in 4 herds based on age. These herds included 277 heifers of positive (POS) FertBV and 251 of negative (NEG) FertBV. The adaptive immune response (IR) was evaluated before puberty at 7.5 mo of age and used to rank animals as high, average, or low for AMIR, CMIR, and overall IR (combined CMIR and AMIR). The animals were studied from 12 wk of age through to the end of their second lactation to measure growth, puberty, and timing and success of fertility phenotypes, including those related to ovulation and pregnancy. Initial analysis indicated no difference in fertility outcomes between cows ranked as high or average for AMIR (n = 55, n = 407, respectively), CMIR (n = 87, n = 354, respectively), and IR (n = 29, n = 470, respectively), so these groups were pooled as HiAv-IR. Proportions of heifers of POS FertBV were similar within HiAv and low categories across AMIR (0.52 and 0.58, respectively), CMIR (0.51 and 0.59, respectively), and IR (0.53 and 0.48, respectively). Heifers with HiAv-IR had a greater average daily weight gain from 13 to 52 wk of age (661 g, 95% confidence interval 652, 669 vs. 619 g, 95% confidence interval 591, 647) and tended to be younger at puberty (371 d, 95% confidence interval 366, 377 vs. 385 d, 95% confidence interval 369, 401) than low-IR heifers. Low-CMIR cows of a NEG FertBV had a >40 d longer calving to first ovulation interval during their first lactation compared with HiAv-CMIR NEG FertBV cows. Low-CMIR cows also had decreased pregnancy rates at both 3 wk (25% ± 7% vs. 42% ± 3%; least squares means ± standard error) and 6 wk (33% ± 7% vs. 54% ± 3%; least squares means ± standard error) into the seasonal breeding period during their first lactation, compared with HiAv-CMIR cows. In summary, although the number of POS and NEG FertBV cows was similar in each immune group; interaction effects between FertBV and immune ranking on reproductive phenotypes are evident when cows were ranked by the overall IR. There were also associations between dairy cows' CMIR ranking and ability to return to estrus and become pregnant early in the breeding period, which will have benefits in seasonal breeding systems.

Relationships between prostaglandin concentrations, a single nucleotide polymorphism in HSD17B12, and reproductive performance in dairy cows.

Prostaglandins are involved in multiple processes important for fertility, with previous work in mice highlighting a potential role for the HSD17B12 gene in prostaglandin biosynthesis. This study aimed to determine the associations among circulating prostaglandin concentrations, a missense SNP in the HSD17B12 gene predicted to disrupt protein function, and fertility traits in first-lactation Holstein-Friesian dairy cows. We used a study population of approximately 500 animals specifically bred to have either a positive (POS, +5%) or negative (NEG, -5%) genetic merit for fertility (FertBV). Genotypes of a previously identified SNP (rs109711583) in HSD17B12 were determined, with 116 animals genotyped as AA, 215 genotyped as AG, and 153 genotyped as GG. Plasma concentrations of prostaglandin E2 and the PGF2α metabolite PGFM were determined at 3 time points (12 mo of age, 4 d postpartum, and 5 wk postpartum during first lactation) in a selection of animals with AA and GG genotypes from both the POS and NEG FertBV groups (n = 33-40 in each genotype for each FertBV group). Binary reproductive traits (yes or no) examined included submission for artificial breeding in the first 3 or 6 wk of the seasonal breeding period; conception to first service; conception during the first 6 wk of the breeding period; and pregnant at the end of the breeding period. Uterine health at 6 wk after calving was examined by evaluating the percentage of polymorphonuclear leukocytes following uterine cytology and by scoring vaginal discharge based on the presence of purulent material. The 3-wk submission rate was increased in animals that carried the G allele of the missense SNP in HSD17B12, but no differences were present among genotypes for 6-wk submission rate. The trait was additive, with each increase of the G allele increasing the 3-wk submission rate by 6 to 7%. We did not observe any consistent associations between SNP alleles and circulating PGE2 concentrations; however, a complex 3-way interaction among time, fertility group, and SNP allele was present for PGFM concentrations. Plasma concentrations of PGE2 were increased approximately 40% at 5 wk postpartum in animals that were submitted for breeding within 3 or 6 wk of the start of the breeding season, and in those that conceived during the first 6 wk of breeding, compared with those that did not. Plasma concentrations of PGFM were decreased approximately 20% in those animals that conceived to their first service and tended to be decreased in animals that were pregnant at the end of the breeding period, compared with those that were not. In summary, associations were observed between the SNP in HSD17B12 and submission rate by d 21 of the breeding season, as well as between circulating prostaglandin concentrations and fertility traits, but the SNP was not consistently linked to changes in prostaglandin concentrations. Thus, the association between submission rate by d 21 of the breeding season and the SNP in HSD17B12 were unlikely driven by changes in prostaglandins.

Response to kisspeptin and gonadotropin-releasing hormone agonist administration in Holstein-Friesian dairy heifers with positive or negative genetic merit for fertility traits.

Previous research has identified that Holstein-Friesian dairy heifers with positive (POS) genetic merit for fertility traits (FertBV) reach puberty earlier than heifers with negative (NEG) FertBV. The hypothalamus-pituitary-gonadal (HPG) axis is functional in heifers before the onset of puberty, with increased LH release evident as heifers progress toward puberty. We investigated the functionality of the HPG axis in peripubertal Holstein-Friesian dairy heifers with divergent POS or NEG FertBV, hypothesizing that the earlier puberty onset of POS heifers is associated with earlier activation of the HPG axis than in NEG heifers. In experiment 1, we tested the dose responsiveness of POS heifers to an intravenous injection of either kisspeptin [Kiss; 2, 4, or 8 µg/kg of body weight (BW); n = 3 per dose] or a GnRH agonist (buserelin; 5, 10, or 20 ng/kg of BW; n = 3 per dose). The use of these 2 agonists investigates the status of the HPG axis in both the hypothalamus (Kiss) and pituitary (buserelin) glands. Doses of 4 µg/kg BW of Kiss and 10 ng/kg BW of buserelin produced submaximal LH responses and were used in experiment 2, in which previously unused POS (n = 22) and NEG (n = 18) FertBV heifers were challenged with both agonists at 10 and 12 mo of age in a partial crossover design. Heifers were randomly allocated to treatment groups, balanced for age and BW. The LH response to buserelin was greater in POS heifers than NEG heifers at 10 mo of age, with no difference in response at 12 mo. The FSH response to buserelin and the LH and FSH responses to Kiss did not differ between the POS and NEG heifers at either age. These results indicate an association between divergent genetic merit for fertility and the LH release to buserelin at 10 mo of age, supporting the hypothesis that gonadotropin responsiveness to a GnRH agonist is more advanced in POS heifers than in NEG heifers.

Temporal profiles describing markers of inflammation and metabolism during the transition period of pasture-based, seasonal-calving dairy cows.

The physiology of the dairy cow while transitioning from pregnancy to lactation is complex, with multifactorial processes studied extensively for the role they play in manifestation of disease along with associated economic losses and compromised animal welfare. Manuscripts outlining associations among nutrition, production, physiology, and genetics variables and transition cow disorders are common in literature, with blood analytes that are central to energy metabolism (e.g., nonesterified fatty acids; NEFA, ß-hydroxybutyrate; BHB) often reported. Immunity and inflammation have increasingly been explored in the pathogenesis and persistence of disorders, with cytokines and acute phase proteins well documented. However, most of these studies have involved cows fed total mixed rations, which may not always reflect profiles of blood analytes and other physiological indicators of transition cow health in grazing cows consuming fresh pasture. Considering the comparatively lesser characterization of these analytes and markers in pasture-based, seasonal-calving dairy cows, we compiled a database consisting of 2,610 cow lactations that span 20 yr of transition cow research in New Zealand. Using this database, analyte profiles from approximately 28 d precalving to 35 d postcalving were identified in dairy cows with a range of genetics, milk production potentials, and pasture-based farm management systems. These profiles characterize changes in energy reserves and metabolism (NEFA, BHB, glucose, insulin, growth hormone, insulin-like growth factor-1, leptin, body condition score, body weight), liver function (globulin, aspartate aminotransferase, glutamate dehydrogenase, gamma-glutamyl transpeptidase, bilirubin, cholesterol, liver triacylglycerides), protein metabolism (albumin, total protein, albumin:globulin ratio, creatinine, urea, creatine kinase), mineral balance (calcium, magnesium, phosphate, potassium, sodium, chloride, bicarbonate), inflammation (IL-1ß, IL-6, haptoglobin, reactive oxygen species, total antioxidant capacity), and uterine health (polymorphonuclear cells, macrophage cells, vaginal discharge score). Temporal changes are generally consistent with previously characterized homeorhetic changes experienced by the dairy cow during the transition from pregnancy to lactation in both pastoral and housed systems. Some of the profiles had not previously been presented for pastoral systems, or in some cases, presented for either system. Our results indicate that moderate-yielding dairy cows undergo similar homeorhetic changes to high-yielding housed cows; however, differences in diet composition result in greater BHB concentrations than expected, based on their milk production and NEFA concentrations. In addition, most cows were able to transition to a state of higher energy requirement following calving, albeit with an increased metabolic challenge in the liver, and only a small percentage of cows were classified with severe hepatic lipidosis or severe hyperketonemia. Increases in metabolic function of the liver were accompanied by changes in indicators of the immune system and changes in mineral balance that, combined, probably reflect the innate response to the transition from gestation to lactation.

Development of an intrauterine infection model in the postpartum dairy cow.

AIMS: To develop an intrauterine infection model for Trueperella pyogenes in postpartum dairy cows and to assess the effect of this infection on the degree of intrauterine inflammation and concentrations of progesterone in serum. METHODS: The oestrous cycles of 36 healthy, non-pregnant, postpartum dairy cows were synchronised. They were then treated by intrauterine infusion of 0.5 g cephapirin before being blocked by age and randomly assigned to treatment with intrauterine infusion of saline (n = 18), 107 (n = 9) or 109 (n = 9) cfu of T. pyogenes, approximately 4 days after the expected time of ovulation (Day 0). Prior to intrauterine infusion on Day 0 and again on Days 3, 7, 10, and 15, cytobrush samples were collected from the uterus of each cow for microbiology and assessment of the percentage of polymorphonuclear neutrophils (PMN%). Blood samples were collected on the same days for measurement of concentrations of progesterone in serum, and uterine lumen diameter was assessed daily using transrectal ultrasonography. RESULTS: Trueperella pyogenes was isolated from 5/18 (28%), 7/9 (78%) and 8/9 (89%) cows infused with saline, 107 or 109 cfu of T. pyogenes, respectively (p < 0.001). Mean PMN% in the control cows did not change over time (p > 0.05), whereas it was higher on Days 7 and 10 than Day 0 in the 107 cfu group, and higher on Days 3 and 10 than Day 0 in the 109 cfu group (p < 0.05). The percentage of observations with uterine lumen diameters >2 mm was higher in cows infused with 107 (29.3 (95% CI = 14.5-44.2)%) or 109 cfu (19.2 (95% CI = 7.0-31.5)%) than in control cows (3.1 (95% CI = 0.1-6.0)%) (p < 0.001). Mean concentrations of progesterone in serum were higher in cows infused with 107 cfu (2.01 (SE 0.19) ng/mL) than cows infused with 109 cfu (1.01 (SE 0.27) ng/mL), with the control group intermediate (1.41 (SE 0.19) ng/mL) (p = 0.03). CONCLUSIONS: Infusion of 107 or 109 cfu of T. pyogenes resulted in the establishment of intrauterine infection in 83% of cows. Infection resulted in increased uterine lumen diameter, and an inflammatory response, i.e. elevated PMN%. CLINICAL RELEVANCE: This intrauterine infection model may be useful for future research on, for example, the pathogenesis of intrauterine infection in postpartum dairy cows.

The microenvironment of ovarian follicles in fertile dairy cows is associated with high oocyte quality.

We hypothesised that heifers and cows with positive genetic merit for fertility would have a follicular microenvironment that resulted in better quality oocytes. To test this, we compared cumulus cell-oocyte complexes (COC) and follicular fluid from preovulatory follicles of 36 Holstein-Friesian nulliparous heifers and 50 primiparous lactating cows with either positive (POS, +5%) or negative (NEG, -5%) fertility breeding values (FertBV). Established gene markers of oocyte quality were measured in individual cumulus cell masses and oocytes, and concentrations of amino acids, steroids, and metabolites were quantified in corresponding follicular fluid and plasma. The timing of visually detectable oestrus in NEG FertBV heifers was inconsistent with their stage of COC maturation. Retrospective analyses of oestrous activity data indicated that NEG FertBV heifers were sampled earlier. Their recovered COC were morphologically less mature and exhibited differential expression of genes that are associated with follicular maturation (lower levels of BMPR2) and protein processing (higher levels of HSP90B1). Despite consistent sampling times being achieved in the lactating cows, lower concentrations of serine, proline, methionine, isoleucine, and non-esterified fatty acids were present in follicular fluid from POS FertBV cows. This was associated with higher expression of gene biomarkers of good oocyte quality (VCAN, PDE8A) in COC recovered from POS FertBV cows. This study supports our hypothesis that the follicular microenvironment in lactating dairy cows with high genetic merit leads to COC with higher metabolic rates and oocytes of superior quality. Moreover, an additional stressor such as lactation is required for this difference to be pronounced.

Changes in plasma electrolytes, minerals, and hepatic markers of health across the transition period in dairy cows divergent in genetic merit for fertility traits and postpartum anovulatory intervals.

Peripartum metabolism and subsequent reproductive performance of dairy cows are linked, with maladaptation over the transition period associated with poor reproductive success. A herd of seasonal calving, grazing dairy cows was established that differed in their genetic merit for fertility traits. The heifers were produced by a customized mating program to achieve a 10-percentage point divergence in the New Zealand fertility breeding value (FertBV) as follows: +5 FertBV (POS) versus -5% FertBV (NEG), while also limiting divergence in other breeding values, including body weight, body condition score, and milk production. In this study, we aimed to characterize differences in metabolic, mineral, and metabolic stress marker profiles during their first postpartum transition period as primiparous heifers and to examine if animals with longer postpartum anestrous intervals (PPAI; more than 66 d compared with less than 35 d) had greater metabolic dysfunction. Blood was sampled at -21, -14, -7, 0, 4, 7, 10, 14, 17, 21, 28, and 35 d relative to calving in 455 primiparous cows and plasma analyzed. The NEG cows had lower concentrations of both plasma nonesterified fatty acids and ß-hydroxybutyrate at d 7 compared with POS cows. Detailed temporal profiling of various metabolic, mineral, and metabolic stress markers was undertaken in a subset of cows (n = 70). Cows were selected retrospectively to create 4 groups in a 2 × 2 factorial design with either a POS or NEG FertBV and either a short (19-35 d) or long (66-131 d) PPAI. The NEG cows tended, on average, to have lower nonesterified fatty acids and ß-hydroxybutyrate concentrations compared with POS cows across the transition period. Mean body weight and body condition score was greatest in NEG cows when compared with the POS cows and an interaction with day demonstrated this only occurred precalving. They also had indications of improved liver health precalving, with higher albumin-to-globulin ratios and lower bilirubin concentrations. Concentrations of aspartate aminotransferase were lower, and the Na-to-Cl ratio was greater in cows with a long versus a short PPAI at d 28 and d 35 after calving, potentially because of cows with a short PPAI (19-35 d) returning to estrous during this time. Magnesium concentrations were lower in NEG cows with a short PPAI from d 21 onwards, indicating NEG cows may metabolically respond to estrous differently than POS cows. The NEG-long PPAI cows had greater gamma-glutamyl transferase concentrations from calving until d 28 and lower bilirubin concentrations throughout the transition period. Together, the results demonstrate significant effects of FertBV on peripartum metabolic status. However, most of the markers tested returned to reference intervals within 4 d after calving or remained within those intervals for the whole transition period, indicating relatively minor biological effects of FertBV on transition period adaptation. The profound differences in reproductive performance among the groups was not explained by underlying differences in metabolic responses during the transition period.

Estrous activity in lactating cows with divergent genetic merit for fertility traits.

This observational study aimed to determine the effect of genetic merit for fertility traits on estrous expression and estrous cycle duration in grazing dairy cows, as measured by an activity monitoring device. A secondary aim was to describe changes in expression of estrus that occur during successive estrous cycles postpartum. Neck-mounted, activity-monitoring devices (Heatime, SCR Engineers Ltd.) were fitted to nulliparous Holstein-Friesian heifers with positive (POS FertBV) or negative genetic merit for fertility traits (NEG FertBV) to capture activity data during their first and second lactations (POS FertBV: n = 242, n = 188; NEG FertBV: n = 159, n = 87 in lactation 1 and 2, respectively). An estrous event was identified when the activity change index exceeded 26 activity units (AU) for 4 h. A total of 1,254 and 892 estrous events were identified in lactation 1 and 2, respectively. Estrous duration was defined as the interval between when the threshold was first exceeded and when activity dropped below the threshold, with no new event starting within 24 h of the end of the previous event. This definition of estrus included cows in which activity crossed the threshold multiple times in a day and were classified as a single estrous event. A second measure, high activity duration, was defined as the total hours that activity exceeded the threshold. To characterize estrous activity, peak activity (above baseline) and total activity (area under the curve of activity above baseline) were measured. Compared with NEG FertBV cows, POS FertBV cows had more active, longer estrous events. In lactation 1, the POS FertBV group had a mean estrous duration and a high activity duration of 12.5 and 12.4 h compared with 11.4 and 11.3 h for the NEG FertBV group [standard error of the difference (SED) = 0.5 and 0.4 h, respectively]. This significant difference also occurred in lactation 2, with a mean estrous duration of 13.1 versus 11.8 h (SED = 0.5 h) and a high activity duration of 13.0 versus 11.8 h (SED = 0.4 h) in the POS and NEG FertBV groups, respectively. Total activity and peak activity were greater in the POS compared with the NEG FertBV group in lactation 1 (peak activity: 65.5 vs. 55.8 AU, SED = 2.4 AU; total activity: 588 vs. 494 AU, SED = 25 AU) and lactation 2 (peak activity: 72.5 vs. 61.2 AU, SED = 2.9 AU; total activity: 648 vs. 541 AU, SED = 30 AU). Estrous cycle duration did not differ between the POS and NEG FertBV groups (lactation 1: 20.4 vs. 20.6 d, SED = 0.25; lactation 2: 20.8 vs. 21.0 d, SED = 0.28). Less estrous activity of the cow was associated with the first postpartum estrus. In contrast, the number of previous estrous events did not consistently affect the duration of the subsequent estrous cycle. The outcomes of this study provide evidence that positive genetic merit for fertility traits is associated with more overt estrous expression. Selection for these traits may improve estrous expression and thus estrous detection in commercial herds.

Investigating anogenital distance and antral follicle count as novel markers of fertility within a herd of cows with positive or negative genetic merit for fertility traits.

Using early-in-life markers of reproductive characteristics may enhance the speed and success of genetic improvement in fertility. We investigated 2 phenotypes that can be measured early in life and are moderately heritable to determine their association with traditional measures of reproductive success and genetic variation under a seasonal-calving, pasture-based system. Cows were bred to be divergent in the New Zealand Fertility Breeding Value, which estimates genetic merit for fertility. Cows consisted of 2 groups with an average positive (+5%) or negative (-5%) genetic merit for fertility traits and were expected to have large diversity in reproductive outcomes. Calves were genotyped at 41 ± 3.1 d of age (mean ± SD; n = 538), and antral follicle counts (AFC) were measured when they were postpubertal heifers before their first breeding (416 ± 15 d old; 92-d range; n = 520). The anogenital distance (AGD) was measured in 478 primiparous cows of this same population 50 to 60 d after the breeding start date when they were 881 ± 25 d old (145-d range). The AGD was shorter in animals with a positive genetic merit for fertility traits (based on parent averages). An indicator of herd reproductive success in a seasonal-calving system (recalving by 6 wk in lactation 2) was chosen for logistic regression with cross-validation, and if significant, a cut-off was calculated that categorized animals into groups. Both linear and quadratic regression was undertaken, and the model with the greatest sensitivity for detection of nonpregnant cows used. The AGD linear model was significant with a sensitivity of 64% and a specificity of 48%. This model resulted in a cut-off of 102 mm, which we used to classify cows as short (≤102 mm) or long (>102 mm) AGD animals. Primiparous cows with a short AGD were more likely to be pregnant within the first 3 and 6 wk of mating, and become pregnant as a primiparous cow, than those with a long AGD. The time from calving to conception was 20 d earlier in short AGD compared with long-AGD cows. None of the models tested for AFC were significant; therefore cows were categorized into 3 groups based on previous work in seasonal systems. However, associations between fertility phenotypes and AFC group were limited. Genomic regions of interest for AGD and AFC did not overlap, indicating phenotypes were genetically independent. Overall, AGD appears as a promising early marker of fertility in seasonal grazing systems.

Positive genetic merit for fertility traits is associated with superior reproductive performance in pasture-based dairy cows with seasonal calving.

New Zealand's Fertility Breeding Value (FertBV) is reported as the percentage of a sire's daughters that calve in the first 42 d of the seasonal calving period and is an estimate of genetic merit for fertility for dairy cattle. Reproductive physiology, milk production, and changes in body weight and body condition score of 2 groups of cows divergent in FertBV (+5.0%: POS; -5.1%: NEG) were characterized during their first 2 lactations. Cows grazed fresh pasture and were managed in a seasonal calving system; they were bred by artificial insemination on observed estrus for the entire breeding period (98 d in lactation 1 and 76 d in lactation 2). During lactation 1, all animals were primiparous and were randomly allocated to 1 of 2 herds, ensuring each herd was balanced for FertBV and expected calving date. During lactation 2, cows that became pregnant during lactation 1 were managed as 1 herd. Cows not inseminated in the first 42 d of the breeding season were examined for the presence of a corpus luteum and treated with an anestrus program. On average, the interval from calving to ovulation was 19 d longer in lactation 1 and 10 d longer in lactation 2 for NEG FertBV cows. The percent of cows submitted for artificial insemination after 21 d (i.e., submission rate) was 38 and 25 percentage points greater in the POS FertBV cows during lactations 1 and 2, respectively. Pregnancy rate from 42 d of breeding was 33 and 30 percentage points greater, respectively. There was no effect of FertBV on vaginal discharge score postcalving; however, POS FertBV cows had a 50% lower risk of having subclinical endometritis (polymorphonuclear leukocytes >7%) 42 d postcalving. Interactions between FertBV and month relative to calving identified that NEG FertBV cows were fatter (greater body condition score) in the month before calving, but thinner between 3 and 5 mo postcalving. There was no effect of FertBV on lactation length, estimated 270-d milk yields, or daily milk, fat, or protein yields, and only small effects on milk fat and protein percentage across the lactations. In summary, the POS FertBV cows had superior uterine health, a shorter calving to ovulation interval, a greater submission rate, and a greater pregnancy rate earlier in the breeding season when compared with the NEG FertBV cohort. Based on these results, these may be useful phenotypes to include in genetic selection indices.

Heifers with positive genetic merit for fertility traits reach puberty earlier and have a greater pregnancy rate than heifers with negative genetic merit for fertility traits.

This study investigated the hypothesis that dairy heifers divergent in genetic merit for fertility traits differ in the age of puberty and reproductive performance. New Zealand's fertility breeding value (FertBV) is the proportion of a sire's daughters expected to calve in the first 42 d of the seasonal calving period. We used the New Zealand national dairy database to identify and select Holstein-Friesian dams with either positive (POS, +5 FertBV, n = 1,334) or negative FertBV (NEG, -5% FertBV, n = 1,662) for insemination with semen from POS or NEG FertBV sires, respectively. The resulting POS and NEG heifers were predicted to have a difference in average FertBV of 10 percentage points. We enrolled 640 heifer calves (POS, n = 324; NEG, n = 316) at 9 d ± 5.4 d (± standard deviation; SD) for the POS calves and 8 d ± 4.4 d old for the NEG calves. Of these, 275 POS and 248 NEG heifers were DNA parent verified and retained for further study. The average FertBV was +5.0% (SD = 0.74) and -5.1% (SD = 1.36) for POS and NEG groups, respectively. Heifers were reared at 2 successive facilities as follows: (1) calf rearing (enrollment to ∼13 wk of age) and (2) grazier, after 13 wk until 22 mo of age. All heifers wore a collar with an activity sensor to monitor estrus events starting at 8 mo of age, and we collected weekly blood samples when individual heifers reached 190 kg of body weight (BW) to measure plasma progesterone concentrations. Puberty was characterized by plasma progesterone concentrations >1 ng/mL in at least 2 of 3 successive weeks. Date of puberty was defined when the first of these samples was >1 ng/mL. Heifers were seasonally bred for 98 d starting at ∼14 mo of age. Transrectal ultrasound was used to confirm pregnancy and combined with activity data to estimate breeding and pregnancy dates. We measured BW every 2 wk, and body condition and stature at 6, 9, 12, and 15 mo of age. The significant FertBV by day interaction for BW was such that the NEG heifers had increasingly greater BW with age. This difference was mirrored with the significant FertBV by month interaction for average daily gain, with the NEG heifers having a greater average daily gain between 9 and 18 mo of age. There was no difference in heifer stature between the POS and NEG heifers. The POS heifers were younger and lighter at puberty, and were at a lesser mature BW, compared with the NEG heifers. As a result, 94 ± 1.6% of the POS and 82 ± 3.2% of the NEG heifers had reached puberty at the start of breeding. The POS heifers were 20% and 11% more likely to be pregnant after 21 d and 42 d of breeding than NEG heifers (relative risk = 1.20, 95% confidence interval of 1.03-1.34; relative risk = 1.11, 95% confidence interval of 1.01-1.16). Results from this experiment support an association between extremes in genetic merit for fertility base on cow traits and heifer reproduction. Our results indicate that heifer puberty and pregnancy rates are affected by genetic merit for fertility traits, and these may be useful phenotypes for genetic selection.

Prevalence of endometritis diagnosed by vaginal discharge scoring or uterine cytology in dairy cows and herds.

Endometritis, diagnosed either by assessing the proportion of nucleated cells that are neutrophils (PMN%) following cytology of the endometrium or by assessing the degree of purulent material within the vagina (purulent vaginal discharge or PVD score), is prevalent among dairy cows. However, limited data exist as to the degree of variation among herds in the prevalence of endometritis diagnosed by these 2 methods. Thus, we undertook a study involving uterine cytological and vaginal sampling at a median of 41 d in milk of 1,807 cows from 100 seasonally breeding dairy herds in New Zealand. The optimal cut-point for PMN% was determined by receiver operator characteristic (ROC) curve analysis using conception to first artificial insemination (AI) as the outcome variable. The prevalence of disease was then calculated at the cow and herd levels, and an estimate of the effect of clustering of cow within a herd was calculated. Reproductive outcomes were collated and associations between endometritis and reproductive outcomes assessed using multivariable models. The optimal PMN% cut-point was ≥2%. The correlation of results for cows within a herd (the intraclass correlation) was 0.03, which was significant; hence, subsequent modeling accounted for this clustering. The cow-level prevalence of PMN% ≥2% was 27.0% [95% confidence interval (CI): 25.0 to 29.1%], whereas the mean within-herd prevalence of PMN% ≥2% was 27.1% (95% CI: 24.7 to 29.6%; range: 5.0 to 63.6%), and the prevalence among herds varied significantly. An elevated PMN% (≥2%) was significantly associated with a reduction in the proportion of cows conceiving to first AI (45.8 vs. 54.5%), a reduced proportion of cows submitted for AI in the first 3 wk of the seasonal breeding program (83.7 vs. 89.3%), and a lower proportion pregnant in the first 3 (44.4 vs. 55.4%) and 6 wk (67.5 vs. 76.4%) of the breeding program relative to cows with a low PMN% (i.e., <2%). A total of 24.6% of cows had a PVD score ≥2. The herd mean prevalence of PVD score ≥2 was 25.1% (95% CI: 22.5 to 27.7%; range: 5.0 to 65.0%) and varied significantly among herds. The level of agreement (kappa) between the PVD score and PMN% was low (16.8%) and nonsignificant. The effects of PVD score and PMN% on reproductive outcomes were independent. The within-herd median prevalence of endometritis based on combining both diagnostic tests and using a Bayesian latent class model was 22.9% (Bayesian 95% CI: 10.4 to 40.1%). We conclude that more than one-fifth of dairy cows have endometritis diagnosed either by PMN% or PVD in seasonal breeding herds when assessed at an average of 41 DIM, which was, on average, 30 d before the start of the seasonal breeding program. There is large and unexplained variation in prevalence of endometritis among herds. The 2 diagnostic methods were both associated with reproductive outcomes but have low levels of agreement between them and their effects appear to be independent.

Feeding synthetic zeolite to transition dairy cows alters neutrophil gene expression.

Synthetic zeolites are used to control the availability of dietary minerals (e.g., Ca, Mg, and P) in dairy cows. Due to calcium demand increasing with lactation onset, most cows become hypocalcemic immediately postpartum, which likely contributes to poorer immune function because calcium is important for immune cell signaling. To overcome postpartum hypocalcemia, we fed transition cows synthetic zeolite A (sodium aluminosilicate) precalving and hypothesized that it would alter calcium and thus neutrophil function during the transition period. Multiparous Holstein-Friesian cows in late gestation were randomly allocated to an untreated control group (n = 10) or a treatment group in which each cow received 500 g of zeolite A daily (n = 10) for 14 d prior to the expected calving date (actual duration = 17 ± 3 d prepartum). The cows grazed pasture, and each was supplemented with 2 kg/d of maize silage (dry matter basis), with or without zeolite, until calving. Blood samples for neutrophil isolation and analysis of plasma indicators of mineral status, energy status, liver function, and inflammation were collected pretreatment (covariate; d -19); on d -14 and -7 precalving; on the day of calving (d 0); and on d 1, 4, 7, and 28 postcalving. Neutrophils were isolated and gene expression was analyzed using microfluidic gene expression arrays. Neutrophil respiratory burst was assessed using stimulation with phorbol 12-myristate 13-acetate and flow cytometry. Plasma calcium and phosphorus revealed a treatment by time interaction; cows offered zeolite had greater plasma calcium concentrations at d 0, 1, and 4 postcalving and plasma phosphorus concentrations were lower in zeolite-treated cows during the precalving period until d 1 postcalving compared with control animals. Zeolite treatment downregulated neutrophil gene expression of CXCR4 and S100A8 and tended to lower gene expression for other immune mediators (CXCR1, IFNG, S100A12, and S100A9) compared with the control. Zeolite treatment did not affect neutrophil respiratory burst or expression of the other genes investigated. Plasma concentrations of cytokine IL-6 were reduced with zeolite treatment, which was most evident immediately postcalving (d 0, 1, and 7). Overall, feeding zeolite precalving had few effects on neutrophil gene expression and function; however, the lower gene expression of neutrophil inflammatory mediators may be due to altered availability of dietary minerals prepartum and indicates that zeolite A may control inflammation during the transition period.

Genotype by environment interactions in fertility traits in New Zealand dairy cows.

New Zealand's seasonal dairy farming system entails a condensed calving pattern with cows required to conceive within approximately 12 wk of the planned start of calving. This has resulted in strong selection for fertility through culling of nonpregnant cows and relatively strong emphasis on fertility in Breeding Worth, the national breeding objective that drives sire selection. Despite this, average herd-level fertility is highly variable across New Zealand dairy farms. We studied genotype by environment interaction in fertility-related traits, with the goal of improving selection decisions in different fertility environments. We used data from the New Zealand national dairy database, which contains records on 3,743,862 animals. Herds were classified into high-, mid-, or low-fertility categories or environments based on herd average fertility performance, and data were analyzed in 2 different ways. First, we estimated genetic parameters when the fertility trait was defined specifically for each fertility environment to determine the extent to which genetic correlations between high- and low-fertility environments differed from 1 and the extent of changes in genetic variance across environments. Second, we used simple regression to evaluate the impact of ancestral genetic merit for fertility on cow fertility phenotypes to compare the effect of changes in genetic merit on phenotypic performance between fertility environments. The genetic standard deviations of fertility-related traits were 1.5 to 3.6 times higher in low-fertility herds than in high-fertility herds, and the genetic correlations between the same fertility-related traits between the high- and low-fertility environments were moderate to high, albeit with high standard errors. The high standard errors of the correlations reflected the low heritabilities of the traits and potential problems of culling bias, particularly for traits expressed in later parities. Regression analysis revealed that the bottom 30% of herds (in terms of fertility) could achieve more than twice the benefit from selection for fertility than the top 30% of herds. Although our analyses do not support separate genetic evaluations of fertility in the different environments, they indicate that low-fertility herds could benefit more from targeted selection of sires with higher fertility estimated breeding values than from selection based solely on the multitrait national index. Conversely, high-fertility herds could focus their sire selection on traits other than fertility, provided they avoid very low fertility sires.