Besnoitia besnoiti lytic cycle in vitro and differences in invasion and intracellular proliferation among isolates
© Frey et al. 2016
Received: 15 December 2015
Accepted: 23 February 2016
Published: 29 February 2016
Bovine besnoitiosis, caused by the protozoan Besnoitia besnoiti, reduces productivity and fertility of affected herds. Besnoitiosis continues to expand in Europe and no effective control tools are currently available. Experimental models are urgently needed. Herein, we describe for the first time the kinetics of standardised in vitro models for the B. besnoiti lytic cycle. This will aid to study the pathogenesis of the disease, in the screening for vaccine targets and drugs potentially useful for the treatment of besnoitiosis.
We compared invasion and proliferation of one B. tarandi (from Finland) and seven B. besnoiti isolates (Bb-Spain1, Bb-Spain2, Bb-Israel, Bb-Evora03, Bb-Ger1, Bb-France, Bb-Italy2) in MARC-145 cell culture. Host cell invasion was studied at 4, 6, 8 and 24 h post infection (hpi), and proliferation characteristics were compared at 24, 48, 72, 96, 120, and 144 hpi.
In Besnoitia spp., the key parameters that determine the sequential adhesion-invasion, proliferation and egress steps are clearly distinct from those in the related apicomplexans Toxoplasma gondii and Neospora caninum. Besnoitia spp. host cell invasion is a rather slow process, since only 50 % of parasites were found intracellular after 3–6 h of exposure to host cells, and invasion still took place after 24 h. Invasion efficacy was significantly higher for Bb-France, Bb-Evora03 and Bb-Israel. In addition, the time span for endodyogeny to take place was as long as 18–35 h. Bb-Israel and B. tarandi isolates were most prolific, as determined by the tachyzoite yield at 72 hpi. The total tachyzoite yield could not be predicted neither by invasion-related parameters (velocity and half time invasion) nor by proliferation parameters (lag phase and doubling time (dT)). The lytic cycle of Besnoitia was asynchronous as evidenced by the presence of three different plaque-forming tachyzoite categories (lysis plaques, large and small parasitophorous vacuoles).
This study provides first insights into the lytic cycle of B. besnoiti isolates and a standardised in vitro model that allows screening of drug candidates for the treatment of besnoitiosis.
Besnoitia besnoiti, a tissue cyst-forming apicomplexan protozoan parasite, is the etiological agent of bovine besnoitiosis. Bovine besnoitiosis is a chronic and debilitating disease of economic importance that affects fertility and can lead to reduced production efficacy (reviewed in ). B. besnoiti has gained the attention of Veterinary Health officials in Europe since bovine besnoitiosis has been identified in previously unaffected regions of this continent during the last years . From the initially recognised foci in Alentejo in Portugal and the Pyrenees, the disease has spread almost all over Spain and France, and has now been observed in Italy, Germany, Switzerland, and most recently in Hungary, Croatia, and Belgium [3–5]. Currently, there are no effective drugs, and no vaccine is licensed in Europe. Closely related to B. besnoiti are B. tarandi that mainly infects reindeer and other wild ungulates , B. caprae that mainly affects goats , and B. bennetti, a parasite infectious for equids .
The life-cycles of these species mentioned above have not yet been elucidated to date. Cats are definitive hosts for several other Besnoitia species of small mammals . A heteroxenous life-cycle is suspected for B. besnoiti, but the definitive host still remains elusive [10, 11]. Cattle represent an important intermediate host for the parasite , and harbour tachyzoites and the cyst-forming bradyzoites. Upon infection, tachyzoites undergo proliferation in endothelial cells and fibroblasts [12, 13]. During this acute phase of the disease, animals develop fever, nasal and ocular discharge, photophobia, edema, and lameness [12–16]. The chronic form of besnoitiosis is characterised by the presence of thick-walled tissue cysts containing bradyzoites. Predilection sites for the cysts are the mucosae of the upper respiratory tract, of the lower genital tract, the scleral conjunctiva of the eyes, and the dermis and tendons of the lower limbs [12, 15–18]. In endemic situations, most animals remain subclinically infected with only about 10 % of animals developing the characteristic disease [19, 20]. The reasons why only a small percentage of animals develop disease are not yet understood. There are parasite- and host-related factors that may influence the outcome of the infection (reviewed in ). Therefore, standardised in vitro and in vivo models for B. besnoiti are urgently needed to study the pathogenesis of the infection and to assess novel drug or vaccine candidates. Several in vitro models have been employed to isolate and propagate Besnoitia parasites, even to test a few drug candidates [22, 23]. However, the lytic cycle has not been described in detail yet, and time points of adhesion-invasion, proliferation and egress may differ from the ones reported in other apicomplexans, or even between isolates. This knowledge is essential to enable comparable studies. It is well known that in other closely related tissue cyst-forming apicomplexan parasites, such as Toxoplasma gondii and Neospora caninum, parasite invasion and proliferation are significantly influenced by parasite strain, host cell line, and activation state, and these phenotypic traits have been associated with virulence [24–26]. Intra-specific variability has been addressed in molecular and virulence studies in T. gondii and N. caninum [24, 26–33]. For T. gondii in Europe and North America, a population structure was demonstrated that is dominated by three major clonal lineages (Types I, II, and III) displaying different behaviour in vitro, contrasting virulence in the mouse model, and distinct genetic hallmarks (reviewed in ). In other parts of the world, especially South America and Asia, highly variable atypical strains of T. gondii circulate [30, 34], some associated with high virulence for humans . For N. caninum, genetic differences among various isolates could be demonstrated [31–33], and interestingly, the in vivo virulence in mice and in cattle could be linked to distinct in vitro behavior in a standardised cell culture assay [25, 36]. So far, similar approaches have not been carried out in B. besnoiti. One study analysing four B. besnoiti isolates of different countries found genetic differences in at least four microsatellites , whereas more recently Gutiérrez-Exposito et al. reported that out of 11 identical B. besnoiti isolates one differed at one microsatellite . However, it is unknown whether different B. besnoiti isolates differ in their in vitro characteristics and finally also in their virulence.
We here report on the first detailed study on the kinetics of the lytic cycle of different Besnoitia spp. isolates in MARC-145 cell culture. The characteristics of host cell invasion and proliferation of one B. tarandi and seven B. besnoiti isolates from six countries and two continents were studied and compared. Knowledge on the lytic cycle of Besnoitia spp. is essential to understand infection kinetics of this parasite. Furthermore, this study demonstrates proof-of-concept for a standardised in vitro model that now allows to compare isolates and to perform meaningful screening of drug candidates.
Besnoitia isolates and cell culture
Besnoitia spp. isolates used in in vitro assays
Year of isolation
Passage numbers (after HFF) used in experiments
Brown Swiss cow
6 to 31
6 to 18
Emilia Romagna Appenines
3 to 11
2 to 14
3 to 18
2 to 15
3 to 13
6 to 11
All isolates were maintained in monolayer cultures of a monkey kidney cell line (MARC-145 cells) employing a parasite-host cell ratio of 2–5:1, in Dulbecco’s modified Eagle medium (DMEM) supplemented with 10 % fetal bovine serum (FBS), 15 mM HEPES (pH 7.2), 2 mM glutamine, penicillin (100 U/ml), streptomycin (100 μg/ml) in 5 % CO2/37 °C. FBS was previously checked for absence of anti-Besnoitia, anti-Toxoplasma, and anti-Neospora IgG by IFAT .
Tachyzoites were harvested 3 days post-infection (p.i.), when most tachyzoites were still intracellular, and parasites were liberated from their host cells by passing them repeatedly through a G25 needle, followed by purification through disposable PD-10 columns (Sephadex G-25; GE-Healthcare, Buckinghamshire, UK) as described earlier for Neospora . Tachyzoite viability was confirmed by trypan blue exclusion and they were counted in a Neubauer chamber. Purified tachyzoites were employed in invasion and proliferation assays.
In order to minimise artefacts due to prolonged maintenance in MARC- 145 cells, all isolates were passaged 4 times through Human Foreskin Fibroblast (HFF) cells, a cell-line they had not been cultured in previously, before being used in the experiments. Afterwards, they were maintained in MARC-145 and used in low passages (Table 1) for the invasion and proliferation experiments.
Invasion assays were carried out using confluent MARC-145 cell monolayers grown in 24 well plates (105 cells/well in DMEM supplemented with 10 % FCS). For infection, the medium supernatant was discarded and monolayers were exposed to tachyzoites suspended in 1 ml medium. In a preliminary experiment, infections were carried out using 100, 200, 500, 1000 and 2000 tachyzoites per well in order to establish a reliable and optimised protocol. Finally, 1000 tachyzoites suspended in DMEM supplemented with 5 % FCS were used for infection, resulting in a parasite: host cell ratio of 1:100. All time points were assessed in triplicates.
Infected cultures were further maintained at 37 °C, 5 % CO2. Initially, invasion was measured after 1 h post infection (hpi), 2, 4, 6, 8, and 24 hpi. For the experimental series, only time points 4, 6, 8, and 24 hpi were retained, because after 1 or 2 hpi, almost no invasion events could be detected. At those time points, the cell monolayer was washed three times with DMEM 5 % FCS to remove non-invaded tachyzoites. Then, 1 ml DMEM 5 % FCS per well was added and the plates were incubated for a total of 72 h. For each isolate, three infected wells were left undisturbed without washing for the entire 72 h. Invasion was measured by counting the number of plaques per well by IFAT assuming that single tachyzoites caused these plaques (see 2.4.). Invasion assays were repeated three times in independent experiments for each isolate. Besnoitia strain Bb-Spain1 was included as an internal control of reproducibility in each series.
For the proliferation assay, 24 well plates with confluent MARC-145 cells (105 cells/well) in DMEM 10 % FCS were used. Monolayers were infected with 106 purified tachyzoites/well suspended in 1 ml DMEM 5 % FCS. After 4 h, the wells were washed three times with DMEM 5 % FCS and infected monolayers were further cultured at 37 °C, 5 % CO2. At time points 24, 48, 72, 96, 120, and 144 hpi, three wells per Besnoitia-isolate were harvested as follows: the supernatants were discarded and the cells were incubated in 100 μl PBS, 100 μl lysis buffer, 10 μl proteinase K (>600 mAU/ml) per well (Qiagen, Hilden, Germany) for 5 min. The lysates were transferred to Eppendorf tubes and stored at −20 °C until further quantitative real-time PCR analysis (see section 2.5). In preliminary experiments performed with Bb-Spain1, no proliferation was observed earlier than 24 hpi, and after 144 hpi all MARC-145 cells were lysed. Proliferation assays were carried out in triplicates, and were repeated in three independent experiments for each isolate. Bb-Spain1 was included as an internal control of reproducibility in each experimental series.
Assessment of invasion assays by immunofluorescence (IFAT)
Assessment of proliferation assays by quantitative real-time PCR
Samples in lysis buffer were thawed and incubated for 10 min at 56 °C. DNA was purified using a spin column protocol for cultured cells according to the manufacturer’s instructions (DNeasy kit, Qiagen). DNA was eluted in 50 μl elution buffer. DNA content and purity of each sample was assessed by UV spectrometry (NanoPhotometer®, Implen GmbH, Munich, Germany).
A modified BbRT2 quantitative real-time PCR (qPCR) assay specific for Besnoitia spp. DNA from ungulates (i.e., B. besnoiti, B. tarandi, B. caprae, and B. bennetti; ) was performed exactly as described earlier . Primers Bb3 (5’-CAA CAA GAG CAT CGC CTT C-3’; 20 μM) and Bb6 (5’-ATT AAC CAA TCC GTG ATA GCA G-3’; 20 μM) were used. As positive control served DNA extracted from a MARC-145 culture infected with Bb-Spain1, while for the negative control the DNA extraction process was performed on water instead of cells. All samples were run in duplicates in parallel in two PCRs: the BbRT2-PCR, and the 28S-PCR  to control for the integrity of the DNA and to detect inhibitory effects. A Besnoitia-standard curve corresponding to 10−1 to 105 Bb-Spain1 tachyzoites was included in each BbRT2-PCR, which allowed expressing the threshold cycle values (Ct-values) obtained for positive samples as tachyzoites. To account for variations in DNA content of the samples, the number of tachyzoites per μl was normalized using the DNA concentration per μl determined using a NanoPhotometer® (Implen, Munich, Germany), and finally expressed as tachyzoites per nanogram DNA.
To assess differences in time points in invasion and proliferation assays for each isolate, and between isolates, ANOVA tests according to Welch (not assuming equal variances between groups) were performed. In case of significant result comparison between time points, pairwise t-test comparison without the use of pooled standard deviation was carried out .
To investigate the invasion and proliferation of Besnoitia isolates as a function of time, non-linear mixed-effect model methods were employed .
Where Asym represents the horizontal asymptote as x tends to infinity and lrc represents the logarithm of the rate constant. The half time, i.e. the time to reach half the value of the asymptote, is then calculated as log 2/e lrc.
Differences in proportion of invasion outcomes between isolates were assessed by χ 2 statistics, followed by Fisher’s exact test with Bonferroni correction for each pairwise comparison. The relationship between proportions of different invasion outcomes across isolates was assessed by Pearson correlation coefficient.
Where Asym represents the horizontal asymptote as x tends to infinity, tmid represents the inflection point i.e. the time value at which y is equal to Asym/2, and scal represents the distance between the inflection point and the point for which y is equal to Asym ∕ (1 + e− 1), i.e. approximately 0.73 x Asym.
Where Y0 is the value of Y when x is equal to zero and k is the rate constant. The doubling time (dT), i.e. the time needed so that Y is doubled, was then calculated as ln(2)/k.
To allow comparison between logistic and asymptotic growth curves, we considered the exponential phase of growth in each curve, as in .
To account for unexplained variation between experiments, we included the assay as a random factor to the non-linear models. To account for heteroscedasticity, we included a variance structure to each model . Model fit was assessed by inspection of the residuals. Doubling time and half time were calculated by delta method . Statistical analysis was performed with the R statistic software (version 3.0.3; ) with the “nlme” package (version 3.1-113, ).
B. besnoiti isolates differ in their invasive capacity
Parameters estimated based on asymptotic growth model and influence of time point on total number of events in invasion assay
Comparison of IRs at the different time points for each isolate showed that for Bb-Spain1 and Bb-Ger1 there were no differences between time points, i.e. after 4 hpi, no statistically significant increase in the number of parasites that had invaded was observed any more. For Bb-Spain2, Bb-Israel, Bb-Evora03, and B. tarandi, an increase until 6 hpi was observed, while Bb-France and Bb-Italy2 continued the invasion process until 24 hpi (Fig. 2). Although not statistically significant, an increase in absolute numbers of invaded parasites was observed for all isolates until 24 hpi. This was also reflected by the half time, i.e. the time necessary for 50 % of invasion events to occur, that varied between 2.8 and 6.5 h among the different isolates. The high invaders Bb-France, Bb-Evora03, and Bb-Israel exhibited longer half times (4.3–6.5 h), compared to the other isolates (2.8–3.6 h) (Table 2).
B. besnoiti isolates exhibit distinctly different proliferation kinetics
Analysis of the tachyzoite yields after 3 and 6 days of culture showed that Bb-Israel and B. tarandi exhibited the highest proliferation rates (Fig. 3b). At day 6 p.i., four groups that statistically differed from each other could be identified (Fig. 3b). Highest values were found for Bb-Israel and B. tarandi (group 1) (p = 0.021), followed by Bb-Ger1 and Bb-France (group 2) and group 3 comprised of Bb-Italy2, Bb-Spain1, and Bb-Evora03. The lowest parasite load was found for Bb-Spain2 (Fig. 3b). Daily microscopical observations every 24 h over 6 days confirmed the values obtained by qPCR: Whereas in Bb-Israel from day 5 on the wells were filled with released tachyzoites, Bb-Spain2 infected cultures still exhibited intact host cell monolayers even after 6 days, alongside lysis plaques with free tachyzoites and intracellular parasites (Additional file 5: Figure S3).
B. besnoiti isolates display different invasion outcomes
This is the first study that describes the kinetics of the B. besnoiti lytic cycle, in order to obtain a standardised in vitro model for the tachyzoite stage. B. besnoiti tachyzoites have been successfully propagated in vitro similarly to the closely related T. gondii and N. caninum (reviewed in ). Furthermore, B. besnoiti tachyzoite cultures have been used to assess the effects of a limited number of potentially interesting compound classes such as thiazolides and arylimidamines [22, 23] and to determine whether these compounds act against invasion and/or proliferation of B. besnoiti tachyzoites. However, the B. besnoiti lytic cycle has not been characterised in depth, which is crucial to understand the pathogenesis of this parasitic infection. The currently available data suggest that a large part of the pathogenesis of the acute stage of besnoitiosis is due to repeated host cell invasion and lysis, which leads to the destruction of vascular endothelial cells . These events, combined with a potential direct toxic effect mediated by the parasite, and associated immunopathology, may lead to vascular lesions and increased vascular permeability (reviewed in ).
In this study, we have covered seven distinct B. besnoiti isolates originating from six different countries, which represent a good part of the isolates available from Europe and Israel.
Until now 15 attempts to grow B. besnoiti have been reported (reviewed in ). However, only a few studies have assessed the suitability of different mammalian primary cell cultures and permanent cell lines for the isolation and propagation of both B. besnoiti and B. tarandi species [8, 41, 50–52]. However, data reported in these studies did not allow to identify a preferred cell type. Apart from variations in the protocols employed, parasite growth may also vary in different cell lines depending on the parasite stage (tachyzoites versus bradyzoites) inoculated in the cell cultures . Thus, we decided to study the in vitro characteristics of tachyzoites of Besnoitia spp. isolates in MARC-145 cells. This cell line was successfully used to isolate B. besnoiti from affected cows ; this study and allowed long-term cultivation of different Besnoitia species [39, 41]; this study. Furthermore, MARC-145 cells were used to assess the in vitro characteristics of different N. caninum isolates [25, 36]. Importantly, the invasion efficiencies and growth rates of the different N. caninum isolates in MARC-145 cells were positively correlated with the virulence of these isolates in mice . This supports the notion that in vitro cultures represent adequate models to study characteristics of apicomplexan parasites. However, prolonged in vitro culture may alter the features, and most notably the virulence, of isolates, as has been shown for N. caninum and for T. gondii [53, 54]. Ideally, isolates with limited number of passages after their original isolation should be used for characterisation, and this was the case for Bb-Spain1, Bb-Spain2, Bb-Italy2 and B. tarandi in this study. Alternatively, isolates can be passaged through an animal host to overcome the effects of culture adaptation . This approach has been used for Bb-Ger1 and Bb-Evora03, both of which were passaged through interferon-gamma knockout (KO) mice that rapidly succumb to the acute infection [6, 8, 41]. Unfortunately, this was not possible for Bb-Israel and Bb-France. However, the passage number of the Bb-France isolate used in this study was low (20 passages in Vero cells) since this isolate was obtained only in 2013 . In order to overcome this bias of differential adaptation to culture in different cell types, all isolates used here were passaged through the same cell type, namely HFF, prior to the detailed characterization of the lytic cycle in MARC-145 cells.
Regarding the data analysis, the use of mixed-models in parasitological studies is recommended . Such models allow controlling for pseudo-replication where data are grouped, for instance when multiple measurements are taken on the same individual, but also controlling for inevitable variances between repetitions of the same experiment, as applies to this study.
The classical sequential steps of the lytic cycle, such as adhesion-invasion, intracellular proliferation and egress, all of which were previously described for N. caninum and T. gondii, are also present in Besnoitia spp. However, some key parameters of the lytic cycle of Besnoitia sp. are clearly distinct from those of T. gondii and N. caninum. The most striking differences were observed in the IRs of Besnoitia sp. as compared with e.g. N. caninum: a 5-times higher multiplicity of infection (MOI) in both the proliferation and the invasion assays had to be used for Besnoitia compared to similar experiments with N. caninum . The values for IRs (5–22 % of tachyzoites led to invasion outcomes) were lower compared to N. caninum that had displayed invasion rates of 20–90 %  as assessed by the same approach of counting formed vacuoles per added tachyzoites. Besnoitia spp. isolates also required more time to invade their host cells compared to N. caninum [25, 36]. At 6 h p.i., B. besnoiti isolates had completed only 50 % of all invasion events, whereas at this time point N. caninum invasion had already reached a plateau (100 %) [25, 36, 57, 58]. A reason for this could be the restricted time of survival in an extracellular environment of N. caninum that loses infectivity rapidly within a few hours . For T. gondii invasion has been shown to be even more efficient, with 75–80 % of tachyzoites invading their host cells within 30 min p.i. . It should be taken into account that isolates can exhibit a limited invasion and proliferation capacity depending on the host cell type . Indeed, the target cells of Besnoitia during natural infection are endothelial cells during the acute phase, and fibroblast and myofibroblasts during the chronic phase of the infection [15, 49]. Accordingly, it would be expected that low invasion rates were correlated with a short extracellular survival period. However, our study clearly showed that Besnoitia spp. isolates were capable to actively invade MARC-145 cells during a time-span of up to 24 h, which indicates that extracellular Besnoitia tachyzoites can survive and retain their infectivity for extended periods of time. This is in stark contrast to N. caninum and T. gondii, which invade their host cells within minutes. N. caninum tachyzoites loose the capacity to infect their host cells within 2–6 h of extracellular maintenance [25, 57]. Extended extracellular survival of Besnoitia spp. tachyzoites might be favourable for mechanical transmission by vectors and in particular tabanids. Experiments investigating survival of B. besnoiti in blood sucking insects reported survival of less than 1 h in Stomoxys calcitrans, less than 3 h in tsetse flies, and less than 24 h in tabanids (reviewed in ).
Following host cell invasion, proliferation of intracellular tachyzoites is initiated. The proliferation assays carried out in this study showed that there is a lag phase of approximately 24 h for Bb-Spain1 (data not shown), and Bb-Italy2 exhibited a lag phase of 48 h. For different N. caninum isolates, lag phases of 8–44 h [25, 36], and for the N. caninum Nc-1 isolate and T. gondii lag phases of 10–12 and 8–10 h have been reported, respectively . Moreover, B. besnoiti isolates exhibit a rather extended dT (17.9–35.2 h depending on the isolate). For N. caninum isolates, dTs of 9.8–14.1 h , and 14–15 h (Nc-1; ) were reported. For the T. gondii RH strain, an even shorter dT of 8–10 h has been observed . After the lag phase, logistic growth has been described for T. gondii and N. caninum . Over a period of 144 h, three isolates of B. besnoiti also exhibited logistic growth, whereas the other five isolates showed an exponential growth pattern. Exponential growth was also observed by Regidor-Cerillo et al.  in different N. caninum isolates over a time span of 68 h. However, in this study the observation phase has been 144 h and represented two lytic cycles for most of our isolates; when only the first 72 h are evaluated (e.g. one lytic cycle), all Besnoitia isolates also followed an exponential growth pattern.
The significant differences between the isolates in both invasion and proliferation parameters observed in this study suggest that there is a considerable degree of intra-species variability in B. besnoiti. Based on the IRs, isolates can be grouped into high invaders (Bb-France, Bb-Evora03 and Bb-Israel), medium invaders (Bb-Ger1) and low invaders (Bb-Spain1, B. tarandi, Bb-Spain2 and Bb-Italy2). IRs in Besnoitia species appeared to be influenced by invasion velocity as shown by the finding that those invading cells with high velocity also exhibited increased IRs. Interestingly, high IR isolates also showed a prolonged extracellular survival evidenced by higher invasion rates observed at 24 hpi compared to 4 and 6 hpi and, consequently, a higher half time of invasion. Thus, a long invasion period, combined with increased invasion velocity, finally led to highest IRs, whereas those isolates that had a short invasion period, e.g. Bb-Spain1 and B. tarandi, did not show a high overall IR.
Another intra-species classification can be carried out when assessing the results of the proliferation assays: according to the tachyzoite yield measured at 72 h p.i., Besnoitia isolates can be grouped into three categories: high prolific (Bb-Israel and B. tarandi), medium prolific (Bb-Ger1 and Bb-France) and low prolific isolates (Bb-Italy2, Bb-Spain1, Bb-Evora03 and Bb-Spain2). The differences observed among isolates increased after two lytic cycles (at 144 hpi). Isolates displayed two different growth patterns (logistic or exponential) although these were not associated to a prolific category.
Another interesting finding obtained after 72 h of culture of infected cells was the simultaneous presence of small vacuoles, large vacuoles and lysis plaques as evidenced by immunofluorescence. Based on the assumption that small vacuoles transform into large vacuoles and eventually into lysis plaques during the lytic cycle, it is conceivable that those isolates producing predominantly lysis plaques would have a more rapid lytic cycle than isolates displaying predominantly large vacuoles or even small vacuoles at the same time point. This indicates that, in comparison to T. gondii and N. caninum, the Besnoitia spp. lytic cycle is rather more asynchronous. This is likely due to an extended time period during which Besnoitia spp. invade their host cells, and the rather extended dTs. After 72 hpi, most isolates had predominantly formed lysis plaques. This was more pronounced for Bb-Spain1, Bb-Israel and B. tarandi, the two latter isolates showed high proliferation rates (i.e. tachyzoite yields at 144 hpi). The limited representation of large vacuoles suggests a prompt egression for Bb-Spain1 tachyzoites. In Bb-France, Bb-Ger1, and Bb-Evora03, about 50 % of invaded tachyzoites had already created a lysis plaque; in those isolates, large vacuoles were more often observed than small vacuoles, and these isolates showed lower proliferation rates. In Bb-Spain2, the picture was dominated by large vacuoles, followed by small vacuoles and only few lysis plaques which correlated with the lowest proliferation rate as determined for this isolate. For Bb-Italy2, however, the interpretation of the invasion outcome was less clear; in this isolate, small vacuoles dominated the picture, but there were more lysis plaques than large vacuoles. In part, this may be explained by the lag phase of 48 h. A possible explanation for the invasion pattern of this isolate could be the presence of faster and slower tachyzoites in the inoculum. We might speculate that during repeated in vitro passage of isolates, the tachyzoites will be exposed to a selection pressure that results in the synchronization of the lytic cycle. However, our results suggest otherwise, since we showed that the speed of the lytic cycle seems to become an intrinsic in vitro characteristic of a given isolate. However, we must keep in mind that the results presented here have been obtained in one specific cell line and that different host cells could influence the proliferation kinetics, as e.g. Schares et al. observed different proliferation rates of the same B. besnoiti isolate in different cell lines .
Interestingly, the tachyzoite yield could not be predicted neither by invasion related parameters (velocity and half time invasion) nor by proliferation parameters (lag phase and dT). In contrast, the IR in N. caninum was shown to impact on the in vitro proliferation of N. caninum isolates of bovine and canine origin .
B. tarandi did not differ from the assessed B. besnoiti isolates, neither with respect to invasion nor proliferation characteristics. Although one isolate is a very limited sample, the close genetic , immunogenic , and proteomic  relationship between B. tarandi and B. besnoiti seems to be reflected in comparable in vitro characteristics at least in the model used here.
Differences in in vitro behavior could be attributed to the biological diversity of these Besnoitia spp. isolates. However, especially the results obtained with Bb-Israel should be interpreted with caution, since this isolate has been maintained for extended periods in cell culture, which may have led to attenuation of virulence, and adaptation to efficient replication in vitro. For N. caninum and T. gondii parasites, this effect had been reported earlier [34, 54]. For N. caninum, a positive correlation between invasion capacity, tachyzoite yield and in vivo virulence in the mouse model could be shown [25, 36]. If this holds true for B. besnoiti as well, Bb-France and Bb-Evora03 would be expected to be the most pathogenic isolates, followed by Bb-Israel and Bb-Ger1. To test this hypothesis, animal experiments would be necessary. Liénard et al. succeeded in provoking clinical besnoitiosis in rabbits upon infection with B. besnoiti bradyzoites freshly isolated from a cow . For testing in vivo virulence of isolates kept in culture, however, no standardis ed animal models are yet available. Furthermore, all isolates used in this study had been obtained from clinically affected cows or bulls with noticeable skin lesions, i.e. no indication on differences in virulence could be derived from this. Thus, there is currently no possibility to address virulence in B. besnoiti isolates. Genetic heterogeneity of B. besnoiti isolates could be another explanation for the different in vitro characteristics observed. In this respect, only very few and somewhat contradictory studies have been published. While Madubata et al.  reported genetic heterogeneity between four B. besnoiti isolates when performing microsatellite analysis, Gutiérrez-Exposito et al.  using the same methods could only find one isolate (Bb-Italy2) that differed from 10 other, homogenous, B. besnoiti isolates. This leads to the conclusion that further molecular markers are needed and that those markers need to be applied to a broad range of isolates to characterize B. besnoiti population structure.
In summary, this study provided the basis for studies on the lytic cycle of Besnoitia tachyzoites. In vitro characteristics of Besnoitia isolates are markedly different from those of the related apicomplexan parasites N. caninum and T. gondii since lower invasion rates, longer extracellular survival, longer dT and longer lag phase were observed. Moreover, this study was the first comparing in vitro characteristics of different Besnoitia isolates showing different invasion efficiency and proliferation rates. Whether in vitro models using primary bovine endothelial cells or fibroblasts i.e. in cells that are natural target cells [15, 41] would reveal the same differences between the isolates observed here will have to be shown by future research. Similarly, further molecular and in vitro studies should aim at comparing all Besnoitia species affecting ungulates, i.e. B. besnoiti, B. tarandi, B. caprae, and B. bennetti, preferably each species being represented by various isolates. We cannot tell at this time to what extent the in vitro characteristics are transferable to in vivo traits as it has been shown for N. caninum [25, 36]. Unfortunately, there is no standardised animal model, preferably even a bovine model, available. Therefore, the standardised in vitro model currently is the best available option for studies on isolate-specific characteristics. Development of control tools against bovine besnoitiosis is urgently needed in order to combat more efficiently the spread of the disease. In particular, drugs are promising options to fight against the acute stage of the disease and avoid spread of the parasite by prophylactic treatment of newly purchased animals. In this context, this study represents also a proof-of-concept standardised in vitro model to test drug candidates.
This study was funded by the Spanish Ministry of Economy and Competitiveness (AGL2013-46442-R) and CYTED (Thematic Network 113RT0469 Protozoovac). CFF and NM are supported by the Swiss National Science Foundation (grants no. PBBEP3_141435 and PBBEP3_139398, respectively). We gratefully acknowledge Vanessa Navarro for her excellent technical assistance, and Andrew Hemphill for critical reading of the manuscript.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Alvarez-García G, Frey CF, Ortega-Mora LM, Schares G. A century of bovine besnoitiosis: an unknown disease re-emerging in Europe. Trends Parasitol. 2013;29:407–15.View ArticlePubMedGoogle Scholar
- EFSA European Food Safety Authority. Bovine Besnoitiosis: an emerging disease in Europe. EFSA J. 2010;8:1499.Google Scholar
- Beck R, Stokovic I, Pleadin J, Beck A. Bovine besnoitiosis in Croatia. In: Proceedings of the 2nd International Meeting on Apicomplexan Parasites in Farm Animals. Kusadasi, Turkey: ApiCOWplexa; 2013. p. 64.Google Scholar
- Hornok S, Fedák A, Baska F, Hofmann-Lehmann R, Basso W. Bovine besnoitiosis emerging in Central-Eastern Europe, Hungary. Parasit Vectors. 2014;7:20.View ArticlePubMedPubMed CentralGoogle Scholar
- Vanhoudt A, Pardon B, De Schutter P, Bosseler L, Sarre C, Vercruysse J, Deprez P. First confirmed case of bovine besnoitiosis in an imported bull in Belgium. Vlaams Diergeneeskundig Tijdschrift, 2015;84:205–11.Google Scholar
- Dubey JP, Sreekumar C, Rosenthal BM, Vianna MC, Nylund M, Nikander S, Oksanen A. Redescription of Besnoitia tarandi (Protozoa: Apicomplexa) from the reindeer (Rangifer tarandus). Int J Parasitol. 2004;34:1273–87.Google Scholar
- Njenga JM, Bwangamoi O, Kangethe EK, Mugera GM, Mutiga ER. Comparative ultrastructural studies on Besnoitia besnoiti and Besnoitia caprae. Vet Res Commun. 1995;19:295–308.View ArticlePubMedGoogle Scholar
- Dubey JP, Sreekumar C, Donovan T, Rozmanec M, Rosenthal BM, Vianna MC, Davis WP, Belden JS. Redescription of Besnoitia bennetti (Protozoa: Apicomplexa) from the donkey (Equus asinus). Int J Parasitol. 2005;35:659–72.Google Scholar
- Dubey JP, Yabsley MJ. Besnoitia neotomofelis n. sp. (Protozoa: Apicomplexa) from the southern plains woodrat (Neotoma micropus). Parasitol. 2010;137:1731–47.View ArticleGoogle Scholar
- Diesing I, Heydorn AO, Matuschka FR, Bauer C, Pipano E, de Waal DT, Potgieter FT. Besnotia besnoiti: studies on the definitive host and experimental infection in cattle. Parasitol Res. 1988;75:114–7.Google Scholar
- Basso W, Schares G, Gollnick NS, Rütten M, Deplazes P. Exploring the life cycle of Besnoitia besnoiti – experimental infection of putative definitive and intermediate host species. Vet Parasitol. 2011;178:223–34.View ArticlePubMedGoogle Scholar
- Pols J. Studies on bovine besnoitiosis with special reference to the aetiology. Onderstepoort J Vet Res. 1960;28:265–356.Google Scholar
- Basson PA, McCully RM, Bigalke RD. Observations on the pathogenesis of bovine antelope strains of Besnoitia besnoiti (Marotel, 1912) infection in cattle and rabbits. Onderstepoort J Vet Res. 1970;37:105–26.PubMedGoogle Scholar
- Jacquiet P, Liénard E, Franc M. Bovine besnoitiosis: epidemiological and clinical aspects. Vet Parasitol. 2010;174:30–6.View ArticlePubMedGoogle Scholar
- Dubey JP, van Wilpe E, Blignaut DJC, Schares G, Williams JH. Development of early tissue cysts and associated pathology of Besnoitia besnoiti in a naturally infected bull (Bos taurus) from South Africa. J Parasitol. 2013;99:459–66.View ArticlePubMedGoogle Scholar
- Gollnick NS, Scharr JC, Schares G, Langenmayer MC. Natural Besnoitia besnoiti infections in cattle: chronology of disease progression. BMC Vet Res. 2015;11:35.View ArticlePubMedPubMed CentralGoogle Scholar
- Basso W, Lesser M, Grimm F, Hilbe M, Sydler T, Trösch L, Ochs H, Braun U Deplazes P. Bovine besnoitiosis in Switzerland: imported cases and local transmission. Vet Parasitol. 2013;198:265–73.Google Scholar
- Frey CF, Gutiérrez-Expósito D, Ortega-Mora LM, Benavides J, Marcén JM, Castillo JA, Casasús I, Sanz A, García-Lunar P, Esteban-Gil A, Álvarez-García G. Chronic bovine besnoitiosis: intra-organ parasite distribution, parasite loads and parasite-associated lesions in subclinical cases. Vet Parasitol. 2013;197:95–103.Google Scholar
- Bigalke RD. New concepts on the epidemiological features of bovine besnoitiosis as determined by laboratory and field investigations. Onderstepoort J Vet Res. 1968;35:3–137.PubMedGoogle Scholar
- Fernández-Garcia A, Álvarez-García G, Risco-Castillo V, Aguado-Martínez A, Marcen JM, Rojo-Montejo S, Castillo JA, Ortega-MoraLM. Development and use of an indirect ELISA in an outbreak of bovine besnoitiosis in Spain. Vet Rec. 2010;166:818–22.Google Scholar
- Alvarez-García G, García-Lunar P, Gutiérrez-Expósito D, Shkap V, Ortega-Mora LM. Dynamics of Besnoitia besnoiti infection in cattle. Parasitol. 2014;141:1419–35.View ArticleGoogle Scholar
- Cortes HCE, Mueller N, Esposito M, Letao A, Naguleswaran A, Hemphill A. In vitro efficacy of nitro- and bromo-thiazolyl-salicylamide (thiazolides) compounds against Besnoiti besnoiti infection in Vero cells. Parasitol. 2007;134:975–85.View ArticleGoogle Scholar
- Cortes HCE, Mueller N, Boykin D, Stephens CE, Hemphill A. In vitro effects of arylimidamides against Besnoitia besnoiti infection in Vero cells. Parasitol. 2011;138:583–92.View ArticleGoogle Scholar
- Saeij JP, Boyle JP, Boothroyd JC. Differences among the three major strains of Toxoplasma gondii and their specific interactions with the infected host. Trends Parasitol. 2005;21:478–81.Google Scholar
- Regidor-Cerillo J, Gómez-Bautista M, Sodupe I, Aduriz G, Alvarez-García G, Del Pozo I, Ortega-Mora LM. In vitro invasion efficiency and intracellular proliferation rate comprise virulence-related phenotypic traits of Neospora caninum. Vet Res. 2011;42:41.Google Scholar
- Regidor-Cerrillo J, Arranz-Solís D, Benavides J, Gómez-Bautista M, Castro-Hermida JA, Mezo M, Pérez V, Ortega-Mora LM, González-Warleta M. Neospora caninum infection during early pregnancy in cattle: how the isolate influences infection dynamics, clinical outcome and peripheral and local immune responses. Vet Res. 2014;45:10.Google Scholar
- Ajzenberg D, Banuls AL, Tibayrenc M, Darde ML. Microsatellite analysis of Toxoplasma gondii shows considerable polymorphism structured into two main clonal groups. Int J Parasitol. 2002;32:27–38.View ArticlePubMedGoogle Scholar
- Ajzenberg D, Banuls AL, Su C, Dumetre A, Demar M, Carme B, Darde ML. Genetic diversity, clonality and sexuality in Toxoplasma gondii. Int J Parasitol. 2004;34:1185–96.Google Scholar
- Grigg ME, Suzuki Y. Sexual recombination and clonal evolution of virulence in Toxoplasma. Microbes Infect. 2003;5:685–90.View ArticlePubMedGoogle Scholar
- Lehmann T, Marcet PL, Graham DH, Dahl ER, Dubey JP. Globalization and the population structure of Toxoplasma gondii. Proc Natl Acad Sci U S A. 2006;103:11423–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Basso W, Schares S, Bärwald A, Herrmann DC, Conraths FJ, Pantchev N, Vrhovec MG, Schares G. Molecular comparison of Neospora caninum oocyst isolates from naturally infected dogs with cell culture-derived tachyzoites of the same isolates using nested polymerase chain reaction to amplify microsatellite markers. Vet Parasitol. 2009;160:43–50.Google Scholar
- Regidor-Cerrillo J, Pedraza-Díaz S, Gómez-Bautista M, Ortega-Mora LM. Multilocus microsatellite analysis reveals extensive genetic diversity in Neospora caninum. J Parasitol. 2006;92:517–24.View ArticlePubMedGoogle Scholar
- Regidor-Cerrillo J, Díez-Fuertes F, García-Culebras A, Moore DP, González-Warleta M, Cuevas C, Schares G, Katzer F, Pedraza-Díaz S, Mezo M, Ortega-Mora LM. Genetic diversity and geographic population structure of bovine Neospora caninum determined by microsatellite genotyping analysis. PLoS One. 2013;8:e72678.Google Scholar
- Su C, Shwab EK, Zhou P, Zhu XQ, Dubey JP. Moving towards an integrated approach to molecular detection and identification of Toxoplasma gondii. Parasitol. 2010;137:1–11.View ArticleGoogle Scholar
- Carme B, Demar M, Ajzenberg D, Dardé ML. Severe acquired toxoplasmosis caused by wild cycle of Toxoplasma gondii, French Guiana. Emerg Infect Dis. 2009;15:656–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Dellarupe A, Regidor-Cerrillo J, Jimenez-Ruiz E, Schares G, Unzaga JM, Venturini MC, Ortega-Mora LM. Comparison of host cell invasion and proliferation among Neospora caninum isolates obtained from oocysts and from clinical cases of naturally infected dogs. Exp Parasitol. 2014;145:22–8.Google Scholar
- Madubata C, Dunams-Morel DB, Elkin B, Oksanen A, Rosenthal BM. Evidence for a recent population bottleneck in an Apicomplexan parasite of caribou and reindeer, B. tarandi. Infect Genet Evol. 2012;12:1605–13.View ArticlePubMedGoogle Scholar
- Gutiérrez-Expósito D, De Luco DF, Arnal M, Regidor-Cerrillo J, Jiménez-Meléndez A, Ortega-Mora LM, Alvarez-García G. Wild ruminants as sentinels of endemic bovine besnoitiosis and first report of Besnoitia besnoiti infection in roe deer. In: Proceedings of the 3rd International Meeting on Apicomplexan Parasites in Farm Animals. Edinburgh,Scotland: ApiCOWplexa; 2015. p. 70.Google Scholar
- Fernández-García A, Risco-Castillo V, Pedraza-Díaz S, Aguado-Martínez A, Álvarez-García G, Gómez-Bautista M, Collantes-Fernández E, Ortega-Mora LM. First isolation of Besnoitia besnoiti from a chronically infected cow in Spain. J Parasitol. 2009;95:474–6.Google Scholar
- Gutiérrez-Expósito D, Ortega-Mora LM, Gajadhar AA, García-Lunar P, Dubey JP, Alvarez-García G. Serological evidence of Besnoitia spp. infection in Canadian wild ruminants and strong cross-reaction between Besnoitia besnoiti and Besnoitia tarandi. Vet Parasitol. 2012;190:19–28.View ArticlePubMedGoogle Scholar
- Schares G, Basso W, Majzoub M, Cortes HC, Rostaher A, Selmair J, Hermanns W, Conraths FJ, Gollnick NS. First in vitro isolation of Besnoitia besnoiti from chronically infected cattle in Germany. Vet Parasitol. 2009;163:315–22.Google Scholar
- Collantes-Fernández E, Zaballos A, Alvarez-García G, Ortega-Mora LM. Quantitative detection of Neospora caninum in bovine aborted fetuses and experimentally infected mice by real-time PCR. J Clin Microbiol. 2002;40:1194–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Dalgaard P. Introductory Statistics with R. 2nd ed. New-York: Springer; 2008.View ArticleGoogle Scholar
- Pinheiro J, Bates D. Mixed effects models in S and S-Plus. New-York: Springer; 2000.View ArticleGoogle Scholar
- Ritz C, Streibig JC. Nonlinear regression with R. New-York: Springer; 2008.Google Scholar
- R core team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2014. http://www.r-project.org. Accessed 16 Oct 2015.
- Pinheiro J, Bates D, DebRoy S, Sarkar D. R core team. nlme: Linear and nonlinear mixed effects models. R package version 3.1-113. 2013.Google Scholar
- Müller J, Hemphill A. In vitro culture systems for the study of apicomplexan parasites in farm animals. Int J Parasitol. 2013;43:115–24.View ArticlePubMedGoogle Scholar
- Langenmayer MC, Gollnick NS, Majzoub-Altweck M, Scharr JC, Schares G, Hermanns W. Naturally acquired bovine besnoitiosis: histological and immunohistochemical findings in acute, subacute, and chronic disease. Vet Pathol. 2015;52:476–88.View ArticlePubMedGoogle Scholar
- Neuman M. Cultivation of Besnoitia besnoiti Marotel, 1912, in cell culture. Tropenmed Parasitol. 1974;25:243–9.PubMedGoogle Scholar
- Göbel E, Widauer R, Reimann M, Munz E. Ultrastructure of the asexual multiplication of Besnoitia besnoiti (Marotel, 1912) in Vero and CRFK-cell cultures. Zbl Vet Med. 1985;B32:202–12.Google Scholar
- Shkap V, Pipano E, Greenblatt C. Cultivation of Besnoitia besnoiti and evaluation of susceptibility of laboratory animals to cultured parasites. Vet Parasitol. 1987;23:169–78.View ArticlePubMedGoogle Scholar
- Nischik N, Schade B, Dytnerska K, Długońska H, Reichmann G, Fischer HG. Attenuation of mouse-virulent Toxoplasma gondii parasites is associated with a decrease in interleukin-12-inducing tachyzoite activity and reduced expression of actin, catalase and excretory proteins. Microbes Infect. 2001;3:689–99.View ArticlePubMedGoogle Scholar
- Bartley PM, Wright S, Sales J, Chianini F, Buxton D, Innes EA. Long-term passage of tachyzoites in tissue culture can attenuate virulence of Neospora caninum in vivo. Parasitol. 2006;133:421–32.View ArticleGoogle Scholar
- Liénard E, Salem A, Jacquiet P, Grisez C, Prevot F, Blanchard B, Bouhsira E, Franc M. Development of a protocol testing the ability of Stomoxys calcitrans (Linnaeus, 1785) (Diptera: Muscidae) to transmit Besnoitia bensoiti (Henry, 1913) (Apicomplexa: Sarcocystidae). Parasitol Res. 2013;112:479–86.Google Scholar
- Paterson S, Lello J. Mixed models: getting the best use of parasitological data. Trends Parasitol. 2003;19:370–5.View ArticlePubMedGoogle Scholar
- Hemphill A, Gottstein B, Kaufmann H. Adhesion and invasion of bovine endothelial cells by Neospora caninum. Parasitol. 1996;112:183–97.View ArticleGoogle Scholar
- Naguleswaran A, Müller N, Hemphill A. Neospora caninum and Toxoplasma gondii: a novel adhesion/invasion assay reveals distinct differences in tachyzoite-host cell interactions. Exp Parasitol. 2003;104:149–58.View ArticlePubMedGoogle Scholar
- Vianna MC, Sreekumar C, Miska KB, Hill DE, Dubey JP. Isolation of Neospora caninum from naturally infected white-tailed deer (Odocoileus virginianus). Vet Parasitol. 2005;129:253–7.View ArticlePubMedGoogle Scholar
- Sundermann CA, Estridge BH. Growth of and competition between Neospora caninum and Toxoplasma gondii in vitro. Int J Parasitol. 1999;29:1725–32.View ArticlePubMedGoogle Scholar
- García-Lunar P, Ortega-Mora LM, Schares G, Gollnick NS, Jacquiet P,Grisez C, Prevot F, Frey CF, Gottstein B, Álvarez-García G. An inter-laboratory comparative study of serological tools employed in the diagnosis of Besnoitia besnoiti infection in bovines. Transbound Emerg Dis. 2012; doi:10.1111/j. 1865–1682.2012.01318.x.Google Scholar
- Liénard E, Pop L, Prevot F, Grisez C, Mallet V, Raymond-Letron I, Bouhsira É, Franc M, Jacquiet P. Experimental infections of rabbits with proliferative and latent stages of Besnoitia besnoiti. Parasitol Res. 2015;114:3815–26.Google Scholar
- Cortes HC, Reis Y, Waap H, Vidal R, Soares H, Marques I, Pereira da Fonseca I, Fazendeiro I, Ferreira ML, Caeiro V, Shkap V, Hemphill A, Leitão A. Isolation of Besnoitia besnoiti from infected cattle in Portugal. Vet Parasitol. 2006;141:226–33.Google Scholar