2.1 Experiment 1
2.1.1Herd management and reproductive activities
In this experiment, a total of 232 cows with genetic composition ¾ Holstein × ¼ Gir, from two farms were used. The animals were on the third lactation and with a body condition score (BCS) between 3 and 3.5 (Edmonson et al., 1989). The first farm had 90 lactating cows with a total of 2400 L/day (27 kg/cow/day). The second farm had 215 lactating animals with 8000 L/day (37 kg/cow/day). On both farms, cows were milked three times a day and had a sprinkler and ventilation system in the milking waiting barn. Drying therapy was performed 60 days before the expected date of calving and those dry cows were kept in a separate compost barn with a specific diet.
Three experimental groups were designed as follows: Group 1 (healthy control): cows without any puerperal disease; Group 2 (retained placenta – RP): cows that did not release fetal membranes until 12 h after calving (Ferreira Ademir de Moraes, 2010); Group 3 (clinical mastitis – CM): cows that presented Grade 1 CM according to the definition proposed by Roberson (2012) until the end of the voluntary waiting period (30 DPP). Cows diagnosed with mastitis after this period were not included in the study.
A vaccination calendar was followed regularly and reproductive vaccines against bovine herpes virus 1 (BoHV-1), bovine viral diarrhea virus (BVDV) and leptospirosis, in addition to rabies and clostridiosis were performed. Cows were treated with bovine recombinant somatotropin (bST) from 60 DPP with an interval of 14 days between applications until reaching 190 days of pregnancy.
For the reproductive management, cows were subjected to the following timed artificial insemination (TAI) protocol: —day 0: insertion of an intravaginal device containing 1.9 g of progesterone (P4) (CIDR®, Zoetis), application of 2.0 mg of estradiol benzoate i.m. (Sincrodiol®, Ouro Fino), and 0.01 mg of GnRH analog i.m. (Sincroforte®, Ouro Fino); day 7: application of 25 mg of dinoprost tromethamine i.m. (Lutalyse®, Zoetis); day 9: intravaginal device withdrawal, 25 mg of dinoprost tromethamine i.m. (Lutalyse®, Zoetis) and 1.2 mg of estradiol cypionate, i.m. (ECP®, Zoetis); day 11: timed AI for all treated animals.
Reproductive evaluations were made weekly. An ultrasound device equipped with a 7.5-MHz linear rectal transducer (Mindray® model DP-2200 VET) was used. Pregnancy diagnosis was performed 32 ± 4 days after AI. Reproduction performance was evaluated by: conception rate at 1st AI; pregnancy rate at 150 DPP, calving-to-1st AI interval and days open.
2.2. Experiment 2
2.2.1. Herd management
In experiment 2, 30 crossbred ¾ Holstein × ¼ Gir dairy cows from three commercial dairy farms were used. Division into groups and inclusion criteria were similar to Experiment 1 as Group 1: healthy control; (n = 10); Group 2: cows with RP were selected (RP-group; n = 10) and Group 3: cows with CM during the postpartum were selected (CM group; n = 10). Only animals in the third lactation and with a BCS between 3 and 3.5 (Edmonson et al., 1989) were included.
The animals received supplementation based on corn silage in addition to a concentrate. Water and mineral salt were available ad libitum. Farm 1 had 280 lactating cows (7500 L/day of milk production; 27 kg/cow/day), Farm 2 had 270 lactating (6000 L/day of milk production; 22 kg/cow/day) and Farm 3 had 195 lactating cows (4200 L/day of milk production; 22 kg/cow/day). In all farms, the cows were milked twice a day. Drying therapy was performed 60 days before the expected date of calving, and the dry cows received a specific diet.
A vaccination calendar was followed regularly complying with all the State law for dairy herd. Cows were treated with bovine recombinant somatotropin (bST) from 60 DPP with an interval of 14 days between applications until 190 days of pregnancy.
2.2.2. Clinical and gynecological examination
Between 40 and 50 DPP, clinical and gynecological examinations were performed. Cows were examined clinically for general health and vital parameters according to Feitosa (2008). The specific examination of the genital tract was performed by rectal palpation and an ultrasound evaluation of the uterine horns to check the uterus and intrauterine fluid was performed; in addition, the characteristics of vaginal discharge was assessed by Metricheck® device (McDougall et al., 2007; Sheldon et al., 2019).
2.2.3. Endometrial samples
After general examination, two endometrial cell samples were collected between 40 and 50 DPP using the cytobrush technique (Kasimanickam et al., 2004). The first sample was used to perform qRT-PCR. After collection, the cytobrush was stored in a sterile cryogenic tube with a screw cap free of DNAse, RNAse and heavy metals. After samples were collected, they were transferred in a liquid nitrogen bottle (1 L) and subsequently stored in a vertical ultra-freezer at − 80 ºC.
The second cytobrush sample was used for endometrial cytological slides for subclinical endometritis diagnosis. The slides were stained using the Rapid Panotic kit (Laborclin®, Pinhais, Brazil), according to the manufacturer's recommendations. Readings were made under an optical microscope at 40X magnification. A total of 100 cells were counted and the percentages of polymorphonuclear (PMN) cells and epithelial cells recorded. Subclinical endometritis was considered when the percentage of PMN cells (neutrophils) was ≥ 5% (Gilbert et al., 2005) compared to epithelial cells.
2.2.4. Real-time polymerase chain reaction (qRT-PCR)
The RNA extraction from the first endometrial cytobrush sample was performed using a Maxwell® 16 Instrument (Promega, Madison, WI, USA) with the extraction kit Maxwell® 16 LEV simply RNA Tissue Kit (AS1280, Promega, Madison, WI, USA). All steps were followed the manufacturer's recommendation. A total of 200 µL of homogenization solution were added to each cryogenic tube containing a cytobrush, then vortexed to recover the cytobrush endometrial secretion. After that, lysis buffer (200 µL) was added and tubes were vortexed again. The supernatants were collected and placed in the cartridges (disposable rack available in the kit). The cartridges were prepared using all reagents and subsequently inserted into the Maxwell® 16 Instrument for RNA extraction, after which the total RNA was quantified on a Quantus™ Fluorometer using the QuantiFluor® RNA System kit (E3310, Promega, Madison, WI, USA) according to the manufacturer's recommendations.
Reverse transcription (RT) was performed using the GoScript™ Reverse Transcription System kit (A5001, Promega, Madison, WI, USA). All steps were performed according to the manufacturer's recommendations. Initially, 5 µL of Mix qDNA removal were previously prepared with 10 µL of the respective sample and this mixture placed in the thermocycler for the first incubation. Meanwhile, RT Mix was prepared to be mixed in the sample tubes and incubated in a thermocycler again for RT enzyme denaturation. Finally, the cDNA of each sample was quantified on a Quantus™ Fluorometer using the QuantiFluor™ ssDNA System kit (E3190, Promega, Madison, WI, USA), and stored in a freezer at − 20°C.
Quantification of genes that encode the pro-inflammatory mediators such L-1β, IL-6, CCL5, and ESR1 and PGR receptors were performed with an Applied Biosystems 7300 Real-Time PCR System (Applied Biosystems, Foster, CA, USA), using the threshold cycle (Ct) comparative method. The Glyceraldehyde-3-Phosphate Dehydrogenase (GAPDH) gene was used as a control. An Axygen® PCR-96 AB-C 96-well plate (Axygen 214 Scientific, Union City, CA, USA) was used and the amplifications were performed in duplicate. The reaction conditions were determined previously by standardized cycles of 50 ºC for 2 min, 95 ºC for 2 min, 40 cycles of 95 ºC for 15 s and 60 ºC for 1 min.
For the reaction, a Quantinova™ Sybr® Green PCR kit was used. The proportions of the reagents used in the mix for each sample were: 5 µL of sybr green PCR Master MIX, 1.0 µL of QN ROX Reference Dye, 0.2 µL of forward primer, 0.2 µL of reverse primer, 1.6 µL of RNAse-free water and 2.0 µL of cDNA at a concentration of 35 ng/µL. The primers used to amplify IL-1β, IL-6, CCL5, ESR1 and PGR were shown in Table 1. The endogenous gene selected was GAPDH, as it is expressed in uterine tissue (Wathes et al., 2009).
Table 1 – Primer sequence for mRNA amplification by real-time polymerase chain reaction (qRT-PCR)
Gene
|
Primer direction
|
Primer sequence
|
IL1β
|
Forward
|
ACGAGTTTCTGTGTGACGCA
|
Reverse
|
TGCAGAACACCACTTCTCGG
|
IL6
|
Forward
|
GGGCTCCCATGATTGTGGTA
|
Reverse
|
GTGTGCCCAGTGGACAGGTT
|
CCL5
|
Forward
|
CATGGCAGCAGTTGTCTTTATCA
|
Reverse
|
CTCTCGCACCCACTTCTTCTCT
|
ESR1
|
Forward
|
CAGGCACATGAGCAACAAAG
|
Reverse
|
TCCAGCAGCAGGTCGTAGAG
|
PGR
|
Forward
|
TCCCCCCACTGATCAACTTG
|
Reverse
|
TCCGAAAACCTGGCAGTGA
|
GAPDH
|
Forward
|
GGTCACCAGGGCTGCTTTTA
|
Reverse
|
TTCCCGTTCTCTGCCTTGAC
|
Note: IL-1β – Interleukin 1β; IL-6 – Interleukin 6; CCL5 – Chemokine 5; ESR1 – α estrogen receptor; PGR – Progesterone receptor; GAPDH – Glyceraldehyde-3-phosphate dehydrogenase.
2.2.5. Statistical analysis
All data were recorded in a Microsoft Excel spreadsheet (Microsoft Corp., Redmond, WA) before statistical analysis. Statistical analyses were performed using SAS Studio (SAS Institute Inc., Cary, NC). Cows were considered as an experimental unit.
Distributions and normality tests were performed using the Univariate procedure. Normality was visually assessed via the Kolmogorov Smirnov method. The number of calving-1st AI interval and number of days open were utilized as a dependent variable and assessed with ANOVA using a generalized linear mixed model, fitted using the GLIMMIX procedure of SAS with a Gaussian distribution and an identity link function, with farm as the random effect and groups (healthy vs. RP vs. CM) as an independent variable. A generalized linear mixed model with the GLIMMIX procedure was used to assess the conception rate at calving-1st AI interval and pregnancy rate at 150 days postpartum as binary dependent variables, using logistic regression with farm as the random effect and groups (healthy vs. RP vs. CM) as an independent variable. Only variables with a P-value of ≤ 0.15 were retained in the final model. Additionally, the effects of CM and RP on calving to calving-1st AI interval and days open were analyzed using univariate models assessed by log-rank tests with PROC LIFETEST, where survival analyses curves were made.
For quantitative variables Kolmogorov–Smirnov test was selected to verify whether or not the data presented parametric distribution. The variables with parametric distribution were subjected to analysis of variance (parametric ANOVA) and Bonferroni's Multiple Comparison post-test. Variables with non-parametric distribution were analyzed using the Kruskal–Wallis test (non-parametric ANOVA) and Dunn’s Multiple Comparison post-test. Statistical significance was defined as P ≤ 0.05 for the main effects and interactions.