Vírus da diarreia viral bovina (bvdv) em rebanhos leiteiros: um estudo de caso-controle pareado e estratégias de construção de modelos

Abstract

O vírus da diarreia viral bovina (BVDV) causa uma das doenças mais importantes de bovinos em termos de custos econômicos e sociais, uma vez que é largamente disseminado na população de gado leiteiro. Os objetivos do trabalho foram estimar a prevalência em nível de rebanho e investigar fatores associados aos níveis de anticorpos em leite de tanque através de um estudo de caso-controle pareado, bem como discutir estratégias de construção de modelos. Para estimar a prevalência de rebanho, amostras de tanque de leite foram selecionadas aleatoriamente (n = 314) de uma população (N = 1604). A prevalência real de BVDV foi de 24,3% (IC95% = 20,1-29,3%). Para o estudo de caso-controle, rebanhos positivos para BVDV (altos níveis de anticorpos) foram classificados como casos (n = 21) e pareados (n = 63) por produção de leite com rebanhos que apresentaram baixos títulos de anticorpos (razão 1:3). Para análise, três modelos multivariáveis foram construídos: 1) modelo completo, onde todas as 21 variáveis independentes foram oferecidas, e dois modelos foram criados de acordo com conhecimento empírico e similaridades entre as variáveis independentes, 2) modelo de fatores animais e 3) modelo de biossegurança. Um questionário foi aplicado (n = 84) para obtenção de informações a respeito de possíveis fatores de risco para BVDV. O modelo completo (Modelo 1) identificou as seguintes variáveis: idade com critério de eliminação (OR = 0,10; IC95% = 0,02 – 0,39; P < 0,01); propriedades que forneceram leite a outras cooperativas anteriormente (OR = 4,13; IC95% = 1,17 - 14,49; P = 0,02) e a presença de piquete de isolamento para animais doentes (OR = 0,14; IC95% = 0,01 – 0,26; P = 0,02). O modelo de biossegurança (Modelo 3) revelou uma associação significativa com o uso de monta natural (OR = 9,03; IC95% = 2,14 – 38,03; P < 0,01); presença de piquete de isolamento para animais doentes (OR = 0,06; IC95% = 0,05 – 0,83; P = 0,03); anos fornecendo leite para a mesma cooperativa (OR = 0,94; IC95% = 0,91 – 0,97; P < 0,01) e contato direto pela cerca entre bovinos de propriedades vizinhas (OR = 5,78; IC95% = 1,41 – 23,67; P = 0,04). O modelo de biossegurança pode ser considerado o “melhor”, pois obteve AIC = 43,880 e BIC = 48,058 menores quando comparado com o modelo completo que obteve AIC = 50,445 e BIC = 53,779. Esta dissertação tem a intenção de promover a discussão sobre as estratégias de construção de modelos especialmente quando se trata de saúde animal. Recomenda-se a aplicação de agrupamento de variáveis independentes como uma boa alternativa na construção de modelos, uma vez que este processo pode levar a uma melhor compreensão a respeito da associação entre cause e efeito de doenças.

Bovine viral diarrhea virus (BVDV) causes one of the most important diseases of cattle in terms of economic costs and welfare, since it is widespread in the dairy cattle population. The aims were to estimate herd prevalence and investigate factors associated with antibodies in bulk tank milk (BTM) in dairy herds through a matched case-control design as well as discuss model-building strategies. To estimate herd prevalence, BTM samples were randomly selected (n = 314) from a population (N = 1604). The true prevalence of BVDV was 24.3% (CI95% = 20.1 - 29.3%). For the case-control study, BVDV antibody-positive herds (high antibody titers) were classified as cases (n = 21) and matched (n = 63) by milk production with herds presenting low antibody titers (ratio of 1:3). For analysis, three multivariable models were built: 1) full model, holding all 21 independent variables; and two models divided according to empirical knowledge and similarity among independent variables, i.e., 2) animal factor model and 3) biosecurity model. A questionnaire was applied (n = 84) to get information about possible BVDV risk factors. The full model (model 1) identified the following variables: age as a culling criteria (OR = 0.10; IC95% = 0.02 – 0.39; P < 0.01); farms that provided milk to other industries previously (OR = 4.13; IC95% = 1.17 – 14.49; P = 0.02); and isolation paddocks for ill animals (OR = 0.14; IC95% = 0.01 – 0.26; P = 0.02). The biosecurity model (model 3) revealed a significant association with the use of natural mating (OR = 9.03; IC95% = 2.14 – 38.03; P < 0.01); isolation paddocks for ill animals (OR = 0.06; IC95% = 0.05 – 0.83; P = 0.03); years providing milk for the same industry (OR = 0.94; CI95% = 0.89 - 0.99; P = 0.02); and direct contact over fences among cattle of neighboring farms (OR = 5.78; IC95% = 1.41 – 23.67; P = 0.04). The biosecurity model could be considered the “best” since AIC = 43,880 and BIC = 48,058 when compared with the full model where AIC = 50,445 and BIC = 53,579. This paper intends to promote discussion about the model-building strategy when animal-health-modeling is on the line. We recommend the application of grouping predictors as a good choice for model building since it could lead to a better understanding of disease-exposure associations.