Abstract

A cross-sectional epidemiological study was conducted to determine the seroprevalence and to identify risk factors for bovine brucellosis seropositivity in traditional and smallholder dairy cattle production systems in the Tanga region of North-eastern Tanzania. The study populations comprised 246 indigenous and 409 crossbred cattle, randomly selected from 105 smallholder dairy and 25 traditional managed herds, respectively. Individual animal and herd-level data were collected using a structured questionnaire. Serum samples were screened for Brucella antibodies using the Rose Bengal Plate Test The overall seroprevalence of Brucella antibodies in the smallholder dairy and traditional managed cattle was 4.1% and 7.3% respectively. The corresponding overall herd prevalence was 10.5% and 20% respectively. Using multivariate logistic regression analysis, closeness to stock route, access to surface drinking water and location were identified as the major risk factors for individual herd seroprevalence. Older animals ( 6 years) were associated with increased risk of sero-positivity compared to animals of age category of 6 years. The results showed that brucellosis is prevalent and widely distributed locally, underscoring the need for further studies including surveillance and institution of preventive and control measures particularly among female young-stock and the general public who are at high risk of contracting brucellosis.

1. Introduction

Tanzania has a wide variety of livestock production systems influenced by climate, feed availability, culture, and farming systems [1]. Smallholder farmers in pastoral and agropastoral production systems account for about 9 9 % of the total livestock population and support the livelihoods of approximately 8 0 % of the population [2]. Under these production systems, livestock grazing range from communal, extensive grazing to stall fed or tethered grazing on smallholdings [3]. Pastoral livestock production system is dominated by indigenous traditional herds whereas agro-pastoralism comprises a range of combination of low-scale crop cultivation and improved or graded livestock rearing. Smallholder dairy farming is concentrated in the high-potential rural areas (highlands) and in many urban and periurban areas of major cities where milk marketing opportunities is high. The predominant production system in Tanga region is the traditional livestock keeping and small-scale dairy production [2].

Brucellosis is considered by the Food and Agriculture Organisation (FAO), the World Health Organisation (WHO), and the Office International des Epizooties (OIE) as one of the most widespread zoonoses in the world [4]. In cattle, the disease is usually caused by Brucella abortus (a gram-negative, facultative intracellular coccobacilli bacterium) and occasionally by B. melitensis and B. suis. Brucellosis is characterised by late term abortion; infertility and reduced milk production as a result of retained placenta, endometritis, and a varying degree of sterility in the males and cows [5].

The history of brucellosis in Tanzania began in 1927 when an outbreak of abortion in cows was reported in Arusha region [6, 7]. Surveys have shown the disease to occur in cattle in various production systems, regions, and zones, with seroprevalence varying considerably [811]. The majority of these studies, which were often conducted purposely, were carried out in parastate farms and in indigenous traditional cattle herds. Surveys carried out in the Eastern zone and in the dairy sector in Tanzania revealed prevalence ranges from 2 . 2 % in small scale (1–10 animals) to 7 . 6 % in large-scale ( > 30 animals) and 1 2 . 3 % in pastoral sector, respectively [8, 9]. Reports from individual dairy cattle in Northern Tanzania showed low prevalence ( < 4%) compared to traditional cattle ( 1 5 % ) in the same zone [6]. In Tanga region, the disease has been insufficiently investigated and information relating to its magnitude, distribution, and risk factors is scant. Such information is important when designing appropriate strategies that would help reduce its prevalence and effects. The present study was carried out in order to gather additional information likely to contribute to data that may be used to devise appropriate national strategies for the control of the disease.

The objective of the serological survey was to determine the prevalence of bovine brucellosis and to identify the associated risk factors under traditional and smallholder dairy production system in the Tanga region, North-eastern Tanzania.

2. Materials and Methods

2.1. Area and Study Population

The study was conducted on both smallholder dairy herds ( 10 graded animals of all ages, breed, and sex and intensively managed) and traditional herds ( 30 indigenous cattle of all ages and sex and extensively managed) in and around Tanga municipality. Characteristics of the study area are described elsewhere [12, 13]. The location of each herd was also classified as urban (within the official boundary of Tanga municipality), periurban (beyond the official boundary of Tanga municipality but within 1 5 k m of the town centre), or rural and grazing system (grazing versus no-grazing). Such classification of herds by location and grazing system could influence, for example, availability of veterinary input supplies, extension services, affecting access to land for forage production, veterinary input suppliers, or animal health service provisions. Herd classifications (by grazing and location) were considered as explanatory variables during data analysis. The type of animals kept under smallholder herds includes taurine breed (Friesian, Ayrshire, Jersey, Simmental) and crosses of these breeds with Bos indicus breeds (Tanzania shorthorn zebu, boran, and Sahiwal). The level of taurine blood varies from 50%–85%. Traditional herds comprise mainly Tanzania short horn zebu (TSHZ).

2.2. Study Design and Farms Selection

A cross-sectional study was undertaken between May 2003 and January 2004. The data bases of farmers under the district livestock department and Tanga dairy development programme [14] were used as the sampling frame. A sampling frame of 1 , 7 3 0 smallholder dairy and traditional herds, comprising 1 2 , 5 0 0 cattle, was used to select 1 3 0 herds to participate in the study. Owing to the practical consideration of logistics and funds, the sampling frame was limited to herds within a radius of 4 0 k m around Tanga municipality.

Simple random sampling using the Excel software (Microsoft Inc., 1999) was used to select the 1 3 0 herds which resulted in the recruitment of 1 0 5 smallholder dairy and 2 5 traditional herds. The sample size of herds and animals was selected to estimate the disease prevalence and association between outcome and explanatory variables at 9 0 percent confidence, with 7 percent precision based on formulae given by Noordhuizen et al. [15].

2.3. Questionnaire Design and Data Collection

Information about each herd and the animals kept was collected by means of a structured questionnaire, which was completed at all the selected herds on a single visit. The questionnaire was designed to comprise mostly closed ended (categorical) questions to ease data processing, minimize variation, and improve precision of responses [16]. The questionnaire was administered using the national Swahili dialect by a veterinary department staff member, who was trained in participatory research methodologies. Important herd and animal level data recorded included cattle location, source of forages (home established, road side, or bought-in hay), sex, breeding method used (natural or artificial insemination), source of drinking water (tap, rain water, shallow well, river, pond), contact with other animals and place (between farms, during grazing, at water source, during mating), herd owner education level (binary variable: illiterate and primary education), and the number of years in livestock farming. Other information sought included presence of other animals on the farm (pigs, sheep, goat), history of abortion and S.19 vaccination, disposal of afterbirth, age determined from birth records and dentition characteristics, categorized as (binary variable; 6 and 6 years), herd-size (binary variable; large 9 and low 9 animals), and the type of floor in the animal house (concrete, dirt) as well as whether or not a system of grazing or zero-grazing was practiced.

2.4. Sample Collection and Handling

Approximately 1 0 m L of blood sample was collected from the jugular vein of each animal in all selected herds using plain vacutainer tube (Becton Dickson, UK). Each sample was labelled using codes describing the specific animal and herd. The tube was set tilted on a table over night at a room temperature to allow clotting. Next morning, the clotted blood in the tubes was centrifuged (at 3 0 0 0 g for 2 0 m i n ) to obtain clear serum. The obtained serum was stored at 2 0 C until tested by Rose Bengal Plate Test (RBPT).

2.5. Rose Bengal Plate Test

All sera samples were screened using RBPT antigen (VLA Weybridge, UK). The test procedure recommended by Alton et al. [17] was followed. Briefly, 3 0 𝜇 L of RBPT antigen and 3 0 𝜇 L of the test serum were placed alongside on the plate, and then mixed thoroughly. The plate was shaken for 4 m i n and the degree of agglutination reactions was recorded. The sample was classified positive if any agglutination was observed and negative if no agglutination. The RBPT, when compared to complement fixation test (CFT), has shown a sensitivity of 9 4 . 2 % and a specificity of 8 7 % on field sera and has been described by other researchers [18, 19]. Confirmation of positives samples with tests of higher sensitivities and specificities such as a CFT or enzyme-linked immunosorbent assay (ELISA) was not done due to the lack of resources (funds) to buy the required kits.

2.6. Statistical Analysis

Herds and individual animal-derived data were stored in Microsoft Access. Descriptive statistics for the animal and herd level explanatory variables examined in the study were developed using Epi-Info version 6 . 0 4 d [20] and Statistix version 8 . 0 . [21]. Relationships between explanatory (independent) variables (herd and animal-level) and outcome variables (Brucella serostatus: negative or positive) were investigated in two steps by logistic regression. In the first step, relationships between each independent and outcome variable were individually investigated. In the second step, any variables that were significantly associated at the 𝑃 < . 1 0 level were included in multivariable models producing, by forwards and backwards substitution and elimination, the most parsimonious models in which all independent variables remained significant at the 𝑃 < . 0 5 level.

3. Results

3.1. Participating Herds Characteristics

All 1 3 0 of the selected herds were visited and the farmers were interviewed, resulting in a 1 0 0 % response rate for participation in the study. Smallholder dairy herds made up 8 0 . 7 % of the herds sampled and traditional herds 1 9 . 3 % . Samples were collected from 6 5 5 animals of which 7 % were males and 9 3 % females. The average age of the animals in the traditional sample was 5 . 7 years compared to 4 . 6 years in the smallholder dairy sample ( 𝑃 < . 0 0 1 ) . The distribution of herds amongst categories of each variable investigated is summarised in Table 1. All traditional herds were grazing animals whereas the majority ( 7 2 % ) of the smallholder dairy herds was kept under zero-grazing conditions. Smallholder dairy herds had an average herd size of six animals, with larger herds being kept under grazing than under zero-grazing conditions ( 𝑃 < . 0 0 1 ) . Traditional herds had an average herd size of 4 6 animals (median 3 0 ). Goats and sheep were significantly kept more commonly in traditional herds than in smallholder dairy herds ( 𝑃 < . 0 0 1 ) . Traditional livestock keepers keep cattle for longer than small-holder farmers, with an average of 1 9 . 1 (median 1 4 ) compared with 9 . 4 (median 9 ) years, respectively. Forty-percent ( 𝑛 = 4 2 ) of the smallholders dairy farmers had more than primary education compared to 1 2 % in the traditional sector. Two-thirds ( 6 6 % ) of the smallholder dairy farmers had at least received training in cattle management through workshop or seminars whereas traditional farmers learned cattle management alone or from parents ( < 0.001). No history of S.19 vaccination was recorded during sampling. Abortion history was recorded in 1 2 % and 3 . 7 % of the herds in traditional and smallholders dairies, respectively.

3.2. Serological Responses to Brucella Infection

Of the 6 5 5 sampled animals, serology results were available from 6 5 4 ( 9 9 . 8 % ) animals. The missing result arose due to loss of labels during storage and transport to laboratories. Over all, animal seroprevalence for Brucella antibodies was 5 . 3 % ( 9 5 % Confidence interval [CI], 3.1–7.8). The animal and herd-level seroprevalence of antibodies to Brucella by production system are shown in Table 2. The corresponding overall herd seroprevalence (at least one animal or herd seropositive) for Brucella antibodies was 1 2 % ( 9 5 % CI, 10.1–13.9). The relationship between age and seroprevalence is shown in Figure 1. In the traditional herds, the prevalence of cattle seropositive for Brucella antibodies increased from 4 % in animals < 4 year-old to 1 4 % in cattle > 6 years of age, while in smallholder dairy herds, the prevalence of cattle seropositive for Brucella antibodies increased from 2 % in young stock < 2 years of age to 1 6 % in stock of > 6 years.

4. Factors Associated with Seropositivity to Brucellosis

4.1. Univariable Analysis

The association of the animal and herd-level categorical explanatory variables at 𝑃 . 1 0 and brucellosis is shown in Table 3. Animal-level categorical variables that qualified at 𝑃 . 1 0 during univariable analysis were age and the introduction of animals during the last two years prior to the current study. Herd-level variables that were significant at 𝑃 . 1 0 during univariable analysis included grazing, herd location, access to surface water, and closeness to cattle stock route or holding ground.

4.2. Multivariable Analysis

Explanatory variables that remained significant for the seropositivity to brucellosis, defined by a sero-positive animal and/or herd in the final multivariable regression model are shown in Table 4. The likelihood that a herd was sero-positive for Brucella increased significantly for being near to the slaughter cattle stock route ( 𝑃 < . 0 0 1 ) . Herds that were accessed to surface water were significantly associated with increased risk compared to the herds which depended on piped water ( O R = 2 . 6 6 , 𝑃 < . 0 0 1 ). Seroprevalence to Brucella varied significantly amongst herd locations, with herds located in urban areas being significantly less likely to be sero-positive compared to rural or periurban located herds ( O R = 0 . 3 7 , 𝑃 < . 0 0 1 for urban). At an animal-level, older animals ( 6 years) were more likely to be seropositive compared to animals 6 years old ( O R = 4 . 0 , 𝑃 < . 0 0 1 ). Seroprevalence varied significantly amongst herd source of drinking water and grazing system but these were strongly confounding significant factors with access to surface drinking water providing the minimal best model fit ( L i k e l i h o o d r a t i o s t a t i s t i c = 1 4 . 1 , 𝑃 = . 0 4 9 ) in the multivariate analysis. Results of Fisher’s Exact Test showed that history of previous abortion ( 𝑃 > . 0 5 ) and stillbirths ( 𝑃 > . 0 5 ) in the individual animals were not significantly associated with Brucella seropositivity. None of the other investigated variables were associated with differences in prevalence values.

5. Discussion

The overall animal seroprevalence of brucellosis in this study was 5 . 3 % . The seroprevalence could reflect a past or present exposure to Brucella organisms, because vaccination against brucellosis has never been practiced in Tanga region, during the last 1 0 years [14]. The prevalence was higher in animals in traditional herds than smallholder dairy herds, being 7 . 3 % and 4 . 1 % , respectively. Also, at herd level, brucellosis was more prevalent in traditional herds, with 2 0 % showing at least one positive animal. Although animal prevalence is low, the high herd prevalence, especially in traditional herds, shows that brucellosis is common and locally widespread.

The low prevalence of brucellosis observed in the smallholder dairy animals when compared to traditional herds agrees with the observations made in other studies [8, 22, 23]. The difference is also likely to be explained by the small-size units and stall feeding that minimises contacts between herds and animals. However, the high proportion of seropositive animals in traditional herds, despite differences between herds, conformed to the results of a recent study in Tanzania which also reported higher prevalence of brucellosis in indigenous TSHZ cattle than in crossbred kept by smallholder dairy farmers [11]. The differences between traditional and crossbred animals are possibly attributed to increased contacts of infected herds/animals and noninfected ones in the indigenous traditional production system, as a result of communal grazing and watering, which become more apparent and acute during the dry period. Moreover, mixed herding and frequent contact with small ruminants and cattle could also be contributing factors to the occurrence of brucellosis.

At herd level, farming close to the stock route and access to surface water emerge as the most important risk factors. Slaughter cattle from outside the Tanga region are trekked along the stock route to be slaughtered at the abattoir. Animals are kept for some days before slaughter at a holding ground where they partly make use of the same communal grazing areas. The brucellosis prevalence in these animals is high [24] and direct, or indirect, contact via pasture exposes grazing animals to diseases carried by these animals. Such practices enhance the exposure potential, especially following abortions, through increased contact and common grazing field and watering points promoting transmission of Brucella organism [25].

The relevant data from East and Central Africa, and from other environments similar to those found in Tanzania, also indicate that Brucella antibodies can be quite widespread among the cattle production systems. For example, in the pastoral animal and herds RBPT-based investigations, seroprevalences of 0 . 7 7 % and 4 6 . 1 % , 9 . 3 % and 3 1 % , 1 4 . 1 % and 4 6 . 2 % have been reported in Ethiopia, Ghana, and Zambia, respectively, [2629]. The apparent geographical variation in the seroprevalence may reflect differences in the levels of natural immunity, management and husbandry practices employed, and sensitivities and specificities of the diagnostic methods used among researchers as well as genetic variation in disease resistance among the breeds maintained in the systems [27, 30].

In both production systems, the seroprevalence for Brucella antibodies increased with age, consistent with previous reports [31, 32]; younger animals are more resistant to primary infection and frequently clear infections, although latent infection occur [33]. The increased likelihood of Brucella force of infection in older animals ( 6 years) could be explained by increased susceptibility to infection due to under nutritional stress combined by the lowered immunity that developed following acute infection [5, 28].

Compared to traditional herds, animals in the smallholder dairy herds become positive at a slightly younger age. It has been described for brucellosis that some of the infected animals, do not become sero-positive until pregnancy [33]. This could partly explain this difference as the age at first calving for crossbred dairy cattle is lower than for traditional zebu cattle [34].

A history of previous abortions or stillbirths was not significantly associated with seropositivity, a finding which is in contrast with reports from other researchers [24, 26]. Fifteen percent of the RBPT positive animals reported a history of abortion, compared to 9 . 7 % of the RBPT-negative animals. This could be explained by the fact that causes of abortion such as vector and nonvector borne diseases which are prevalent in the study area were important factors rather than brucellosis [13, 35]. Similar observations were made by other investigators [5, 36]. Although RBPT is generally accepted and accredited by OIE [37] as the definite test for Brucella screening in many countries, the test has several limitations [38]. Apart from being highly sensitive especially in vaccinated animals, the test has the major drawback of producing false-negative reactors. False-negatives are usually obtained during early stages of incubation or immediate after abortion whereas false-positives occur due to the presence of I g M as a result of S.19 vaccination and colostral antibodies in young stock. The lower specificity is the reason that normally RBPT-positive animals are either confirmed with the complement fixation test (CFT) or enzyme-linked immunosorbent assay (ELISA).

The study further revealed that sex did not show significant association with the seroprevalence for Brucella antibodies. The absence of association with the seroprevalence is in agreement with that of other workers [39] and possibly due to the small sample size of the male cattle sampled.

The possible sources of bias in this study include the fact that herds fell into aggregates with large spaces between them. Furthermore, several of the brucellosis positive herds showed clear clustering. This might have lead to unreliable estimates of seroprevalence. This was overcome during study designing stage by incorporating random sampling procedures that would give a 9 0 % to 9 5 % probability of detection of a disease in a herd with a minimum intraherd prevalence of 1 0 % [40]. Sampling an equal number of herds from each production system that had large sampling frame differences might also have introduced bias. The sampled herds in smallholder dairy production represented 9 . 3 % ( 4 0 9 / 4 3 7 5 ) of the eligible population while in traditional indigenous cattle production it was 3 % ( 2 4 6 / 8 1 2 5 ) . This was overcome by random selection. Another possible source of bias in this study is the problem of some variables influencing and being associated with Brucella seropositivity. Confounding was addressed by using multivariable model. However, other confounders may potentially have resulted in residual confounding. These include grazing and animal access to surface drinking water.

6. Conclusion

This study reports that bovine brucellosis is prevalent and widely distributed locally. Seroprevalence was higher in traditional than smallholder dairy herds suggesting traditional herds to be at a higher risk of brucellosis, most likely due to the commingling with small ruminants and frequent contact at communal grazing and watering points. Animals were more likely to be seropositive if they were in traditional herds, accessed to surface running water and graze close to a trade/ slaughter animal stock route. This may reflect the increased likelihood of contacts with contaminated pasture and water from aborted materials, although history of abortions was not a significant risk factor. Consistent with this, the likelihood of encounter increased significantly and logically with age. Further studies are needed to understand the dynamics of transmission cycles and institution of preventive and control measures particularly among female young-stock, and to identify alternative management practices to replace those that are risk factors for animal and human infections.

Conflict of Interests Statement

None of the authors of this paper has a financial or personal relationship with other people or organization that could inappropriately influence or bias the content of the paper.

Acknowledgments

Field work and data collection were funded as part of study leave grant awarded to L. Schoonman. The authors are grateful to all participating farmers for their cooperation. The livestock field officers and technician at Tanga are thanked for their field and technical assistance. Thanks are extended to the Director of Veterinary Service, Tanzania, for permission to publish these studies.