Quantifying transmission of Mycobacterium avium subsp. paratuberculosis among group-housed dairy calves.

Johne's disease (JD) is a chronic enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), with control primarily aimed at preventing new infections among calves. The aim of the current study was to quantify calf-to-calf transmission of MAP among penmates in an experimental trial. Newborn Holstein bull calves (n = 32) were allocated into pens of 4, with 2 inoculated (IN) calves and 2 calves that were contact exposed (CE). Calves were group-housed for 3 months, with frequent collection of fecal and blood samples and tissue collection after euthanasia. The basic reproduction ratio (R0) was estimated using a final size (FS) model with a susceptible-infected model, based on INF-γ ELISA and tissue culture followed by qPCR. In addition, the transmission rate parameter (β) for new shedding events was estimated using a general linearized method (GLM) model with a susceptible-infected-susceptible model based on culture, followed by qPCR, of fecal samples collected during group housing. The R0 was derived for IN and CE calves separately, due to a difference in susceptibility, as well as differences in duration of shedding events. Based on the FS model, interferon-γ results from blood samples resulted in a R 0 IG of 0.90 (0.24, 2.59) and tissue culture resulted in a R 0 T of 1.36 (0.45, 3.94). Based on the GLM model, the R0 for CE calves to begin shedding (R 0 CE ) was 3.24 (1.14, 7.41). We concluded that transmission of MAP infection between penmates occurred and that transmission among calves may be an important cause of persistent MAP infection on dairy farms that is currently uncontrolled for in current JD control programs.

Introduction

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne’s disease (JD), a chronic enteritis primarily affecting ruminants and causing substantial losses to dairy industries worldwide [1-3]. There is currently no treatment, cure, or vaccine for the prevention of MAP infection in cattle; therefore, control is primarily based on preventing transmission of MAP and reducing new infections within the herd [2, 4–6].

The primary route of MAP transmission is fecal–oral through contaminated feed, milk, water and the environment caused by infectious animals intermittingly shedding MAP in their feces [7-9]. To decrease transmission from cows to calves and limit exposure of young stock to MAP, calves are removed from adult cows as soon as possible after birth and placed in calf barns or pens [10, 11]; however, calves up to 1 year of age have demonstrated susceptibility to MAP infection and calves can begin fecal shedding of MAP bacteria as early as 2 weeks after exposure [12-15]. Therefore, separating calves from cows and subsequent group housing may not be an effective method for prevention of new infections in young stock. Despite evidence that calf-to-calf transmission can occur [14, 16], findings regarding implications and impact of this route of transmission on control within herds were inconsistent [16-20]. Additionally, statistical analyses of data from 21 MAP infected farms suggested that transmission among calves was necessary to explain observed patterns of transmission within a herd [21].

Susceptible-infected-susceptible (SIS) models are modified from susceptible-infected-recovered (SIR) models; both have been used to model transmission of pathogens in a population [22-25]. Transmission of pathogens is quantified using the basic reproduction ratio (R0), the average number of new cases caused by 1 typically infectious individual introduced into a completely susceptible population [26, 27]. The threshold at which an outbreak can occur is when the R0 value is > 1, whereas an infection is certain to fade out in a population if R0 is < 1 [26]. Transmission dynamics for MAP infection are notoriously difficult to determine, due to long incubation, latently infected animals, variability of diagnostic tests and lack of long-term randomized control studies [28]. Deterministic and stochastic mathematical modelling techniques have been used to investigate MAP spread and transmission parameters [21, 29, 30], assess impact of varying infectious animals, as well as varying levels of environmental contamination [9, 31, 32], economic impacts of disease [33] and to determine effectiveness of control programs and interventions on spread and control of this disease [34, 35]. Due to complications regarding MAP infection, these models all rely on educated guesses regarding infectivity and susceptibility of animals in a herd, as well as impacts of various transmission routes. Increasing knowledge regarding transmission within a herd will enhance understanding of disease maintenance and spread within a herd and enable control programs to be optimised to better manage spread of infection [36]. Although transmission among calves is possible, sufficient quantitative information regarding amount of transmission is lacking. The aim of the current study was, therefore, to quantify calf-to-calf transmission of MAP infection among penmates and estimate R0 for calves identified as infected, or infectious, based on various diagnostic tests.

Materials and methods

Experimental design

Study design, calf collection and sample collection were as described [14]. Briefly, 32 newborn Holstein–Friesian bull calves were purchased from 13 Alberta (Canada) dairy farms that had tested negative for MAP for at least 4 years based on culture of environmental samples and milk ELISA or individual fecal sampling of the herd. Calves were assigned to infection status and pen based on birth order and entry into the Biosecurity Level-2 facility, with 7 experimental group pens each consisting of 2 inoculated calves (IN) and 2 recipient contact-exposed calves (CE) (Figure 1). The last 4 calves to enter the barn were designated controls. At time of group housing, IN calves were on average 4 weeks of age, whereas CE calves were 2 weeks of age in each experimental pen. Calves were group-housed in pens of 4 (2 IN, 2 CE) for 3 months after inoculation. After 3 months of group housing, IN calves were euthanized and remaining calves in experimental pens were individually housed (87 days after 1st day of group housing) for an additional 3 mo. At the end of the trial, all remaining calves were euthanized for tissue sampling. Control calves were group-housed for the entirety of the study. All protocols and the experimental design were approved by the University of Calgary Veterinary Sciences Animal Care Committee (protocol AC14-0168).

Figure 1

Experimental design for quantification of subsp. among group-housed calves. One of 7 experimental pens represented, containing 2 inoculated (IN) calves (grey circles) and 2 contact exposed (CE) calves (white circles). a) 2 newborn calves entered the barn and were housed in a clean environment for 2 weeks; b) IN calves were inoculated with 5 × 108 CFU’s of MAP at 2 weeks of age; c) 2 weeks of IN housing to allow pass-through of inoculum. During this time, 2 newborn calves entered the barn and were designated as CE and housed apart from IN calves; d) first day of group housing. IN and CE calves were co-mingled in a clean pen and remained group-housed for 3 months; e) IN calves removed for necropsy, CE individually housed for additional 3 months.

Inoculum

Inoculum preparation was as described [14]. A virulent strain from a clinical case of JD (cow 69) in Alberta was used for inoculation. Two calves in each experimental pen were inoculated with an oral dose of 2.5 × 108 CFU on 2 consecutive days at 2 weeks of age. To ensure shedding, this dose was higher than the minimum dose of 106 recommended for experimental infection and detection in tissues [37, 38]. Concurrently, this dose was also lower than previous experimental trials that used “high” doses greater than 109 [39], which are unlikely to be comparable to natural infectious doses and ensuing shedding. Following inoculation, calves remained individually housed for 2 weeks to allow passive shedding of the inoculum to cease before being group-housed with recipient calves.

Sampling

Fecal, blood and tissue samples were all collected as described [14]. Briefly, individual fecal samples were collected three times/week for 3 months, starting 1 day after the onset of group housing. Samples were stored at 4 °C and processed within 7 days after collection. All fecal samples were processed using modified TREK ESP II culture media (TREK para-JEM®; TREK Diagnostic Systems, Cleveland, OH, USA) as described [12, 14]. Following processing, fecal samples were incubated at 37 °C for 49 days.

Blood samples were collected weekly from the jugular vein, alternating between sides. Within 2 h after collection, samples were transported to the laboratory (in an insulated box with hot n class="Chemical">water bottles at 25–35 °C).

After 3 months of group housing (i.e., 4 months of age), IN calves were euthanized by intravenous injection of barbiturates (Euthanyl Forte®, DIN 00241326, Bimeda-MTC Animal Health Inc., Cambridge, ON, Canada), necropsied and tissue samples collected. After an additional 3 months of being housed individually, CE calves were euthanized, necropsied and tissues collected. Thirteen tissue samples were collected from each calf, including sections of the ileum, jejunum, duodenum, all associated lymph nodes and spleen. After transportation to the laboratory, samples were processed immediately for culture, as described [14]. Briefly, 2.5 g of tissue was dissociated and disinfected prior to incubation within paraJEM® culture bottles and incubated at 37 °C for 49 days.

Detection of MAP

Following liquid culture of fecal and tissue samples for 49 days, DNA was extracted as described [14, 40]. A duplex qPCR was performed targeting the MAP-specific F57 region and an internal amplification control (IAC) with primers, probes and IAC sequences identical to those described [41]. Samples were considered positive if the cycle threshold (CT) value was < 40.

Detection of an immune response

Following transportation to the laboratory, blood samples were processed for detection of IFN-γ release, as described [14, 42]. Each whole-blood sample was treated with 100 μL avium Purified Protein Derivative (aPPD; 0.3 mg/mL; Canadian Food Inspection Agency, Ottawa, ON, Canada), 100 μL of pokeweed mitogen (positive stimulation control; 0.3 mg/mL; Sigma-Aldrich Canada Co., Oakville, ON, Canada) and 100 μL sterile PBS (negative stimulation control). Following overnight incubation at 37 °C and centrifugation serum was collected and assayed using the BOVIGAM® sandwich ELISA (Prionics, La Vista, NE, USA). The %IFN-γ was calculated as follows: [(PPD Johnin-negative assay control)/(positive–negative assay control)] × 100 [42, 43].

Quantification of MAP transmission

GLM and FS model

The transmission rate parameter β represents the average number of new infections in a fully susceptible population caused by 1 typically infectious animal per unit of time [27]. The basic reproduction ratio (R0) represents the average number of new infections in a totally susceptible population caused by 1 typically infectious calf during its infectious period and was estimated using 2 approaches. The generalized linear model (GLM) was based on an SIS model (infectious status determined during the experiment based on fecal shedding) to quantify transmission parameter β (from which R0 can be derived). The GLM model is based on fecal shedding, which can only identify infectious animals and not infected animals. Therefore, it is a model quantifying the transmission parameter for a calf to shed MAP. The final size (FS) models were based on an SI model (infection status determined at end point, with no animals recovering) to estimate R0 directly [23], in which infection dynamics were based on number of recipient animals, i.e. susceptible for the FS, or not infectious for the GLM (S), infectious for the GLM (shedding), or infected for the FS (I) and total number of animals (N). Thus, in a classical model for the GLM analysis, the rate at which susceptible animals (S) are infected is given by βSI/N and the probability of a single susceptible calf becoming infected during a period is:where is the total transmission rate parameter [27]. Data were analyzed using STATA 11 (StataCorp LP, College Station, TX, USA). All confidence intervals refer to 95% confidence intervals.

General linearized model (GLM)

The GLM used a binomial distribution, with the dependent variable being number of new cases, which in this case are shedding events (C) and total number of susceptible (S) calves as the binomial total. The analysis was done with a complementary log–log (cloglog) link function, a binomial error term and an offset explained below [27, 44].

The expression for the GLM was:where is the intercept and thus the logarithm of the transmission parameter is the (intercept) regression coefficient and is the offset variable;  = expected number of cases (C) during the infectious interval (t, t + ) divided by the number of susceptible individuals (S) at the start of the time interval (i.e. at t); βT = total transmission rate parameter; MAPT = total amount of infectious material (animals and environment), also at the start of the time interval (t); Δt = duration of the time interval; and NT = total number of animals at the start of the time interval (t) as this measures the size of the area given that density is constant [45]. Infectious material in the environment, or shed by an infectious calf (I), were assumed to be equally infectious as the primary route of MAP transmission is fecal–oral. Thus on any given day, transmission occurs due to shedding calves and the amount of bacteria in the environment due to prior shedding in the pen.

Note that the βT is the total transmission rate parameter. Due to several potential routes of transmission, different transmission effects may contribute to βT and hence different transmission rate parameters can be estimated depending on the population composition, as discussed below. Additionally, infectious calves can either cause new infections through direct contact (IIN or ICE) with susceptible calves, or through environmental contamination caused by infectious calves (EIN or ECE), as discussed below. Therefore,

Note the definition of n class="Gene">MAPT implies the following:

The GLM analysis is based on interval data and was therefore based on fecal samples collected over the course of the trial using the SIS model. At any given time interval within each pen, infectious (I) individuals were those with MAP-positive fecal results, susceptible individuals (S) had MAP-negative fecal results and new cases (C) were calves with negative fecal culture results at the previous observation point that became infectious (I) in the current time interval. Changes in individual status were observed for each new time interval within each of the 7 experimental pens, either 2 or 3 days apart.

All time intervals where a pen had 0 susceptible (S) calves were removed from the analysis (n = 9). The 95% confidence interval (CI) of the estimated β parameters was calculated using the standard error of the mean of log β.

Estimation of MAP transmission through the environment: βE

Environmental contamination caused by fecal shedding calves was included in the model as a route of transmission for infection (expressed as parameter βE). On the first day of group housing, calves were introduced to a clean, MAP-free environment. Therefore, environmental contamination was the result of MAP accumulating in the environment from either IN or CE calves shedding. We assumed that environmental contamination on a specific day (Et) resulted from accumulation of MAP in the environment by infectious individuals (either IN or CE) on the previous day (I(t−1)), as well as the remaining exposed MAP in the environment from all anterior shedding events prior to t − 1 (Et−1). Environmental contamination caused by previous shedding events were discounted by the concealment rate (γ) for every day after the shedding event. Concealment and resulting exposure rate of MAP were estimated as described, based on survival rate definitions [22] and accounting for potential differences in contamination caused by IN or CE calves. Briefly, ET was calculated as the sum of environmental contamination caused by IN (EIN) and CE (ECE) calves, .

The same concealment rate (γ) was applied to both EIN and ECE and was calculated based on the assumption that fecal–oral transmission from the environment will be the same as fecal–oral transmission from an infectious pen mate on the same day [22]. Total oral transmission (IT = IIN + ICE) directly to other calves within the same time period (t) as the infectious material is shed, is basically the same as will occur from the equivalent but diminished amount of material in subsequent time periods. MAP remains viable in the environment for extended periods of time and the infectious dose is not known; it can, therefore, be assumed that a calf picking MAP up directly from fecal shedding of a pen mate (I) or from the environmental contamination (E) would have equal probability of getting infected. However, over time, MAP will be concealed (diminished) in the environment and exposure of the other animals decreased, resulting in lower transmission [22]. Using new cases (C) as the result of environmental contamination (new cases that occurred following shedding in the pen at previous sampling point (t − 1)) for all time intervals in all pens, we estimated γ to be equal 0.1422 day−1 as described [22]. A sensitivity analysis in which concealment rate was increased to 0.5 and decreased to 0.001 was performed to ensure the best estimate was made for the model and was in agreement with current literature [46]. Time intervals (t) in the current study were either 2 or 3 days long; therefore, the exposure (survival) rate (σ) was calculated based on the following equation: with ∆t = 2 or 3 days. ET for period t was then calculated as follows:where is the environmental contamination caused by direct shedding of infectious material of IN calves in the time interval prior to t; is the environmental contamination caused by direct shedding of infectious (I) CE calves in the time interval period prior to time t; and is the amount of environmental accumulation resulting from all previous shedding events leading up to a time interval i.e. up to t-1 (accounting for the concealment rate over time). Total environmental contamination (ET) for a specific period (t) was calculated based on the number of calves shedding in the time interval immediately prior to t (t−1) and any environmental contamination already in the prior interval was added to account for the accumulation effect.

Transmission rate parameters

The transmission rate parameter βT is dependent on the population composition in a pen and represents all transmission that could occur between the pen mates in a group housing pen. Susceptibility and infectivity among calves can be variable, with part of that variation attributed to observable differences, e.g. inoculation status. These differences in shedding patterns between penmates influence the transmission rate parameter β, as not all calves can be assumed to have the same susceptibility and infectivity. All calves were either inoculated (IN) or contact exposed (CE) and this dichotomy can influence susceptibility to begin and to continue shedding, as well as impact how infectious a calf is i.e. propensity to cause shedding in other calves by the shedding calf. To differentiate susceptibility status of IN and CE calves, we used an explanatory dummy variable (INO) classified as 0 or 1 for the recipient animal to be either inoculated (INO = 1) or contact infected (INO = 0). Additionally, the fraction of infectivity for either IN or CE can be calculated based on direct fecal shedding and contamination in the environment. For example, the expression for the fraction of infectivity caused by infectious IN calves in one time interval is:where is the total infectivity caused by IN calves (MAPIN) and is total infectivity caused by all calves (MAPT).

Additionally, we can test whether transmission by direct shedding in the current period, or in subsequent periods due to environmental contamination, are the same. These comparisons are based on calculating the fraction of the total environmental infectivity where and is equal to total infectivity in the environment (ET).

Yet another explanatory variable we added was the day since the start of group housing (day). This variable has the same value for all calves in the pen for each observed interval and thus the corresponding regression coefficient contained both susceptibility and infectivity effects. When main effects were significant, potential interactions were also examined.

Thus, based on the above, the following transmission route parameters were quantified using GLM (Figure 2), due to transmission caused by: direct contact between IN calves (), direct contact from IN calves to CE calves (), environmental contamination between IN calves (), environmental contamination from INto CE calves (), direct contact between CE calves (), direct contact from CE calves to IN calves (), environmental contamination between CE calves () and environmental contamination from CE calves to IN calves (), where and all are based on regression coefficients in the following equation:The regression coefficients can be interpreted as follows: C0 is the intercept and is the log of the transmission rate parameter of the reference that is direct transmission from contact-exposed calves (that have begun fecal shedding) to other contact-exposed calves at day 0 of the experiment; C1 is the extra effect of the recipient being of the inoculated type (susceptibility to begin shedding); C2 is the extra effect on transmission with days since the start of infection (susceptibility and infectivity); C3 is the extra effect on infectivity if the infectious calf is of the inoculated type; C4 is the extra infectivity of a unit of infectious material arriving at the recipient through the environment; and C5 is the effect of the interaction of susceptibility of an inoculated recipient and the day since the start of the experiment.

Figure 2

The GLM model for 1 time interval. The transmission rate parameters depend on the number of infectious IN and CE calves (IIN and ICE, respectively) and/or on the amount of MAP in the environment caused by IN or CE shedding (EIN and ECE, respectively). E depends on MAP secretion in current time interval from the infectious calves (either IN or CE) and amount of MAP remaining in the environment from excretion by the infected calves on previous days (t-1) weighted by σ.

From this, transmission rate parameters can be derived as follows (assuming that all regression coefficients were significantly different from zero):

is the transmission rate parameter from an infectious IN calf to a susceptible In class="Chemical">N calf through direct contact;

is the transmission rate parameter from an infectious IN calf to a susceptible In class="Chemical">N calf through the environment;

is the transmission rate parameter from an infectious CE calf to a susceptible CE n class="Species">calf through direct contact;

is the transmission rate parameter from an infectious CE calf to a susceptible CE n class="Species">calf through the environment;

is the transmission rate from an infectious IN calf to a susceptible CE n class="Species">calf through direct contact;

is the transmission rate parameter from an infectious INto a susceptible CE n class="Species">calf through the environment;

is the transmission rate parameter from an infectious CE to a susceptible IN n class="Species">calf through direct transmission; and

is the transmission rate parameter from an infectious CE to a susceptible CE calf through the environment.

Regression coefficients were tested for being different from zero. For regression coefficients not significantly different from zero (p > 0.05), the corresponding explanatory variable was dropped from the model unless it led to confounding (> 25% change in the regression coefficient of the other explanatory variables). Interaction terms were only calculated when both main effects were statistically significant.

Infectious periods: TIN and TCE

Infectious periods were calculated based on duration of shedding events for IN and CE calves as described [14]. The first time a calf tested MAP-positive in fecal samples was considered day 1 of its infectious period. The infectious period then would end if the next sample was negative, or continue until a negative sample was detected. A shedding event was defined as continuous streak of positive fecal samples and ranged from 2 days (1 positive fecal sample) up to 54 days (23 consecutive positive fecal samples) during the 90+ days of group housing.

Mean duration of infectious period for IN n class="Species">calves (TIN) and CE calves (TCE) was calculated as the average length of shedding events for IN calves or CE calves, respectively. The 95% confidence intervals (CI) were calculated using the standard variation of the mean of log TIN and log TCE for IN calves and CE calves, respectively.

Using the GLM: R0CE and R0IN

The R0 can be estimated by multiplying the rate of transmission (β) by the infectious period (T). The concealment rate by the environment due to new bedding being added and infectious material being covered was also taken into consideration using the variable , thus represents the total exposure rate of MAP in the environment for each day of the infectious period T. The equation for estimating R0 is:

The 95% CI of R0 was calculated using the variance and regression constant of the GLM results (log β) and the variance and the average of the logarithm of the infectious period T. R0 is calculated separately for CE and IN, due to differences in infectiousness and susceptibility to begin shedding. Therefore, the estimated basic reproduction value for CE calves is R0CE and R0IN for inoculated calves.

Final size (FS) model: R0T, R0IG

The FS model was based on the SI model, which uses the total number of susceptible animals remaining at the end of the experiment [23]. Two FS models were used to account for varying definitions of MAP infection. In the first FS model, infected (I) calves were counted as those with at least one culture-positive tissue sample by the end of the experiment, and remaining susceptible (S) calves as those with no positive tissue sample. R0T was the estimate of the number of new tissue-positive infections that would result from the introduction of one tissue-positive calf to the susceptible population. In the second FS model, infected (I) calves were those with at least 1 positive INF-γ sample by the end of the experiment, whereas remaining susceptible (S) calves had no positive INF-γ samples. R0IG was the estimate of the number of new INF-γ positive infections that would result from the introduction of 1 INF-γ positive calf to the susceptible population. For both definitions, IN calves started the trial with infected status (I), and CE calves were considered susceptible (S), as all IN calves met the INF-γ and tissue culture criteria for categorization as infected calves. By the end of the trial, outcomes for the number of cases (C) in a pen could be 0, 1 or 2 where the number of susceptible calves will be 2, 1, or 0 (respectively) depending on the extent of transmission, and the formula for the probability of each outcome is as follows [27]:

These probabilities for each single experiment can be combined assuming that the experiments are independent by using the multinomial distribution with the 3 p values and the multinomial total being 17.

Thus, R0 can be derived analytically by the maximum likelihood estimate (MLE) given these probabilities for the definition of tissue infection (R0T) or In class="Chemical">NF-γ infection (R0IG) [27, 47].

Results

Detection of infection

Descriptive results were presented more extensively in the previous report [14]. All IN and CE calves had MAP-culture positive fecal samples at various times during group housing; however, CE calves ceased shedding after individual housing (Table 1). All IN calves and 5/14 CE calves had positive INF-γ samples, whereas all IN calves and 7/14 CE calves had at least 1 MAP culture-positive tissue sample (Table 1). All control calves tested positive for MAP in feces on day 21 of group housing, with 1 calf having an additional positive fecal sample on day 56. Positive fecal samples among control calves were likely due to contamination during laboratory processing and not fecal shedding, as all control animals tested positive on the same day. All control calves had negative tissue results and negative INF-γ results.

Table 1

MAP culture results detected by F57-specific qPCR for fecal and tissue samples and INF-γ detection from whole blood samples

Pen	Calf ID	IN:CE	MAP detection with F57 qPCR following culture of individual fecal samples
			Days following group housing
			0	3	5	7	10	12	14	17	19	21	24	26	28	31	33
1	1	IN	I		C	I		C		C			C	I	I		C
1	2	IN				C	I	I		C	I		C	I	I	I	I
1	15	CE								C	I
1	16	CE									C
2	3	IN		C	I	I			C	I	I	I	I	I	I		C
2	4	IN	I		C			C	I	I	I	I	I	I	I		C
2	17	CE				C		C	I							C	I
2	18	CE						C						C			C
3	5	IN					C		C	I		C	I	I	I	I	I
3	6	IN							C	I	I	I	I	I	I	I	I
3	19	CE						C	I
3	20	CE							C						C		C
4	7	IN	I			C	I				C	I	I		C
4	8	IN	I	I		C	I	I	I	I	I	I	I	I	I	I	I
4	21	CE												C		C	I
4	22	CE			C				C					C		C
5	9	IN				C	I			C	I	I			C	I	I
5	10	IN				C		C	I	I	I	I	I	I	I	I	I
5	23	CE														C
5	24	CE				C							C			C
6	11	IN						C		C		C				C
6	12	IN	I						C	C		C	I				C
6	25	CE									C					C
6	26	CE													C
7	13	IN		C				C			C	I	I		C		C
7	14	IN			C	I			C			C	I	I	I		C
7	27	CE					C			C			C			C
7	28	CE									C		C

“I” indicate an infectious animal and “C” indicates a new shedding event. CE calves were individually housed 87 days following group housing (italic).

CE contact exposed, IN inoculated.

MAP culture results detected by F57-specific qPCR for fecal and tissue samples and In class="Chemical">NF-γ detection from whole blood samples

GLM

There was no difference in infectiousness (ability to cause shedding in contact exposed calves) between CE and IN calves (C3 = 0), environmental or direct oral transmission (C4 = 0) and thus also not for the interaction of the 2 variables. Susceptibility to begin shedding was different for IN and CE calves (C1 ≠ 0), likely due to the difference in exposure (IN vs CE); therefore, the contrast was included as an explanatory variable in the final GLM model (p < 0.001). Explanatory variables also included in the final model were number of days following group housing date (day), the interaction between housing date and inoculation status and the intercept. All estimates of transmission rate parameters were calculated for the beginning of the experiment at “day” = 0; therefore, the interaction between “day” and susceptibility status was not included in the final model (Table 2). When the interaction was included in the model, there was no significant difference in how the model fit the data and the difference in the coefficients are shown (Table 2). No coefficients of the explanatory variables changed more than 11%; therefore, it was concluded that no confounding was present.

Table 2

Estimates for coefficients and their 95% confidence for fecal shedding of subsp. in final General Linearized Model (GLM) with and without the interaction term for start day*susceptibility status

Model	Model coefficient (code)	Estimate (day−1)	p value	95% confidence interval	AIC value for model fit
Final model with no interaction between susceptibility and day					530.67
	Susceptibility status (C1)	1.41	< 0.001	1.05, 1.77
	Time since group housing (C2)	− 0.027	< 0.001	− 0.035, − 0.019
	Infectiousness of contact exposed calves (C0)	− 1.84	< 0.001	− 2.21, − 1.47
Model including interaction between susceptibility and day					527.86
	Susceptibility status (C1)	2.08	< 0.001	1.37,2.79
	Time since group housing (C2)	− 0.014	0.033	− 0.028, − 0.001
	Interaction (C5)	− 0.019	0.028	− 0.036, − 0.002
	Infectiousness of contact exposed calves (C0)	− 2.26	< 0.001	− 2.82, − 1.71

As IN and CE n class="Species">calves had different infectious intervals, R0 had to be estimated separately for these 2 groups of calves. The transmission rate parameter βCE_CE was estimated as 0.158 per day (0.109, 0.230). The average infectious period for CE calves was 2.91 (0.98, 6.08) days. Therefore, the equation to calculate R0CE is as follows:

Based on the GLM, the estimated basic reproduction value for CE calves (R0CE) was 3.24 (1.41, 7.41). Thus, the transmission of fecal shedding events among calves by themselves without the effect of inoculation is significantly above one, i.e. MAP is transmitted among calves.

For IN n class="Species">calves, the transmission rate parameter βIN_IN was 0.649 per day (0.437, 0.965) and the average infectious period was 7.49 (0, 25.3) days; therefore, the equation to calculate the reproduction value for IN calves was as follows:

The estimated reproduction value for IN n class="Species">calves R0IN was 24.6 (4.57, 133.3). Note that as inoculated calves were only inoculated at the beginning of the trial, the estimate is not in itself a reproduction ratio of a population of calves, but rather reflects the initial susceptibility to begin shedding of calves exposed in a different way (inoculation) then by calf-to-calf transmission.

FS model

When quantifying transmission based on the INF-γ definition of MAP infection, in 3 pens no susceptible (S) calves became infected (I), in 2 pens only 1 of the 2 susceptible (S) calves became infected (I) and in 1 pen, all susceptible (S) calves became infected (I). Therefore, the estimated reproduction value for defining infected calves as infected based on INF-γ (R0IG) was 0.90 (0.24, 2.59).

When quantifying transmission based on being tissue culture-positive as a definition of MAP infection, in 2 pens none of the susceptible (S) calves became infected (I), in 3 pens, 1 susceptible calf (S) became infected and in 2 pens all susceptible (S) calves became infected (I). Therefore, the estimated reproduction value for defining infected calves based on tissue (R0T) was 1.36 (0.45, 3.94).

Discussion

Quantification of transmission of MAP infection based on fecal shedding and accumulation of environmental contamination indicated that one CE calf that has a new shedding event within a completely susceptible population of calves in a clean environment could cause approximately 3 calves to begin shedding (R0CE = 3.24; 1.41, 7.41). This estimate represents MAP transmission for calves in a herd entering a clean environment as it starts from a “naturally” infected calf. Using MAP-positive tissue culture, or cellular immune response to MAP, one infectious calf would cause approximately one CE calf to become infected in group housing (R0T = 1.36; and R0IG = 0.90, respectively). Thus, the R0 value for MAP transmission among naturally infected calves lies between 0.9 and 3.24 depending on the definition of infection (fecal shedding, INF-γ response, or tissue infection). Our study indicated that transmission of MAP between group-housed dairy calves occurred (R > 1) and that potential shedding outbreaks may occur. When including the interaction between “day” and susceptibility based on exposure status, the transmission estimate for CE calves decreases, indicating that over time infectiousness or susceptibility to begin shedding decreases. However, this change may be due to several extraneous variables specific to the trial and the environmental contamination of MAP (bed pack, addition of clean material, cleaning of environment, etc.) as well as the possibility that CE calves are already infected and shedding events are due to intermittent active shedding rather than the result of penmates shedding. Therefore, the model without the interaction provides a better estimate of the transmission estimate for group-housed calves entering a clean environment (corresponding to the beginning of the trial, day = 0).

An important factor to take into consideration is that duration of fecal shedding was determined by the number of days that fecal samples tested positive for MAP. It is possible that due to intermittent shedding, detection of a positive shedding event may be influenced by these patterns [14] and affect the length of shedding event (or infectious period) in the current model. However, fecal samples were collected more frequently than previously reported [12, 48] and cultured, we expected to have captured more shedding events and more accurately estimated the infectious period. Additionally, the estimate of the transmission rate parameter beta is based on the same definition of infectiousness for calves. Therefore, although the interval of shedding may be underestimated, there is no effect on the reproduction ratio.

That the R0CE value for fecal shedding among calves was higher than previously described [16, 17], this may be due to the fact that new shedding events were measured and the difference between active shedding and passive shedding cannot be differentiated. Although these types of shedding are not discernable, evidence of fecal shedding is typically assumed to be an indicator of infection and any shedding event is used to identify infected animals [16, 17]. Additionally, environmental contamination was taken into consideration in the current study. All calves shed MAP into the group pen and due to the resilience of MAP, the presence of MAP in this environment remained a source of infection for a prolonged interval. Presence of MAP in the environment, especially resulting from a high herd prevalence of MAP, has a large impact on fecal–oral transmission from cows to calves [49-51]. In the analysis of our experiments, only when environmental contamination with MAP was taken into consideration did the model better capture transmission of fecal shedding between calves, based on AIC values. Both IN and CE calves shed MAP into the environment with no difference between infectivity rates (βIN = βCE); however, there was a difference in the IN or CE calves’ susceptibilities to begin shedding and duration of shedding events. The IN calves were more likely to begin shedding and had longer shedding events, leading to a higher point estimate R0 value (24.6) than CE calves (3.2). Indirect transmission of infection (indicated by positive tissue and immune responses) that may occur through the environment is likely to have an impact on calf-to-calf transmission, especially considering the large amount of intermittent shedding that occurs within these animals and opportunities to miss detection of shedding events. It should be taken into consideration that the R0CE may be an overestimation of new cases of shedding, as intermittent shedding could not be taken into consideration and each new shedding event was considered a new case of shedding rather than a result from the original exposure. However, it should be considered that all CE calves ceased to shed following individual housing in a clean environment. With both IN and CE calves, it is possible that individuals were already infected and shedding events were due to intermittent shedding, especially among IN calves. However, it is not possible to differentiate which events were due to intermittent shedding, new infections and active shedding, or passive shedding. That all IN and CE animals shed at least twice during group-housing, FS estimation based on fecal samples would result in an infinite R0 value; therefore, the GLM model based on new shedding events was deemed most appropriate.

One of the difficulties associated with MAP is that an animal can be “infectious” and/or “infected”; therefore a GLM model was used to identify “infectious” calves, whereas the FS model was used to identify “infected” calves and transmission of infected status [52]. Although fecal shedding may not be a true indication of MAP infection and could be the result of either passive or active shedding, the evidence that 1 CE calf might cause 3 additional calves in the population to have shedding events has important implications for JD control programs. Any pass-through of infected material (pass-through shedding) has potential to cause infection within the shedding calf; therefore, new shedding events may be of importance to control programs, as these calves may develop an active infection. Although all calves had positive fecal shedding detected throughout the study, not all had positive tissue samples. Transmission of fecal shedding among calves may be one explanation for the low, consistent prevalence and maintenance of MAP infections on farm. Of all CE calves that were fecal shedding, half (7/14) had culture-positive tissue samples and the estimated R0T value for “transmission” of MAP-positive tissue samples from 1 CE calf to a susceptible calf was 1.36. However, this was not significantly larger than 1. Calves that are tissue-positive, may cause 1 other calf in the group pen to become tissue-positive through fecal shedding, potentially leading to infection and symptoms later in life; however, this is only part of all MAP transmission dynamics that occur on farm. Additionally, it should be taken into consideration that testing of tissue samples in non-clinical calves may lead to a high number of false-negative samples, due to low MAP concentrations in tissue and few samples sites that are positive, which would lead to an underestimation of the number of infected calves [13, 53–55]. The tendency for calves to shed MAP, but only some to have detectable positive tissue samples indicating infection, may be one explanation for the low, consistent prevalence and maintenance of MAP infections on farm.

Epidemiological models are a useful tool when testing a hypothesis in the field is difficult or improbable. This is the case with MAP, as infections are slow progressing with a long latent period, susceptibilities and infectivities may vary between age groups, there are low sensitivities and specificities of diagnostic tests and not all transmission routes (including the environmental impact) are clearly understood [28]. These models are simplified representations based on current knowledge designed to investigate specific hypotheses such as: MAP transmission dynamics in a herd [32, 56], impact of control programs [33, 35], predicting fadeout and persistence of disease [19, 29] and effects of infectious young stock on disease control [20]. Outcomes of each of these models rely on the availability of knowledge regarding the underlying mechanisms of the disease, as well as assumptions made to fill in knowledge gaps and the modelling objective to answer the question at hand [28]. Due to contradictory findings regarding the impact of calf-to-calf transmission [16, 19, 20, 29, 56], inconsistent values are used to model the impact of infection and transmission rates of young stock and subsequent herd transmission dynamics. Note that in the earlier experimental quantification of transmission [16] the estimate of R was below 1 and we now know that this could have been caused by not taking the environment into account in the estimation, because the time to first contact infection underestimates the transmission rate parameter as transmission increases after being placed in a clean environment. The current study provides R0 values for calf-to-calf transmission that can be used to better model transmission dynamics on farm, with potential to examine different infectious abilities. The fecal shedding R0 value demonstrates occurrence of new shedding events caused by infectious calves and the contaminated environment in group-housed pens; however, this value did not necessarily represent infections that will persist over time. Based on INF-γ or tissue results, a more conservative estimate for persistent infections in calves could be obtained. Although these different R0 values correspond to different parameters of a model and different definitions of infection, these results may enable better, more accurate modelling of infection and transmission of MAP within dairy farms, leading to more effective control and prevention programs.

Currently, control programs focus on decreasing transmission to susceptible calves from infectious dams by separating young stock soon after birth, based on the assumption calves will not transmit infections to each other [10]. However, results from this study indicated that shedding of MAP may be more easily transmitted to other pen mates than previously assumed, not only from direct contact, but potentially more importantly, due to environmental contamination caused by shedding animals. Individual calf housing may be one solution to the risk of calf-to-calf transmission; however, this method of calf housing may not be possible for all herds, due to increased labour costs [57], transition to automated feeding systems [58], decreased calf welfare [59-61] and the need for effective cleaning between successive calves. A second solution would be regular and vigorous maintenance of a clean group-housed environment, in which contaminated material is removed and replaced with new bedding and base, and sides of the enclosure are disinfected. Additionally, young stock should be considered for inclusion of MAP testing in a herd testing program, as up to 2% of young stock on positive farms have been identified positively for MAP shedding [15]. Animals with positive fecal samples, or whose environment tests are positive for MAP, should be immediately removed from group-housing and monitored/tested in the future to determine infection status. Although not all infected or infectious animals may be detected due to intermittent shedding and lack of sensitivity in immune tests, this would be an added measure that may decrease MAP prevalence on dairy farms.

Results from this research indicate a potential for group-housed calves to cause shedding among penmates, leading sometimes to more extensive infection as evidenced by tissue and immune responses; however, consideration should be made for parameters of the experimental situation when applying results to the field. Future research is required regarding calf-transmission that occurs on farms and effects of a contaminated environment, as positive young stock have been detected in positive herds [15]. Ideally, a longitudinal study on infected farms should be conducted in which fecal samples, environmental samples and blood samples are collected regularly from all penmates that are group-housed, beginning at birth, for the entire duration of group housing. These longitudinal studies will lead to a better understanding of transmission rates and infection status that may occur due to direct contact and/or environmental contamination that occurs within a commercial dairy herd.

In conclusion, in this study, transmission of MAP among group-housed calves was quantified using GLM and FS modelling, based on various definitions of infectiousness or infection. The GLM model was based on the propensity of new shedding events to occur, whereas the FS model was based on tissue or immune responses as a definition of infection. Based on the definition of infection or infectious, the R0 value for MAP was between 0.9 and 3.24. Using the GLM model of fecal shedding over time and an SIS model, the R0CE for an infectious contact-exposed animal was estimated to be 3.24. However, this model only considered fecal shedding and therefore was not a perfect indication of infection among calves. Alternatively, FS models using positive tissue and immune responses as a definition of infection had R0 values equal to 1.36 and 0.9, respectively. Together, this study provided evidence that transmission of MAP infection among group-housed calves was not only present, but may be a contributing factor to maintenance of infections within a herd and therefore, should be seriously considered in future Johne’s disease prevention and control programs.