Associations between Heart Rate Variability Parameters and Housing- and Individual-Related Variables in Dairy Cows Using Canonical Correspondence Analysis

Noninfectious Causes of Pregnancy Loss at the Late Embryonic/Early Fetal Stage in Dairy Cattle

In cattle, initial pregnancy diagnosis takes place during the late embryonic/early fetal stage of gestation. From this point onward, pregnancy loss may occur in up to one fifth of pregnancies before the initial pregnancy diagnosis is confirmed. This means the early identification of risk factors is a key part of pregnancy diagnosis and herd management. The various factors responsible for pregnancy losses are classified into infectious and noninfectious. Among the noninfectious causes, several dam-related (circumstances of the individual pregnancy or milk production) and herd-related factors causing stress have been well established. In this review, we summarize the impacts of these noninfectious factors and predict associated risks of pregnancy loss.

Nonlinear canonical correspondence analysis and its application

The canonical correspondence analysis (CCA) is a multivariate direct gradient analysis method performing well in many fields, however, when it comes to approximating the unimodal response of species to an environmental gradient, which still assumes that the relationship between the environment and the weighted species score is linear. In this work, we propose a nonlinear canonical correspondence analysis method (NCCA), which first determines the most appropriate nonlinear explanatory factor through two screenings by correlation and LASSO regression, and successively uses the linear regression method and the improved heuristic optimal quadratic approximation method to fit the chi-square transformation values of the response variables. Thus, our method effectively reflects the nonlinear relationship between the species and the environment factors, and a biplot is employed to visualize the effects of the later on the distribution of species. The results from applying this method over a real dataset show that the NCCA method not only maintains the advantages of the polynomial canonical correspondence analysis (PCCA) proposed by Makarenkov (2002), but also outperforms Makarenkov's method in explaining the variance of response variables.

Altered oral microbiome in Sudanese Toombak smokeless tobacco users carries a newly emerging risk of squamous cell carcinoma development and progression

There are an estimated 6-10 million smokeless tobacco (Toombak) users in Sudan, the majority being males. Toombak is known to be a carcinogenic product that is likely to modify the oral microbiome spatiality into a high-risk potential for the development and progression of oral cancer, but previous studies are lacking in this field. Here, we endeavour for the first time the exploration of the oral microbiome in key mucosal areas of the oral cavity and assess the microbiome variations in premalignant and oral squamous cell carcinoma (OSCC) samples from both users and non-users of Toombak. 16S rRNA sequencing was performed on DNA obtained from pooled saliva, oral mucosa and supragingival plaque from 78 Sudanese users and non-users of Toombak, aged between 20 and 70 years. In 32 of the pooled saliva samples, the mycobiome (fungal) environment was analysed through ITS sequencing. Then, 46 formalin-fixed paraffin-embedded samples of premalignant and OSCC samples were collected, and their associated microbiomes sequenced. The oral Sudanese microbiome was found to be enriched in Streptococcaceae, but Staphylococcaceae were significantly more abundant amongst Toombak users. Genera enriched in the oral cavity of Toombak users included Corynebacterium_1 and Cardiobacterium while in non-users, Prevotella, Lactobacillus and Bifidobacterium were prominent. Aspergillus was the most abundant fungus in the mouths of Toombak users with a marked loss of Candida. The genus Corynebacterium_1 was abundant in the buccal, floor of the mouth and saliva microbiomes as well as in oral cancer samples from Toombak users indicating a possible role for this genus in the early stages of oral cancer development. An oral cancer microbiome that favours poor survival and metastasis in those who use Toombak also emerged that includes the genera Stenotrophomonas and Schlegelella. Those utilising Toombak carry an altered oral microbiome that may be an additional risk factor for this products carcinogenicity to the oral structures. These significant microbiome modulations are a newly emerging key driving factor in oral cancer development and progression in Toombak users while it is also shown that Toombak users carry an oral cancer microbiome that may increase the potential for a poorer prognosis.

Associations between cow-level parameters and heart rate variability as a marker of the physiological stress response in dairy cows

To maintain and enhance cow productivity and welfare, it is important that we can accurately assess and understand how cows respond to the physiological demands of gestation and lactation. Several methods have been developed for assessing the physiological responses to stressors and for detecting distress in cattle. Heart rate (HR) variability (HRV) is a non-invasive measure of autonomic nervous system activity and consequently a component of the physiological response to stress. In cattle, HRV has been successfully used to measure autonomic responses to a variety of health conditions and management procedures. The objectives of this study were to determine whether, among commercial Holstein Friesian cows and across farms, relationships exist between cow-level factors, HR and HRV. HRV parameters were compared with production records for 170 randomly selected, Holstein-Friesian-cows on 3 commercial dairy farms. Production data included parity, days in milk (DIM), milk yield, somatic cell count (SCC), % butterfat and protein, body condition score (BCS) and genetic indices. Fixed-effect, multivariable linear regression models were constructed to examine the association between cow-level variables and HRV parameters. Statistically significant relationships were found between HR and farm, temperature and BCS, and between HRV parameters and farm, rectal temperature, BCS, DIM, and percentage butterfat. Given the significant association between farms and several of the indices measured, it is recommended that care must be taken in the interpretation of HRV studies that are conducted on animals from a single farm. The current study indicated that within clinically normal dairy cattle HRV differed with the percentage of butterfat and BCS. Based on the relationships reported previously between HRV and stress in dairy cattle these results suggest that stress may be increased early in lactation, in cows with BCS <2.75 that are producing a high percentage of butterfat milk. Future work could focus on the physiological mechanisms through which these factors and their interactions alter HRV and how such physiological stress may be managed within a commercial farm setting.

Effects of Heat Stress on Heart Rate Variability in Free-Moving Sheep and Goats Assessed With Correction for Physical Activity

Heart rate variability (HRV) is the heart beat-to-beat variation under control of the cardiovascular function of animals. Under stressed conditions, cardiac activity is generally regulated with an upregulated sympathetic tone and withdrawal of vagal tone; thus, HRV monitoring can be a non-invasive technique to assess stress level in animals especially related to animal welfare. Among several stress-induced factors, heat stress is one of the most serious causes of physiological damage to animals. The aim of this study was to assess the effects of heat stress on HRV in small ruminants under free-moving conditions. In three experimental periods (June, August, and October), inter-beat intervals in sheep and goats (three for each) in two consecutive days were measured. HRV parameters were calculated from the inter-beat interval data by three types of analyses: time domain, frequency domain, and non-linear analyses. The temperature–humidity index (THI) was used as an indicator of heat stress, and vectorial dynamic body acceleration (VeDBA) was calculated to quantify the physical activity of the animals tested. First, we investigated correlations of THI and VeDBA with HRV parameters; subsequently, THI was divided into five categories according to the values obtained (≤ 65, 65–70, 70–75, 75–80, and >80), and the effects of the THI categories on HRV parameters were investigated with and without correcting for the effects of physical activity based on the VeDBA. The results indicated that HRV significantly decreased with increasing THI and VeDBA. For non-linear HRV parameters that were corrected for the effects of physical activity, it was suggested that there would be a threshold of THI around 80 that strongly affected HRV; high heat stress can affect the autonomic balance of animals non-linearly by inducing the sympathetic nervous system. In conclusion, to assess psychophysiological conditions of unrestrained animals by HRV analysis, the confounding effect of physical activity on HRV should be minimized for a more precise interpretation of the results.

Rearing in female-only groups and dietary mineral supplementation improves sow welfare in the early parities and lifetime performance

The lifetime performance of commercial sows relies on longevity, which is dependent on good health and reproductive performance. However, there is a high rate of wastage of sows in the early parities, which is influenced by the way they are managed and housed during rearing. This study investigated the carry-over effect of gilt rearing strategy on the measures of welfare and performance. Eighty sows were reared using a two by two factorial design: rearing group composition [GC; female-only (FEM) or mixed-sex (MIX) from weaning] with or without supplementary minerals (CON = control diet; SUPP = control + Cu, Zn, and Mn) from 5 wk into the finisher stage. Once served, gilts were managed in a dynamic group gestation pen and fed a standard gestating sow diet. Locomotory ability was scored (0 to 5) and salivary cortisol measured five times during the first gestation, and human approach tests were carried out on day 108. Hooves were scored for injuries and legs for bursas at day 70 of the first gestation, at first weaning, and at the second farrowing. Sow behavior in the hoof scoring crate (movement, vocalization, and handling ease) was also recorded. The number of piglets born alive and dead during the first five parities was recorded as was the performance of the first litter to finish. Data were analyzed using general or generalized linear mixed models, as appropriate, using SAS (v 9.4). There was no effect (P > 0.05) of rearing treatment on locomotory ability, bursa score, the total number of piglets born, or on offspring growth. However, there was an interaction between GC and supplementary minerals (P < 0.05) on salivary cortisol levels with MIX × SUPP sows having the highest levels. Total hoof scores and heel erosion scores were higher in sows reared in MIX groups (P < 0.01), and CON sows tended to have higher horizontal crack scores (P = 0.06). Sows from MIX kicked more at weaning than FEM (P < 0.05) and tended to be more fearful in the forced human approach test (P = 0.1) where they are scored on their reaction to being approached. They also had more stillborn piglets across all five parities than FEM (P < 0.05). Overall, rearing replacement sows in FEM groups and dietary mineral supplementation had minimal but beneficial effects on their subsequent welfare and performance.

Clinical Significance of Fractional Anisotropy Measured in Peritumoral Edema as a Biomarker of Overall Survival in Glioblastoma: Evidence Using Correspondence Analysis

Introduction

Fractional anisotropy (FA), a diffusion tensor image (DTI) derived biomarker is related to invasion, infiltration, and extension of glioblastoma (GB). We aimed to evaluate FA values and their association with intervals of overall survival (OS).

Materials and Methods

Retrospective study conducted in 36 patients with GB included 23 (63.9%) males, 46 ± 14 y; and 13 (36.1%) females, 53 ± 13; followed up for 36 months. We measured FA at edema, enhancing rim, and necrosis. We created two categorical variables using levels of FA and intervals of OS to evaluate their relationships. Kaplan-Meier method and correspondence analysis evaluated the association between OS (grouped in 7 six-month intervals) and FA measurements.

Results

Median FA values were higher in healthy brain regions (0.351), followed by peritumoral edema (0.190), enhancing ring (0.116), and necrosis (0.071). Pair-wise comparisons among tumor regions showed a significant difference, P < 0.001. The median OS for all patients was 19.3 months; variations in the OS curves among subgroups was significant χ² (3) = 8.48, P = 0.037. Correspondence analysis showed a significant association between FA values in the edema region and the survival intervals χ² (18) = 30.996, P = 0.029.

Conclusions

Alternative multivariate assessment using correspondence analysis might supplement the traditional survival analysis in patients with GB. A close follow-up of the variability of FA in the peritumoral edema region is predictive of the OS within specific six-month interval subgroup. Further studies should focus on predictive models combining surgical and DTI biomarkers.

Expression of heat-shock protein 72 mRNA in relation to heart rate variability of Sahiwal and Karan-Fries in different temperature-humidity indices

AIM

To investigate the effect of temperature-humidity index (THI) on the expression pattern of heat-shock protein 72 (HSP72) mRNA of Sahiwal and Karan-Fries (KF) cattle in different THIs.

MATERIALS AND METHODS

Five different periods of a year were selected based on combinations of T_max/T_min, viz., P1: <20°C/<10°C; P2: >20°C/<10°C, P3: <30°C/<15°C; P4: >35°C/<20°C, and P5: >35°C/>20°C. The THI was calculated from the records of temperature and relative humidity in different periods. Heart rate variability (HRV) was calculated from electrocardiogram records in different periods. HSP72 mRNA expression was estimated by reverse transcription polymerase chain reaction.

RESULTS

The THI recorded during P1, P2, P3, P4, and P5 were 55.5, 60.3, 70.1, 74.5, and 79.0, respectively. THI in P4 and P5 were stressful to animals. HSP72 mRNA expression increased during cold stress in P1 in Sahiwal and heat stress in P4 and P5 in both Sahiwal and KF. Sahiwal maintained increased HSP72 mRNA expression longer than KF without causing a significant change in HRV.

CONCLUSION

Both low THI in winter and high THI in summer increased HSP72 mRNA of Sahiwal and KF without significant change in HRV. Thermotolerance of Sahiwal could be due to the maintenance of higher HSP72 expression longer than KF in prolonged heat stress in summer.

Timing of obstetrical assistance affects peripartal cardiac autonomic function and early maternal behavior of dairy cows

Peripartal autonomic nervous system function and early maternal behavior were investigated in 79 multiparous Holstein-Friesian cows. Animals were allocated into four groups based on the technology of calving management: 1) unassisted calving in a group pen (UCG; N=19), 2) unassisted calving in an individual pen (UCI; N=21), 3) assisted calving with appropriately timed obstetrical assistance (ACA; N=20), and 4) assisted calving with premature obstetrical assistance (ACP; N=19). Heart rate, the high frequency (HF) component of heart rate variability (HRV) as a measure of vagal activity and the ratio between the low frequency (LF) and HF components (LF/HF ratio) as a parameter of sympathetic nervous system activity were calculated. Heart rate and HRV parameters were presented as areas under the curves (AUC) for the following periods: 1) prepartum period (between 96h before the onset of calving restlessness and the onset of restlessness), 2) parturition (between the onset calving restlessness and delivery), and 3) postpartum period (during a 48-h period after delivery). Pain-related behaviors were recorded during parturition (i.e., the occurrence of vocalization and stretching the neck towards the abdomen) and during a 2-h observation period after calving (i.e., the occurrence of vocalization, stretching the neck towards the abdomen and the duration of standing with an arched back). Early maternal behavior was observed during the first 2h following calving as follows: 1) latency and duration of sniffing calf's head/body, and 2) latency and duration of licking calf's head/body. No difference was found across groups in autonomic function before the onset of calving restlessness. Area under the heart rate curve was higher in ACP cows during parturition (39.6±2.5beats/min×h) compared to UCG, UCI and ACA animals (AUC=13.1±0.9beats/min×h, AUC=22.3±1.4beats/min×h and AUC=25.0±2.1beats/min×h, respectively). Area under the heart rate curve did not differ across the UCG, UCI and ACA groups during the postpartum period (AUC=65.2±16.7beats/min×h, AUC=58.0±14.2beats/min×h and AUC=62.9±12.1beats/min×h, respectively) but it was higher in ACP cows compared to the former groups (AUC=269.1±36.3beats/min×h). During parturition, area under the HF curve reflected a lower vagal tone (AUC=-30.5±1.6n.u.×h) in cows with premature obstetrical assistance than in animals that calved individually without farmer assistance (AUC=2.7±0.4n.u.×h) or with appropriately timed assistance (AUC=3.2±1.2n.u.×h). During parturition, LF/HF ratio showed greater sympathetic activity in ACP cows than in animals from any other group. Area under the HF curve was similar across UCG, UCI and ACA cows (AUC=-232.1±42.0n.u.×h, AUC=-163.4±35.6n.u.×h and AUC=-331.4±56.2n.u.×h, respectively) during the postpartum period and was the lowest in ACP cows (AUC=-1025.6±44.2n.u.×h) reflecting a long-term stress load in the latter group. During parturition, both vocalization and stretching the neck towards the abdomen occurred more often in UCG cows than in cows from any other groups, and the incidence of both behaviors was statistically higher in ACP cows than in UCI and ACA animals. There were no significant differences across groups in these behaviors during the 2-h postpartum observation. UCG cows had a shorter latency and a longer duration of maternal grooming during the first 2h following delivery compared to any other groups. UCI and ACA dams spent more time with licking the calf within the 2-h period after calving and had a shorter latency to sniff and lick the offspring compared to cows that received premature assistance. Group calving is less stressful for cows than calving in an individual pen either with or without obstetrical assistance. Calving in a group or with appropriately timed farmer assistance supports the expression of early maternal behavior and lead to a rapid postpartum recovery of the autonomic nervous system. Premature obstetrical assistance means stress for cows during parturition, leads to a prolonged postpartum recovery of the autonomic nervous system and inhibits the expression of early maternal behavior.

Hypothalamic-pituitary-adrenal and cardiac autonomic responses to transrectal examination differ with behavioral reactivity in dairy cows

Behavior, hypothalamic-pituitary-adrenal axis, and cardiac autonomic nervous system (ANS) activity were evaluated in response to transrectal examination in nonlactating Holstein-Friesian cows with different behavioral reactivity. According to behavioral reactions shown to the procedure of fixing the heart rate (HR) monitors, the 20 cows with the highest and the 20 cows with the lowest behavioral reactivity were involved in the study (high responder, n=20; and low responder, n=20, respectively). Activity of the ANS was assessed by HR and HR variability parameters. Blood and saliva were collected at 5 min before (baseline) and 0, 5 10, 15, 20, 30, 40, 60, and 120 min after the examination to determine cortisol concentrations. The examination lasted for 5 min. Cardiac parameters included HR, the root mean square of successive differences between the consecutive interbeat intervals, the high frequency (HF) component of heart rate variability, and the ratio between the low frequency (LF) and HF parameter (LF/HF). Following the examination, peak plasma and saliva cortisol levels and the amplitude of the plasma and saliva cortisol response were higher in high responder cows than in low responders. Areas under the plasma and saliva cortisol response curves were greater in high responder cows. Plasma and salivary cortisol levels correlated significantly at baseline (r=0.91), right after examination (r=0.98), and at peak levels (r=0.96). Area under the HR response curve was higher in low responder cows; however, maximum HR and the amplitude of the HR response showed no differences between groups. Minimum values of both parameters calculated for the examination were higher in high responders. Following the examination, response parameters of root mean square of successive differences and HF did not differ between groups. The maximum and the amplitude of LF/HF response and area under the LF/HF response curve were lower in low responder cows, suggesting a lower sympathetic activation of the ANS. Although changes in behaviors indicated that the procedure was painful for the animals, no differences were observed either in vocalization or in attendant behavior between groups during the examination. Our results demonstrate that behaviorally more reactive animals exhibit increased plasma and salivary cortisol concentrations and higher cardiac autonomic responsiveness to transrectal examination than less reactive cows. Salivary cortisol may substitute for plasma cortisol when assessing response of cattle to stress.