DOI: 10.1590/1809-6891v20e-46898
ZOOTECNIA


QUALIDADE FÍSICO-QUÍMICA E MICROBIOLÓGICA DO LEITE BOVINO EM DIFERENTES SISTEMAS DE PRODUÇÃO E ESTAÇÕES DO ANO

PHYSICAL-CHEMICAL AND MICROBIOLOGICAL QUALITY OF BOVINE MILK IN DIFFERENT PRODUCTION SYSTEMS AND SEASONS



Diego Prado de Vargas1* - ORCID http://orcid.org/0000-0003-0927-9469
José Laerte Nörnberg2 - ORCID http://orcid.org/0000-0002-8366-4480
Rudolf Brand Scheibler3 - ORCID http://orcid.org/0000-0002-8885-9009
Fábio Antunes Rizzo3 - ORCID http://orcid.org/0000-0001-8573-8845
Luciano Antônio Ritt4 - ORCID http://orcid.org/0000-0003-2062-4622
Marceli Pazini Milani5 - ORCID https://orcid.org/0000-0001-5948-3017



1Universidade de Santa Cruz do Sul, Santa Cruz do Sul, RS, Brasil.
2Universidade Federal de Santa Maria, Santa Maria, RS, Brasil.
3Universidade Federal de Pelotas, Pelotas, RS, Brasil.
4Universidade Federal do Rio Grande do Sul, Porto Alegre, RS, Brasil.
5Instituto Federal Farroupilha, Santa Maria, RS, Brasil.
*Autor para correspondência - diegodevargas@hotmail.com



Resumo
O objetivo deste trabalho foi avaliar a interação entre sistemas produtivos e estações do ano sobre a qualidade do leite bovino. Três propriedades foram avaliadas mensalmente, durante as estações, em quatro graus de especialização dos sistemas de produção: altamente especializado, especializado, semiespecializado e não especializado. Em termos de qualidade do leite, consideraram-se os teores de gordura, proteína, sólidos totais, contagem de células somáticas, contagem bacteriana total, acidez titulável em graus Dornic, densidade, índice crioscópico, estabilidade ao alizarol, potencial de hidrogênio (pH) e teste de redução do azul de metileno. Além destes, devido à sua importância para o rendimento industrial dos produtos lácteos de base proteica, foi realizado o fracionamento da proteína total do leite em equivalente proteico do nitrogênio não proteico, da proteína verdadeira e da caseína. A acidez Dornic foi maior na primavera e no inverno, enquanto os valores mais elevados de contagem de células somáticas foram no verão. Os teores de proteína total e proteína verdadeira foram menores durante o verão e elevaram-se os sólidos totais no outono. Com a diminuição dos níveis de especialização das propriedades rurais, observaram-se aumento na contagem bacteriana total e diminuição do nitrogênio não proteico. A densidade apresentou comportamento inverso ao da gordura, enquanto o potencial de hidrogênio (pH) se expressou de modo antagônico à acidez titulável. Na análise multivariada dos dados, foi possível identificar as combinações entre estratos produtivos e estações do ano que apresentavam a melhor qualidade da proteína do leite. Nesse sentido, todos os arranjos de sistemas de produção com o outono, o sistema altamente especializado no inverno, o semiespecializado no verão e o especializado na primavera tiveram qualidade proteica positiva, destacando-se pelas menores contagens de células somáticas a combinação entre a estação primavera e o sistema semiespecializado de produção.
Palavras-chave: análise multivariada, análise de componentes principais, qualidade do leite, tipologia, unidades produtoras de leite.

Abstract
The objective of this work was to evaluate the effect of the interaction between production systems and seasons on physical-chemical and microbiological qualities of bovine milk. Three farms were assessed monthly, along the seasons, for four specialization levels of the production systems: highly specialized, specialized, partially specialized, and non-specialized. In terms of milk quality, the levels of fat, protein, total solids, somatic cell count, total bacterial count, titratable acidity in degrees Dornic, density, cryoscopic index, alizarol stability, potential of hydrogen (pH) and methylene blue reduction test. In addition, due to their importance for the industrial yield of protein-based dairy products, total milk protein fractionation was performed in protein equivalent of non-protein nitrogen, true protein and casein.The milk Dornic acidity was higher in spring and winter, while the highest somatic cell count values were in the summer. Total and true protein levels were lower during the summer and total solids increased in the autumn. With the decrease in the levels of specialization of farms, there was an increase in total bacterial count and decrease in non-protein nitrogen. The density presented an inverse behavior to that of the fat while the hydrogen potential (pH) was expressed in a way antagonistic to the acidity. In the multivariate analysis of the data it was possible to identify the combinations between production strata and stations that presented the best milk protein quality. In this sense, all the arrangements of production systems with autumn, winter highly specialized, summer partially specialized, spring specialized had positive protein quality, being highlighted by the smaller of somatic cells counts in the spring partially specialized combination.
Keywords: multivariate analysis, principal components analysis, milk quality, typology, milk production units.


Recebido em 15 de maio de 2017.
Aceito em 04 de fevereiro de 2019


Introdução



O leite brasileiro é produzido com expressiva heterogeneidade de sistemas de produção, incluindo um reduzido número de produtores mais especializados e um elevado número de produtores com intermediária ou nenhuma especialização, diferindo-se principalmente pela produtividade leiteira e pelo manejo nutricional e higiênico-sanitário adotado (1).

É de extrema importância para a cadeia láctea mundial identificar os sistemas de exploração e manejos adotados em rebanhos leiteiros ao longo das estações do ano, que favoreçam a obtenção de um leite de maior qualidade, englobando-se a composição centesimal (gordura, proteína e sólidos totais) e aspectos higiênico-sanitários (contagem de células somáticas e contagem bacteriana total) (2-4).

Da mesma maneira, o leite que chega à plataforma de recepção da indústria de lácteos também deve atingir requisitos estabelecidos pelo setor normativo quanto a acidez titulável em graus Dornic, densidade, índice crioscópico, estabilidade ao alizarol, potencial de hidrogênio (pH) e teste de redução do azul de metileno (TRAM). Além destes, o fracionamento da proteína do leite pode ser de grande importância para o processo de industrialização, já que as caseínas têm relação direta com o rendimento industrial de produtos lácteos de base proteica.

Sistemas de produção e estações do ano interferem na qualidade físico-química e microbiológica do leite (5), porém existem poucas informações referentes ao efeito conjunto desses fatores na qualidade dessa matéria-prima. Assim, objetivou-se verificar as variações na qualidade físico-química e microbiológica do leite bovino nas diferentes estações climáticas e dos sistemas de produção atualmente empregados na região Sul do Brasil.


Material e métodos



Durante o ano de 2015, foram avaliadas mensalmente 12 unidades produtoras de leite (UPL), três de cada uma das quatro classes estudadas, nas estações primavera (PRIM), verão (VER), outono (OUT) e inverno (INV), totalizando 144 amostras analisadas.

De acordo com suas características de produção, os distintos estratos produtivos foram classificados em sistema altamente especializado (AE), especializado (ES), semiespecializado (SE) e não especializado (NE) (Tabela 1).


Os teores de gordura, proteína, lactose e sólidos totais foram determinados em equipamento Bentley(r) 2000, enquanto a contagem de células somáticas (CCS) e a contagem bacteriana total (CBT), por citometria de fluxo, utilizando equipamento Somacount(r) 300 e Bactocount(r) IBC, respectivamente. Essas determinações foram obtidas a partir de laudos oficiais emitidos pelo Laboratório de Serviço de Análises de Rebanhos Leiteiros (SARLE), órgão credenciado junto ao Ministério da Agricultura, Pecuária e Abastecimento (MAPA).

A avaliação das amostras quanto ao índice crioscópico (realizada em crioscópio eletrônico digital ITRMK 540), densidade a 15°C (realizada por meio de termolactodensímetro de Quevenne), acidez titulável em graus Dornic e estabilidade ao alizarol, foi realizada segundo metodologia convencional recomendada pelo Ministério da Agricultura, Pecuária e Abastecimento(6). A estimativa da carga microbiológica do leite foi realizada pelo teste de redução do azul de metileno (TRAM), conforme BRASIL(7), enquanto o pH determinado por aparelho Testo - modelo 206, seguindo manual de instrução do equipamento.

O fracionamento do nitrogênio total (NT) em nitrogênio não proteico (NNP), nitrogênio da proteína verdadeira (N-PV) e nitrogênio caseínico (N-cas), bem como o percentual equivalente proteico do nitrogênio não proteico (P-NNP), proteína verdadeira (P-PV) e caseína (P-cas) foram determinados segundo Silva et al.(8).

Com o intuito de linearizar os dados, a CCS foi transformada em escore linear de células somáticas ECS = [log2(CCS/100)]+3(9) e a CBT transformada pelo logaritmo natural da CBT normal (CBTt).

Os resultados foram analisados em delineamento inteiramente casualizado com medida repetida no tempo, utilizando o procedimento MIXED e empregando estrutura paramétrica especial nas matrizes de (co)variância, conforme o modelo estatístico:
Yijk = μ + αi +γk+ (αγ)ik + εijk,
em que Yijk é o valor observado pela k-ésima estação e j-ésima UPL do i-ésimo sistema produtivo; μ + αi +γk+ (αγ)ik é a média para o sistema produtivo i na estação k, contendo os efeitos fixos para sistema produtivo i, estação k, e interação sistema produtivo × estação ik; εijk é o erro aleatório associado com a mensuração na estação k na j-ésima UPL designada ao sistema produtivo i. Assumiu-se uma estrutura de medidas repetidas para os resíduos, ou seja, pressupostos que εijk ~ N (0;∑), em que ∑ representa a matriz de (co)variância residual que mais adequou-se aos dados, em que as testadas foram: AR (1), autorregressiva de primeira ordem; VC, componentes de variância; CS, simetria composta; ARH(1), autorregressiva heterogênea de primeira ordem; ARMA(1,1), média móvel autorregressiva de primeira ordem, e TOEP, Toeplitz. As medidas foram estimadas com o comando LSMEANS (Least Squares Means) e comparadas pelo teste de Tukey ao nível de 5% de significância.

No conjunto original das variáveis que representam a qualidade físico-química e microbiológica do leite, procedeu-se a análise de variância multivariada pelo procedimento GLM e o comando MANOVA, em que foi realizado o teste de Wilks, para testar a hipótese de que os vetores de médias das distintas combinações entre sistemas produtivos e estações fossem nulos.

Após, foi utilizada a análise de componentes principais (Principal component analysis - PCA), seguida da análise de agrupamento (Cluster analysis), utilizando o algoritmo hierárquico aglomerativo de Ward como método de agrupamento e a distância euclidiana quadrática como medida de dissimilaridade.

As análises estatísticas foram executadas no aplicativo SAS(r) System for Windows(tm) versão 9.0 (SAS Institute Inc., Cary - NC, USA).


Resultados



As estações do ano tiveram efeito significativo sobre a CCS, a acidez titulável em graus Dornic, o pH, a densidade, os teores de proteína total, a proteína verdadeira e os sólidos totais (Tabela 2 e Tabela 3). A acidez titulável do leite foi maior na primavera e no inverno, enquanto os mais elevados valores de CCS foram encontrados no verão. Os teores de proteína total e verdadeira foram menores durante o verão (Tabela 3) e elevaram-se os sólidos totais no outono (Tabela 3).



Houve efeito significativo de níveis de especialização nas variáveis CBT, acidez titulável em graus Dornic, pH, densidade e P-NNP. Com a diminuição dos níveis de especialização das UPL, observou-se aumento na CBT (Tabela 2). Por outro lado, com a elevação da especialização dos estratos produtivos, observou-se um aumento do P-NNP (Tabela 3).

A densidade apresentou comportamento inverso ao da gordura, enquanto o pH se expressou de modo antagônico à acidez titulável em graus Dornic (Tabela 2). As únicas variáveis que apresentaram interação significativa (P<0,05) entre níveis de especialização e estações do ano foram o índice crioscópico do leite e a caseína.

Na análise de variância multivariada, verificou-se diferença (P<0,05) entre os vetores de médias para as combinações de sistemas produtivos e estações, tanto pelo teste Wilks quanto pelos testes Pillai, Hotelling-Lawley e Roy. Os três primeiros autovalores foram significativos pelo teste da razão de verossimilhança, sendo que o primeiro autovalor explicou 45,28%, o segundo 22,71% e o terceiro 14,04%, ou seja, os três primeiros autovalores explicaram 82,04% da variação total dos dados relacionados à qualidade do leite (Figura 1).


As variáveis ECS (r=-0,85; P<0,001), CBTt (r=0,77; P<0,001), acidez titulável em graus Dornic (r=0,90; P<0,001) e alizarol (r=-0,85; P<0,001) explicaram a variabilidade do CP1 (componente principal 1), ao passo que a proteína verdadeira (r=0,86; P<0,001) e a caseína (r=0,88; P<0,001) explicaram a variabilidade do CP2 (componente principal 2) e o nitrogênio não proteico (r=0,85; P<0,001) do CP3 (componente principal 3) (Figura 1).


Discussão



A acidez titulável em graus Dornic do leite foi maior na primavera e no inverno (Tabela 2). Além disso, no verão e no outono observaram-se os maiores valores de CCS (Tabela 2), o que pode ser indicativo de mastite subclínica, contribuindo para o maior teor de sódio e menores teores de cálcio, fósforo e potássio no leite (11), portanto, para a menor acidez titulável nessas estações climáticas (Tabela 2).

Assim, apesar de as estações do ano não influenciarem na CBT do leite, os maiores valores desse indicador higiênico-sanitário são observados na primavera e no inverno, devido à ação dos microrganismos sobre a lactose, o que resultaria no aumento do ácido láctico no leite (11) e, em conjunto com os possíveis efeitos que a mastite subclínica exerce sobre a sua composição, poderia resultar nos maiores valores de acidez titulável observados na primavera e no inverno (Tabela 2).

Nesse sentido, os maiores valores de CCS encontrados no verão podem ser resultantes das altas temperaturas dessa estação na região Sul do Brasil, o que poderia resultar em estresse térmico e limitação da resposta imunológica dos animais a doenças e sua susceptibilidade à invasão do úbere por microrganismos, resultando em aumento na CCS do leite (5).

O teor de proteína total e verdadeira foi menor durante o verão (Tabela 3). Isto se deve à baixa qualidade nutricional das principais pastagens tropicais (Sorghum sudanense e Pennisetum glaucum) em comparação às temperadas (Avena strigosa e Lolium multiflorum) encontradas na região de realização do estudo (12). Nesse sentido, as pastagens tropicais apresentam menores teores de proteína degradável no rúmen e carboidratos fermentáveis disponíveis que as pastagens temperadas (13), o que desfavorece a formação de proteína microbiana em nível ruminal e diminui os valores de proteína total e verdadeira no leite.

A elevação no conteúdo de sólidos totais no outono deveu-se principalmente à elevação no conteúdo de proteína e gordura do leite na estação referida (Tabela 3). Isto pode ser resultante dos efeitos diminutivos que o vazio forrageiro outonal exerce sobre a produção leiteira, o que pode ter causado uma concentração relativa desses constituintes, aumentando-se os sólidos totais do leite.

O aumento na CBT, com a diminuição dos níveis de especialização das UPL (Tabela 2), deveu-se, possivelmente, ao tipo e ao local de ordenha, bem como ao emprego ou não do resfriamento (Tabela 1). Assim, o sistema altamente especializado, por contar com sala de ordenha própria com circuito fechado e tanque de refrigeração a granel (Tabela 1), apresentou os menores valores de CBT (Tabela 2), enquanto o semiespecializado e o não especializado, por realizar ordenha balde ao pé em local não apropriado (Tabela 1), evidenciaram os maiores valores desse indicador higiênico- sanitário (Tabela 2). Nesse sentido, possivelmente com o aumento da CBT, aumentaram-se as bactérias mesófilas no leite, elevando-se a utilização da lactose como principal substrato desses microrganismos, o que poderia estar contribuindo para o aumento da acidez titulável em graus Dornic do leite à medida que diminuíram os níveis de especialização dos estratos produtivos (Tabela 2).

Observou-se aumento do NNP com a elevação dos níveis de especialização (Tabela 3), o que pode ser resultante do maior aporte dietético de proteína (13), já que essas UPL utilizavam maiores quantidade de concentrado na dieta total (Tabela 1). Somando-se a isso, apesar da ausência de significância dos estratos produtivos sobre a CCS, com o aumento das células somáticas no leite, eleva-se a atividade proteolítica de origem leucocitária (neutrófilos e macrófagos), o que poderia estar contribuindo para o maior valor de NNP observado no nível altamente especializado, já que estes, além de fornecerem mais de 10 Kg de concentrado/vaca/dia, têm os maiores valores de CCS do estudo (Tabela 2).

Além disso, as características dietéticas dos diferentes estratos produtivos podem estar contribuindo para a interação significativa (P<0,05) observada entre os níveis de especialização e estações do ano para a variável caseína (Tabela 3). Nesse sentido, o nível não especializado parece estar com um comprometimento no aporte de nutrientes necessários para a síntese de caseína, pois não fornece silagem de milho como principal volumoso (Tabela 1) e é dependente da variação sazonal das pastagens, o que, no verão, devido à baixa qualidade das gramíneas tropicais em relação às temperadas, resultou na redução dos valores da caseína do leite nesse estrato produtivo (Tabela 3).

Evidenciou-se interação (P<0,05) entre níveis de especialização e estações do ano no índice crioscópico do leite (Tabela 2). O alto valor de crioscopia do nível altamente especializado (-0,537 °H) e especializado (-0,508 °H) no outono pode ser resultante do vazio forrageiro outonal, já que esses estratos produtivos aumentam o fornecimento de concentrado nessa época do ano, resultando na elevação dos valores de NNP (14), o que pode ter resultado em aumento dos valores de crioscopia no leite (Tabela 2).

Na análise multivariada dos dados, a partir das correlações dos componentes principais com as variáveis, pode-se inferir que o componente principal 1 (CP1) representa a qualidade microbiológica, enquanto o componente principal 2 (CP2) e o componente principal 3 (CP3), a qualidade proteica do leite. Assim, a localização no plano cartesiano das combinações entre estratos produtivos e estações climáticas pode auxiliar na identificação do leite com maior rendimento industrial em relação a produtos lácteos de base proteica.

Portanto, combinações situadas no primeiro e no segundo quadrantes apresentam melhor qualidade da proteína do leite, o que aumentaria o rendimento industrial na fabricação de queijos. No outono, todos os níveis de especialização apresentaram qualidade proteica positiva (CP1 positivo). Devido à escassez forrageira dessa estação, os estratos altamente especializado, especializado e semiespecializado aumentaram o fornecimento de silagem de milho e concentrado para suprir as necessidades animais, o que aumentaria a proteína verdadeira e a caseína do leite (13). Por outro lado, nas propriedades não especializadas, é possível que esteja ocorrendo efeito de concentração relativa, já que esse estrato produtivo está mais susceptível ao vazio forrageiro outonal, o que diminuiria o seu volume de produção e concentraria essas frações nitrogenadas do leite. No entanto, ressalta-se que, entre estas, as situadas no primeiro quadrante (especialmente o SE-PRIM) apresentaram os menores valores de CCS, o que beneficiaria ainda mais o rendimento desses produtos à base proteica (15).

Ainda sobre a qualidade microbiológica do leite, no inverno, observa-se que os níveis mais especializados estão próximos da origem do plano cartesiano, enquanto os de menor especialização deslocaram-se para o quarto quadrante, pela maior contaminação bacteriana e pela acidez titulável em graus Dornic (Figura 1). Isto se deve provavelmente à maior umidade e à pluviometria nessas épocas do ano, o que aumentou o acúmulo de sujidades e fontes de contaminação do leite nos níveis de menor especialização, já que nesses estratos produtivos faziam-se uso da ordenha balde ao pé em galpões (Tabela 1).

Nesse sentido, ao se observar a segunda projeção bidimensional, nota-se que o SE-INV, além de deslocar-se para o quarto quadrante pela maior contaminação bacteriana e pela acidez titulável em graus Dornic do leite, alocou-se nessa posição pelos elevados valores de NNP. Isto se deve à contaminação por bactérias capazes de fermentar a lactose, formando ácido láctico, aumentando a acidez titulável e a proteólise bacteriana, o que pode evidenciar alto potencial de degradação do leite, resultando em aumento de suas substâncias não proteicas (16).


Conclusões



As estações do ano interferem na proteína total, na proteína verdadeira, nos sólidos totais e na contagem de células somáticas, enquanto os sistemas de produção de leite, na contagem bacteriana total e no nitrogênio não proteico. Esses dois fatores de variação influenciam separadamente a densidade, a acidez e o potencial de hidrogênio e, conjuntamente, a crioscopia e a caseína. A acidez titulável do leite em graus Dornic foi maior na primavera e no inverno, enquanto os valores mais elevados de contagem de células somáticas foram no verão. Os teores de proteína total e verdadeira foram menores durante o verão e elevaram-se os sólidos totais no outono. Com a diminuição dos níveis de especialização das unidades produtoras de leite, observaram-se aumento na contagem bacteriana total e diminuição do nitrogênio não proteico. A densidade apresentou comportamento inverso ao da gordura, enquanto o potencial de hidrogênio (pH) se expressou de modo antagônico à acidez titulável. Todas as combinações de sistemas de produção com o outono, o inverno - altamente especializado, o verão - semiespecializado e a primavera - especializado tiveram qualidade proteica positiva, destacando-se pelas menores contagens de células somáticas a combinação primavera-semiespecializado.


Agradecimentos


À Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) e à Universidade Federal de Santa Maria (UFSM), pelo apoio e pelo auxílio financeiro prestado, bem como à empresa de lácteos estudada, pela cessão do laboratório de controle de qualidade para a realização das análises físico-químicas.


Referências


1. Gabbi AM, McManus CM, Silva AV, Marques LT, Zanela MB, Stumpf MP, Fischer V. Typology and physical‑chemical characterization of bovine milk produced with different productions strategies. Agricultural Systems. 2013; 121(1):130‑134. English.
2. Bodenmüller Filho A., Damasceno JC, Previdelli ITS, Santana RG, Ramos CECO, Santos GT. Tipologia de sistemas de produção baseada nas características do leite. Revista. Brasileira de. Zootecnia. 2010; 39(8):1832-1839. Portuguese.
3. Vargas DP, Nörnberg JL, Mello RO, Sheibler RB, Milani MP, Mello FCB. Correlações entre contagem bacteriana total e parâmetros de qualidade do leite. Revista Brasileira de Ciência Veterinária. 2013; 20(4):241-247. Portuguese.
4. Vargas, DP, Nörnberg JL, Mello RO, Sheibler RB, Mello FCB, Milani MP. Correlações entre contagem de células somáticas e parâmetros físico-químicos e microbiológicos de qualidade do leite. Ciência Animal Brasileira. 2014; 5(4):473-483. Portuguese.
5. Vargas DP, Nörnberg JL, Sheibler RB, Júnior JS, Rizzo FA, Wagner R. Qualidade e potencial nutracêutico do leite bovino em diferentes sistemas de produção e estações do ano. Pesquisa Agropecuária Brasileira. 2015; 50(12):1208-1219. Portuguese.
6. Tronco VM. Manual para inspeção da qualidade do leite. Santa Maria: UFSM; 1997. 166p. Portuguese.
7. BRASIL, Ministério da Agricultura Pecuária e do Abastecimento, Secretaria de Defesa Agropecuária, Departamento de Defesa Animal, Manual de métodos microbiológicos para alimentos, Coordenação Geral de Laboratório Animal, 2ª revisão. 136p., 1991/1992. Portuguese.
8. Silva PHF, Pereira DBC, Oliveira LL, Costa Junior LCG. Físico-química do leite e derivados - métodos analíticos. São Paulo: Oficina de Impressão Gráfica e Editora Ltda.; 1997. 190 p. Portuguese.
9. Shook GE. Genetic improvement of mastitis through selection on somatic cell count. The Veterinary Clinics of North America: Food Animal Practice. 1993; 9(3):563-581. English.
10. STATISTICAL ANALYSIS SYSTEM - SAS. The SAS system for windows. v.9.0 Cary: SAS Institute Inc., 2002.
11. Machado PF, Pereira AR, Sarríes GA. Composição do leite de tanques de rebanhos brasileiros distribuídos segundo sua contagem de células somáticas. Revista Brasileira de Zootecnia. 2000; 29(6):1883- 1886.
12. Fanti MGN, Almeida KE, Rodrigues, AM, Silva, RC, Florence ACR, Gioielli LA, Oliveira MN. Contribuição ao estudo das características físico‑químicas e da fração lipídica do leite orgânico. Ciência e Tecnologia de Alimentos. 2008; 28:259‑265. Portuguese.
13. Noro G., González FHD, Campos R, Dürr JW. Fatores ambientais que afetam a produção e a composição do leite em rebanhos assistidos por cooperativas no Rio Grande do Sul. Revista Brasileira de Zootecnia. 2006; 35(3):1129-1135. Portuguese.
14. Kędzierska‐Matysek M, Litwińczuk Z, Florek M, Barłowska J. The effects of breed and other factors on the composition and freezing point of cow's milk in Poland. International Journal of Dairy Technology. 2011; 64(3): 336-342. English.
15. Coelho KO, Mesquita AJ, Machado PF, Lage ME, Meyer PM, Reis A.P. Efeito da contagem de células somáticas sobre o rendimento e a composição físico-química do queijo muçarela. Arquivos Brasileiros de Medicina Veterinária e Zootecnia. 2014; 66(4):1260-1268. Portuguese.
16. Arcuri EF, Silva PDL, Brito MAVP, Brito JRF, Lange CC, Magalhães MMA. Contagem, isolamento e caracterização de bactérias psicrotróficas contaminantes de leite cru refrigerado. Ciência Rural. 2008; 38(8): 2250-2255. Portuguese.