Next Article in Journal
Neutron-Induced Nuclear Cross-Sections Study for Plasma Facing Materials via Machine Learning: Molybdenum Isotopes
Next Article in Special Issue
Characterization of Lactic Acid Bacteria in Pecorino di Farindola Cheese and Manufacturing with a Lacticaseibacillus paracasei Autochthonous Culture
Previous Article in Journal
TickPhone App: A Smartphone Application for Rapid Tick Identification Using Deep Learning
Previous Article in Special Issue
White Brined Cheese Production by Incorporation of a Traditional Milk-Cereal Prebiotic Matrix with a Candidate Probiotic Bacterial Strain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Association of Breed of Sheep or Goats with Somatic Cell Counts and Total Bacterial Counts of Bulk-Tank Milk

by
Daphne T. Lianou
1,
Charalambia K. Michael
1,
Natalia G. C. Vasileiou
2,
Dimitra V. Liagka
1,
Vasia S. Mavrogianni
1,
Mariangela Caroprese
3 and
George C. Fthenakis
1,*
1
Veterinary Faculty, University of Thessaly, 43100 Karditsa, Greece
2
Faculty of Animal Science, University of Thessaly, 41110 Larissa, Greece
3
Department of Agriculture, Food, Natural Resources and Engineering (DAFNE), University of Foggia, 71122 Foggia, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(16), 7356; https://doi.org/10.3390/app11167356
Submission received: 22 June 2021 / Revised: 6 August 2021 / Accepted: 9 August 2021 / Published: 10 August 2021
(This article belongs to the Special Issue Food Microbiology: Dairy Products Microbiota)

Abstract

:
The objective was to describe potential associations of somatic cell counts (SCC) and total bacterial counts (TBC) in bulk-tank milk from sheep and goat farms with breeds of these animals in Greece. In total, 325 dairy sheep flocks and 119 dairy goat herds were visited for the collection of milk; the breed of animals in farms was evaluated for a potential association with SCC or TBC. The most frequently seen sheep breeds were the Lacaune (95 flocks) and the Chios (44 flocks). The most frequently seen goat breeds were the indigenous Greek (Capra prisca) (50 herds) and the Murciano-Granadina (13 herds). In a multivariable analysis, the breed and the application of machine-milking in sheep flocks, and the breed and the management system in goat herds emerged as significant factors for increased SCC (>0.75 × 106 cells mL−1) in bulk-tank milk. Further, the month of lactation at sampling in sheep flocks emerged a significant factor for increased TBC (>1500 × 103 cfu mL−1) in bulk-tank milk.

1. Introduction

There are over 6,000,000 sheep and 3,000,000 goats in Greece. The predominant type of production is dairy, with over 95% of animals in the country being milked. Many reports are available, presenting morphological characteristics of sheep and goat breeds in Greece (e.g., for sheep: [1,2,3,4,5], for goats: [6,7]).
A genetic background to mastitis in small ruminants has been identified [8,9,10,11]. Indeed, in the study of Davies et al. [12], mastitis was considered the top ranking sheep disease amenable for genetic studies. An assessment of mastitis resistance in genetic studies is usually performed by using somatic cell counts in milk; increased somatic cell counts can be employed as an indicator of infection and therefore of subclinical mastitis [13]. Depending on the breed, the heritabilities for somatic cell counts range from 0.09 to 0.30, whilst those for production traits vary between 0.30 and 0.60 [13].
In healthy mammary glands of ewes and does, there are epithelial cells and leucocytes, which predominate. The epithelial cells originate from the mammary gland itself, whilst the origin of leucocytes is the blood. Leucocyte types in the milk of healthy animals include macrophages (40–85%), neutrophils (5–35%) and lymphocytes (10–20%) [14,15,16]. Leucocytes are part of the innate defences of the mammary gland and play a significant protective role against invading bacteria. Macrophages are the first among the leucocytes present in the mammary gland to counteract the invading pathogens and to initiate the defensive leucocytic response [17]. As part of the defence response of the animal during infection, leucocytes (initially neutrophils) enter into the infected mammary gland, which leads to an increase in the numbers of cell counts in the milk. Therefore, an increased number of somatic cells in the milk of individual animals reflects the presence of an inflammatory process in the mammary gland, i.e., mastitis. With regard to the bulk-tank milk, somatic cell counts are employed to indicate the existence of mastitis in the animals of the farm. In sheep, according to Berthelot et al. [18], somatic cell counts in the bulk-tank milk indicate a mastitis prevalence of 15% in the flock. In the European Union, the current legislation does not provide an officially acceptable threshold for somatic cell counts in the raw milk of sheep and goats. In the United States of America, there is relevant legislation only for goat milk [19], allowing up to 0.75 × 106 cells mL−1 in the bulk-tank milk.
Total bacterial counts in the bulk-tank milk express the amount of microbial contamination therein. Bacteria in bulk-tank milk can originate from within the mammary gland, from the udder and teat skin, and from the surface of equipment for milk handling and storage. Hence, total bacterial counts reflect the possible presence of mastitis, udder cleaning procedures during the milking routine and environmental contamination in the farm environment (particularly in the milking parlour). The European Union legislation has set a threshold for total bacterial counts in the bulk-tank milk of small ruminants equal to 1,500,000 colony-forming units (cfu) mL−1 for milk that would undergo thermal processing and 500,000 cfu mL−1 for milk that would be used for direct consumption [20].
The objective of the work presented herein is to describe the potential associations of somatic cell counts (SCC) and total bacterial counts (TBC) in the bulk-tank milk from sheep and goat farms with breeds of these animals in Greece, with special reference to the potential role of animal breeds.

2. Materials and Methods

2.1. Sheep and Goat Farms, Sampling and Laboratory Examinations

An extensive cross-sectional study was performed. In total, 325 dairy sheep flocks and 119 dairy goat herds in all 13 administrative regions of Greece (Figure 1) were included into the study and visited for the collection of samples and information. Further details regarding these farms are presented in Table S1. Farms were selected by collaborating veterinarians on a convenience basis (willingness of farmers to accept a visit by university personnel for an interview and sample collection). The principal investigators (D.T.L. and G.C.F.) visited all of the farms in the study for sample collection.
During each visit, we collected four 20 mL samples from the bulk-tank milk using appropriate techniques for sampling (two samples were for cell counting, and two samples were for bacterial counting). Then, we stored these samples at 0.0 to 4.0 °C using ice packs in portable refrigerators. We performed somatic cell counting on each of the two samples within 4 h after collection, and we used the remaining two samples for bacterial counting, which was performed within 24 h after the collection of samples. At this stage, two sub-samples were created and processed from each of the four samples, so that we performed each separate test four times (each one in different sub-samples). The transportation of samples to the laboratory was by the investigators and by car; the samples collected from farms in the islands were also transported as ice-packed accompanying luggage by airplane or boat.
The somatic cell counts (SCC) (Lactoscan SCC; Milkotronic Ltd., Nova Zagora, Bulgaria) were measured on each of the four relevant sub-samples. For total bacterial counts (TBC) in the milk samples, the procedures detailed by Laird et al. [21] were followed and the total bacterial count in the initial sample was calculated. After the completion of sample processing for TBC, the temperature of the respective samples was measured and in no case was found to exceed 3.8 °C.

2.2. Data Management and Analysis

During cell counting and bacterial counting, for each bulk-tank milk sample, the results of the two sub-samples from each sample were averaged and, then, the two means were again averaged for the final result regarding each bulk-tank milk.
The data were entered into Microsoft Excel and analysed using SPSS v. 21 (IBM Analytics, Armonk, NY, USA). A basic descriptive analysis was performed. Exact binomial confidence intervals (CI) were obtained. For all statistical analyses and according to the standards employed for analyses of SCC, these were transformed to somatic cell scores (SCS) by using the equation SCS = log2(SCC/100) + 3 [22,23]. TBC were transformed to log10, and the transformed data were used in the analyses. At the end, for presentation of the results, the transformed findings were back-transformed into 100 × 2(SCS−3) and 10log data.
Initially, SCC and TBC were compared between farms in accordance with the breed of the animals therein, by using an analysis of variance.
Then, the outcomes for “increased SCC in bulk-tank milk” (i.e., with SCC over 0.75 × 106 cells mL−1) and “increased TBC in bulk-tank milk” (i.e., with TBC over 1500 × 103 cfu mL−1) were studied. The following variables were assessed for the potential association with SCC or TBC in the bulk-tank milk of the farms: breed, management system applied in the farm, application of machine-milking in the farm, month of lactation period at sampling, number of animals in the farms, average per animal yearly milk production in the farm and application of teat disinfection at milking. For each of these variables, categories were created. Exact binomial confidence intervals were obtained. Initially, the importance of predictors was assessed using cross-tabulation with Pearson’s chi-square test univariable analysis without random effects. Subsequently, multivariable models were created using mixed-effects logistic regression with farms as the random effect, as defined by Ŷ = b0 + b1X1 + b2X2 + … + biXi (Y: dependent variable, X: independent variables) and initially offering to the model all variables, which achieved a significance of p < 0.2 in the univariable analysis. The variables were removed from the initial model by backwards elimination. The p value of removal of a variable was assessed by the likelihood ratio test, and for those with a p value of >0.2, the variable with the largest probability was removed. This process was repeated until no variable could be removed with a p value of >0.2. The variables required for the final multivariable model on each occasion are shown in Table S2.
Separate analyses were performed for sheep flocks and goat herds. Initially, the analyses for the outcomes “increased SCC in bulk-tank milk” and “increased TBC in bulk-tank milk” were performed by taking into account all breeds, independent of the number of farms with each breed. Then, the analyses were repeated by taking into account only farms with the six most frequently recorded sheep breeds and the four most frequently recorded goat breeds.
For the evaluation of the zootechnical characteristics, Pearson’s chi-square test or analysis of variance was used as appropriate.
In all analyses, statistical significance was defined at p < 0.05.

3. Results

In sheep farms, the breed most frequently seen was the Lacaune breed (95 flocks), followed by the Chios breed (44 flocks). In goat farms, the breed most frequently seen was the indigenous Greek breed (Capra prisca) (50 herds), followed by the Murciano-Granadina breed (13 herds). With regard to sheep, breeds of Greek origin (including indigenous and crossbred animals) were present in 186 flocks (57.2%) and imported breeds were present in 139 flocks (42.8%). With regard to goats, breeds of Greek origin (including indigenous and crossbred animals) were present in 74 herds (62.2%) and imported breeds were present in 45 herds (37.8%). The various breeds of sheep and goats present in the farms into the study are shown in Table 1; their zootechnical characteristics are summarized according to the origin of the breeds (i.e., Greek or imported) in Table 2.
In farms with breeds of Greek origin, for sheep and goats, semi-extensive management was applied more frequently (40.9% and 59.5% of farms, respectively). In contrast, in farms with imported breeds, semi-intensive management was applied more frequently (44.6% and 44.4% of farms, respectively) (p < 0.0001). In the former farms, average yearly milk production was significantly lower than in the latter: 181 ± 6 L (mean ± standard error of the mean) and 155 ± 9 L per animal for sheep and goat farms, respectively, versus 243 ± 8 and 276 ± 20 L per animal, respectively (p < 0.0001). There was a significant progressive reduction in average yearly milk production from intensive to extensive for both sheep (p < 0.0001) and goats (p = 0.0002) (Figure 2).

3.1. Effect of Breed in Sheep Flocks

The geometric mean SCC and TBC according to animal breed in the sheep flocks of the study are shown in Table 1.

3.1.1. Effect on Somatic Cell Counts

When flocks with all breeds were considered, no significant difference was seen in SCC between the breeds (p = 0.22). However, when only flocks with the six most frequently recorded breeds were taken into account, a clear difference in SCC was seen between the breeds (p = 0.004) (Table 1). There was no significant difference in SCC between flocks with breeds of Greek origin or imported ones (0.501 × 106 versus 0.473 × 106 cells mL−1, respectively; p = 0.24).
For increased SCC (>0.75 × 106 cells mL−1) in bulk-tank milk, when flocks with all breeds were considered, a significant association was seen during the univariable analysis only for the application of machine-milking in the farms (p = 0.001). When only flocks with the six most frequently recorded breeds were taken into account, a significant association was seen during the univariable analysis for the application of machine-milking in the farms (p = 0.006) and the number of animals in the flock (p = 0.033) (Table S3). Among the variables included in the multivariable analysis (Table S2), the following emerged as significant factors for increased SCC in the flocks in both analyses: the breed (p = 0.025 and 0.014, respectively) (Figure 3) and the application of machine-milking in the farms (p = 0.0004 and 0.003, respectively) (Table 3).

3.1.2. Effect on Total Bacterial Counts

No significant difference was seen in TBC between the breeds, either when flocks with all breeds were considered or when only flocks with the six most frequently recorded breeds were taken into account (p > 0.78) (Table 1). There was no significant difference in TBC between flocks with breeds of Greek origin or imported ones (451 × 103 versus 336 × 103 cfu mL−1, respectively; p = 0.11).
For increased TBC (>1500 × 103 cfu mL−1) in bulk-tank milk, when flocks with all breeds were considered, a significant association was seen during the univariable analysis for the application of machine-milking in the farms (p = 0.032) and the month of lactation at sampling (p = 0.006). When only flocks with the six most frequently recorded breeds were taken into account, a significant association was seen during the univariable analysis for the management system applied in flocks (p = 0.045) (Table S4). Among the variables included in the multivariable analysis (Table S2), the following emerged as significant factors for increased TBC in the flocks: the month of lactation at sampling for both analyses (p = 0.0008 and 0.014, respectively) and the average per animal yearly milk production in the farm (p = 0.037, only when all breeds were considered (Table 4).

3.2. Effect of Breed in Goat Herds

The geometric mean SCC and TBC according to animal breed in the goat herds of the study are shown in Table 1.

3.2.1. Effect on Somatic Cell Counts

When herds with all breeds were considered, no significant difference was seen in SCC between the breeds (p = 0.23). Similar findings were seen when only herds with the four most frequently recorded breeds were taken into account (p = 0.48) (Table 1). There was no significant difference in SCC between herds with breeds of Greek origin or imported ones (0.894 × 106 versus 0.754 × 106 cells mL−1, respectively; p = 0.11).
For increased SCC (>0.75 × 106 cells mL−1) in bulk-tank milk, in both analyses (i.e., with herds with all breeds and with herds with the four most frequently recorded breeds), a significant association was seen during the univariable analysis only for the breed (p = 0.012 and 0.001, respectively) (Table S5). In the multivariable analysis (Table S2), the breed emerged as a significant factor for increased SCC in the herds in both analyses (p = 0.012 and 0.042, respectively) (Figure 4), whilst the management system was found to be significant in the analysis with herds with the four most frequently recorded breeds (p = 0.01) (Table 5).

3.2.2. Effect on Total Bacterial Counts

No significant difference was seen in TBC between the breeds, either when flocks with all breeds were considered or when only flocks with the six most frequently recorded breeds were taken into account (p > 0.23) (Table 1). There was no significant difference in TBC between herds with breeds of Greek origin or imported ones (522 × 103 versus 694 × 103 cfu mL−1, respectively; p = 0.29).
For increased TBC (>1500 × 103 cfu mL−1) in bulk-tank milk, no significant association was seen during the univariable analysis with any variable in both analyses (p > 0.06) (Table S5). Among the variables included in the multivariable analysis (Table S2), none emerged as a significant factor for increased TBC (p > 0.06).

4. Discussion

In this work, we studied the sheep and goat breeds present in small ruminant farms in Greece and we evaluated possible differences between the breeds in SCC and TBC. Farms from all regions of the country were included into the study; hence, conditions prevailing throughout Greece had been taken into account and regional factors had a smaller influence. Additionally, during the study, we employed consistent methodologies for sampling and sample processing, whilst the same investigators performed specific tasks and always maintained consistency, which all helped to minimize possible bias. Sheep and goat milk are important agricultural products in Greece, and the country produces approximately 25% of the total European small ruminant milk. However, the possible association of the sheep and goat breeds prevalent in the country in the quality of milk have not been hitherto described.
A significant proportion of the breeds in the farms was imported breeds. This indicates the continuous efforts of farmers to improve the genetic material in their flocks/herds and to achieve better production, with special reference to milk production. It is noteworthy that all imported breeds are characterized for high milk yield. This was reflected in the higher average milk production achieved in those farms than in farms with breeds of Greek origin. Consequently, this could also be seen in the management system applied in the farms; higher milk production requires increased input into farms: infrastructures (e.g., installation of milking systems), reduced grazing, increased staff numbers, etc., all of which are in line with the intensive or semi-intensive management system.
For SCC, as the current EU legislation does not provide a legal threshold, we employed the threshold of 0.75 × 106 cells mL−1, following the relevant US regulations and taking into account that dairy factories in the country use this threshold to regulate payments to farmers according to the quality of milk delivered. A similar practice is also applied in Spain, where payments for milk by dairy factories similarly depend upon the SCC therein [25].
For sheep, genomic selection has been shown to have good accuracy for mastitis resistance [26] but only when predicting resistance in closely related animals [27]. In the present study, it was found that breed was of significance for increased SCC, i.e., above the threshold of 0.75 × 106 cells mL−1. In Lacaune-breed sheep, a long-standing program for improvements in the resistance to mastitis has been running for over 15 years now [8,28]. The development of genomic methods and tools for applications to sheep and the sequencing of the sheep genome [29] have greatly contributed in those efforts [13,30]. In a recent study, Oget et al. [31] have confirmed the QTL regions associated with mastitis resistance in sheep of that breed.
The Greek breed most frequently recorded in the farms was the Chios. These animals are considered highly susceptible to mastitis. For example, Kiossis et al. [32] have reported that Staphylococcus chromogenes can survive in the mammary glands of Chios-breed ewes during the dry-period and cause subclinical infection after the subsequent lambing; in another study, Vasileiou et al. [33] have reported that the prevalence of subclinical mastitis in flocks with Chios-breed animals was 32%, higher than the national average of 26% [33]. Although genomic regions associated with specific mastitis traits in Chios-breed sheep have been identified [34], unfortunately, in Greece, no genetic improvement programmes have been developed for resistance to mastitis despite the significance of the breed for the sheep population of the country. This was reflected in the fact that there was a high proportion of farms with Chios-breed animals that had increased somatic cell counts in the bulk-tank milk; in fact, the results indicated that farms with Chios-breed animals had higher odds for increased SCC. Possible reasons for the increased susceptibility of the breed could be the bad udder conformation in animals of this breed, which adversely affects correct milking and can predispose animals to mammary infections [2]. In Latxa- and Sarda-breed animals, a positive association was found between SCC and udder conformation [35,36], suggesting that udders with what is perceived to be a good shape would be less affected by subclinical mastitis. Moreover, in Chios-breed animals, a peri-parturient immunosuppression associated with macrophage and neutrophil function has been found and can contribute to the development of mastitis [37].
In contrast, there was a low proportion of farms with animals of the Mytilini-breed that had increased SCC in bulk-tank milk. Future studies could be based on breed as a potential model for mastitis resistance among Greek breeds. Unfortunately, however, the breed is native to the island of Lesvos, where bluetongue is endemic [38] and, consequently, the Ministry of Agriculture has imposed a total ban on animal movements outside that island.
In goats, there is significantly less work regarding potential resistance to mastitis. Some work has been carried out in Alpine-breed [11] or Saanen-breed [39] animals. Again, in the Greek breeds, no efforts toward relevant genetic improvement have been made and this was reflected in farms with the local Indigenous Greek breed (Capra prisca) being at higher odds of having increased SCC in bulk-tank milk. However, no systematic efforts have been carried out as in sheep for resistance to the infection.
However, apart from breed, other factors were also found to be important for increased SCC, specifically the use of hand-milking for sheep and the application of semi-extensive or extensive management systems for goats (in which most often hand-milking is used). Similar findings have been reported in Spain by Gonzalo et al. [40], who studied four breeds (Assaf, Awassi, Churra and Castelana) and found that SCC in hand-milked flocks were higher than in machine-milked ones. In hand-milked ewes, staphylococci from the hands of milkers can be transmitted to the teats of ewes [41]. Likely as a consequence of that, hand-milking has been found to result in increased bacterial colonisation within the teat duct [42]. These bacteria can lead to mastitis, when the local defences of the teat are impaired [43,44], during which SCC increases, and reflect higher SCC in the bulk-tank milk of the farm.
In contrast, in goats, no such differences were reported in a study of two breeds (Alpine and Nubian) and SCC did not differ between machine- or hand-milked herds [45]. In goat farms, we found that the application of extensive or semi-extensive management systems can be significant factors for increased SCC. In these conditions, indeed hand-milking is more frequently practiced and certainly this may be a contributing factor, but not to as large an extent as in sheep [45]. Other factors can also play a role in these conditions; for example, incorrect nutrition can lead in increased incidence of mastitis [46] causing increased SCC, but on the other hand, it should be noted that the reduced milk production found in herds under extensive or semi-extensive management can also be a partial determinant for the increased SCC [47].
With regard to TBC, no effect of breed was evident. This was reasonable given that total bacterial counts depend on sources other than the mammary gland. TBC can reflect bacterial populations on the udder skin (transferred to milk at the time of milking) as well as from the surface of milk handling and storage equipment, e.g., teatcups, pipelines of the milking parlour and milk tanks, or from the hands of milkers themselves in flocks/herds where hand-milking is applied. Although there is a correlation between SCC and TBC in milk, the role of animal breed was not found to be important for TBC in bulk-tank milk.

5. Conclusions

In an extensive countrywide study throughout Greece, it was found that the animal breed had an effect on the presence of increased somatic cell counts (>0.75 × 106 cells mL−1) in the bulk-tank milk of sheep and goat farms.
The above also underlines the importance of introducing mastitis resistance as a trait for potential genetic improvement in the sheep and goat industry in Greece. In a country where sheep and goats constitute a large part of the agriculture sector and which has a large proportion of milk production in Europe [48], genetic improvement for mastitis resistance will contribute to better health and welfare of sheep and goats and to increased productivity on farms.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/app11167356/s1, Table S1: Zootechnical characteristics of 444 small ruminant farms in Greece, classified in accordance with the breed of animals therein. Table S2: Details of multivariable models employed for the evaluation of “increased somatic cell counts in bulk-tank milk” (i.e., with SCC above 0.75 × 106 cells mL−1) and ”increased total bacterial counts in bulk-tank milk” (i.e., with TBC above 1500 × 103 cfu mL−1) with the zootechnical characteristics of 444 small ruminant farms in Greece. Table S3: Associations of zootechnical characteristics with increased somatic cell counts in bulk-tank milk (i.e., above 0.75 × 106 cells mL−1) in 325 sheep farms in Greece, as found in a univariable analysis. Table S4: Associations of zootechnical characteristics with increased total bacterial counts in bulk-tank milk (i.e., above 1500 × 103 cfu mL−1) in 325 sheep farms in Greece, as found in a univariable analysis. Table S5: Associations of zootechnical characteristics with increased somatic cell counts in bulk-tank milk (i.e., above 0.75 × 106 cells mL−1) in 119 goat farms in Greece, as found in a univariable analysis. Table S6: Associations of zootechnical characteristics with increased total bacterial counts in bulk-tank milk (i.e., above 1500 × 103 cfu mL−1) in 119 goat farms in Greece, as found in a univariable analysis.

Author Contributions

Conceptualization, D.T.L., N.G.C.V., V.S.M. and G.C.F.; methodology, D.T.L., C.K.M., N.G.C.V., V.S.M. and G.C.F.; validation, G.C.F.; formal analysis, D.T.L. and G.C.F.; investigation, D.T.L., C.K.M., N.G.C.V., D.V.L. and G.C.F.; data curation, D.T.L. and G.C.F.; writing—original draft preparation, D.T.L. and G.C.F.; writing—review and editing, D.T.L., N.G.C.V., V.S.M., M.C. and G.C.F.; visualization, D.T.L. and G.C.F.; supervision, V.S.M., M.C. and G.C.F.; project administration, G.C.F.; funding acquisition, G.C.F. All authors have read and agreed to the published version of the manuscript.

Funding

Partial funding for this project was received by the Research Committee of the University of Thessaly. The funding body 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. The contribution is gratefully acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available on request from the corresponding author. The data are not publicly available as they form part of the Ph.D. thesis of the first author, which has not yet been examined, approved and uploaded in the official depository of Ph.D. thesis from Greek Universities.

Acknowledgments

The help received by E. Albanaki, M. Ali, E. Angelidou, C. Arsenopoulos, E. Charalambidis, M. Choutea, T. Chovolos, P. Dagkas, A. Daousi, V. Delistamatis, G. Dimou, I. Farmakis, E. Fotiadi, C. Gogos, D.A. Gougoulis, A. Hantzara, K.S. Ioannidi, T. Kanavos, T. Karatzinos, T. Kazakos, N.G. Kordalis, A. Koukoufikas, I. Lakafosi, G. Loufardaki, C. Mathiopoulos, I. Michalopoulos, G. Mousourakis, P. Ntais, D.C. Orfanou, S. Pakos, C. Pantopoulos, I. Papadakis, A. Papathanasiou, A. Papoutsi, C. Peleki, N. Petis, A. Pechlivanidis, A.P. Politis, P. Skourfounta, S.A. Spanos, M.C. Statira, E. Tassos, P. Toulatos and S. Tsakalidis is greatly appreciated; these veterinarians contributed by selecting the flocks and by arranging the visits to these flocks.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pourlis, A.F. A review of morphological characteristics relating to the production and reproduction of fat-tailed sheep breeds. Trop. Anim. Health Prod. 2011, 43, 1267–1287. [Google Scholar] [CrossRef]
  2. Gelasakis, A.I.; Arsenos, G.; Valergakis, G.E.; Oikonomou, G.; Kiossis, E.; Fthenakis, G.C. Study of factors affecting udder traits and assessment of their interrelationships with milking efficiency in Chios breed ewes. Small Rumin. Res. 2012, 103, 232–239. [Google Scholar] [CrossRef]
  3. Mastranestasis, I.; Ekateriniadou, L.V.; Ligda, C.; Theodorou, K. Genetic diversity and structure of the Lesvos sheep breed. Small Rumin. Res. 2015, 130, 54–59. [Google Scholar] [CrossRef]
  4. Kominakis, A.; Hager-Theodorides, A.L.; Saridaki, A.; Antonakos, G.; Tsiamis, G. Genome-wide population structure and evolutionary history of the Frizarta dairy sheep. Animal 2017, 11, 1680–1688. [Google Scholar] [CrossRef]
  5. Koutsouli, P.; Simitzis, P.; Theodorou, G.; Massouras, T.; Bizelis, I.; Politis, I. The effect of milking frequency reduction from twice to once daily on mammary physiology and animal welfare of two dairy Greek sheep breeds. Small Rumin. Res. 2017, 147, 18–24. [Google Scholar] [CrossRef]
  6. Michailidou, S.; Tsangaris, G.T.; Tzora, A.; Skoufos, I.; Banos, G.; Argiriou, A.; Arsenos, G. Analysis of genome-wide DNA arrays reveals the genomic population structure and diversity in autochthonous Greek goat breeds. PLoS ONE 2019, 14, e0226179. [Google Scholar]
  7. Argyriadou, A.; Gelasakis, A.I.; Banos, G.; Arsenos, G. Genetic improvement of indigenous Greek sheep and goat breeds. J. Hell. Vet. Med. Soc. 2020, 71, 2063–2072. [Google Scholar] [CrossRef]
  8. Rupp, R.; Boichard, D. Genetics of resistance to mastitis in dairy cattle. Vet. Res. 2003, 34, 671–688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Rupp, R.; Bergonier, D.; Dion, S.; Hygonenq, M.C.; Aurel, M.R.; Robert-Granié, C.; Foucras, G. Response to somatic cell count-based selection for mastitis resistance in a divergent selection experiment in sheep. J. Dairy Sci. 2009, 92, 1203–1219. [Google Scholar] [CrossRef] [Green Version]
  10. Rupp, R.; Lagriffoul, G.; Astruc, J.M.; Barillet, F. Genetic parameters for milk somatic cell scores and relationships with production traits in French Lacaune dairy sheep. J. Dairy Sci. 2003, 86, 1476–1481. [Google Scholar] [CrossRef]
  11. Rupp, R.; Huau, C.; Caillat, H.; Fassier, T.; Bouvier, F.; Pampouille, E.; Clement, V.; Palhiere, 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]
  12. Davies, G.; Genini, S.; Bishop, S.C.; Giuffra, E. An assessment of opportunities to dissect host genetic variation in resistance to infectious diseases in livestock. Animal 2009, 3, 415–436. [Google Scholar] [CrossRef] [Green Version]
  13. Oget, C.; Tosser-Klopp, G.; Rupp, R. Genetic and genomic studies in ovine mastitis. Small Rumin. Res. 2019, 176, 55–64. [Google Scholar] [CrossRef]
  14. Cuccuru, C.; Moroni, P.; Zecconi, A.; Casu, S.; Caria, A.; Contini, A. Milk differential cell counts in relation to total counts in Sardinian ewes. Small Rumin. Res. 1997, 25, 169–173. [Google Scholar] [CrossRef]
  15. Boutinaud, M.; Jammes, H. Potential uses of milk epithelial cells: A review. Reprod. Nutr. Dev. 2002, 42, 133–147. [Google Scholar] [CrossRef] [PubMed]
  16. Paape, M.J.; Wiggans, G.R.; Bannerman, D.D.; Thomas, D.L.; Sanders, A.H.; Contreras, A.; Moroni, P.; Miller, R.H. Monitoring goat and sheep milk somatic cell counts. J. Dairy Sci. 2007, 68, 114–125. [Google Scholar] [CrossRef]
  17. Katsafadou, A.I.; Politis, A.P.; Mavrogianni, V.S.; Barbagianni, M.S.; Vasileiou, N.G.C.; Fthenakis, G.C.; Fragkou, I.A. Mammary defences and immunity against mastitis in sheep. Animals 2019, 9, 726. [Google Scholar] [CrossRef] [Green Version]
  18. Berthelot, X.; Lagriffoul, G.; Concordet, D.; Barillet, F.; Bergonier, D. Physiological and pathological thresholds of somatic cell counts in ewe milk. Small Rumin. Res. 2006, 62, 27–31. [Google Scholar] [CrossRef]
  19. US Department of Health and Human Services. Grade “A” Pasteurized Milk Ordinance; 2017 Revision; U.S. Department of Health and Human Services, Public Health Service, Food and Drug Administration: Washington, DC, USA, 2017.
  20. European Union. Regulation (EC) No. 853/2004 of the European Parliament and of the Council of 29 April 2004 laying down specific hygiene rules foron the hygiene of foodstuffs. Off. J. Eur. Union 2004, 47, 55–205. [Google Scholar]
  21. Laird, D.T.; Gambrel-Lenarz, S.A.; Scher, F.M.; Graham, T.E.; Reddy, R. Microbiological Count Methods. In Standard Methods for the Examination of Dairy Products, 17th ed.; Wehr, H.M., Frank, J.F., Eds.; APHA Press: Washington, DC, USA, 2004; pp. 153–186. [Google Scholar]
  22. Wiggans, G.R.; Shook, G.E. A lactation measure of somatic cell count. J. Dairy Sci. 1987, 70 (Suppl. 13), 2666–2672. [Google Scholar] [CrossRef]
  23. Franzoi, M.; Manuelian, C.L.; Penasa, M.; De Marchi, M. Effects of somatic cell score on milk yield and mid-infrared predicted composition and technological traits of Brown Swiss, Holstein Friesian, and Simmental cattle breeds. J. Dairy Sci. 2020, 103, 791–804. [Google Scholar] [CrossRef] [PubMed]
  24. European Food Safety Authority. Scientific opinion on the welfare risks related to the farming of sheep for wool, meat and milk production. EFSA J. 2014, 12, 3933–4060. [Google Scholar]
  25. Gonzalo, C.; Juarez, M.T.; Garcia-Jimeno, M.C.; De La Fuente, L.F. Bulk tank somatic cell count and total bacterial count are affected by target practices and milking machine features in dairy sheep flocks in Castilla y Leon region, Spain. Small Rumin. Res. 2019, 178, 22–29. [Google Scholar] [CrossRef]
  26. Duchemin, S.I.; Colombani, C.; Legarra, A.; Baloche, G.; Larroque, H.; Astruc, J.M.; Barillet, F.; Robert-Granie, C.; Manfredi, E. Genomic selection in the French Lacaune dairy sheep breed. J. Dairy Sci. 2012, 95, 2723–2733. [Google Scholar] [CrossRef] [Green Version]
  27. Bishop, S.C. Genetic resistance to infections in sheep. Vet. Microbiol. 2015, 181, 2–7. [Google Scholar] [CrossRef]
  28. Detilleux, J.C.; Leroy, P.L. Application of a mixed normal mixture model for the estimation of mastitis-related parameters. J. Dairy Sci. 2000, 83, 2341–2349. [Google Scholar] [CrossRef]
  29. Jiang, Y.; Xie, M.; Chen, W.; Talbot, R.; Maddox, J.F.; Faraut, T.; Wu, C.; Muzny, D.M.; Li, Y.; Zhang, W. The sheep genome illuminates biology of the rumen and lipid metabolism. Science 2014, 344, 1168–1173. [Google Scholar] [CrossRef] [Green Version]
  30. Kijas, J.W.; Townley, D.; Dalrymple, B.P.; Heaton, M.P.; Maddox, J.F.; McGrath, A.; Wilson, P.; Ingersoll, R.G.; McCulloch, R.; McWilliam, S. A genome wide survey of SNP variation reveals the genetic structure of sheep breeds. PLoS ONE 2009, 4, e4668. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Oget, C.; Allain, C.; Portes, D.; Foucras, G.; Stella, A.; Astruc, J.M.; Sarry, J.; Tosser-Klopp, G.; Rupp, R. A validation study of loci associated with mastitis resistance in two French dairy sheep breeds. Gen. Sel. Evol. 2019, 51, 5. [Google Scholar] [CrossRef]
  32. Kiossis, E.; Brozos, C.N.; Petridou, E.; Zdragas, A.; Papadopoulos, T.; Boscos, C. Study on the possible survival of Staphylococcus chromogenes through the dry period in dairy ewes. Small Rumin. Res. 2013, 115, 124–129. [Google Scholar] [CrossRef]
  33. Vasileiou, N.G.C.; Cripps, P.J.; Ioannidi, K.S.; Chatzopoulos, D.C.; Gougoulis, D.A.; Sarrou, S.; Orfanou, D.C.; Politis, A.; Calvo Gonzalez-Valerio, T.; Argyros, S.; et al. Extensive countrywide field investigation of subclinical mastitis in sheep in Greece. J. Dairy Sci. 2018, 101, 7297–7310. [Google Scholar] [CrossRef] [Green Version]
  34. Banos, G.; Bramis, G.; Bush, S.J.; Clark, E.L.; McCulloch, M.E.B.; Smith, J.; Schulze, G.; Arsenos, G.; Hume, D.A.; Psifidi, A. The genomic architecture of mastitis resistance in dairy sheep. BMC Genom. 2017, 18, 624. [Google Scholar] [CrossRef] [PubMed]
  35. Legarra, A.; Ramon, M.; Ugarte, E.; Pérez-Guzmán, M.D.; Arranz, J. Economic weights of somatic cell score in dairy sheep. Animal 2007, 1, 205–212. [Google Scholar] [CrossRef] [Green Version]
  36. Sechi, S.; Salaris, S.; Carta, A.; Casu, S. Relationships between SCC and udder morphology traits in Sardinian sheep. In Proceedings of the 5th International Symposium on the Challenge to Sheep and Goat Milk Sectors, Alghero, Italy, 18–20 April 2007; p. 68. [Google Scholar]
  37. Theodorou, G.; Fragou, S.; Chronopoulou, R.; Kominakis, A.; Rogdakis, E.; Politis, I. Study of immune parameters in three Greek dairy sheep breeds during the periparturient period. J. Dairy Sci. 2007, 90, 5567–5571. [Google Scholar] [CrossRef] [PubMed]
  38. World Organisation for Animal Health. World Animal Health Information System. Available online: https://wahis.oie.int/#/home (accessed on 18 June 2021).
  39. Martin, P.; Palhiere, I.; Maroteau, C.; Clement, V.; David, I.; Tosser-Klopp, G.; Rupp, R. Genome-wide association mapping for type and mammary health traits in French dairy goats identifies a pleiotropic region on chromosome 19 in the Saanen breed. J. Dairy Sci. 2018, 101, 5214–5221. [Google Scholar] [CrossRef]
  40. Gonzalo, C.; Carriedo, J.A.; Blanco, M.A.; Beneitez, E.; Juarez, M.T.; De La Fuente, L.F.; San Primitivo, F. Factors of variation influencing bulk tank somatic cell count in dairy sheep. J. Dairy Sci. 2005, 88, 969–974. [Google Scholar] [CrossRef]
  41. Marco Melero, J.C. Mastitis in Laxta Breed Sheep: Epidemiology, Diagnosis and Control. Ph.D. Thesis, University of Zaragoza, Zaragoza, Spain, 1994. [Google Scholar]
  42. Mavrogianni, V.S.; Cripps, P.J.; Fthenakis, G.C. Bacterial flora and risk of infection of the ovine teat duct and mammary gland throughout lactation. Prev. Vet. Med. 2007, 79, 163–173. [Google Scholar] [CrossRef]
  43. Fragkou, I.A.; Mavrogianni, V.S.; Cripps, P.J.; Gougoulis, D.A.; Fthenakis, G.C. The bacterial flora in the teat duct of ewes can protect against and can cause mastitis. Vet. Res. 2007, 38, 525–545. [Google Scholar] [CrossRef] [Green Version]
  44. Fragkou, I.A.; Papaioannou, N.; Cripps, P.J.; Boscos, C.M.; Fthenakis, G.C. Teat lesions predispose to invasion of the ovine mammary gland by Mannheimia haemolytica. J. Comp. Pathol. 2007, 137, 239–244. [Google Scholar] [CrossRef]
  45. Zeng, S.S.; Escobar, E.N. Effect of breed and milking method on somatic cell count, standard plate count and composition of goat milk. Small Rumin. Res. 1996, 19, 169–175. [Google Scholar] [CrossRef]
  46. Caroprese, M.; Giannenas, I.; Fthenakis, G.C. Interactions between nutritional approaches and defences against microbial diseases in small ruminants. Vet. Microbiol. 2015, 181, 8–14. [Google Scholar] [CrossRef] [PubMed]
  47. Jones, G.M.; Pearson, R.E.; Clabaugh, G.A.; Heald, C.W. Relationships between somatic cell counts and milk production. J. Dairy Sci. 1984, 67, 1823–1831. [Google Scholar] [CrossRef]
  48. Pulina, G.; Milan, M.J.; Lavin, M.P.; Theodoridis, A.; Morin, E.; Capote, J.; Thomas, D.L.; Francesconi, A.H.D.; Caja, G. Current production trends, farm structures, and economics of the dairy sheep and goat sector. J. Dairy Sci. 2018, 101, 6715–6729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Locations of 444 small ruminant farms around Greece that were visited for bulk-tank milk sampling.
Figure 1. Locations of 444 small ruminant farms around Greece that were visited for bulk-tank milk sampling.
Applsci 11 07356 g001
Figure 2. Average yearly milk production per animal in sheep (solid pattern) or goat (motif pattern) farms in Greece, applying the intensive (green), semi-intensive (blue), semi-extensive (yellow) or extensive (red) management system.
Figure 2. Average yearly milk production per animal in sheep (solid pattern) or goat (motif pattern) farms in Greece, applying the intensive (green), semi-intensive (blue), semi-extensive (yellow) or extensive (red) management system.
Applsci 11 07356 g002
Figure 3. Proportion of flocks with increased somatic cell counts (>0.75 × 106 cells mL−1) in the bulk-tank milk of sheep flocks in Greece in accordance with the breed of sheep therein.
Figure 3. Proportion of flocks with increased somatic cell counts (>0.75 × 106 cells mL−1) in the bulk-tank milk of sheep flocks in Greece in accordance with the breed of sheep therein.
Applsci 11 07356 g003
Figure 4. Proportion of herds with increased somatic cell counts (>0.75 × 106 cells mL−1) in the bulk-tank milk of goat herds in Greece in accordance with the breed of goats therein.
Figure 4. Proportion of herds with increased somatic cell counts (>0.75 × 106 cells mL−1) in the bulk-tank milk of goat herds in Greece in accordance with the breed of goats therein.
Applsci 11 07356 g004
Table 1. Geometric mean somatic cell counts and total bacterial counts in bulk-tank milk in small ruminant farms in Greece according to animal breed in the farms.
Table 1. Geometric mean somatic cell counts and total bacterial counts in bulk-tank milk in small ruminant farms in Greece according to animal breed in the farms.
BreedFarms (n)Geometric Mean (95% CI 1)
Somatic Cells Counts (Cells mL−1)Total Bacterial Counts (cfu mL−1)
Sheep farms
Assaf300.528 × 106 (0.427 × 106–0.652 × 106)329 × 103 (191 × 103–562 × 103)
Awassi10.227 × 106495 × 103
Boutsko20.374 × 106 (0.333 × 106—0.418 × 106)425 × 103 (93 × 103—1905 × 103)
Chios440.592 × 106 (0.467 × 106—0.754 × 106)442 × 103 (185 × 103—977 × 103)
Crossbreds430.547 × 106 (0.439 × 106—0.679 × 106)547 × 103 (288 × 103—776 × 103)
Friesarta120.457 × 106 (0.296 × 106—0.708 × 106)461 × 103 (204 × 103—1047 × 103)
Friesian130.530 × 106 (0.421 × 106—0.665 × 106)216 × 103 (135 × 103—347 × 103)
Karagouniko50.642 × 106 (0.252 × 106—1.638 × 106)272 × 103 (50 × 103—1479 × 103)
Kefallinia10.172 × 106677 × 103
Lacaune950.453 × 106 (0.398 × 106—0.515 × 106)360 × 103 (240 × 103—550 × 103)
‘Local’ 2550.516 × 106 (0.412 × 106—0.643 × 106)496 × 103 (389 × 103—891 × 103)
Mytilini180.276 × 1 06 (0.192 × 106—0.396 × 106)408 × 103 (209 × 103—794 × 103)
Sfakia60.597 × 1 06 (0.412 × 106—0.866 × 106)357 × 103 (257 × 103—490 × 103)
Goat farms
Alpine90.647 × 106 (0.477 × 106—0.878 × 106)659 × 103 (309 × 103—1413 × 103)
Crossbreds180.799 × 106 (0.579 × 106—1.104 × 106)305 × 103 (123 × 103—759 × 103)
Damascus180.936 × 106 (0.775 × 106—1.135 × 106)1009 × 103 (550 × 103—1862 × 103)
Kefallinia11.831 × 1067482 × 103
Indigenous Greek
(Capra prisca)
500.932 × 106 (0.802 × 106—1.081 × 106)589 × 103 (437 × 103—871 × 103)
Murciano-Granadina130.657 × 106 (0.471 × 106—0.915 × 106)608 × 103 (309 × 103—1202 × 103)
Saanen50.653 × 106 (0.424 × 106—1.009 × 106)281 × 103 (100 × 103—794 × 103)
Skopelos50.766 × 106 (0.415 × 106—1.407 × 106)625 × 103 (219 × 103—1950 × 103)
1 CI: confidence interval; 2 the term used by the respective farmers to refer to a variety of small-scale breeds, not always related between them, and each one prevailing only in some areas of the country with limited geographic dissemination.
Table 2. Zootechnical characteristics of 444 small ruminant farms in Greece, summarized according to the origin of the breeds of animals therein.
Table 2. Zootechnical characteristics of 444 small ruminant farms in Greece, summarized according to the origin of the breeds of animals therein.
Origin of Breed 1Farms (n)Management System Applied in the FarmsApplication of
Machine-Milking
Average No. of Animals per Farm 4Average Yearly Milk Production per Animal (L) 4
I 2s-Is-EEY 3N
Sheep farms
Greek1861869762312858305 ± 16181 ± 6
Imported139256231112712247 ± 11243 ± 8
Goat farms
Greek742944193539234 ± 22155 ± 9
Imported457201713114244 ± 42276 ± 20
1 Greek: indigenous (i.e., breed native of Greece or created in the country) or crossbred (i.e., cross between imported or indigenous breeds), imported: breed brought into Greece from another country; 2 management system applied in farms classified according to the standards of the European Food Safety Authority [24]: I: intensive, s-I: semi-intensive, s-E: semi-extensive, E: extensive; 3 Y: yes, N: no; 4 mean ± standard error of the mean.
Table 3. Results of multivariable analysis for increased somatic cell counts (>0.75 × 106 cells mL−1) in the bulk-tank milk of sheep flocks in Greece (mixed effects logistic regression and analysis with only flocks with the six most frequently recorded breeds).
Table 3. Results of multivariable analysis for increased somatic cell counts (>0.75 × 106 cells mL−1) in the bulk-tank milk of sheep flocks in Greece (mixed effects logistic regression and analysis with only flocks with the six most frequently recorded breeds).
Variables (n = 2)Odds Ratio 1
(95% Confidence Intervals)
p
Breed 0.014
Assaf (n = 30)3.429 (0.649–18.114)0.15
Chios (n = 44)5.539 (1.132–27.108)0.035
Friesarta (n = 12)2.667 (0.373–19.061)0.33
Friesian (n = 13)1.455 (0.177–11.937)0.73
Lacaune (n = 95)1.870 (0.394–8.873)0.43
Mytilini (n = 18)reference
Application of machine-milking 0.003
Yes (n = 186)reference
No (n = 26)3.129 (1.342–7.296)0.008
1 Odds ratios calculated against the lowest prevalence associations of the variables.
Table 4. Results of the multivariable analysis for increased total bacterial (>1500 × 103 cells mL−1) in the bulk-tank milk of sheep flocks in Greece (mixed effects logistic regression and analysis with only flocks with the six most frequently recorded breeds).
Table 4. Results of the multivariable analysis for increased total bacterial (>1500 × 103 cells mL−1) in the bulk-tank milk of sheep flocks in Greece (mixed effects logistic regression and analysis with only flocks with the six most frequently recorded breeds).
Variable (n = 1)Odds Ratio 1
(95% Confidence Intervals)
p
Month of lactation period at sampling 0.014
Up to 3rd month (n = 64)2.644 (1.067–6.552)0.036
3rd to 6th month (n = 54)2.698 (1.054–6.908)0.039
After 6th month (n = 94)reference
1 Odds ratios calculated against the lowest prevalence associations of the variables.
Table 5. Results of multivariable analysis for increased somatic cell counts (>0.75 × 106 cells mL−1) in the bulk-tank milk of goat herds in Greece (mixed effects logistic regression and analysis with only herds with the four most frequently recorded breeds).
Table 5. Results of multivariable analysis for increased somatic cell counts (>0.75 × 106 cells mL−1) in the bulk-tank milk of goat herds in Greece (mixed effects logistic regression and analysis with only herds with the four most frequently recorded breeds).
Variables (n = 2)Odds Ratio 1
(95% Confidence Intervals)
p
Breed 0.042
Alpine (n = 9)1.125 (0.183–6.935)0.90
Damascus (n = 18)5.850 (1.222–27.995)0.027
Indigenous Greek (Capra prisca) (n = 50)9.000 (2.295–35.296)0.002
Murciana-Granadina (n = 18)reference
Management system 0.01
Intensive (n = 7)reference
Semi-intensive (n = 22)3.611 (0.570–22.898)0.17
Semi-extensive (n = 43)5.769 (0.988–33.677)0.05
Extensive (n = 18)12.500 (1.600–97.650)0.016
1 Odds ratios calculated against the lowest prevalence associations of the variables.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lianou, D.T.; Michael, C.K.; Vasileiou, N.G.C.; Liagka, D.V.; Mavrogianni, V.S.; Caroprese, M.; Fthenakis, G.C. Association of Breed of Sheep or Goats with Somatic Cell Counts and Total Bacterial Counts of Bulk-Tank Milk. Appl. Sci. 2021, 11, 7356. https://doi.org/10.3390/app11167356

AMA Style

Lianou DT, Michael CK, Vasileiou NGC, Liagka DV, Mavrogianni VS, Caroprese M, Fthenakis GC. Association of Breed of Sheep or Goats with Somatic Cell Counts and Total Bacterial Counts of Bulk-Tank Milk. Applied Sciences. 2021; 11(16):7356. https://doi.org/10.3390/app11167356

Chicago/Turabian Style

Lianou, Daphne T., Charalambia K. Michael, Natalia G. C. Vasileiou, Dimitra V. Liagka, Vasia S. Mavrogianni, Mariangela Caroprese, and George C. Fthenakis. 2021. "Association of Breed of Sheep or Goats with Somatic Cell Counts and Total Bacterial Counts of Bulk-Tank Milk" Applied Sciences 11, no. 16: 7356. https://doi.org/10.3390/app11167356

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