Qualidade e potencial nutracêutico do leite bovino em diferentes sistemas de produção e estações do ano

Quality and nutraceutical potential of bovine milk in different production systems and seasons

Resumo:

O objetivo deste trabalho foi avaliar o efeito da interação entre sistemas produtivos e estações do ano sobre a qualidade química e microbiológica do leite bovino, além de propor um índice de potencial nutracêutico do perfil de ácidos graxos do leite. Três propriedades foram avaliadas mensalmente, durante as diferentes estações do ano, em quatro graus de especialização dos sistemas de produção: altamente especializado, especializado, semiespecializado e não especializado. Os sistemas de produção e as estações do ano interferem de forma conjunta no perfil de ácidos graxos e, de forma isolada, na qualidade química e microbiológica do leite. As maiores contagens de células somáticas e os menores conteúdos de proteína foram observados no verão, e o grau de especialização das unidades produtivas esteve indiretamente relacionado à contagem bacteriana total no leite. No inverno, sistemas não especializados produziram leite com o melhor índice nutracêutico, que apresentaram os maiores teores de ácidos graxos poli-insaturados, ácido rumênico (CLA, 18:2n7-c9,t11) e ácido t10,c12-octadecadienoico (CLA, 18:2n6-t10,c12).

Termos para indexação:
ácido linoleico conjugado; alimentos funcionais; perfil de ácidos graxos; qualidade do leite bovino

Abstract:

The objective of this work was to evaluate the effect of the interaction between production systems and seasons on the chemical and microbiological qualities of bovine milk, besides proposing a potential nutraceutical index of milk fatty acid profiles. Three farms were assessed monthly, along the different seasons, for four specialization levels of the production systems: highly specialized, specialized, partially specialized, and nonspecialized. The production systems and the seasons jointly interfere in fatty acid profiles, and singly interfere in the chemical and microbiological quality of milk. The highest somatic cell counts and the lowest protein contents were observed in the summer, and the specialization level of the production units had an inverse relationship with the total bacterial count. During winter, nonspecialized systems produced milk with the best nutraceutical index, with the highest contents of polyunsaturated fatty acids, rumenic acid (CLA, 18:2n7-c9,t11) and t10,c12-octadecadienoic acid (CLA, 18:2n6-t10,c12).

Index terms:
conjugated linoleic acid; functional foods; fatty acid profiles; bovine milk quality

Introdução

O leite é produzido com expressiva heterogeneidade de sistemas de produção, em todo o território brasileiro, com menor número de produtores especializados (sistemas intensivos de produção) e elevado número de produtores com intermediária ou nenhuma especialização (sistemas extrativistas de produção). Esses sistemas diferem, principalmente, pela produtividade leiteira e pelo manejo nutricional e higiênico-sanitário adotado (Gabbi et al., 2013).

Para a cadeia láctea mundial, é de grande relevância, identificar os sistemas de manejo e produção de rebanhos leiteiros que favoreçam a obtenção de leite com maior qualidade, ao longo das diferentes estações do ano (Bodenmüller Filho et al., 2010).

Quanto à qualidade do leite, do ponto de vista legal e prático, consideram-se a composição centesimal - que inclui os teores de gordura, proteína e sólidos totais - e os aspectos higiênico-sanitários, que englobam padrões internacionais relacionados à contagem de células somáticas (CCS) e à contagem bacteriana total (CBT) (Vargas et al., 2013, 2014).

Com o maior acesso à informação sobre o papel dos alimentos na saúde, o conceito de qualidade tem-se tornado mais abrangente. Assim, uma nova classe de alimentos tem surgido, denominada de "alimentos nutracêuticos", capazes de proporcionar benefícios à saúde do consumidor, que incluem prevenção ou tratamento de doenças (Tripathi, 2014).

Com foco no potencial nutracêutico de alguns ácidos graxos, tem-se buscado a elevação dos conteúdos de ácido rumênico (CLA, 18:2n7-c9,t11) e ácido t10,c12-octadecadienoico (CLA, 18:2n6-t10,c12) no leite, que podem contribuir para a prevenção e combate de muitas doenças em humanos, como obesidade, doenças cardiovasculares, diabetes do tipo-2, e alguns tipos de câncer (Kratz et al., 2013; Lahlou et al., 2014).

Sistemas de produção e estações do ano podem interferir na qualidade química e microbiológica do leite (Gabbi et al., 2013); porém, há pouca informação disponível quanto ao efeito conjunto desses fatores sobre a qualidade e o potencial nutracêutico do leite.

O objetivo deste trabalho foi avaliar o efeito da interação entre sistemas de produção e estações do ano sobre a qualidade química e microbiológica do leite bovino, além de propor um índice de potencial nutracêutico do perfil de ácidos graxos.

Material e Métodos

Três propriedades foram avaliadas, mensalmente, durante as diferentes estações do ano, em quatro graus de especialização dos sistemas de produção: altamente especializado, especializado, semiespecializado e não especializado (Tabela 1), tendo-se analisado um total de 144 amostras, de 12 propriedades.

Tabela 1:
Características dos sistemas de produção de leite altamente especializado (AE), especializado (ES), semiespecializado (SE) e não especializado (NE).

• ⁽¹⁾
  Sim, silagem de milho como base do volumoso diário; Não, silagem de milho como parte do volumoso diário.

As amostras foram coletadas em 2012, em fazendas localizadas em municípios compreendidos desde a Região Central do Rio Grande do Sul até o início da Serra Gaúcha (29°15'16" a 29°29'13"S; e 51°31'52" a 51°59'50"W). Paralelamente, foram coletadas amostras dos alimentos destinados aos animais, para sua caracterização (Tabela 2). Os lipídeos foram extraídos das amostras de leite e de alimento - silagem de milho, concentrado, pastagens temperadas (aveia-preta, Avena strigosa, e azevém, Lolium multiflorum) e pastagens tropicais (capim-sudão, Sorghum sudanense, e milheto, Pennisetum glaucum), conforme Bligh & Dyer (1959). A transesterificação e a metilação seguiram a técnica descrita por Christie (1982). O perfil de ácidos graxos foi determinado de acordo com Simionato et al. (2010). A identificação dos ésteres metílicos de ácidos graxos (Emag) foi realizada pela comparação dos tempos de retenção dos constituintes da amostra, tendo-se utilizado mistura composta de 37 Emag (24056, Supelco, 37 Component Fame Mix) e padrão do éster metílico do ácido transvacênico (18:1n7-t11) (46905-U, trans-11-octadecenoic methyl ester), com mistura dos isômeros do éster metílico do ácido rumênico (CLA, 18:2n7-c9,t11) e ácido t10,c12-octadecadienoico (CLA, 18:2n6-t10,c12) (O5632, Linoleic acid, conjugated methyl ester).

Tabela 2:
Caracterização do perfil de ácidos graxos (g 100g^-1 do total de ácidos graxos), dos principais constituintes da dieta recebida pelos animais, nos sistemas de produção altamente especializado (AE), especializado (ES), semiespecializado (SE) e não especializado (NE).

• ⁽¹⁾
  Aveia-preta (Avena strigosa) e azevém (Lolium multiflorum).
• ⁽²⁾
  Capim-sudão (Sorghum sudanense) e milheto (Pennisetum glaucum). AAL, ácido a-linolênico.

As amostras de alimento foram expressas em percentagem do total de ácidos graxos (g 100g^-1), e as de leite em lipídeos totais (mg g^-1), em relação ao padrão interno [tricosanoato de metila (23:0) (T9900, methyl tricosanoate)]. Os cálculos foram realizados conforme Simionato et al. (2010), com a equação: C (mg g^-1) = (A_x × m_23:0 × F_RT)/(A_23:0 × M_A × F_CT), em que A_X é a área dos ésteres metílicos dos ácidos graxos; A_23:0, a área do padrão interno; m_23:0, a massa do padrão interno adicionado à amostra (em mg); M_A, a massa da amostra (em g); F_RT, o fator teórico de resposta dos ésteres metílicos de ácidos graxos; F_CT, o fator de conversão para expressão dos resultados em miligramas de ácidos graxos por grama de lipídeos totais.

Os resultados da determinação dos ácidos graxos do leite foram utilizados para calcular os índices de aterogenicidade e trombogenicidade (Ulbricht & Southgate, 1991), bem como a razão entre ácidos graxos hipocolesterolêmicos e hipercolesterolêmicos (Santos-Silva et al., 2002) e ácidos graxos desejáveis (AGD) (Costa et al., 2008). O índice de AGD leva em consideração o somatório de ácidos graxos monoinsaturados (ΣAGMI), poli-insaturados (ΣAGPI) e o ácido esteárico (C18:0), e não somente os ácidos com reconhecido potencial nutracêutico. Portanto, este estudo propõe uma equação, com ácidos graxos que apresentem esse potencial, denominada de índice nutracêutico do leite bovino (IN): IN = (ácido transvacênico + ácido rumênico + ácido t10,c12-octadecadienoico)/Σn6:n3. De acordo com a equação, o numerador é diretamente proporcional à concentração de ácidos graxos benéficos à saúde, e inversamente proporcional à relação ômega 6 e ômega 3 (Σn6:n3), o que está de acordo com a diminuição recomendada nesta razão (Simopoulos, 2006).

Os teores de gordura, proteína, lactose e sólidos totais do leite foram determinados em equipamento Bentley 2000; enquanto a CCS e a CBT foram obtidas por citometria de fluxo, em equipamento Somacount 300 e Bactocount IBC, respectivamente. Com o intuito de linearizar os dados, a CCS foi transformada em escore linear de células somáticas para ECS = [log2(CCS/100)]+3 (Shook, 1993), e a CBT, transformada pelo logaritmo natural da CBT normal.

Os resultados foram analisados conforme delineamento inteiramente casualizado, com medidas repetidas no tempo, tendo-se utilizado o procedimento Mixed e empregado estrutura paramétrica especial nas matrizes de variância e covariância, conforme o modelo estatístico: Y_ijk = μ + α_i +γ_k+ (αγ)_ik + ε_ijk, em que: Y_ijk é 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, que contém os efeitos fixos para o sistema produtivo i, estação k, bem como a interação sistema produtivo x estação ik; e ε_ijk é o erro aleatório, associado à mensuração na estação k na j-ésima unidade de produção de leite (UPL), designada ao sistema produtivo i. Assumiu-se uma estrutura de medidas repetidas para os resíduos; ou seja, utilizaram-se os pressupostos de que ε_ijk ~ N (0;∑), em que: ∑ representa a matriz de covariância e covariância residual que mais adequou-se aos dados. As matrizes 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 e comparadas pelo teste de Tukey, a 5% de probabilidade.

No conjunto original das variáveis que representam o potencial nutracêutico do leite, procedeu-se à análise de variância multivariada pelo procedimento GLM e o comando Manova, tendo-se realizado o teste de Wilks para testar a hipótese de que os vetores de médias das diferentes combinações entre sistemas produtivos e estações eram nulos.

Em seguida, utilizou-se a análise de componentes principais (ACP), seguida da análise de agrupamento, com uso do 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 realizadas com o aplicativo SAS, versão 9.0 (SAS Institute Inc., Cary, NC, USA).

Resultados e Discussão

As estações do ano produziram efeito significativo sobre os teores de proteína e sólidos totais, e sobre a CCS (Tabela 3). O teor de proteína foi menor durante o verão, em razão da baixa qualidade nutricional das pastagens tropicais (Sorghum sudanense e Pennisetum glaucum), em comparação à das temperadas (Avena strigosa e Lolium multiflorum) (Fanti et al., 2008). Além disso, os animais, possivelmente, sofreram com o estresse térmico, uma vez que a temperatura máxima média da região onde se desenvolveu o estudo é de 31,2°C, o que limita o consumo de matéria seca e diminui o teor de proteína do leite (Fanti et al., 2008). Segundo Vargas et al. (2014), animais em ambientes com temperaturas elevadas apresentam menor capacidade de respostas a doenças, o que aumenta sua suscetibilidade à invasão do úbere por microrganismos e, no presente trabalho, resultou nos maiores valores de CCS no verão (Tabela 3).

Tabela 3:
Média±erro-padrão dos teores de proteína total, gordura, sólidos totais, contagem de células somáticas (CCS) e contagem bacteriana total (CBT) do leite bovino, nos sistemas de produção altamente especializado (AE), especializado (ES), semiespecializado (SE) e não especializado (NE), nas diferentes estações do ano⁽¹⁾.

• ⁽¹⁾
  Médias seguidas por letras iguais, minúsculas nas colunas e maiúsculas nas linhas, não diferem pelo teste Tukey, a 5% de probabilidade.

A elevação do conteúdo de sólidos totais no outono deveu-se, principalmente, à elevação do conteúdo de proteína e gordura do leite.

O aumento da CBT, com a diminuição dos níveis de especialização das UPL (Tabela 3), deveu-se, possivelmente, ao tipo e ao local de ordenha, bem como ao emprego ou não do resfriamento. Assim, o sistema altamente especializado, por possuir sala de ordenha própria, com circuito fechado e tanque de refrigeração a granel (Tabela 1), apresentou os menores valores de CBT.

Houve interação significativa entre os sistemas produtivos e as estações, para todos os ácidos graxos saturados, monoinsaturados e poli-insaturados (Tabelas 4, 5 e 6), bem como para os índices de qualidade da fração lipídica (Tabela 5).

Tabela 4:
Média±erro-padrão dos teores de ácidos graxos (mg g^-1) saturados e somatório dos ácidos graxos saturados (ΣAGS) no leite bovino, nos sistemas de produção altamente especializado (AE), especializado (ES), semiespecializado (SE) e não especializado (NE), nas diferentes estações do ano⁽¹⁾.

• ⁽¹⁾
  Médias seguidas por letras iguais, minúsculas nas colunas e maiúsculas nas linhas, não diferem pelo teste Tukey, a 5% de probabilidade.

Tabela 5:
Média±erro-padrão dos ácidos graxos (mg g^-1) monoinsaturados, somatório dos ácidos graxos monoinsaturados (ΣAGMI), atividade da enzima delta(9)-dessaturase, índice de aterogenicidade, índice de trombogenicidade, razão de ácidos graxos hipocolesterolêmico e hipercolesterolêmico (h/H), ácidos graxos desejáveis e índice nutracêutico do leite bovino, nos sistemas de produção altamente especializado (AE), especializado (ES), semiespecializado (SE) e não especializado (NE), nas diferentes estações do ano⁽¹⁾.

• ⁽¹⁾
  Médias seguidas por letras iguais, minúsculas nas colunas e maiúsculas nas linhas, não diferem, pelo teste Tukey, a 5% de probabilidade.

Tabela 6:
Média±erro-padrão ajustadas dos ácidos graxos (mg g^-1) poli-insaturados, somatório dos ácidos graxos ômega 3 (Σn3), dos ácidos graxos ômega 6 (Σn6), da relação ômega 6 e ômega 3 (Σn6:n3) e do total de ácidos graxos poli-insaturados (ΣAGPI) do leite bovino, nos sistemas produtivos altamente especializado (AE), especializado (ES), semiespecializado (SE) e não especializado (NE), nas distintas estações do ano⁽¹⁾.

• ⁽¹⁾
  Médias seguidas por letras iguais, minúsculas nas colunas e maiúsculas nas linhas, não diferem, pelo teste Tukey, a 5% de probabilidade.

Com a diminuição do grau de especialização dos sistemas produtivos, observou-se menor concentração do somatório dos ácidos graxos saturados (ΣAGS), nas estações primavera, verão e outono (Tabela 4). Isto se atribui à maior participação de forragens frescas (Chilliard et al., 2007) na dieta dos extratos menos especializados.

No inverno, observou-se menor concentração de ácido linoleico (18:2n6-c9,c12) e maiores valores de ácido rumênico e do somatório dos ácidos graxos poli-insaturados (ΣAGPI), nos extratos produtivos menos especializados (Tabela 6). Esses resultados estão relacionados ao predomínio de uso dos extratos menos especializados, de pastagens de aveia-preta (Avena strigosa) e azevém (Lolium multiflorum), em detrimento da silagem de milho e de concentrado, o que teria resultado em maiores concentrações de ácido alfa-linolênico (18:3n3-c9,c12,c15) do que de ácido linoleico e oleico (18:1n9-c9) na dieta (Tabela 2). Conforme Dewhurst et al. (2006), isto aumenta as concentrações de ácido rumênico e ΣAGPI e diminui a concentração de ácido linoleico no leite, o que corrobora os resultados obtidos no presente trabalho (Tabela 6).

O extrato produtivo mais especializado, apesar de não ter significativamente diferido dos dois extratos seguintes, aparentemente favoreceu a formação de ácido rumênico, a ponto de seus teores terem-se igualado estatisticamente aos observados no extrato com a menor especialização, que apresentou o maior valor desse ácido graxo (Tabela 6). Assim, pode-se inferir que a relação volumoso:concentrado, nos três primeiros níveis de especialização, influenciou esse resultado. No sistema produtivo altamente especializado, que utiliza a maior suplementação com concentrado na dieta (Tabela 1), os limites da relação volumoso:concentrado foram comprometidos, o que segundo Chilliard et al. (2007) diminui a bio-hidrogenação ruminal dos ácidos graxos de cadeia longa e inibe a etapa final desse processo, em que o ácido transvacênico (18:1n7-t11) é convertido a ácido esteárico (18:0). Assim, após a dessaturação do ácido transvacênico na glândula mamária pela enzima delta(9)-dessaturase, os teores de ácido rumênico no leite aumentam (Tabela 6).

Os extratos altamente especializado e especializado não diferiram quanto aos teores de ácido rumênico; no entanto, nos sistemas semiespecializado e não especializado, os menores valores para esse ácido foram constatados no outono, e os maiores, no inverno e na primavera. Como nos extratos mais especializados o processo de formação desse ácido graxo depende principalmente da concentração de ácido linoleico na dieta, e as fontes dietéticas nesses manejos que o contêm em sua composição (silagem de milho e concentrado) (Tabela 2) não são tão suscetíveis a variações sazonais, não houve variação nos teores de ácido rumênico (CLA, 18:2n7 - c9, t11) entre esses extratos (Tabela 6).

Nos níveis de menor especialização, a formação de ácido rumênico depende da participação do ácido alfa-linolênico, principal constituinte das pastagens frescas na dieta (Tabela 2). Portanto, esses sistemas produtivos são mais suscetíveis a alterações climáticas, com maiores teores do ácido alfa-linolênico nas pastagens temperadas, que são exploradas até a primavera, e menores nas pastagens tropicais. Além disso, a menor disponibilidade forrageira no outono resulta em menores níveis dietéticos de ácido alfa-linolênico, com reflexo no conteúdo do ácido rumênico no leite (Tabela 6).

Portanto, os extratos produtivos altamente especializado e não especializado apresentam características importantes para o aumento do potencial nutracêutico do leite, uma vez que o ácido linoleico conjugado (CLA) está relacionado à redução de uma série de doenças degenerativas, como aterosclerose, à prevenção e tratamento do diabetes mellitus, além de apresentar efeito anticarcinogênico, ser estimulante do sistema imunológico, reduzir a pressão arterial e atuar no metabolismo lipídico corporal, especialmente na redução da gordura corporal e no aumento da massa magra (Kratz et al., 2013; Lahlou et al., 2014). No entanto, o principal antecessor metabólico do ácido rumênico em ruminantes, o ácido transvacênico, por sua conformação geométrica (trans) poderia exercer ação prejudicial à saúde; contudo, Haug et al. (2007) constataram que, em humanos, este ácido graxo pode ser convertido ao seu sucessor.

As principais diferenças observadas no presente estudo estão relacionadas ao comportamento do ácido linoleico e do ácido alfa-linolênico nos diferentes sistemas de produção e estações do ano. As alterações observadas quanto aos conteúdos destes ácidos graxos no leite estão ligadas a modificações no conteúdo de seus sucessores - respectivamente o ácido araquidônico (20:4n6-c5,c8,c11,c14) e o timnodônico (20:5n3-c5,c8,c11,c14,c17) (Tabela 6). De acordo com Simopoulos (2006), a ação de enzimas dessaturases e elongases converte os sucessores em seus respectivos derivados. Entretanto, é importante ressaltar que essas duas famílias de ácidos graxos competem pelas mesmas enzimas envolvidas na dessaturação e no alongamento de suas cadeias (Simopoulos, 2006). Assim, conforme o autor, o balanço entre ômega 6 e ômega 3 (Σn6:n3), na dieta, assume grande importância na nutrição humana, tendo-se preconizado a diminuição dessa relação.

Os valores de Σn6:n3 diminuíram à medida que a especialização dos sistemas produtivos também diminuiu (Tabela 5). Como Simopoulos (2006) preconiza Σn6:n3 de 2:1 a 4:1, os melhores resultados foram obtidos no sistema não especializado.

O perfil de ácidos graxos mais apropriado para a alimentação humana foi obtido com a diminuição da especialização dos sistemas produtivos e nas estações com maior participação de gramíneas temperadas, ou seja, no inverno e na primavera. Assim, o nível não especializado, no inverno, apresentou o melhor índice nutracêutico do estudo (Tabela 5), pois, além de apresentar maior participação de forrageiras frescas na dieta, o perfil de ácidos graxos nas pastagens temperadas (Tabela 2) favorece a formação dos ácidos transvacênico, rumênico e do ácido t10,c12-octadecadienoico, além de diminuir a Σn6:n3 no leite (Roy et al., 2006; Chilliard et al., 2007).

O índice nutracêutico proposto mostrou maior poder de diferenciação entre os sistemas de produção, pois as diferenças no índice acentuaram-se entre os extratos mais e menos especializados, principalmente nas estações verão e inverno (Tabela 5). Assim, em estudos de manipulação da dieta que visem à obtenção de um perfil de ácidos graxos mais adequado para o leite de consumo, sugere-se o índice nutracêutico como parâmetro de avaliação da qualidade da fração lipídica do leite.

Além disso, o índice nutracêutico pode ser considerado adequado para analisar a capacidade do leite de contribuir na prevenção e combate de doenças cardiovasculares, uma vez que o nutracêutico apresentou correspondência inversa com os índices de aterogenecidade e trombogenicidade, propostos por Ulbricht & Southgate (1991), e relação direta com a razão entre ácidos graxos hipocolesterolêmico e hipercolesterolêmico (Santos-Silva et al., 2002) (Tabela 5). O índice também se mostrou mais eficiente na avaliação do potencial nutracêutico do leite do que o conteúdo de ácidos graxos desejáveis (Costa et al., 2008), que não apresentou correspondência com os índices estudados no presente trabalho.

Com a análise de variância multivariada, verificou-se diferença significativa entre os vetores de médias, para as diferentes combinações de sistemas produtivos e estações do ano, tanto pelo teste de Wilks, quanto pelos testes Pillai, Hotelling-Lawley e Roy. Os três primeiros autovalores foram significativos pelo teste da razão de verossimilhança; o primeiro autovalor explicou 70,72%, e o segundo, 15,51%, ou seja, os dois primeiros autovalores explicaram 86,23% da variação total dos dados relacionados ao potencial nutracêutico do leite bovino (Figura 1). Portanto, a análise de componentes principais permitiu visualizar a similaridade e dissimilaridade entre as diferentes combinações entre sistemas produtivos e estações.

Figura:1:
Projeção bidimensional dos escores das diferentes combinações entre os sistemas de produção - altamente especializado (AE), especializado (ES), semiespecializado (SE) e não especializado (NE) - e as diferentes estações do ano - primavera (PRIM), verão (VER), outono (OUT) e inverno (INV), além das cargas (correlações entre as variáveis e os componentes principais) das variáveis: ácido vacênico, total de ácidos graxos monoinsaturados (ΣAGMI), poli-insaturados (ΣAGPI), relação dos ácidos graxos ômega 6 e ômega 3 (N6/N3) e os isômeros de ácido linoleico conjugado (CLA), ácido rumênico e ácido t10,c12-octadecadienoico (18:2n6-t10,c12).

As variáveis ΣAGPI (r=0,69, p<0,01), ΣAGMI (r=0,78, p<0,001), Σn6:n3 (r=-0,74, p<0,01), ácido rumênico (r=0,97, p<0,001) e ácido t10,c12-octadecadienoico (r=0,93, p<0,001) explicaram a variabilidade do CP1 (componente principal 1), ao passo que as variáveis Σn6:n3 (r=0,55, p<0,05) e ΣAGPI (r=0,58, p<0,05) explicaram a variabilidade do CP2 (Figura 1). A partir das correlações dos componentes principais com as variáveis, pode-se observar que quanto mais à direita (variação explicada pelo CP1) e acima (variação explicada pelo CP2) as combinações se localizarem no plano cartesiano, maior será o potencial nutracêutico do leite, uma vez que essas posições apresentam valores mais elevados de ΣAGPI, ΣAGMI, ácido rumênico e ácido t10,c12-octadecadienoico, e menores de Σn6:n3; exceto pela variação explicada pelo CP2 em relação ao Σn6:n3, que apresentou maiores valores em posições mais elevadas.

A análise de agrupamento mostrou que foi quatro o número ideal de grupos formados entre as 16 possíveis combinações entre sistemas produtivos e estações do ano, em que: o grupo I englobou a classe inverno-não especializado; o grupo II, as classes primavera-semiespecializado, primavera-não especializado, verão-semiespecializado, verão-não especializado, inverno-altamente especializado, inverno-especializado e inverno-semiespecializado; o grupo III, as classes primavera-altamente especiali-zado, primavera-especializado, verão-altamente especializado, verão-especializado e outono-altamente especializado; e o grupo IV, as classes outono-especializado, outono-semiespecializado e outono-não especializado (Figura 1).

Portanto, no outono, todos os sistemas produtivos situaram-se no grupo IV, exceto o altamente especializado, que se localizou no grupo III; esse extrato produtivo foi o único capaz de mitigar os efeitos negativos da carência de forragens no outono sobre a qualidade da fração lipídica do leite, o que culminou com um aumento de ΣAGPI, ΣAGMI, ácido rumênico e ácido t10,c12-octadecadienoico, variação esta explicada pelo CP1 (Figura 1).

Entretanto, é importante ressaltar que, por meio da análise das variáveis que explicaram a variabilidade do CP2, observou-se aumento na Σn6:n3 no nível altamente especializado, no outono, em comparação aos outros sistemas produtivos. Isto se atribui, segundo Chilliard et al. (2007), ao incremento do fornecimento de fontes dietéticas com elevados teores de ácido linoleico, como a silagem de milho e o concentrado (Tabela 2), o que elevaria a absorção intacta desse principal representante da família ômega 6 no leite (Figura 1).

Além da combinação outono-altamente especializado, outras quatro assentaram-se no grupo III, compostas pelos sistemas produtivos altamente especializado e especializado, nas estações verão e primavera (Figura 1). Porém, nota-se que as combinações dessas estações com os níveis de menor especialização migraram para o grupo II. Essa migração pode ser explicada pelas diferenças no manejo nutricional entre os diferentes extratos produtivos. Os sistemas mais especializados necessitaram desocupar áreas antes destinadas ao pastejo, para destiná-las ao plantio do milho, entre a primavera e o verão, o que resultou em diminuição do fornecimento de forrageiras frescas, em razão do aumento da oferta de silagem de milho e concentrado. Isso manteve esses sistemas de produção no grupo III (Figura 1), principalmente em razão dos maiores valores de Σn6:n3 (variação explicada pelo CP1 e CP2) e menores de ΣAGPI, ΣAGMI, ácido rumênico e ácido t10,c12-octadecadienoico, em comparação ao grupo II (variação explicada principalmente pelo CP1).

Esse manejo nutricional, entretanto, não ocorreu em propriedades dos extratos produtivos menos especializados, que apresentaram menor ou nenhum fornecimento de silagem de milho na dieta (Tabela 1). O extrato semiespecializado, que também utiliza silagem de milho e concentrado, mas em menores proporções, apesar de situar-se no grupo II, deslocou-se em direção ao grupo III (Figura 1), o que corrobora a explicação anterior.

No grupo II, além das combinações dos sistemas semiespecializados e não especializados com as estações primavera e verão, ocorreram todas as combinações dos extratos produtivos com a estação inverno, exceto para a combinação inverno-não especializado, que foi o único representante do grupo I, caracterizado pelo leite de maior potencial nutracêutico. Essa variação foi explicada principalmente pelo CP1, já que esses grupos diferiram minimamente quanto ao CP2.

Conclusões

1. Os sistemas de produção e as estações do ano interferem conjuntamente no perfil de ácidos graxos do leite bovino, mas isoladamente em sua qualidade físico-química e microbiológica.

2. As maiores contagens de células somáticas e os menores valores de proteína são obtidos no verão e, com a diminuição do grau de especialização das unidades produtoras, a contagem bacteriana total aumenta.

3. Os leites bovinos de maior potencial nutracêutico são obtidos no inverno, nos sistemas não especializados de produção, que apresentam os maiores teores de ácidos graxos poli-insaturados, ácido rumênico e ácido t10,c12-octadecadienoico.

4. O sistema de produção não especializado, no inverno, apresenta os melhores índices de aterogenicidade, trombogenicidade, razão de ácidos graxos hipocolesterolêmicos e hipercolesterolêmicos e, em especial, o melhor índice nutracêutico.

5. Em razão de sua alta capacidade em diferenciar alterações no manejo nutricional dos animais, o índice nutracêutico deve ser preferido para a avaliação da qualidade da fração lipídica do leite bovino.

Agradecimentos

À Fundação de Amparo à Pesquisa do Estado do Rio Grande do Sul (Fapergs), à Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Capes) e ao Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), pelo apoio e auxílio financeiro prestados.

Referências

• BLIGH, E.G.; DYER, W.J. A rapid method of total lipid extraction and purification. Canadian Journal of Biochemistry and Physiology, v.37, p.911-917, 1959. DOI: 10.1139/o59-099.
  » https://doi.org/10.1139/o59-099.
• BODENMÜLLER FILHO, A.; DAMASCENO, J.C.; PREVIDELLI, I.T.S.; SANTANA, R.G.; RAMOS, C.E.C. de O.; SANTOS, G.T. dos. Tipologia de sistemas de produção baseada nas características do leite. Revista Brasileira de Zootecnia, v.39, p.1832-1839, 2010. DOI: 10.1590/S1516-35982010000800028.
  » https://doi.org/10.1590/S1516-35982010000800028.
• CHILLIARD, Y.; GLASSER, F.; FERLAY, A.; BERNARD, L.; ROUEL, J.; DOREAU, M. Diet, rumen biohydrogenation and nutritional quality of cow and goat milk fat. European Journal of Lipid Science and Technology, v.109, p.828-855, 2007. DOI:10.1002/ejlt.200700080.
  » https://doi.org/10.1002/ejlt.200700080
• CHRISTIE, W.W. A simple procedure for rapid transmethylation of glycerolipids and cholesteryl esters. Journal of Lipid Research, v.23, p.1072-1077, 1982.
• COSTA, R.G.; MESQUITA, Í.V.U.; QUEIROGA, R. de C.R. do E.; MEDEIROS, A.N. de; CARVALHO, F.F.R. de; BELTRÃO FILHO, E. Características químicas e sensoriais do leite de cabras Moxotó alimentadas com silagem de maniçoba. Revista Brasileira de Zootecnia, v.37, p.694-702, 2008. DOI: 10.1590/S1516-35982008000400016.
  » https://doi.org/10.1590/S1516-35982008000400016.
• DEWHURST, R.J.; SHINGFIELD, K.J.; LEE, M.R.F.; SCOLLAN, N.D. Increasing the concentrations of beneficial polyunsaturated fatty acids in milk produced by dairy cows in high-forage systems. Animal Feed Science and Technology, v.131, p.168-206, 2006. DOI: 10.1016/j.anifeedsci.2006.04.016.
  » https://doi.org/10.1016/j.anifeedsci.2006.04.016.
• FANTI, M.G.N.; ALMEIDA, K.E. de; RODRIGUES, A.M.; SILVA, R.C. da; FLORENCE, A.C.R.; GIOIELLI, L.A.; OLIVEIRA, M.N. de. 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, v.28, p.259-265, 2008. DOI: 10.1590/S0101-20612008000500039.
  » https://doi.org/10.1590/S0101-20612008000500039.
• GABBI, A.M.; McMANUS, C.M.; SILVA, A.V.; MARQUES, L.T.; ZANELA, M.B.; STUMPF, M.P.; FISCHER, V. Typology and physical-chemical characterization of bovine milk produced with different productions strategies. Agricultural Systems, v.121, p.130-134, 2013. DOI: 10.1016/j.agsy.2013.07.004.
  » https://doi.org/10.1016/j.agsy.2013.07.004.
• HAUG, A.; HØSTMARK, A.T.; HARSTAD, O.M. Bovine milk in human nutrition - a review. Lipids in Health and Disease, v.25, p.1-16, 2007. DOI: 10.1186/1476-511x-6-25.
  » https://doi.org/10.1186/1476-511x-6-25.
• KRATZ, M.; BAARS, T.; GUYENET, S. The relationship between high-fat dairy consumption and obesity, cardiovascular, and metabolic disease. European Journal of Nutrition, v.52, p.1-24, 2013. DOI: 10.1007/s00394-012-0418-1.
  » https://doi.org/10.1007/s00394-012-0418-1.
• LAHLOU, M.N.; KANNEGANTI, R.; MASSINGILL, L.J.; BRODERICK, G.A.; PARK, Y.; PARIZA, M.W.; FERGUSON, J.D.; WU, Z. Grazing increases the concentration of CLA in dairy cow milk. Animal, v.8, p.1191-1200, 2014. DOI: 10.1017/S1751731114000998.
  » https://doi.org/10.1017/S1751731114000998.
• ROY, A.; FERLAY, A.; SHINGFIELD, K.J.; CHILLIARD, Y. Examination of the persistency of milk fatty acid composition responses to plant oils in cows given different basal diets, with particular emphasis on trans-C18:1 fatty acids and isomers of conjugated linoleic acid. Animal Science, v.82, p.479-492, 2006. DOI: 10.1079/asc200658.
  » https://doi.org/10.1079/asc200658.
• SANTOS-SILVA, J.; BESSA, R.J.B.; SANTOS-SILVA, F. Effect of genotype, feeding system and slaughter weight on the quality of light lambs: II. Fatty acid composition of meat. Livestock Production Science, v.77, p.187-194, 2002. DOI: 10.1016/S0301-6226(02)00059-3.
  » https://doi.org/10.1016/S0301-6226(02)00059-3.
• SHOOK, G.E. Genetic improvement of mastitis through selection on somatic cell count. Veterinary Clinics of North America: Food Animal Practice, v.9, p.563-581, 1993.
• SIMIONATO, J.I.; GARCIA, J.C.; SANTOS, G.T. dos; OLIVEIRA, C.C.; VISENTAINER, J.V.; SOUZA, N.E. de. Validation of the determination of fatty acids in milk by gas chromatography. Journal of the Brazilian Chemical Society, v.21, p.520-524, 2010. DOI: 10.1590/S0103-50532010000300018.
  » https://doi.org/10.1590/S0103-50532010000300018.
• SIMOPOULOS, A.P. Evolutionary aspects of diet, the omega-6/omega-3 ratio and genetic variation: nutritional implications for chronic diseases. Biomedicine & Pharmacotherapy, v.60, p.502-507, 2006. DOI: 10.1016/j.biopha.2006.07.080.
  » https://doi.org/10.1016/j.biopha.2006.07.080.
• TRIPATHI, M.K. Effect of nutrition on production, composition, fatty acids and nutraceutical properties of milk. Advances in Dairy Research, v.2, p.1-11, 2014. DOI: 10.4172/2329-888X.1000115.
  » https://doi.org/10.4172/2329-888X.1000115.
• ULBRICHT, T.L.V; SOUTHGATE, D.A.T. Coronary heart disease: seven dietary factors. Lancet, v.338, p.985-992, 1991. DOI: 10.1016/0140-6736(91)91846-M.
  » https://doi.org/10.1016/0140-6736(91)91846-M.
• VARGAS, D.P. de; NÖRNBERG, J.L.; MELLO, R. de O.; SHEIBLER, R.B.; MILANI, M.P.; MELLO, F.C.B. Correlações entre contagem bacteriana total e parâmetros de qualidade do leite. Revista Brasileira de Ciência Veterinária, v.20, p.241-247, 2013.
• VARGAS, D.P. de; NÖRNBERG, J.L.; MELLO, R. de O.; SHEIBLER, R.B.; BREDA, F.C.; MILANI, M.P. 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, v.15, p.473-483, 2014. DOI: 10.4322/rbcv.2014.009.
  » https://doi.org/10.4322/rbcv.2014.009.

Datas de Publicação

• Publicação nesta coleção
  Dez 2015