Next Article in Journal
Innovative Use of Olive, Winery and Cheese Waste By-Products as Functional Ingredients in Broiler Nutrition
Previous Article in Journal
Describing and Characterizing the Literature Regarding Umbilical Health in Intensively Raised Cattle: A Scoping Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Repeatability of Health and Welfare Traits and Correlation with Performance Traits in Dairy Goats Reared under Low-Input Farming Systems

by
Sotiria Vouraki
1,*,
Athanasios I. Gelasakis
2,
Vasileia Fotiadou
1,
Georgios Banos
1,3 and
Georgios Arsenos
1
1
Laboratory of Animal Husbandry, School of Veterinary Medicine, Faculty of Health Sciences, Aristotle University, 54124 Thessaloniki, Greece
2
Department of Animal Science, School of Animal Biosciences, Agricultural University of Athens, 11855 Athens, Greece
3
Department of Animal and Veterinary Sciences, Scotland’s Rural College (SRUC), Easter Bush, Midlothian EH25 9RG, UK
*
Author to whom correspondence should be addressed.
Submission received: 12 May 2022 / Revised: 1 June 2022 / Accepted: 9 June 2022 / Published: 11 June 2022
(This article belongs to the Section Veterinary Biomedical Sciences)

Abstract

:
The objectives of the study were to estimate the repeatability of health and welfare traits and investigate their association with performance in three breeds of dairy goats reared under low-input farming systems in Greece. A total of 1210 goats of Eghoria (n = 418), Skopelos (n = 429), and Damascus (n = 363) breeds were assessed. Udder health, parasitic resistance, welfare, milk yield and quality, and body condition score were recorded monthly for two milking periods. Udder health records included somatic cell count (SCC) and total viable count (TVC). Based on combinations of SCC and TVC and thresholds set at >106 cells/mL and >2 × 104 cfu/mL, respectively, additional udder health phenotypes were defined. Parasitism included myiasis, tick infestation, gastrointestinal nematode (GIN) and cestode faecal egg count (FEC), and lungworm faecal larval count (FLC). Infection with each of the endoparasites was defined based on FEC/FLC. Welfare assessment parameters included the presence of ear and horn injuries, ocular and nasal discharge, body and udder abscesses, injury and lesions on the skin of different regions, diarrhoea, hernias, overgrown hooves, arthritis, lameness, and udder asymmetry. Trait repeatability and animal correlations were estimated. Significant (p < 0.05) repeatability was reported for all udder health and most welfare traits in all breeds, GIN and cestode FEC, and GIN and lungworm infection in Eghoria, and myiasis in Skopelos. Correlations of health and most of welfare traits with performance were non-significant or favourable. Overall, results demonstrate potential to improve health and welfare of the studied breeds without compromising performance.

1. Introduction

Dairy goat farming is a significant animal husbandry activity with an important socioeconomic and environmental role in disadvantaged areas [1]. Given the ability of goats, especially of indigenous breeds, to utilise diverse vegetation in mountainous and marginal regions, they are traditionally reared under low-input pastoral farming systems [2,3]. The notion is that the sustainability of such systems depends on animal productivity [2]. However, the latter is defined by the status of animal health and welfare that is directly associated with environmental conditions, housing conditions, and herd health management [4,5].
Grazing at natural pastures favours internal and external parasitism, which affects goat health and welfare status [4,5,6]. Moreover, traditional goat shelters are frequently characterised by high stocking density, inadequate bedding, and poor hygiene, which predispose animals to bacterial infections, particularly subclinical mastitis, with a cumulative incidence of 24–32% recently reported in low-input pastoral farming systems [7]. Additionally, poor husbandry conditions and hygiene status can lead to welfare issues including injuries, abscesses and skin lesions, udder problems, and limb disorders [4,5,8].
Preventive and treatment measures to mitigate the occurrence and consequences of the above health and welfare issues can be challenging. Specifically, the repeated use of anthelmintics for controlling gastrointestinal nematodes (GIN) has led to increasing anthelmintic resistance [9,10,11]. At the same time, chemical control of ectoparasites with zero withdrawal time is perceived as rather expensive by the farmers [12]. Moreover, preventive strategies against subclinical mastitis can be cumbersome and expensive, especially for large-size herds [13]. Genetic selection for disease-resistant animals is a promising complementary approach to the above measures. Previous studies have shown the possibility for selective breeding towards enhanced udder health [14,15,16,17] and GIN resistance [18,19,20] in goats.
Genetic selection, however, is often limited by the lack of reliable pedigree records in low-input pastoral dairy goat farms [5,20]. Nevertheless, investigating between-animal (co)variation of performance, health, and welfare traits, even without pedigree data, could still inform decision-making for improved management and selection practices at the animal level towards enhanced overall productivity. Such an investigation is feasible by estimating the repeatability of health and welfare traits and the animal correlations with performance traits. This could help identify any potential for improving goat health and welfare status without compromising performance.
In goats, relevant studies are scarce. Previous research has investigated the repeatability of somatic cell count (SCC) [14,16,17,21]. To the best of our knowledge, there are no references regarding other indicators of udder health status such as total viable count (TVC). Moreover, studies on goat parasitism have reported trait repeatability for GIN based on faecal egg count (FEC) [9,10,22,23], but there has been no research on resistance to other endoparasites or ectoparasites. There have been no relevant studies on welfare-related traits and their association with goat performance.
The objectives of the present study were to (i) estimate the repeatability of a wide range of health and welfare traits measured on goats of three breeds reared under low-input pastoral farming systems in Greece, and (ii) investigate their association with milk production and body condition score (BCS).

2. Materials and Methods

2.1. Farms and Animals

A total of 1210 dairy goats of two indigenous Greek breeds (Eghoria, n = 418 and Skopelos, n = 429) and one imported breed (Damascus, n = 363) were used. Animals were randomly selected from seven farms (two with Eghoria, two with Skopelos, and three with Damascus goats) located in northern and central Greece (Figure 1). All herds were selected as representatives of the studied breeds and the low-input pastoral farming system, the typical goat production system in Greece [24]. This system is characterised by grazing shrubland and woodland areas throughout the year, traditional management practices, and random mating of animals.

2.2. Phenotypic Data and Recording Protocols

Data collected during FP7-SOLID project (Sustainable Organic Low Input Dairying, 266367, 2011–2015) were used for this study. The available dataset included individual monthly records for two consecutive milking periods (corresponding to kiddings between November 2011 and January 2012, and November 2012 and January 2013; five records per milking period) of animal health, welfare, and production traits.
Specifically, a standard protocol for welfare trait assessment was developed based on Anzuino et al. [25]. Prior to milking, each goat was restrained and assessed through observation and palpation of the head, body, udder, and limbs for ear and horn injuries, head skin lesions (hair loss and/or dermatitis), ocular and nasal discharge, body abscesses and injuries, diarrhoea, hernia, overgrown claws, arthritis, udder asymmetry, udder abscess, and udder skin lesions. Moreover, a locomotion score was assigned using the 5-point scale (1: normal gait; 2: no obvious lameness when standing, occasional limping when walking; 3: lifting foot while standing and moderate lameness when walking; 4: shifting stance and severe lameness when walking; 5: unwilling to bear weight on one foot when standing or walking) of Ley et al. [26]. All the above traits were consistently assessed by the same veterinarian. A detailed description of the welfare traits assessed is available in the study by Gelasakis et al. [5].
All goats were also thoroughly examined individually for the presence of hard ticks throughout the body and warbles beneath the skin on the back caused by Przhevalskiana silenus larvae. Additionally, a faecal sample was collected to measure FEC of GIN and cestodes, and faecal larval counts (FLC) of lungworms using the modified McMaster technique [27].
The BCS of each goat was assessed by palpation of the dorsal lumbar region and recorded on the 5-point scale (1: emaciated to 5: obese, in 0.5-point increments) of Russel et al. [28]. Subsequently, a milk sample was collected aseptically from both udder halves to be tested for TVC. Each goat was then hand-milked into a bucket, the milk was weighed, and a milk sample was collected from the milking bucket to assess SCC and milk composition (fat, protein lactose, and solids-non-fat (SNF) content). Sampling and analyses of milk samples have been described in detail in previous studies [7,29,30].

2.3. Phenotypic Data Handling

The official AT method of the International Committee of Animal Recording [31] was used to calculate individual daily milk yield of the studied goats by doubling the recorded yield of the morning or evening milking. Daily milk fat, protein, lactose, and SNF yields were also calculated from milk and the corresponding content records.
Moreover, some additional health and welfare phenotypes were defined based on the measurements described above. Specifically, a subclinical mastitis index (SMI) was derived using the following formula:
SMI = 0.6 SCC     Mean SD   +   0.4 TVC     Mean SD
where Mean and SD correspond to the mean and standard deviation, respectively, of each trait (SCC and TVC) within each studied breed.
Thresholds were also set for SCC and TVC at >106 cells/mL and >2 × 104 cfu/mL [7], respectively, and used to define three additional udder health phenotypes: (i) UHP1, scored as 0 or 1 if at least one of the traits was below or both were above the thresholds, respectively, (ii) UHP2, scored as 0–2 if both traits were below, only one was above, or both exceeded the thresholds, and (iii) UHP3 scored as 0–3 if both traits were below, only TVC was above, only SCC was above, or both were above the thresholds.
Furthermore, FEC and FLC records were used to determine the presence or absence of infection with GIN and cestodes, and lungworms, respectively. In all cases, records of ≥50 FEC or FLC were indicative of infection [32]. We also defined an overall endoparasite infection index, where at least one of these parasites was detected. Moreover, the presence of at least 10 ticks attached to the body of goats was considered as tick infestation [6,33]. A myiasis case was defined by the presence of at least two warbles beneath the skin on the back of the animal.
Regarding welfare traits, the locomotion score assigned to each animal was used to define the presence or absence of lameness. In addition, measurements on each body part were used to define the presence or absence of head, body, limb and udder problems, and skin lesions and injuries.
Descriptive statistics of all animal traits considered in the study are shown in Table 1, Table 2, Table 3 and Table 4. The final dataset used for the analyses is presented in Dataset S1.

2.4. Data Analysis

Milk production traits SCC and TVC were logarithmically transformed to ensure normality of distributions. FEC and FLC records, which were also significantly skewed, were transformed using Tukey’s Ladder of Powers with R package “rcompanion”.
Preliminary analyses were performed to identify environmental factors with statistically significant (p < 0.05) effects on the studied traits. The effects of farm, period of kidding, age at kidding, days from kidding, and interactions between them were tested.
(Co)variance components of health and welfare phenotypes with milk production traits and BCS were estimated within breed in a series of bivariate statistical analyses implemented with Bayesian Markov Chain Monte Carlo methods and the R software package “MCMCglmm” [34]. All analyses were based on the following model:
Y ijmn = μ + F i P j + b 1 * A + b 2 * D + A m + e ijmn
where Yijmn = studied trait (nth measurement on animal m); μ = overall population mean; FiPj = fixed effect of the interaction between farm (I = two levels for Skopelos and Eghoria and three for Damascus) and period of kidding (j =two levels); b1 = regression coefficient on age at kidding A (months); b2 = regression coefficient on days from kidding D; Am = random effect of the animal m; eijmn = random residual effect.
A probit link function for binomial distribution was fitted to model 2 for the analyses of binary traits.
Weekly informative priors were used for the random animal and residual effects, whereas for the fixed effects the default normal prior distribution with null mean and a large variance (1010 was used. For binary traits, residual variance was fixed to one [34]. For each bivariate model, three chains of 13,000 to 2,300,000 iterations with a burn-in period of 3000 iterations and a thinning interval of 10 to 2000 samples were used (depending on each model’s convergence and autocorrelation diagnostics). Convergence of the models was tested with visual inspection of estimate plots and the Gelman and Rubin’s convergence diagnostic, where values above 1 indicate lack of convergence [35]. Moreover, autocorrelation across chains was tested for all lag values greater than zero (values below 0.1 were considered acceptable).
Repeatability and animal correlations between studied traits were estimated from the posterior means of corresponding variance and covariance values after convergence. In the case of binary traits, results were expressed on the latent (liability) scale. Repeatability estimates of all studied traits were derived as the mean of estimates from multiple bivariate analyses.

3. Results

3.1. Trait Repeatability Estimates

Estimates of repeatability for health and welfare traits within breed (Eghoria, Skopelos, and Damascus) are shown in Table 5 and Table 6, respectively; for comparison, estimates of milk production traits and BCS are listed in Supplement Table S1. Statistically significant (p < 0.05) estimates were derived for all udder health traits ranging from 0.08 to 0.50, 0.09 to 0.59, and 0.08 to 0.43 in Eghoria, Skopelos, and Damascus goats, respectively. In all breeds, the most repeatable udder health trait was UHP1 (Table 5). Of the parasitism traits, significant repeatability estimates were found for GIN and cestode FEC, and GIN, lungworm and overall endoparasite infection in Eghoria goats (0.09–0.32) and for cestode infection (0.27) and myiasis (0.34) in Skopelos goats (Table 5). Regarding welfare, significant repeatability estimates were found for all traits except for ocular discharge, diarrhoea, injury and hernia in all breeds, lameness and arthritis in Eghoria and Skopelos goats, and nasal discharge in Skopelos goats. Significant repeatability estimates for welfare traits were 0.16–0.99, 0.19–0.97, and 0.20–0.98 in Eghoria, Skopelos, and Damascus goats, respectively; the most repeatable traits were ear, horn, and total injuries (Table 6). All studied milk production traits were significantly repeatable with estimates ranging from 0.26 to 0.41, 0.35 to 0.47, and 0.21 to 0.52 in Eghoria, Skopelos, and Damascus goats, respectively (Supplementary Table S1). Significant repeatability estimates were also found for BCS in all three breeds (0.31–0.47; Supplement Table S1).

3.2. Animal Correlations of Health and Welfare Traits with Milk Production Traits and Body Condition Score

Statistically significant animal correlations of health and welfare traits with milk production traits and BCS within breed are shown in Table 7 and Table 8, respectively. All animal correlations estimated are presented in detail in Supplement Tables S2–S4. Notably, a negative sign in the correlation is in the desirable direction signifying reduced incidence of health and welfare issues as performance improves. Significant negative correlations were reported between most udder health and milk production traits in all breeds (Table 7). Moreover, GIN and endoparasite infections were also negatively correlated with milk yield, lactose yield, and SNF yield in Eghoria goats (Table 7). Regarding welfare traits, negative correlations were reported between head skin lesions and protein yield, and total skin lesions and fat yield in Eghoria goats, and between udder asymmetry and all milk production traits in Damascus goats. In Damascus goats, arthritis and limb problems were also negatively correlated with BCS (Table 8). In contrast, udder abscess was positively correlated with protein yield in Eghoria goats and with most milk production traits in Skopelos goats (Table 8). In the latter breed, positive correlations were also identified between udder problems and protein yield, and of overgrown claws and limb problems with all milk production traits (Table 8).

4. Discussion

This study investigated the repeatability of a wide range of animal health and welfare traits and their association with milk production and BCS in three dairy goat breeds reared under low-input pastoral farming systems in Greece. Overall, significant between-animal variation was revealed for udder health, resistance to parasitism, and most studied welfare traits in Eghoria, Skopelos, and Damascus goats. Moreover, animal correlations of health and welfare with performance traits were mostly either non-significantly different from zero or favourable, implying that improving one would also benefit the other set of traits.
Few previous studies have estimated the repeatability of milk SCC, an indicator of subclinical mastitis, and significant estimates (0.31–0.59) have been reported in Saanen, Alpine, and mixed populations of goats [14,16,17,21]. Repeatability estimates of SCC reported in the present study for Eghoria and Skopelos goats were within the aforementioned range. The value found in Damascus goats was lower (0.22), which could be possibly attributed to differences in the studied breeds and/or farming systems and practices.
The ability of milk SCC to predict subclinical mastitis is reportedly lower in goats compared to dairy sheep and cows since it may be influenced by many more physiological and environmental factors in the former [36,37,38]. Nevertheless, there is no literature on the repeatability of other udder health indicator traits. The present study is the first to investigate TVC as well as combinations of SCC and TVC in this regard. The latter included a weighted SMI in which a greater emphasis (60%) was placed in SCC due to ease of measurement and cost-effectiveness of the trait [37], and three additional udder health phenotypes (UHP1, UHP2, and UHP3) defined according to thresholds for SCC and TVC (>103 cells/mL and >2 × 104 cfu/mL, respectively). These thresholds were based on the study by Gelasakis et al. [7] in which ca. 80% of the goat milk samples had exceeded both thresholds and found to have positive milk microbiological cultures. According to our results, significant between-animal variation exists for all the above udder health traits in Eghoria, Skopelos, and Damascus goats. UHP1 seems to be the best candidate since it was the most repeatable trait (0.43–0.59). However, it should be noted that including UHP1 in management and selection practices would require routine recording of TVC, hence implying an additional cost for the farmer. In this regard, decisions could be based only on SCC given that significant repeatability estimates were found for this trait in all three breeds. Moreover, despite previous research suggesting a lower association with intramammary infections in goats compared to other species, Rupp et al. [13] recently provided evidence of SCC being a valuable predictor of subclinical mastitis in Alpine goats that could be efficiently used for selection purposes and informed culling decisions.
Regarding endoparasites, the repeatability of GIN FEC has been widely studied in extensively reared dairy goats and significant medium to high estimates (0.25–0.84) have been reported [9,10,22] in cases of natural infection. These findings are in broad accordance with those obtained in the present study in Eghoria goats, although much lower estimates (0.11) are reported here. This could be associated with the small sample size (n = 84) and number of records (n = 584) available in our study. There are no previous reports on the repeatability of cestode FEC and lungworm FLC. In our study, significant estimates were found only for cestode FEC in Eghoria goats. However, cestode and lungworm load manifested in FEC and FLC was low in our data, which may have led to an underestimation of between-animal variation [20]. Therefore, in a separate series of analyses, infection with each endoparasite was defined as a binary trait (presence or absence of infection). These analyses revealed increased between-animal variation compared to the corresponding traits based on FEC and FLC. Higher repeatability estimates (026–0.32) were derived for GIN, lungworm, and overall endoparasite infection in Eghoria goats, which were within the range of previous reports [9,22]. Since data distribution in our study indicates that lungworm and cestode infections may not be a serious problem in the three studied dairy goat breeds, future research with a higher sample size and parasitic load could help to further investigate between-animal variation for these traits.
To our knowledge, there are no previous reports on the repeatability of infestation with ectoparasites in goats. In beef cattle, significant estimates (ca. 0.50) have been reported for myiasis resulting from buffalo and horn flies [39,40]. These findings are in accordance with those of the present study in Skopelos goats for myiasis with Przhevalskiana silenus larvae. Moreover, previous reports in cattle [41,42,43,44] and sheep [45,46] have indicated significant repeatability estimates for tick infestation (0.21–0.45 and 0.39–0.44, respectively) under natural or experimental challenge. However, other estimates were found to not significantly differ from zero in Angus cattle [47]. The latter results are in general agreement with those of our study in Eghoria and Skopelos goats. According to Giglioti et al. [44], animal age and frequency of examination and disease recording could influence the repeatability estimates of tick infestation. In the present study, animals of 2–5 years of age were included, and records were collected in one-month intervals. Further research within different age groups and with shorter examination intervals is warranted. Moreover, additional data with a higher frequency of ectoparasite infestation would be needed to properly estimate between-animal variance of these traits.
All studied welfare traits were significantly repeatable, except for those associated with problems in very low frequency (ocular discharge, diarrhoea, injuries, and hernias in all breeds, lameness and arthritis in Eghoria and Skopelos goats, and nasal discharge in Skopelos goats). According to available literature, our study is the first to evaluate the repeatability of welfare traits in goats. Moreover, there is a notable shortage of relevant literature in other farm animal species for most welfare-related traits. Currently, there are several genetic studies of lameness in dairy cows and significant heritability estimates (0.10–0.54) have been reported [48,49,50,51]. Heritability can indicate the lower limit of repeatability. In this regard, the moderate to high repeatability estimate of lameness found in Damascus goats is in general agreement with the results of the aforementioned studies, although the latter were on another species. Moderate to high repeatability estimates were also found for other studied limb problems (arthritis and overgrown claws) in the Damascus breed. Moderate estimates were found for udder problems with the most repeatable trait being udder asymmetry in all breeds. Udder asymmetry is often the result of a former intramammary infection and is usually irreversible [25], thus explaining its high repeatability estimate. The very high repeatability (>0.90) of ear and horn injuries could also be possibly attributed to records pertaining to the same injuries rather than new incidents. Therefore, further research on distinct incidents of injuries is warranted.
Overall, according to our results, between-animal variance accounts for a significant proportion of the total phenotypic variance of most of the studied traits in Eghoria, Skopelos, and Damascus goats. Such results suggest that there is potential to improve goat health and welfare through management and selection practices. Specifically, results indicate that it is possible to predict the future health and welfare status of individual animals and identify those most probable to (i) be repeatedly infected with GIN and, consequently, produce excretions responsible for pasture contamination in Eghoria breed, (ii) contribute to high incidence of myiasis in Skopelos breed, and (iii) be prone to udder health and welfare issues in all breeds. Based on the above, preventive and control strategies could be applied selectively to individual animals rather than in the whole herd. This may contribute to reducing anthelmintic and antibiotic resistance, while being cost-effective for the farmer. These findings are consistent with outcomes from previous studies on other dairy goat breeds [9,10,22]. Moreover, selection to reduce subclinical mastitis, GIN infections, and myiasis and improve the welfare status of animals could be applied based on goat records early in life, thereby enabling timely and informed culling decisions.
Furthermore, given that between-animal variation is partly genetic, the studied traits with significant repeatability could also be heritable. In our study, it was not possible to disentangle the additive genetic variance from the permanent environmental variance due to lack of pedigree data. Previous research in goats supports the possibility of selective breeding towards resistance to subclinical mastitis based on SCC [14,15,16,17] and resistance to GIN [11,18,19,52]. Myiasis resistance has also been shown to be moderately heritable in beef cattle (0.47) [39]. Finally, low but significant heritability estimates (≤0.15) have been generally reported for lameness in cattle and sheep [48,49,50,53] and for overall foot health in dairy cows [54].
According to the findings of the present study, implementing management and selection practices to improve goat health is not expected to compromise performance of the studied breeds. All animal correlations between udder health and milk production traits were either non-significant or favourable. Given that part of these correlations has a genetic basis, our findings are consistent with previous genetic correlations of SCC with milk, fat and protein yield reported in Alpine and Saanen goats [15]. However, Scholtens et al. [21] reported slightly positive and unfavourable correlations of SCC with milk and protein yield in New Zealand goats. Nevertheless, the generally favourable correlations between udder health and milk production traits suggest that it is possible to simultaneously improve udder health and milk production of the three studied breeds. Moreover, this is not expected to compromise goat BCS since all correlations with udder health traits did not significantly differ from zero.
In Eghoria goats, animal correlations of GIN infection with milk production traits were favourable, in agreement with the results of Morris et al. [23] in Saanen goats. However, unfavourable genetic correlations between these traits have been reported in Saanen goats by Heckendron et al. [11]. Contradictory results could be possibly explained by differences in GIN species, pasture infection pressure, and/or immunological responses of animals in the different studies [11,55]. Moreover, the present is the first study to estimate correlations between GIN FEC and BCS in goats. The non-significant estimates found in our study agree with the genetic correlations reported by Boareki et al. [56] in sheep. There is no literature regarding the association between myiasis and goat performance. In beef cattle, unfavourable genetic correlations have been reported of myasis with growth traits [40]. These findings were not supported by our study, where no significant animal correlations were found with BCS and milk production traits. The discrepancy may be attributed to species differences, as well as the type and severity of myiasis.
Our results suggest that the welfare status of the studied dairy goat breeds could also be improved without compromising their overall performance since, in most cases, no significant relevant associations were reported. Moreover, we found that skin lesions and udder asymmetry could be reduced by improving milk quality, based on the favourable animal correlations with some milk composition traits in Eghoria and Damascus goats, respectively. However, it should be noted that some welfare and milk production traits were unfavourably associated. Specifically, Skopelos goats with higher milk production and better milk quality seem to be more vulnerable to limb problems (mostly overgrown claws) and udder problems (mostly udder abscesses). Such results could be partly attributed to environmental factors such as grazing on soft and wet soil, inadequate flooring, and/or high supplementation with concentrate feedstuffs [5,25]. Nevertheless, between-animal covariance also includes genetic components. Relevant literature in goats is scarce. In dairy cows, an unfavourable genetic relationship between milk production and welfare indicators, including lameness and udder-type traits, has been reported [57]. Uribe et al. [54] also reported an antagonistic association between milk production traits and culling for limb problems. Such results are in agreement with the animal correlations of the present study for Skopelos goats. Therefore, selecting for higher milk production could possibly increase limb and udder issues in Skopelos goats. However, it might be possible to improve milk production while also reducing udder abscess through practices targeting minimisation of intramammary infections and udder injuries. Udder abscess usually results from intramammary infections and/or injuries, which as shown in the present study, are not unfavourably associated with milk production traits. In addition, such antagonistic traits can be combined in a suitably constructed phenotypic selection index towards overall improvement. This index would be based on adjusted animal phenotypes, each weighed according to relative economic importance, between-animal trait variance, and co-variance with the other traits of interest [58,59,60,61], and applied to inform culling decisions. Finally, results also suggest that arthritis and overall limb issues could be reduced in Damascus goats by improving BCS or maintaining it at desirable levels (score of 2.5 to 3). This finding is in accordance with previous research in dairy cows that reported favourable genetic correlations between foot health and BCS [51,62,63].
Future research aiming to estimate genetic correlations between the studied traits would provide results to underpin comprehensive selective breeding and genetic improvement programmes. However, this would require the collection of accurate pedigree data. As previously indicated, in low-input pastoral farming systems, pedigree recording is hindered by poor infrastructure, uncontrolled natural mating, and absence of artificial insemination [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]. Moreover, farmers are not well informed of the benefits of genetic improvement and reluctant to invest time and resources in record keeping [20]. However, there are possibilities to overcome some of these challenges through genomic evaluation and selection practices. Specifically, genotyping key animals with genome-wide DNA arrays could be a promising alternative to pedigree data availability in low-input dairy goat systems. Customised arrays for genomic markers identified to be associated with the studied traits could be developed for this purpose [20]. However, a complete feasibility and economic analysis would be needed to determine the optimal use of genomic technologies in low-input systems.
Finally, when practices to improve goat performance, health, and welfare are systematically implemented, it is advisable to periodically carry out welfare and performance assessments to assess improvements. Selection, even at the phenotypic level, could result in allelic frequency changes in the populations over time [64,65], which might impact on the estimation of (co)variance components. Moreover, ongoing changes in environmental conditions, including climate change, may also affect between-animal (co)variance estimates for many traits of interest [66].

5. Conclusions

Results of the present study indicate that there is significant between-animal variation for udder health, resistance to parasitism, and welfare traits to support management and selection practices aiming to improve the health and welfare status of Eghoria, Skopelos, and Damascus goats reared under low-input pastoral systems. Such practices could be implemented without compromising animal performance. Importantly, in most cases, subclinical mastitis, GIN infection, udder asymmetry, and skin lesions could be reduced when selecting for higher milk yield and quality. In some cases, arthritis and overall limb issues could be reduced by improving animal BCS. Few antagonistic trait associations revealed in our study, exemplified by the unfavourable correlation of certain welfare issues with milk production, could be addressed with appropriate management practices and possibly the implementation of a phenotypic selection index to underpin culling decisions. In all cases, accurate and systematic record keeping is essential for improving overall goat performance in low-input pastoral systems.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/vetsci9060289/s1, Dataset S1: Total dataset of goat traits used for the analyses; Table S1: Repeatability estimates (standard error in parentheses) for performance traits in Eghoria, Skopelos, and Damascus goats; all estimates were statistically greater than zero (p < 0.05); Table S2: Animal correlations (standard error in parentheses) of health and welfare traits with performance traits in Eghoria goats; Table S3: Animal correlations (standard error in parentheses) of health and welfare traits with performance traits in Skopelos goats; Table S4: Animal correlations (standard error in parentheses) of health and welfare traits with performance traits in Damascus goats.

Author Contributions

Conceptualisation, S.V., A.I.G., V.F., G.B. and G.A.; methodology, S.V., V.F., G.B. and G.A.; validation, S.V. and G.B.; formal analysis, S.V. and V.F.; investigation, A.I.G. and G.A.; resources, A.I.G. and G.A.; data curation, S.V.; writing—original draft preparation, S.V.; writing—review and editing, S.V., A.I.G., V.F., G.B. and G.A.; visualisation, S.V.; supervision, A.I.G., G.B. and G.A.; project administration, A.I.G., G.B. and G.A.; funding acquisition, G.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Union’s H2020 research and innovation program for the project SMAll RuminanTs breeding for Efficiency and Resilience (SMARTER), grant number 772787. Collection of data was funded by the European Commission Seventh Framework Programme (FP7 Food, Agriculture and Fisheries, Biotechnology) for the project Sustainable Organic and Low Input Dairying (SOLID), grant number FP7-KBBE-266367.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Research Committee of the Aristotle University of Thessaloniki (26362/03.05.2011) within the framework of SOLID project.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data presented in this study are contained within the article and supplementary material (Dataset S1 and Supplement 1).

Acknowledgments

Results presented in this work were produced using the AUTH Compute Infrastructure and Resources.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Kosgey, I.S.; Okeyo, A.M. Genetic improvement of small ruminants in low-input, smallholder production systems: Technical and infrastructural issues. Small Rumin. Res. 2007, 70, 76–88. [Google Scholar] [CrossRef]
  2. Barillet, F. Genetic improvement for dairy production in sheep and goats. Small Rumin. Res. 2007, 70, 60–75. [Google Scholar] [CrossRef]
  3. Simões, J.; Abecia, J.A.; Cannas, A.; Delgadillo, J.A.; Lacasta, D.; Voigt, K.; Chemineau, P. Managing sheep and goats for sustainable high yield production. Animal 2021, 15, 100293. [Google Scholar] [CrossRef] [PubMed]
  4. Sevi, A.; Casamassima, D.; Pulina, G.; Pazzona, A. Factors of welfare reduction in dairy sheep and goats. Ital. J. Anim. Sci. 2009, 8, 81–101. [Google Scholar] [CrossRef]
  5. Gelasakis, A.I.; Valergakis, G.E.; Arsenos, G. Health and welfare of indigenous goat breeds from dairy farms in Greece. In Sustainable Goat Production in Adverse Environments: Volume I; Simões, J., Gutiérrez, C., Eds.; Springer: Cham, Switzerland, 2017; pp. 223–246. [Google Scholar] [CrossRef]
  6. Vouraki, S.; Gelasakis, A.I.; Papanikolopoulou, V.; Papadopoulos, E.; Arsenos, G. Association of Hard Ticks (Ixodidae) Infestation with Milk Production and Udder Health of Extensively Reared Dairy Goats. Animals 2022, 12, 354. [Google Scholar] [CrossRef] [PubMed]
  7. Gelasakis, A.I.; Angelidis, A.S.; Giannakou, R.; Filioussis, G.; Kalamaki, M.S.; Arsenos, G. Bacterial subclinical mastitis and its effect on milk yield in low-input dairy goat herds. J. Dairy Sci. 2016, 99, 3698–3708. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Battini, M.; Vieira, A.; Barbieri, S.; Ajuda, I.; Stilwell, G.; Mattiello, S. Invited review: Animal-based indicators for on-farm welfare assessment for dairy goats. J. Dairy Sci. 2014, 97, 6625–6648. [Google Scholar] [CrossRef] [Green Version]
  9. Hoste, H.; Le Frileux, Y.; Pommaret, A. Distribution and repeatability of faecal egg counts and blood parameters in dairy goats naturally infected with gastrointestinal nematodes. Res. Vet. Sci. 2001, 70, 57–60. [Google Scholar] [CrossRef]
  10. Hoste, H.; Chartier, C.; Le Frileux, Y. Control of gastrointestinal parasitism with nematodes in dairy goats by treating the host category at risk. Vet. Res. 2002, 33, 531–545. [Google Scholar] [CrossRef] [Green Version]
  11. Heckendorn, F.; Bieber, A.; Werne, S.; Saratsis, A.; Maurer, V.; Stricker, C. The genetic basis for the selection of dairy goats with enhanced resistance to gastrointestinal nematodes. Parasite 2017, 24, 32. [Google Scholar] [CrossRef] [Green Version]
  12. Rostang, A.; Devos, J.; Chartier, C. Review of the Eprinomectin effective doses required for dairy goats: Where do we go from here? Vet. Parasitol. 2020, 277, 108992. [Google Scholar] [CrossRef] [PubMed]
  13. Rupp, R.; Huau, C.; Caillat, H.; Fassier, T.; Bouvier, F.; Pampouille, E.; Clément, V.; Palhière, I.; Larroque, H.; Tosser-Klopp, G.; et al. Divergent selection on milk somatic cell count in goats improves udder health and milk quality with no effect on nematode resistance. J. Dairy Sci. 2019, 102, 5242–5253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Apodaca-Sarabia, C.A.; Lopez-Villalobos, N.; Blair, H.T.; Prosser, G.G. Genetic parameters for somatic cell score in dairy goats estimated by random regression. Proc. N. Z. Soc. Anim. Prod. 2009, 69, 206–209. [Google Scholar]
  15. Rupp, R.; Clément, V.; Piacere, A.; Robert-Granié, C.; Manfredi, E. Genetic parameters for milk somatic cell score and relationship with production and udder type traits in dairy Alpine and Saanen primiparous goats. J. Dairy Sci. 2011, 94, 3629–3634. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Maroteau, C.; Palhière, I.; Larroque, H.; Clément, V.; Ferrand, M.; Tosser-Klopp, G.; Rupp, R. Genetic parameter estimation for major milk fatty acids in Alpine and Saanen primiparous goats. J. Dairy Sci. 2014, 97, 3142–3155. [Google Scholar] [CrossRef] [PubMed]
  17. Bagnicka, E.; Lukaszewicz, M.; Ådnøy, T. Genetic parameters of somatic cell score and lactose content in goat’s milk. J. Anim. Feed Sci. 2016, 25, 210–215. [Google Scholar] [CrossRef]
  18. Vagenas, D.; Jackson, F.; Russel, A.J.F.; Merchant, M.; Wright, I.A.; Bishop, S.C. Genetic control of resistance to gastro-intestinal parasites in crossbred cashmere-producing goats: Responses to selection, genetic parameters and relationships with production traits. Anim. Sci. 2002, 74, 199–208. [Google Scholar] [CrossRef] [Green Version]
  19. Mandonnet, N.; Menendez-Buxadera, A.; Arquet, R.; Mahieu, M.; Bachand, M.; Aumont, G. Genetic variability in resistance to gastro-intestinal strongyles during early lactation in Creole goats. Anim. Sci. 2006, 82, 283–287. [Google Scholar] [CrossRef]
  20. Zvinorova, P.I.; Halimani, T.E.; Muchadeyi, F.C.; Matika, O.; Riggio, V.; Dzama, K. Breeding for resistance to gastrointestinal nematodes–the potential in low-input/output small ruminant production systems. Vet. Parasitol. 2016, 225, 19–28. [Google Scholar] [CrossRef] [Green Version]
  21. Scholtens, M.R.; Lopez-Villalobos, N.; Garrick, D.; Blair, H.; Lehnert, K.; Snell, R. Genetic parameters for total lactation yields of milk, fat, protein, and somatic cell score in New Zealand dairy goats. Anim. Sci. J. 2020, 91, e13310. [Google Scholar] [CrossRef]
  22. Vlassoff, A.; Bisset, S.A.; McMurtry, L.W. Faecal egg counts in Angora goats following natural or experimental challenge with nematode parasites: Within-flock variability and repeatabilities. Vet. Parasitol. 1999, 84, 113–123. [Google Scholar] [CrossRef]
  23. Morris, C.A.; Wheeler, M.; Hosking, B.C.; Watson, T.G.; Hurford, A.P.; Foote, B.J.; Foote, J.F. Genetic parameters for milk yield and faecal nematode egg count in Saanen does. N. Z. J. Agric. Res. 1997, 40, 523–528. [Google Scholar] [CrossRef]
  24. Gelasakis, A.I.; Rose, G.; Giannakou, R.; Valergakis, G.E.; Theodoridis, A.; Fortomaris, P.; Arsenos, G. Typology and characteristics of dairy goat production systems in Greece. Livest. Sci. 2017, 197, 22–29. [Google Scholar] [CrossRef]
  25. Anzuino, K.; Bell, N.J.; Bazeley, K.J.; Nicol, C.J. Assessment of welfare on 24 commercial UK dairy goat farms based on direct observations. Vet. Rec. 2010, 167, 774–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Ley, S.J.; Waterman, A.E.; Livingston, A.; Parkinson, T.J. Effect of chronic pain associated with lameness on plasma cortisol concentrations in sheep: A field study. Res. Vet. Sci. 1994, 57, 332–335. [Google Scholar] [CrossRef]
  27. Ministry of Agriculture, Fisheries and Food (MAFF). Manual of Veterinary Parasitological Laboratory Techniques, 3rd ed.; H.M.S.O.: London, UK, 1986; ISBN 0112427243.
  28. Russel, A.J.F.; Doney, J.M.; Gunn, R.G. Subjective assessment of body fat in live sheep. J. Agric. Sci. 1969, 72, 451–454. [Google Scholar] [CrossRef]
  29. Gelasakis, A.I.; Angelidis, A.S.; Giannakou, R.; Arsenos, G. Bacterial subclinical mastitis and its effect on milk quality traits in low-input dairy goat herds. Vet. Rec. 2018, 183, 449. [Google Scholar] [CrossRef]
  30. Vouraki, S.; Gelasakis, A.I.; Ekateriniadou, L.V.; Banos, G.; Arsenos, G. Impact of polymorphisms at the PRNP locus on the performance of dairy goats reared under low-input pastoral farming systems. Small Rumin. Res. 2019, 174, 77–82. [Google Scholar] [CrossRef]
  31. International Committee for Animal Recording. ICAR Recording Guidelines. International Agreement of Recording Practices. Available online: http://www.icar.org/wp-content/uploads/2016/Guidelines-Edition-2016.pdf/ (accessed on 7 June 2018).
  32. Papadopoulos, E.; Arsenos, G.; Sotiraki, S.; Deligiannis, C.; Lainas, T.; Zygoyiannis, D. The epizootiology of gastrointestinal nematode parasites in Greek dairy breeds of sheep and goats. Small Rumin. Res. 2003, 47, 193–202. [Google Scholar] [CrossRef]
  33. Papadopoulos, B.; Morel, P.C.; Aeschlimann, A. Ticks of domestic animals in the Macedonia region of Greece. Vet. Parasitol. 1996, 63, 25–40. [Google Scholar] [CrossRef]
  34. Hadfield, J.D. MCMC methods for multi-response generalised linear mixed models: The MCMCglmm R package. J. Stat. Softw. 2010, 33, 1–22. [Google Scholar] [CrossRef] [Green Version]
  35. Gelman, A.; Rubin, D.B. Inference from iterative simulation using multiple sequences. Stat. Sci. 1992, 7, 457–472. [Google Scholar] [CrossRef]
  36. Boettcher, P.J.; Moroni, P.; Pisoni, G.; Gianola, D. Application of a finite mixture model to somatic cell scores of Italian goats. J. Dairy Sci. 2005, 88, 2209–2216. [Google Scholar] [CrossRef] [Green Version]
  37. Stuhr, T.; Aulrich, K. Intramammary infections in dairy goats: Recent knowledge and indicators for detection of subclinical mastitis. Landbauforschung 2010, 60, 267–279. [Google Scholar]
  38. Rainard, P.; Foucras, G.; Boichard, D.; Rupp, R. Invited review: Low milk somatic cell count and susceptibility to mastitis. J. Dairy Sci. 2018, 101, 6703–6714. [Google Scholar] [CrossRef] [Green Version]
  39. Brown, A.H., Jr.; Steelman, C.D.; Johnson, Z.B.; Rosenkrans, C.F., Jr.; Brasuell, T.M. Estimates of repeatability and heritability of horn fly resistance in beef cattle. J. Anim. Sci. 1992, 70, 1375–1381. [Google Scholar] [CrossRef] [Green Version]
  40. Burrow, H.M. Variances and covariances between productive and adaptive traits and temperament in a composite breed of tropical beef cattle. Livest. Prod. Sci. 2001, 70, 213–233. [Google Scholar] [CrossRef]
  41. Mackinnon, M.J.; Meyer, K.; Hetzel, D.J.S. Genetic variation and covariation for growth, parasite resistance and heat tolerance in tropical cattle. Livest. Prod. Sci. 1991, 27, 105–122. [Google Scholar] [CrossRef]
  42. Silva, A.M.D.; Alencar, M.M.D.; Regitano, L.C.D.A.; Oliveira, M.C.D.S.; Júnior, W.B. Artificial infestation of Boophilus microplus in beef cattle heifers of four genetic groups. Genet. Mol. Biol. 2007, 30, 1150–1155. [Google Scholar] [CrossRef]
  43. Cardoso, F.F.; Gomes, C.C.G.; Sollero, B.P.; Oliveira, M.M.; Roso, V.M.; Piccoli, M.L.; Higa, R.H.; Yokoo, M.J.; Caetano, A.R.; Aguilar, I. Genomic prediction for tick resistance in Braford and Hereford cattle. J. Anim. Sci. 2015, 93, 2693–2705. [Google Scholar] [CrossRef]
  44. Giglioti, R.; de Oliveira, H.N.; Bilhassi, T.B.; Portilho, A.I.; Okino, C.H.; Marcondes, C.R.; de Sena Oliveira, M.C. Estimates of repeatability and correlations of hemoparasites infection levels for cattle reared in endemic areas for Rhipicephalus microplus. Vet. Parasitol. 2018, 250, 78–84. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Grøva, L.; Sae-Lim, P.; Olesen, I. Genetic parameters of tick–infestation on lambs of the Norwegian White sheep breed. In Proceedings of the 10th World Congress on Genetics Applied to Livestock Production, Vancouver, BC, Canada, 17–22 August 2014. [Google Scholar]
  46. Cloete, S.W.P.; Cloete, J.J.E.; Scholtz, A.J. Genetic parameters for tick count and udder health in commercial and indigenous ewes in South Africa. Vet. Parasitol. 2016, 230, 33–42. [Google Scholar] [CrossRef] [PubMed]
  47. Giglioti, R.; Oliveira, H.N.; Santana, C.H.; Ibelli, A.M.G.; Néo, T.A.; Bilhassi, T.B.; Rabelo, M.D.; Machado, R.Z.; Brito, L.G.; Oliveira, M.D.S. Babesia bovis and Babesia bigemina infection levels estimated by qPCR in Angus cattle from an endemic area of São Paulo state, Brazil. Ticks Tick Borne Dis. 2016, 7, 657–662. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Kadarmideen, H.N.; Thompson, R.; Simm, G. Linear and threshold model genetic parameters for disease, fertility and milk production in dairy cattle. Anim. Sci. 2000, 71, 411–419. [Google Scholar] [CrossRef] [Green Version]
  49. Berry, D.P.; Bermingham, M.L.; Good, M.; More, S.J. Genetics of animal health and disease in cattle. Ir. Vet. J. 2011, 64, 5. [Google Scholar] [CrossRef] [Green Version]
  50. Weber, A.; Stamer, E.; Junge, W.; Thaller, G. Genetic parameters for lameness and claw and leg diseases in dairy cows. J. Dairy Sci. 2013, 96, 3310–3318. [Google Scholar] [CrossRef] [Green Version]
  51. Kougioumtzis, A.; Valergakis, G.E.; Oikonomou, G.; Arsenos, G.; Banos, G. Profile and genetic parameters of dairy cattle locomotion score and lameness across lactation. Animal 2014, 8, 20–27. [Google Scholar] [CrossRef]
  52. Mandal, A.; Sharma, D.K. Inheritance of faecal nematode egg count in Barbari goats following natural Haemonchus contortus infection. Vet. Parasitol. 2008, 155, 89–94. [Google Scholar] [CrossRef]
  53. O’Brien, A.C.; McHugh, N.; Wall, E.; Pabiou, T.; McDermott, K.; Randles, S.; Fair, S.; Berry, D.P. Genetic parameters for lameness, mastitis and dagginess in a multi-breed sheep population. Animal 2017, 11, 911–919. [Google Scholar] [CrossRef]
  54. Uribe, H.A.; Kennedy, B.W.; Martin, S.W.; Kelton, D.F. Genetic parameters for common health disorders of Holstein cows. J. Dairy Sci. 1995, 78, 421–430. [Google Scholar] [CrossRef]
  55. Doeschl-Wilson, A.B.; Vagenas, D.; Kyriazakis, I.; Bishop, S.C. Exploring the assumptions underlying genetic variation in host nematode resistance. Genet. Sel. Evol. 2008, 40, 241–264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Boareki, M.N.; Schenkel, F.S.; Willoughby, O.; Suarez-Vega, A.; Kennedy, D.; Cánovas, A. Comparison between methods for measuring fecal egg count and estimating genetic parameters for gastrointestinal parasite resistance traits in sheep. J. Anim. Sci. 2021, 99, skab341. [Google Scholar] [CrossRef] [PubMed]
  57. Oltenacu, P.A.; Broom, D.M. The impact of genetic selection for increased milk yield on the welfare of dairy cows. Anim. Welf. 2010, 19, 39–49. [Google Scholar]
  58. Hazel, L.N. The genetic basis for constructing selection indexes. Genetics 1943, 28, 476–490. [Google Scholar] [CrossRef] [PubMed]
  59. Hazel, L.N.; Dickerson, G.E.; Freeman, A.E. The selection index—Then, now, and for the future. J. Dairy Sci. 1994, 77, 3236–3251. [Google Scholar] [CrossRef]
  60. Henderson, C.R. Selection index and expected genetic advance. In Statistical Genetics and Plant Breeding, National Academy of Sciences; National Research Council Publication: Washington, DC, USA, 1963; pp. 141–163. [Google Scholar]
  61. Céron-Rojas, J.J.; Crossa, J. The linear phenotypic selection index theory. In Linear Selection Indices in Modern Plant Breeding; Springer: Cham, Switzerland, 2018; pp. 15–42. [Google Scholar] [CrossRef] [Green Version]
  62. Boettcher, P.J.; Dekkers, J.C.M.; Warnick, L.D.; Wells, S.J. Genetic analysis of clinical lameness in dairy cattle. J. Dairy Sci. 1998, 81, 1148–1156. [Google Scholar] [CrossRef]
  63. Van Dorp, T.E.; Dekkers, J.C.M.; Martin, S.W.; Noordhuizen, J.P.T.M. Genetic parameters of health disorders and relationships with 305-day milk yield and conformation traits of registered Holstein cows. J. Dairy Sci. 1998, 81, 2264–2270. [Google Scholar] [CrossRef]
  64. Alatalo, R.V.; Gustafsson, L.; Lundberg, A. Phenotypic selection on heritable size traits: Environmental variance and genetic response. Am. Nat. 1990, 135, 464–471. [Google Scholar] [CrossRef]
  65. Kruuk, L.E.; Merilä, J.; Sheldon, B.C. Phenotypic selection on a heritable size trait revisited. Am. Nat. 2001, 158, 557–571. [Google Scholar] [CrossRef]
  66. Killen, S.S.; Adriaenssens, B.; Marras, S.; Claireaux, G.; Cooke, S.J. Context dependency of trait repeatability and its relevance for management and conservation of fish populations. Conserv. Physiol. 2016, 4, cow007. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Map of Greece illustrating the regions in which the studied goat herds were located.
Figure 1. Map of Greece illustrating the regions in which the studied goat herds were located.
Vetsci 09 00289 g001
Table 1. Descriptive statistics of milk production, body condition score, udder health (SCC, TVC, SMI) and parasitism (faecal egg and larval counts) traits in the studied goat breeds.
Table 1. Descriptive statistics of milk production, body condition score, udder health (SCC, TVC, SMI) and parasitism (faecal egg and larval counts) traits in the studied goat breeds.
EghoriaSkopelosDamascus
TraitN 1Mean (±SD 2)N 1Mean (±SD 2)N 1Mean (±SD 2)
Daily milk yield (g)2597829.32 (379.64)27661359.53 (738.40)21241827.88 (1053.49)
Daily fat yield (g)256139.46 (17.67)268563.06 (31.85)211076.49 (40.77)
Daily protein yield (g)256330.82 (13.59)269350.39 (26.69)209565.00 (33.32)
Daily lactose yield (g)256236.91 (17.96)268859.77 (32.55)208978.29 (45.47)
Daily SNF 3 yield (g)256775.03 (34.65)2695122.40 (65.45)2111158.90 (87.48)
BCS 4 (1–5)29202.31 (0.37)28302.54 (0.32)22982.44 (0.41)
SCC 5 (×103 cells/mL)26732050.43 (3590.82)21671999.01 (1990.09)20012816.54 (4125.64)
TVC 6 (×103 cfu/mL)2858157.26 (415.80)265797.6 (319.21)1969201.00 (463.86)
SMI 72667−0.01 (0.91)20870.00 (0.93)1734−0.04 (0.87)
GIN 8 FEC 9 (eggs/g)584267.29 (426.67)5468.52 (28.18)81493.30 (360.23)
Cestode FEC 9 (eggs/g)5841.80 (16.07)5465.59 (32.92)8140.98 (16.97)
Lungworm FLC 10 (larvae/g)5846.51 (33.01)54510.18 (59.19)8140.06 (1.75)
1 N = number of records; 2 SD = standard deviation; 3 SNF = solids-non-fat; 4 BCS = body condition score (1 = emaciated, 5 = obese); 5 SCC = somatic cell count; 6 TVC = total viable count; 7 SMI = subclinical mastitis index; 8 GIN = gastrointestinal nematodes; 9 FEC = faecal egg count; 10 FLC = faecal larval count.
Table 2. Frequency (%) of udder health phenotypes derived from milk somatic cell count and total viable count in the studied goat breeds.
Table 2. Frequency (%) of udder health phenotypes derived from milk somatic cell count and total viable count in the studied goat breeds.
TraitLevelsEghoriaSkopelosDamascus
UHP1 1073.6876.9566.49
126.3223.0533.51
UHP2 2050.3654.8640.60
123.3222.0925.89
226.3223.0533.51
UHP3 3050.3654.8640.60
18.776.376.46
214.5515.7219.43
326.3223.0533.51
1 UHP1 = scored as 0 if somatic cell count ≤106 cells/mL and/or total viable count ≤2 × 104 cfu/mL, or 1 if somatic cell count >106 cells/mL and total viable count >20 × 103 cfu/mL; 2 UHP2 = scored as 0 if somatic cell count ≤106 cells/mL and total viable count ≤2 × 104 cfu/mL, 1 if somatic cell count >106 cells/mL or total viable count >2 × 104 cfu/mL, respectively, or 2 if somatic cell count >106 cells/mL and total viable count >2 × 104 cfu/mL; 3 UHP3 = scored as 0 if somatic cell count ≤106 cells/mL and total viable count ≤2 × 104 cfu/mL, 1 if only total viable count >2 × 104 cfu/mL, 2 if only somatic cell count ≤106 cells/mL or 3 if somatic cell count >106 cells/mL and total viable count >2 × 104 cfu/mL.
Table 3. Frequency (%) of parasitism in the studied goat breeds.
Table 3. Frequency (%) of parasitism in the studied goat breeds.
EghoriaSkopelosDamascus
TraitN 1Frequency (%)N 1Frequency (%)N 1Frequency (%)
Tick infestation292113.6328308.8022980.26
Myiasis29210.2428304.5522980.00
GIN 2 infection58466.4454611.7281427.40
Cestode infection5842.055464.768140.61
Lungworm infection5845.145457.348140.12
Endoparasite infection58470.0354519.2381427.64
1 N = number of records; 2 GIN = gastrointestinal nematodes.
Table 4. Frequency (%) of welfare issues in the studied goat breeds.
Table 4. Frequency (%) of welfare issues in the studied goat breeds.
EghoriaSkopelosDamascus
Body PartTraitN 1Frequency (%)N 1Frequency (%)N 1Frequency (%)
HeadEar injuries29208.2928301.8022984.09
Horn injuries292011.2328308.90229818.28
Head skin lesions29207.6028302.3722988.66
Nasal discharge292010.1428300.8822982.39
Ocular discharge29200.4128300.0722980.57
Head problems292031.43283013.36229831.38
BodyAbscess29211.51283013.3922989.88
Diarrhoea29200.5128300.0422980.52
Injury29210.4528300.2522980.35
Hernia29210.0028300.2822980.00
Body problems29213.15283014.42229810.92
LegsLameness29210.1028300.0422981.39
Overgrown claws29204.3528304.38229837.21
Arthritis29200.0328300.1122984.79
Leg problems29214.4528304.52229839.21
UdderUdder asymmetry292038.18283020.99229832.46
Udder abscess292028.08283025.51229829.77
Udder skin lesions29201.5128301.1022982.05
Udder problems292054.25283039.01229851.61
TotalSkin lesions29208.8028303.43229810.40
Injuries292019.25283010.71229821.76
1 N = number of records.
Table 5. Repeatability estimates (standard error in parentheses) for health (udder health and parasitism) traits in Eghoria, Skopelos, and Damascus goats.
Table 5. Repeatability estimates (standard error in parentheses) for health (udder health and parasitism) traits in Eghoria, Skopelos, and Damascus goats.
TraitEghoriaSkopelosDamascus
SCC 1 (×103 cells/mL)0.45 (0.02) *0.47 (0.02) *0.22 (0.02) *
TVC 2 (×103 cfu/mL)0.22 (0.02) *0.30 (0.02) *0.15 (0.02) *
SMI 30.20 (0.02) *0.28 (0.03) *0.10 (0.02) *
UHP1 4 (0–1)0.50 (0.04) *0.59 (0.04) *0.43 (0.05) *
UHP2 5 (0–2)0.08 (0.01) *0.09 (0.02) *0.08 (0.01) *
UHP3 6 (0–3)0.14 (0.02) *0.15 (0.02) *0.14 (0.02) *
Tick infestation (0–1)0.03 (0.02)0.03 (0.01)NE 10
Myiasis (0–1)NE 100.34 (0.10) *NE 10
GIN 7 FEC 8 (eggs/g, Tukey)0.11 (0.04) *0.001 (0.002)0.001 (0.002)
Cestode FEC 8 (eggs/g, Tukey)0.09 (0.03) *0.04 (0.03)NE 10
Lungworm FLC 9 (larvae/g, Tukey)0.04 (0.02)0.001 (0.002)NE 10
GIN 7 infection (0–1)0.26 (0.07) *0.06 (0.04)0.01 (0.02)
Cestode infection (0–1)0.23 (0.19)0.27 (0.15)NE 10
Lungworm infection (0–1)0.32 (0.11) *0.03 (0.03)NE 10
Endoparasite infection (0–1)0.29 (0.07) *0.04 (0.03)0.002 (0.005)
1 SCC = somatic cell count; 2 TVC = total viable count; 3 SMI = subclinical mastitis index; 4 UHP1 = scored as 0 if somatic cell count ≤106 cells/mL and/or total viable count ≤2 × 104 cfu/mL, or 1 if somatic cell count >106 cells/mL and total viable count >2 × 104 cfu/mL; 5 UHP2 = scored as 0 if somatic cell count ≤106 cells/mL and total viable count ≤2 × 104 cfu/mL, 1 if somatic cell count >106 cells/mL or total viable count >2 × 104 cfu/mL, respectively, or 2 if somatic cell count >106 cells/mL and total viable count >2 × 104 cfu/mL; 6 UHP3 = scored as 0 if somatic cell count ≤106 cells/mL and total viable count ≤2 × 104 cfu/mL, 1 if only total viable count >2 × 104 cfu/mL, 2 if only somatic cell count ≤106 cells/mL or 3 if somatic cell count >106 cells/mL and total viable count >2 × 104 cfu/mL; 7 GIN = gastrointestinal nematodes; 8 FEC = faecal egg count; 9 FLC = faecal larval count; 10 NE = not estimable. * Indicates statistically significant repeatability estimates (p < 0.05).
Table 6. Repeatability estimates (standard error in parentheses) for welfare traits in Eghoria, Skopelos, and Damascus goats.
Table 6. Repeatability estimates (standard error in parentheses) for welfare traits in Eghoria, Skopelos, and Damascus goats.
TraitEghoriaSkopelosDamascus
Ear injuries (0–1)0.95 (0.01) *0.96 (0.02) *0.91 (0.03) *
Horn injuries (0–1)0.99 (0.004) *0.97 (0.01) *0.98 (0.01) *
Head skin lesions (0–1)0.18 (0.04) *0.23 (0.07) *0.43 (0.06) *
Nasal discharge (0–1)0.21 (0.04) *NE 10.25 (0.07) *
Ocular discharge (0–1)NENE 1NE 1
Head problems (0–1)0.61 (0.03) *0.78 (0.03) *0.74 (0.03) *
Abscess (0–1)0.41 (0.07) *0.49 (0.04) *0.46 (0.05) *
Diarrhoea (0–1)NE 1NE 1NE 1
Injury (0–1)NE 1NE 1NE 1
Hernia (0–1)NE 1NE 1NE 1
Body problems (0–1)0.22 (0.05) *0.45 (0.04) *0.42 (0.05) *
Lameness (0–1)NE 1NE 10.51 (0.08) *
Overgrown claws (0–1)0.17 (0.05) *0.46 (0.07) *0.21 (0.04) *
Arthritis (0–1)NE 1NE 10.87 (0.04) *
Limb problems (0–1)0.16 (0.05) *0.45 (0.07) *0.20 (0.04) *
Udder asymmetry (0–1)0.52 (0.03) *0.52 (0.04) *0.63 (0.03) *
Udder abscess (0–1)0.37 (0.03) *0.50 (0.03) *0.39 (0.04) *
Udder skin lesions (0–1)0.43 (0.08) *0.54 (0.08) *0.29 (0.08) *
Udder problems (0–1)0.36 (0.03) *0.41 (0.03) *0.46 (0.03) *
Total skin lesions (0–1)0.16 (0.04) *0.19 (0.05) *0.35 (0.05) *
Total injuries (0–1)0.95 (0.01) *0.95 (0.01) *0.93 (0.01) *
1 NE = not estimable. * Indicates statistically significant repeatability estimates (p < 0.05).
Table 7. Statistically significant (p < 0.05) animal correlations (standard error in parentheses) between health and milk production traits in Eghoria, Skopelos, and Damascus goats.
Table 7. Statistically significant (p < 0.05) animal correlations (standard error in parentheses) between health and milk production traits in Eghoria, Skopelos, and Damascus goats.
Milk Production Traits
BreedHealth TraitsMilk Yield
(g, ln)
Fat Yield
(g, ln)
Protein Yield
(g, ln)
Lactose Yield
(g, ln)
SNF Yield
(g, ln)
EghoriaSCC 1 (cells/ mL, ln)−0.15 (0.06)−0.22 (0.06)NS 8−0.20 (0.06)−0.15 (0.07)
TVC 2 (cfu/ mL, ln)−0.14 (0.07)−0.20 (0.07)NS 8−0.19 (0.07)−0.12 (0.07)
SMI 3−0.21 (0.07)−0.26 (0.08)NS 8−0.26 (0.07)−0.19 (0.07)
UHP1 4 (0–1)−0.21 (0.07)−0.26 (0.08)NS 8−0.28 (0.07)−0.20 (0.07)
UHP2 5 (0–2)−0.22 (0.08)−0.31 (0.09)NS 8−0.30 (0.08)−0.21 (0.09)
UHP3 6 (0–3)−0.18 (0.07)−0.27 (0.07)NS 8−0.25 (0.07)−0.18 (0.07)
GIN 7 infection (0–1)−0.39 (0.17)NS 8NS 8−0.39 (0.18)−0.37 (0.17)
Endoparasite infection (0–1)−0.39 (0.17)NS 8NS 8−0.39 (0.17)−0.36 (0.18)
SkopelosSCC 1 (cells/mL, ln)−0.22 (0.06)−0.29 (0.06)−0.18 (0.06)−0.27 (0.06)−0.23 (0.06)
TVC 2 (cfu/mL, ln)−0.23 (0.06)−0.30 (0.07)−0.18 (0.07)−0.28 (0.06)−0.24 (0.06)
SMI 3−0.24 (0.07)−0.33 (0.07)−0.20 (0.07)−0.29 (0.07)−0.25 (0.07)
UHP1 4 (0–1)−0.30 (0.07)−0.39 (0.07)−0.24 (0.07)−0.35 (0.07)−0.31 (0.07)
UHP2 5 (0–2)−0.31 (0.09)−0.39 (0.09)−0.22 (0.09)−0.37 (0.09)−0.31 (0.09)
UHP3 6 (0–3)−0.22 (0.07)−0.30 (0.07)−0.15 (0.07)−0.28 (0.07)−0.22 (0.07)
DamascusSCC 1 (cells/mL, ln)−0.25 (0.09)−0.24 (0.09)NS 8−0.27 (0.09)−0.20 (0.09)
UHP1 4 (0–1)−0.25 (0.09)−0.21 (0.10)NS 8−0.33 (0.10)NS 8
UHP2 5 (0–2)−0.23 (0.09)−0.34 (0.09)NS 8−0.33 (0.09)−0.22 (0.09)
UHP3 6 (0–3)−0.20 (0.08)−0.30 (0.08)NS 8−0.27 (0.08)−0.19 (0.08)
1 SCC = somatic cell count; 2 TVC = total viable count; 3 SMI = subclinical mastitis index; 4 UHP1 = scored as 0 if somatic cell count ≤106 cells/mL and/or total viable count ≤2 × 104 cfu/mL, or 1 if somatic cell count >106 cells/mL and total viable count >2 × 104 cfu/mL; 5 UHP2 = scored as 0 if somatic cell count ≤106 cells/mL and total viable count ≤2 × 104 cfu/mL, 1 if somatic cell count >106 cells/mL or total viable count >2 × 104 cfu/mL, respectively, or 2 if somatic cell count >106 cells/mL and total viable count >2 × 104 cfu/mL; 6 UHP3 = scored as 0 if somatic cell count ≤106 cells/mL and total viable count ≤2 × 104 cfu/mL, 1 if only total viable count >2 × 104 cfu/mL, 2 if only somatic cell count ≤106 cells/mL or 3 if somatic cell count >106 cells/mL and total viable count >2 × 104 cfu/mL; 7 GIN = gastrointestinal nematodes; 8 NS = non-significant.
Table 8. Statistically significant (p < 0.05) animal correlations (standard error in parentheses) between welfare and performance traits in Eghoria, Skopelos, and Damascus goats.
Table 8. Statistically significant (p < 0.05) animal correlations (standard error in parentheses) between welfare and performance traits in Eghoria, Skopelos, and Damascus goats.
Performance Traits
BreedHealth TraitsMilk Yield
(g, ln)
Fat Yield
(g, ln)
Protein Yield
(g, ln)
Lactose Yield
(g, ln)
SNF Yield
(g, ln)
BCS 1
(1–5)
EghoriaUdder abscess (0–1)NS 2NS 20.17 (0.07)NS 2NS 2NS 2
Head skin lesions (0–1)NS 2NS 2−0.28 (0.12)NS 2NS 2NS 2
Total skin lesions (0–1)NS 2−0.27 (0.13)NS 2NS 2NS 2NS 2
SkopelosUdder abscess (0–1)0.23 (0.06)NS 20.24 (0.07)0.18 (0.07)0.21 (0.07)NS 2
Udder problems (0–1)NS 2NS 20.15 (0.07)NS 2NS 2NS 2
Overgrown claws (0–1)0.36 (0.10)0.44 (0.10)0.38 (0.10)0.37 (0.09)0.40 (0.10)NS 2
Limb problems (0–1)0.34 (0.09)0.45 (0.10)0.37 (0.10)0.36 (0.10)0.39 (0.10)NS 2
DamascusUdder asymmetry (0–1)−0.24 (0.08)−0.24 (0.08)−0.20 (0.08)−0.23 (0.08)−0.23 (0.08)NS 2
Arthritis (0–1)NS 2NS 2NS 2NS 2NS 2−0.22 (0.09)
Limb problems (0–1)NS 2NS 2NS 2NS 2NS 2−0.28 (0.10)
1 BCS = body condition score (1 = emaciated, 5 = obese); 2 NS = non-significant.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vouraki, S.; Gelasakis, A.I.; Fotiadou, V.; Banos, G.; Arsenos, G. Repeatability of Health and Welfare Traits and Correlation with Performance Traits in Dairy Goats Reared under Low-Input Farming Systems. Vet. Sci. 2022, 9, 289. https://0-doi-org.brum.beds.ac.uk/10.3390/vetsci9060289

AMA Style

Vouraki S, Gelasakis AI, Fotiadou V, Banos G, Arsenos G. Repeatability of Health and Welfare Traits and Correlation with Performance Traits in Dairy Goats Reared under Low-Input Farming Systems. Veterinary Sciences. 2022; 9(6):289. https://0-doi-org.brum.beds.ac.uk/10.3390/vetsci9060289

Chicago/Turabian Style

Vouraki, Sotiria, Athanasios I. Gelasakis, Vasileia Fotiadou, Georgios Banos, and Georgios Arsenos. 2022. "Repeatability of Health and Welfare Traits and Correlation with Performance Traits in Dairy Goats Reared under Low-Input Farming Systems" Veterinary Sciences 9, no. 6: 289. https://0-doi-org.brum.beds.ac.uk/10.3390/vetsci9060289

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop