Skip to main content

Immune response profile of caruncular and trophoblast cell lines infected by high- (Nc-Spain7) and low-virulence (Nc-Spain1H) isolates of Neospora caninum

Abstract

Background

Bovine neosporosis, one of the main causes of reproductive failure in cattle worldwide, poses a challenge for the immune system of pregnant cows. Changes in the Th-1/Th-2 balance in the placenta during gestation have been associated with abortion. Cotyledon and caruncle cell layers form the maternal-foetal interface in the placenta and are able to recognize and induce immune responses against Neospora caninum among other pathogens. The objective of the present work was to elucidate the immunomodulation produced by high- (Nc-Spain7) and low-virulence (Nc-Spain1H) isolates of N. caninum in bovine trophoblast (F3) and caruncular cells (BCEC-1) at early and late points after infection. Variations in the mRNA expression levels of toll-like receptor-2 (TLR-2), Th1 and Th2 cytokines (IL-4, IL-10, IL-8, IL-6, IL-12p40, IL-17, IFN-γ, TGF-β1, TNF-α), and endothelial adhesion molecules (ICAM-1 and VCAM-1) were investigated by RT-qPCR, and protein variations in culture supernatants were investigated by ELISA.

Results

A similar pattern of modulation was found in both cell lines. The most upregulated cytokines in infected cells were pro-inflammatory TNF-α (P < 0.05–0.0001) and IL-8 (P < 0.05–0.001) whereas regulatory IL-6 (P < 0.05–0.001) and TGF-β1 (P < 0.05–0.001) were downregulated in both cell lines. The measurement of secreted IL-6, IL-8 and TNF-α confirmed the mRNA expression level results. Differences between isolates were found in the mRNA expression levels of TLR-2 (P < 0.05) in both cell lines and in the mRNA expression levels (P < 0.05) and protein secretion of TNF-α (P < 0.05), which were higher in the trophoblast cell line (F3) infected with the low-virulence isolate Nc-Spain1H.

Conclusions

Neospora caninum infection is shown to favor a pro-inflammatory response in placental target cells in vitro. In addition, significant immunomodulation differences were observed between high- and low-virulence isolates, which would partially explain the differences in virulence.

Background

Bovine neosporosis is one of the main transmissible causes of abortion in cattle worldwide [1,2,3]. The etiological agent of bovine neosporosis is Neospora caninum, an obligate intracellular parasite closely related to the zoonotic agent Toxoplasma gondii. Transplacental transmission is the main route of transmission in cattle [4] and the placenta can play a key role in the pathogenesis of bovine neosporosis [5, 6]. The direct damage produced by the multiplication of the parasite in placental and foetal tissues has been proposed as one of the possible causes of abortion observed during N. caninum infections. Importantly, the placenta is considered to be an immune regulatory organ since it acts as a modulator of foetal and maternal immune responses. In fact, an immune-mediated pathogenesis has also been suggested as a possible cause of abortion [7]. It has been shown that the multiplication of the parasite in the placenta alters the immunological balance at the maternal-foetal interface with an increase of local pro-inflammatory IFN-γ, IL-12p40 and TNF-α cytokines which could compromise the gestation, together with an increase in IL-4 and IL-10 levels [8, 9], which avoids the immunological rejection of the foetus but favours the multiplication and vertical transmission of the parasite [5, 10]. Trophoblast and caruncular cells are able to recognize pathogens and secrete cytokines and chemokines that recruit immune cells in the damaged area [11,12,13]. Thus, both cell types play a fundamental role in the initiation of innate immune responses at the placental level as well as in the development of an adaptative immune response for the pregnant dam and foetus.

Previous in vivo studies have shown the influence of the isolate on the dynamics and outcome of the infection in pregnant bovine models and in the cytokine profiles induced during the infection ([9, 14,15,16], Jiménez-Pelayo et al. unpublished data). To date, only one recent study has utilized an in vitro model consisting of immortalized bovine trophoblasts (F3) from the foetus and caruncular cells (BCEC-1) from the dam. The aim of the study was to elucidate the interactions between tachyzoites and the host cells that resemble the maternal-foetal interface of the bovine placentome while also taking into account the influence of the isolate. Maternal cells, where both isolates showed similar phenotypic traits, presented higher resistance to the infection than trophoblast cells, where the high- (Nc-Spain7) and the low-virulence (Nc-Spain1H) isolates showed marked differences in proliferation [17].

However, the interactions between the parasite and the placental target cells from an immunological point of view have not been investigated in vitro until now. Thus, the objective of the present study was to compare the immune response profiles of the bovine placental cells in vitro after the infection with two N. caninum isolates of different virulence. Messenger RNA expression levels of TLR-2, pro-inflammatory cytokines IL-8, IL-12p40, IL-17, IFN-γ, TNF-α, anti-inflammatory/regulatory cytokines TGF-β1, IL-4, IL-6 and IL-10 as well as ICAM-1 and VCAM-1 endothelial adhesion molecules were determined at 4 and 24 hours post-infection (hpi) in maternal caruncular (BCEC-1) and foetal trophoblast (F3) cell cultures and protein secretion was assessed in culture supernatants by ELISA.

Results

Expression profile of TLR-2

Our results showed that N. caninum infection for 4 h in BCEC-1 cells resulted in a significant upregulation of TLR-2 expression in Nc-Spain1H-infected cells compared with that of negative control cells (Kruskal–Wallis H-test followed by Dunnʼs multiple comparison test: χ2 = 16.2, df = 2, P = 0.0001) and with that of BCEC-1 cells infected with the high-virulence isolate Nc-Spain7 (Mann–Whitney U-test: U(8) = 8, Z = 2.591, P = 0.0007). In F3 cultures, statistical significance was not found at either 4 or 24 hpi between infected groups and the control group. However, comparing both isolates, lower expression of TLR-2 was found in the F3 cultures infected with Nc-Spain7 than in the F3 cultures infected with the low-virulence isolate Nc-Spain1H at 4 hpi (Mann–Whitney U-test: U(8) = 16, Z = 2.287, P = 0.0315) (Fig. 1a).

Fig. 1
figure 1

TLR-2, IL-8, TNF-α, IL-6, IL-12p40, TGF-β1, ICAM-1 and VCAM-1 transcript expression. Scatter-plot graphs of relative mRNA expression levels (as x-fold change) of TLR-2 (a), IL-8 (b), TNF-α (c), IL-6 (d), IL-12p40 (e), TGF-β1 (f), ICAM-1 (g) and VCAM-1 (h) in F3 and BCEC-1 cell cultures at 4 and 24 hpi with Nc-Spain7 and Nc-Spain1H isolates. Data are represented as individual points. Horizontal lines represent median values for each group. ****P < 0.0001, ***P < 0.001, **P < 0.01, *P < 0.05. Unbracketed symbols represent differences with respect to the control group, while significant differences between isolates are denoted by horizontal square brackets

Pro-inflammatory and regulatory cytokine modulation

The pro-inflammatory cytokines IL-8 (Kruskal–Wallis H-test: χ2 = 19.52, df = 2, P < 0.0001 in BCEC-1 and χ2 = 17.56, df = 2, P = 0.0002 in F3) and TNF-α (Kruskal–Wallis H-test: χ2 = 19.73, df = 2, P < 0.0001 in BCEC-1 and χ2 = 19.4, df = 2, P < 0.0001 in F3) were upregulated in both cell types at 4 hpi compared to the respective control groups (Fig. 1b, c). At 24 hpi, IL-8 expression was still increased in BCEC-1 cells infected by both isolates (Kruskal–Wallis H-test followed by Dunnʼs multiple comparison test: χ2 = 16.63, df = 2, P = 0.0003 and χ2 = 10.19, df = 2, P = 0.0117 for Nc-Spain7 and Nc-Spain1H, respectively); however, the increment of IL-8 had disappeared at 24 hpi in F3-infected cells with respect to the control group. When both isolates were compared, Nc-Spain1H induced a higher expression of TNF-α than the high-virulence isolate Nc-Spain7 at 4 hpi in F3 cells (Mann–Whitney U-test: U(8) = 0, Z = 2.579, P < 0.0001). Protein levels of the pro-inflammatory cytokines IL-8 and TNF-α were also investigated in the supernatant of control and infected cultures at different time-points post-infection. A higher secretion of IL-8 was found for both isolates in BCEC-1 cells at 24 hpi (Kruskal–Wallis H-test: χ2 = 15.87, df = 2, P = 0.0004) and in F3 cells at 56 hpi (Kruskal–Wallis H-test: χ2 = 13.74, df = 2, P = 0.001) with respect to the control group (Fig. 2a, b). Secretion of TNF-α was higher in BCEC-1 cells infected by both isolates (Kruskal–Wallis H-test: χ2 = 18.9, df = 2, P < 0.0001; Fig. 2c) and in F3 cells infected by Nc-Spain1H (Kruskal–Wallis H-test followed by Dunnʼs multiple comparison test: χ2 = 16, df = 2, P < 0.0001) at 8 hpi, although an earlier secretion of TNF-α was also found in F3 cells infected by Nc-Spain1H (χ2 = 14, df = 2, P = 0.0003) at 4 hpi (Fig. 2d, e). As observed with the TNF-α mRNA expression, Nc-Spain1H induced a higher secretion of TNF-α than did the high-virulence isolate Nc-Spain7 at 4 hpi (Mann–Whitney U-test: U(8) = 0, Z = 2.736, P = 0.0002) and at 8 hpi (U(8) = 0, Z = 2.305, P = 0.0002) in the F3 cultures (Fig. 2d, e).

Fig. 2
figure 2

IL-8, TNF-α and IL-6 secretion levels in culture supernatants. Scatter-plot graphs representing the concentration of IL-8 (pg/ml) in BCEC-1 (a) and F3 (b) supernatants infected with Nc-Spain7 and Nc-Spain1H at 24 and 56 hpi, respectively, the concentration of TNF-α (pg/ml) in the BCEC-1 supernatants at 8 hpi (c) and in the F3 supernatants at 4 hpi (d) and 8 hpi (e), and the concentration of IL-6 (pg/ml) in the BCEC-1 supernatants at 4 hpi (f). Data are represented as individual points. Horizontal lines represent median values for each group. ****P < 0.0001, ***P < 0.001, **P < 0.01, *P < 0.05. Unbracketed symbols represent differences with respect to the control group, while significant differences between isolates are denoted by horizontal square brackets

The expression levels of other important cytokines associated with N. caninum infection, such as IL-12p40 and IL-6 (Fig. 1d, e), were modified in placental cells after parasite infection. Specifically, IL-6 levels were downregulated in BCEC-1 infected by Nc-Spain1H and Nc-Spain7 at 4 hpi (Kruskal–Wallis H-test: χ2 = 16.08, df = 2, P = 0.0003) and F3 cultures infected by Nc-Spain1H at 24 hpi (χ2 = 10.5, df = 2, P = 0.0052). IL-12p40 was also downregulated but only in infected F3 cultures at 4 hpi (χ2 = 12.99, df = 2, P = 0.0015). Differences between isolates in the modulation of IL-6 and IL-12p40 were not found. In addition, we observed that the caruncular cell layer did not express IL-12p40 mRNA at any time point. The decrease in the expression of IL-6 observed in infected BCEC-1 cells was confirmed by the decrease in the secretion levels of that protein found in the supernatants from BCEC-1 cultures infected with both isolates at 4 hpi, although statistically significant differences were not found (Kruskal–Wallis H-test: χ2 = 2.765, df = 2, P = 0.251), probably because of the high deviation between samples (Fig. 2f).

Finally, pro-inflammatory IL-17 and IFN-γ responses were not detected in any cell lines at 4 nor at 24 hpi.

We also studied the mRNA levels of the anti-inflammatory cytokines TGF-β1, IL-4 and IL-10. Remarkably, we observed a decrease in the expression levels of TGF-β1 in both cell lines infected with both isolates. Specifically, a decrease was observed at 4 hpi in F3 cultures (Kruskal–Wallis H-test: χ2 = 18.44, df = 2, P < 0.0001) and at 24 hpi in BCEC-1 cultures (χ2 = 12.02, df = 2, P = 0.0025) (Fig. 1f). No differences between isolates were observed in the mRNA expression levels of TGF-β1. There was not detection in bovine placental cells of the anti-inflammatory cytokine IL-4 or the regulatory cytokine IL-10 at 4 and 24 hpi.

Endothelial adhesion molecule (ICAM-1 and VCAM-1) expression

The adhesion molecule ICAM-1 was expressed by both cell lines at 4 and 24 hpi. However, only a slight decrease in the mRNA expression levels of ICAM-1 was observed in the BCEC-1 cultures infected with Nc-Spain7 at 24 hpi compared to the control group although statistical significance was not found (Kruskal–Wallis H-test: χ2 = 5.894, df = 2, P = 0.0525) (Fig. 1g). VCAM-1 expression was detected only in the F3 cultures at 24 hpi, but differences between the infected and the control groups were not found in this culture at this time point (Fig. 1h).

Results of mRNA expression levels and protein secretion from statistical tests are reported in Additional file 1: Tables S1 and S2, respectively.

Discussion

Transmission of N. caninum across the placenta makes this organ key in the pathogenesis of bovine neosporosis. Innate immune signalling is crucial at the maternal-foetal interface, where vertical transmission of pathogens to the foetus can have profound pathological outcomes. Trophoblasts and other cell types within the placenta may also be involved in the physiological protection of the placenta [18]. Trophoblast cells have been shown to respond to some infections by producing pro-inflammatory cytokines and chemokines and endometrial or decidual cells can produce and secrete a variety of cytokines, participating in the attraction and activation of immune effector cells [19, 20]. However, these innate immune mechanisms are unexplored at the maternal-foetal interface during N. caninum infection in pregnant cattle [21].

The expression of TLRs has been described in trophoblasts and other cell types within the placenta [22]. Specifically, TLR-2 was overexpressed in bovine trophoblast cell cultures at 8 hpi [23] and TLR-3, 7 and 8 have been implicated in N. caninum recognition in the bovine placenta [21, 24]. In our study, differential activation of TLR-2 in the F3 and BCEC-1 cultures was observed. An upregulation of TLR-2 was found in BCEC-1-infected cultures, especially in those infected with the low-virulence isolate Nc-Spain1H. The caruncular part of the placentome showed a higher expression of several TLRs, suggesting that the initial recognition of N. caninum at the placental level would occur in the maternal side of placenta [21]. Taking into account the data shown in Jiménez-Pelayo et al. [17], which confirm the higher proliferation of N. caninum in trophoblast cells, an important role of placental TLR-2 in the immune response against N. caninum seems plausible. TLR activation is crucial for initiating the innate immune responses responsible for the elimination of intracellular parasites such as N. caninum, and the signalling pathway activated by TLR-2 leads to an increase in the transcription factors NF-κβ and AP-1, which trigger the synthesis of pro-inflammatory cytokines (TNF-α, IL-6, IL-12 and IL-1β) and chemokines (IL-8, RANTES) [25].

Despite differential modulation of host TLR-2, both cell types presented a similar variation in the IL-6, TNF-α and IL-8 expression levels in infected cultures. The pro-inflammatory IL-8 and TNF-α cytokines were upregulated, and secretion of the proteins in the supernatants of both cell lines was also detected by ELISA. IL-8, a cytokine with neutrophil chemotactic/activating and T-cell chemotactic activity both in vivo and in vitro, is important in the recruitment of leukocytes to the endometrium and may be a potential mediator of placental macrophage infiltration [26], which might help to eliminate the parasite. IL-8 upregulation has already been observed in bovine umbilical vein endothelial cells (BUVECs) infected by T. gondii and N. caninum [27] as well as in bovine trophoblastic cells and placentomes from cows infected with Brucella abortus [28]. TNF-α is an inflammatory cytokine whose expression has also been described for epithelial cells [29]. TNF-α is expressed in all cell types of the trophoblastic lineage and provokes a variety of biological effects on placental and endometrial cell types [29]. In addition, TNF-α, through the NF-κβ signalling pathway, coordinates the inflammatory response via the induction of other cytokines (IL-1 and IL-6) and chemokines (IL-8) and via the upregulation of adhesion molecules (ICAM-1 and VCAM-1) [30, 31], playing a role favouring protective immunity in infectious diseases [32]. There are several lines of experimental evidence indicating that TNF-α plays a role not only in immunity to N. caninum but also in the immunopathology of neosporosis. TNF-α expression and secretion may reduce the parasite presence in the placenta by inhibiting the intracellular multiplication of the parasite [33] and participating in parasite proliferation control mechanisms [34]; however, TNF-α expression is detrimental to pregnancy maintenance [8].

IL-6 expression levels were diminished in infected cultures of BCEC-1 and F3 at 4 and 24 hpi, respectively. The classification of IL-6 as Th1 or Th2 has been considered controversial since it can have characteristics of both depending on the dose, the cellular source and the gestational stage studied [35]. Currently, the presence of IL-6 displaces the Th1/Th2 balance towards a Th2 response [36]. However, in vivo models of N. caninum infection have shown IL-6 upregulation [37, 38, Jiménez-Pelayo et al. unpublished data]. This response pattern may be related to a protective action that protects the foetus and allows gestation even if the animals are born infected [10]. The decrease in IL-6 observed could be explained by the following: (i) IL-6 expression levels were affected by the high antigenic dose administered (MOI 10), resulting in downregulation [39]; (ii) the time points were not adequate for detecting the peak of IL-6 expression, and the observed decrease may be the consequence of the rapid reduction in IL-6 expression after a peak of expression [40,41,42]; or (iii) other cell types are implicated in the upregulation of IL-6 that was observed in vivo.

The anti-inflammatory cytokine TGF-β1 was also found to be downregulated in F3 and BCEC-1 cultures at 4 and 24 hpi, respectively. Several members of the TGF-β superfamily have been suggested to regulate trophoblast cell functions, and their dysregulation has been implicated in pregnancy-associated diseases. TGF-β1 is crucial in neutralizing the inflammatory responses induced by Th1-type cytokines [5]. This effect has already been observed in previous works where the reduction of TGF-β1 was shown to be beneficial for controlling N. caninum growth but detrimental for the adequate maintenance of pregnancy [38, 43].

The reduction of pro-inflammatory IL-12p40 observed in trophoblast cells, together with the lack of expression of IL-12p40 in BCEC-1 cultures, disagrees with the results of previous experimental infections [8, 9, 38, 44, 45]. Similarly, the lack of expression of pro-inflammatory IFN-γ, essential for controlling parasite infection [1, 46], and anti-inflammatory IL-4 and IL-10, related to placental protection during N. caninum infections [8, 9], lead us to hypothesize that the upregulation of IL-12p40, IFN-γ, IL-4 and IL-10 observed in vivo could be attributed to immune cells present in the placenta, such as dendritic cells, NK cells or macrophages. Therefore, trophoblast and/or caruncular cells would not be responsible for the direct production of these cytokines, although the assayed time points may not have been appropriate for their detection.

Finally, ICAM-1 and VCAM-1 expression were not modulated by the parasite infection. These adhesion molecules participate in the recruitment of inflammatory immune cells [47] and promote the adherence of monocytes to endothelial cells [48]. The upregulation of ICAM-1 has been observed in in vitro infections with apicomplexan parasites [27, 49, 50]. The absence of modulation observed in this work may be explained by differences in the timing of the expression of ICAM-1 and VCAM-1 [27, 49, 50] or by the lack of stimuli such as the acute-phase protein C-reactive protein (CRP), which is produced by the liver in response to IL-6 [51].

As mentioned above, the parasite isolate is a key factor in the outcome of the infection. In general, differences in the modulation between high- and low-virulence isolates were not remarkable in trophoblast or caruncular cells, with the exception of the mRNA expression levels of TLR-2 and TNF-α. TLR-2 levels were more upregulated by Nc-Spain1H infection than by Nc-Spain7 infection in both cell lines, which led us to hypothesize that the high-virulence isolate would activate less of the TLR recognition system, reducing the immune responses triggered by TLR-2. The inhibition of the TLRs implicated in the recognition of Trypanosoma cruzi and T. gondii in HPCVE increased the parasite burden and, importantly, TLR-2 inhibition prevented the secretion of IL-6 and IL-10, increasing parasite damage [52, 53]. The low-virulence isolate Nc-Spain1H activates the expression of TLR-2, starting an inflammatory response, which may be the cause of the lower proliferation of this isolate [17, 54] in addition to being one of the causes that explains the higher levels of TNF-α in Nc-Spain1H-infected cells, especially in trophoblast cells. Our results suggest that differential activation of the TLRs by the isolates of differing virulences should be subject to future research since they may be responsible for the biological differences observed both in vitro and in vivo.

The low-virulence isolate Nc-Spain1H also induced higher expression of TNF-α in F3. A higher TNF-α response may more efficiently control the proliferation of Nc-Spain1H in F3 cultures, which could explain the observations made by Jiménez-Pelayo et al. [17] where lower proliferation of Nc-Spain1H was observed in these cells. The lower expression of TNF-α observed during the early stage of infection of trophoblasts with the high-virulence isolate Nc-Spain7 supports the hypothesis that this isolate may modify by yet unknown mechanisms the pro-inflammatory response by trophoblast cells. However, how Nc-Spain7 is able to evade the immune response and maintain lower levels of TNF-α expression in F3 remains unknown. On the other hand, these results suggest that pro-inflammatory cytokines such as TNF-α could have a minor impact in placental damage than postulated in previous works [8, 55], but other mechanisms should be implicated in placental damage in vivo and the occurrence of abortion, such as the high multiplication ability showed by the high-virulence isolate Nc-Spain7 [17].

Conclusions

The results presented in this manuscript suggest that placental cells participate in the innate immune response at the maternal-foetal interface via a rapid pro-inflammatory response characterized by the overexpression of IL-8 and TNF-α and the downregulation of TGF-β1 and IL-6. Slight differences were detected when the immunomodulatory response induced by the high and low virulent N. caninum isolates was compared. The higher expression of TLR-2 in the F3 and BCEC-1 cells and the TNF-α in F3 cells infected with the low-virulence isolate Nc-Spain1H may indicate a higher stimulation of the immune response by this isolate or a higher immunomodulation of Nc-Spain7, which could explain the biological differences observed in vitro and in vivo. F3 and BCEC-1 cultures seem to be a good tool for the study of the TLR activation mechanisms by N. caninum. Finally, we observed that cytokines such as IFN-γ, IL-4 or IL-10, which are commonly upregulated in the placenta after N. caninum infection, are not expressed in F3 and BCEC-1 cells; we conclude that the trophoblast and caruncular epithelial cells are not implicated in the production of these cytokines in the placenta or that other pathways/cells/molecules are needed for their production.

Methods

Parasites and cell cultures

A full description of the Nc-Spain1H and Nc-Spain7 parasites and cell cultures of bovine caruncular epithelial (BCEC-1) and bovine trophoblast cells (F3) is provided in a previous report [17]. Briefly, Nc-Spain7 and Nc-Spain1H isolates were obtained from healthy, congenitally infected calves [56, 57] and tachyzoites were maintained in a MARC-145 culture as described previously [54]. The number of culture passages of both N. caninum isolates was limited (passages from 9 to 11) to maintain their biological in vivo behavior [58].

The BCEC-1 and F3 cell lines were kindly donated by Dr C. Pfarrer from the University of Veterinary Medicine Hannover and maintained following the protocols described in the literature [59, 60].

Infection of the cultures, collection and preservation of the samples

BCEC-1 and F3 cells were seeded in 25 cm2 culture flasks adjusting the number of cells in order to obtain a confluent monolayer after 24 h of culture. F3 was seeded at 106 cells per flask, whereas BCEC-1 was seeded at a concentration of 2 × 106 cells per flask. Tachyzoites were recovered from MARC-145 cultures when most of the parasites were still inside parasitophorous vacuoles; tachyzoites were purified using disposable PD-10 Desalting Columns (G.E. Healthcare, Amersham, UK) as previously described [54]. The parasite viability was checked by trypan blue exclusion, and the tachyzoites were counted. Multiplicity of infection (MOI) of 8 (8 × 106 tachyzoites in F3 and 16 × 106 tachyzoites in BCEC-1) and 10 (107 tachyzoites in F3 and 2 × 107 tachyzoites in BCEC-1) from the Nc-Spain7 and Nc-Spain1H isolates, respectively, were inoculated into confluent monolayers of F3 and BCEC-1 quickly after collection. Due to the differences observed in the infection rate between isolates [17], different MOIs of each isolate were selected with the aim of obtaining cultures infected with the same quantity of each parasite at 4 and 24 hpi. This way possible differences in the modulation of the mRNA expression levels between isolates could be attributed to differences in their biological behavior and not to the differences in the parasite burden. In addition, cultures were infected with high doses of both parasites to get a high infection of the cultures at 4 and 24 hpi so that the RNA from uninfected cells did not mask possible differences in RNA expression levels induced by the infection. The flasks were incubated at 37 °C until collection of the samples. The supernatants were collected at different time points (4, 8, 24 and 56 hpi) and stored at −80 °C for the detection of proteins by ELISA. The cultures were harvested at 4 or 24 hpi by scraping, centrifugation at 1350×g for 15 min at 4 °C and resuspending the pellet in 300 µl of RNAlater® (Qiagen, Hilden, Germany). The samples were stored at −80 °C prior to RNA extraction.

Two independent experiments were carried out and four replicates were obtained in each experiment.

RNA extraction, reverse transcription and quantitative real-time PCR

The mRNA expression levels of TLR-2, pro-inflammatory cytokines IL-6, IL-8, IL-12p40, IL-17, IFN-γ, TNF-α, anti-inflammatory/regulatory cytokines TGF-β1, IL-4 and IL-10 as well as ICAM-1 and VCAM-1 endothelial adhesion molecules were determined by real-time RT-PCR in the F3 and BCEC-1 cell layers infected with the high-virulence isolate (Nc-Spain7) and the low-virulence isolate (Nc-Spain1H) of N. caninum at an early (4 hpi) and a late (24 hpi) time point.

RNA was extracted using a commercial Maxwell® 16 LEV simplyRNA Purification kit (Promega, Madison, WI, USA) following the manufacturer’s recommendations. RNA integrity was checked by 1% agarose gel and RNA concentrations were determined using a NanoPhotometer® spectrophotometer (Implen, Munich, Germany). cDNA was obtained by reverse transcription of 2.5 µg of RNA using the master mix SuperScript® VILO™ cDNA Synthesis kit (Invitrogen, Paisley, UK), which was diluted 1:20 in molecular grade water for the qPCR assays.

The PCRs were performed using 12.5 µl of Power SYBR® Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA), 10 pmol of each primer (except for TLR-2 primers which were used at a concentration of 22.5 pmol) and 5 μl of diluted cDNA samples in an ABI 7300 Real Time PCR System (Applied Biosystems). The primers used for the qPCR reactions are shown in Table 1. β-Actin and GAPDH were used as housekeeping genes, obtaining comparable Ct values for all the samples. For each target gene, a seven-point standard curve was included in each batch of amplifications based on 10-fold serial dilutions starting at 10 ng/µl of plasmid DNA. The relative quantification of the mRNA expression levels (x-fold change in expression) was carried out by the comparative 2−ΔΔCt method [63].

Table 1 Sequences of primers used for cytokine real-time PCR (qPCR) and standard curve data

Measurement of cytokines in supernatants of BCEC-1 and F3 cell cultures by ELISA

Protein concentrations of the cytokines that showed variations in the mRNA expression levels were determined in the culture supernatants at 4, 8, 24 and 56 hpi using commercial ELISA kits. The levels of IL-6, IL-8 and TNF-α cytokines were measured in the supernatants of the BCEC-1 and F3 cells by sandwich ELISAs using a Bovine IL-6 ELISA Reagent kit (ESS0029; Thermo Fisher Scientific, Waltham, MA, USA), Bovine IL-8 (CXCL8) ELISA Development kit (3114-1A-6; Mabtech AB, Stockholm, Sweden) and Bovine TNF-α ELISA kit (EBTNF; Thermo Fisher Scientific) following the manufacturers’ instructions. The sensitivity limits of these assays were 78 pg/ml for IL-6, 25 pg/ml for IL-8 and 100 pg/ml for TNF-α.

Statistical analysis

TLR, cytokine and endothelial adhesion molecule mRNA expression levels, as well as differences in the protein secretion between infected and control groups, were analysed using the non-parametric Kruskal–Wallis test, followed by Dunn’s multiple comparison test for all pairwise comparisons. In addition, to assess differences between both infected groups a Mann–Whitney test was performed for each molecule analysed. The statistical significance for all the analyses was established with P < 0.05. GraphPad Prism v.5.01 software (GraphPad Software, San Diego, CA, USA) was used to perform all statistical analyses and create all the graphical illustrations.

Abbreviations

BCEC-1:

bovine caruncular epithelial cell line 1

F3:

bovine placental trophoblast cell line

qPCR:

real-time polymerase chain reaction

hpi:

hours post-infection

BVD:

bovine viral diarrhea

DMEM:

Dulbecco’s modified Eagle medium

FCS:

foetal calf serum

PBS:

phosphate buffered saline

TLR-2:

Toll-like Receptor 2

IL:

interleukin

IFN:

interferon

TGF:

transforming growth factor

TNF:

tumoral necrosis factor

ICAM:

intercellular adhesion molecule

VCAM:

vascular cell adhesion molecule

ELISA:

enzyme-linked immunosorbent assay

BUVECs:

bovine umbilical vein endothelial cells

HPCVE:

human placental chorionic villi explants

MOI:

multiplicity of infection

References

  1. Innes EA, Wright S, Bartley P, Maley S, Macaldowie C, Esteban-Redondo I, et al. The host-parasite relationship in bovine neosporosis. Vet Immunol Immunopathol. 2005;108:29–36.

    Article  Google Scholar 

  2. Dubey JP, Buxton D, Wouda W. Pathogenesis of bovine neosporosis. J Comp Pathol. 2006;134:267–89.

    Article  CAS  Google Scholar 

  3. Dubey JP, Schares G, Ortega-Mora LM. Epidemiology and control of neosporosis and Neospora caninum. Clin Microbiol Rev. 2007;20:323–67.

    Article  CAS  Google Scholar 

  4. Williams DJ, Hartley CS, Bjorkman C, Trees AJ. Endogenous and exogenous transplacental transmission of Neospora caninum—how the route of transmission impacts on epidemiology and control of disease. Parasitology. 2009;136:1895–900.

    Article  CAS  Google Scholar 

  5. Entrican G. Immune regulation during pregnancy and host-pathogen interactions in infectious abortion. J Comp Pathol. 2002;126:79–94.

    Article  CAS  Google Scholar 

  6. Innes EA. The host-parasite relationship in pregnant cattle infected with Neospora caninum. Parasitology. 2007;134:1903–10.

    Article  Google Scholar 

  7. Quinn HE, Ellis JT, Smith NC. Neospora caninum: a cause of immune-mediated failure of pregnancy? Trends Parasitol. 2002;18:391–4.

    Article  CAS  Google Scholar 

  8. Rosbottom A, Gibney EH, Guy CS, Kipar A, Smith RF, Kaiser P, et al. Upregulation of cytokines is detected in the placentas of cattle infected with Neospora caninum and is more marked early in gestation when fetal death is observed. Infect Immun. 2008;76:2352–61.

    Article  CAS  Google Scholar 

  9. Regidor-Cerrillo J, Arranz-Solis D, Benavides J, Gomez-Bautista M, Castro-Hermida JA, Mezo M, et al. 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.

    Article  Google Scholar 

  10. Innes EA, Andrianarivo AG, Bjorkman C, Williams DJ, Conrad PA. Immune responses to Neospora caninum and prospects for vaccination. Trends Parasitol. 2002;18:497–504.

    Article  CAS  Google Scholar 

  11. Montes MJ, Tortosa CG, Borja C, Abadia AC, González-Gómez F, Ruiz C, et al. Constitutive secretion of interleukin-6 by human decidual stromal cells in culture. Regulatory effect of progesterone. Am J Reprod Immunol. 1995;34:188–94.

    Article  CAS  Google Scholar 

  12. Steinborn A, Von Gall C, Hildenbrand R, Stutte H, Kaufmann M. Identification of placental cytokine-producing cells in term and preterm labor. Obstet Gynecol. 1998;91:329–35.

    Article  CAS  Google Scholar 

  13. Steinborn A, Geisse M, Kaufmann M. Expression of cytokine receptors in the placenta in term and preterm labour. Placenta. 1998;19:165–70.

    Article  CAS  Google Scholar 

  14. Rojo-Montejo S, Collantes-Fernández E, Blanco-Murcia J, Rodríguez-Bertos A, Risco-Castillo V, Ortega-Mora LM. Experimental infection with a low virulence isolate of Neospora caninum at 70 days gestation in cattle did not result in foetopathy. Vet Res. 2009;40:49.

    Article  Google Scholar 

  15. Caspe SG, Moore DP, Leunda MR, Cano DB, Lischinsky L, Regidor-Cerrillo J, et al. The Neospora caninum-Spain 7 isolate induces placental damage, fetal death and abortion in cattle when inoculated in early gestation. Vet Parasitol. 2012;189:171–81.

    Article  CAS  Google Scholar 

  16. Dellarupe A, Regidor-Cerrillo J, Jimenez-Ruiz E, Schares G, Unzaga JM, Venturini MC, et al. 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.

    Article  CAS  Google Scholar 

  17. Jiménez-Pelayo L, García-Sánchez M, Regidor-Cerrillo J, Horcajo P, Collantes-Fernández E, Gómez-Bautista M, et al. Differential susceptibility of bovine caruncular and trophoblast cell lines to infection with high and low virulence isolates of Neospora caninum. Parasit Vectors. 2017;10:463.

    Article  Google Scholar 

  18. Bevilacqua E, Hoshida MS, Amarante-Paffaro A, Albieri-Borges A, Zago Gomes S. Trophoblast phagocytic program: roles in different placental systems. Int J Dev Biol. 2010;54:495–505.

    Article  Google Scholar 

  19. Mitsunari M, Yoshida S, Shoji T, Tsukihara S, Iwabe T, Harada T, et al. Macrophage-activating lipopeptide-2 induces cyclooxygenase-2 and prostaglandin E2 via toll-like receptor 2 in human placental trophoblast cells. J Reprod Immunol. 2006;72:46–59.

    Article  CAS  Google Scholar 

  20. Gillaux C, Mehats C, Vaiman D, Cabrol D, Breuiller-Fouche M. Functional screening of TLRs in human amniotic epithelial cells. J Immunol. 2011;187:2766–74.

    Article  CAS  Google Scholar 

  21. Marin MS, Hecker YP, Quintana S, Pérez S, Leunda MR, Cantón G, et al. Toll-like receptors 3, 7 and 8 are upregulated in the placental caruncle and fetal spleen of Neospora caninum experimentally infected cattle. Vet Parasitol. 2017;236:58–61.

    Article  CAS  Google Scholar 

  22. Koga K, Mor G. expression and function of Toll-like receptors at the maternal-fetal interface. Reprod Sci. 2008;15:231–42.

    Article  CAS  Google Scholar 

  23. Horcajo P, Jimenez-Pelayo L, Garcia-Sanchez M, Regidor-Cerrillo J, Collantes-Fernandez E, Rozas D, et al. Transcriptome modulation of bovine trophoblast cells in vitro by Neospora caninum. Int J Parasitol. 2017;47:791–9.

    Article  CAS  Google Scholar 

  24. Marin MS, Hecker YP, Quintana S, Pérez S, Leunda MR, Cantón G, et al. Immunization with inactivated antigens of Neospora caninum induces toll-like receptors 3, 7, 8 and 9 in maternal-fetal interface of infected pregnant heifers. Vet Parasitol. 2017;243:12–7.

    Article  CAS  Google Scholar 

  25. Liu J, Cao X. Cellular and molecular regulation of innate inflammatory responses. Cell Mol Immunol. 2016;13:711.

    Article  CAS  Google Scholar 

  26. Adams DH, Rlloyd A. Chemokines: leucocyte recruitment and activation cytokines. Lancet. 1997;349:490–5.

    Article  CAS  Google Scholar 

  27. Taubert A, Krull M, Zahner H, Hermosilla C. Toxoplasma gondii and Neospora caninum infections of bovine endothelial cells induce endothelial adhesion molecule gene transcription and subsequent PMN adhesion. Vet Immunol Immunopathol. 2006;112:272–83.

    Article  CAS  Google Scholar 

  28. Carvalho Neta AV, Stynen AP, Paixao TA, Miranda KL, Silva FL, Roux CM, et al. Modulation of the bovine trophoblastic innate immune response by Brucella abortus. Infect Immun. 2008;76:1897–907.

    Article  Google Scholar 

  29. Haider S, Knöfler M. Human tumour necrosis factor: physiological and pathological roles in placenta and endometrium. Placenta. 2009;30:111–23.

    Article  CAS  Google Scholar 

  30. Haraldsen G, Kvale D, Lien B, Farstad IN, Brandtzaeg P. Cytokine-regulated expression of E-selectin, intercellular adhesion molecule-1 (ICAM-1), and vascular cell adhesion molecule-1 (VCAM-1) in human microvascular endothelial cells. J Immunol. 1996;156:2558–65.

    CAS  PubMed  Google Scholar 

  31. Cavalcanti YV, Brelaz MC, Neves JK, Ferraz JC, Pereira VR. Role of TNF-alpha, IFN-gamma, and IL-10 in the development of pulmonary tuberculosis. Pulm Med. 2012;2012:745483.

    Article  Google Scholar 

  32. Robbins JR, Zeldovich VB, Poukchanski A, Boothroyd JC, Bakardjiev AI. Tissue barriers of the human placenta to infection with Toxoplasma gondii. Infect Immun. 2012;80:418–28.

    Article  CAS  Google Scholar 

  33. Yamane I, Kitani H, Kokuho T, Shibahara T, Haritani M, Hamaoka T, et al. The inhibitory effect of interferon gamma and tumor necrosis factor alpha on intracellular multiplication of Neospora caninum in primary bovine brain cells. J Vet Med Sci. 2000;62:347–51.

    Article  CAS  Google Scholar 

  34. Jesus EE, Pinheiro AM, Santos AB, Freire SM, Tardy MB, El-Bacha RS, et al. Effects of IFN-gamma, TNF-alpha, IL-10 and TGF-beta on Neospora caninum infection in rat glial cells. Exp Parasitol. 2013;133:269–74.

    Article  CAS  Google Scholar 

  35. Jauniaux E, Gulbis B, Schandene L, Collette J, Hustin J. Molecular interactions during pregnancy: distribution of interleukin-6 in maternal and embryonic tissues during the first trimester. Mol Hum Reprod. 1996;2:239–43.

    Article  CAS  Google Scholar 

  36. Diehl S, Rincón M. The two faces of IL-6 on Th1/Th2 differentiation. Mol Immunol. 2002;39:531–6.

    Article  CAS  Google Scholar 

  37. Pinheiro AM, Costa SL, Freire SM, Ribeiro CS, Tardy M, El-Bacha RS, et al. Neospora caninum: early immune response of rat mixed glial cultures after tachyzoites infection. Exp Parasitol. 2010;124:442–7.

    Article  CAS  Google Scholar 

  38. Almería S, Araujo RN, Darwich L, Dubey JP, Gasbarre LC. Cytokine gene expression at the materno-foetal interface after experimental Neospora caninum infection of heifers at 110 days of gestation. Parasite Immunol. 2011;33:517–23.

    Article  Google Scholar 

  39. Liao Y, Zhang Y, Liu X, Lu Y, Zhang L, Xi T, et al. Maternal murine cytomegalovirus infection during pregnancy up-regulates the gene expression of Toll-like receptor 2 and 4 in placenta. Curr Med Sci. 2018;38:632–9.

    Article  CAS  Google Scholar 

  40. Gayle DA, Beloosesky R, Desai M, Amidi F, Nuñez SE, Ross MG. Maternal LPS induces cytokines in the amniotic fluid and corticotropin releasing hormone in the fetal rat brain. Am J Physiol Regul Integr Comp Physiol. 2004;286:1024–9.

    Article  Google Scholar 

  41. Ashdown H, Dumont Y, Ng M, Poole S, Boksa P, Luheshi G. The role of cytokines in mediating effects of prenatal infection on the fetus: implications for schizophrenia. Mol Psychiatry. 2006;11:47.

    Article  CAS  Google Scholar 

  42. Beloosesky R, Gayle DA, Amidi F, Nunez SE, Babu J, Desai M, et al. N-Acetyl-cysteine suppresses amniotic fluid and placenta inflammatory cytokine responses to lipopolysaccharide in rats. Obstet Gynecol. 2006;194:268–73.

    CAS  Google Scholar 

  43. Arranz-Solís D, Benavides J, Regidor-Cerrillo J, Horcajo P, Castaño P, del Carmen Ferreras M, et al. Systemic and local immune responses in sheep after Neospora caninum experimental infection at early, mid and late gestation. Vet Res. 2016;47:2.

    Article  Google Scholar 

  44. Rosbottom A, Gibney H, Kaiser P, Hartley C, Smith RF, Robinson R, et al. Up regulation of the maternal immune response in the placenta of cattle naturally infected with Neospora caninum. PLoS ONE. 2011;6:e15799.

    Article  CAS  Google Scholar 

  45. Almería S, Serrano-Perez B, Darwich L, Domingo M, Mur-Novales R, Regidor-Cerrillo J, et al. Foetal death in naive heifers inoculated with Neospora caninum isolate Nc-Spain7 at 110 days of pregnancy. Exp Parasitol. 2016;168:62–9.

    Article  Google Scholar 

  46. Baszler TV, Long MT, McElwain TF, Mathison BA. Interferon-gamma and interleukin-12 mediate protection to acute Neospora caninum infection in BALB/c mice. Int J Parasitol. 1999;29:1635–46.

    Article  CAS  Google Scholar 

  47. Etienne-Manneville S, Chaverot N, Strosberg AD, Couraud PO. ICAM-1-coupled signaling pathways in astrocytes converge to cyclic AMP response element-binding protein phosphorylation and TNF-alpha secretion. J Immunol. 1999;163:66874.

    Google Scholar 

  48. Deisher TA, Haddix TL, Montgomery KF, Pohlman TH, Kaushansky K, Harlan JM. The role of protein kinase C in the induction of VCAM-1 expression on human umbilical vein endothelial cells. FEBS Lett. 1993;331:285–90.

    Article  CAS  Google Scholar 

  49. Silva LM, Vila-Viçosa MJ, Cortes HC, Taubert A, Hermosilla C. Suitable in vitro Eimeria arloingi macromeront formation in host endothelial cells and modulation of adhesion molecule, cytokine and chemokine gene transcription. Parasitol Res. 2015;114:113–24.

    Article  Google Scholar 

  50. Maksimov P, Hermosilla C, Kleinertz S, Hirzmann J, Taubert A. Besnoitia besnoiti infections activate primary bovine endothelial cells and promote PMN adhesion and NET formation under physiological flow condition. Parasitol Res. 2016;115:1991–2001.

    Article  CAS  Google Scholar 

  51. Zhang D, Chen L, Li S, Gu Z, Yan J. Lipopolysaccharide (LPS) of Porphyromonas gingivalis induces IL-1β, TNF-α and IL-6 production by THP-1 cells in a way different from that of Escherichia coli LPS. Innate Immun. 2008;14:99–107.

    Article  Google Scholar 

  52. Castillo C, Muñoz L, Carrillo I, Liempi A, Medina L, Galanti N, et al. Ex vivo infection of human placental chorionic villi explants with Trypanosoma cruzi and Toxoplasma gondii induces different Toll-like receptor expression and cytokine/chemokine profiles. Am J Reprod Immunol. 2017;78:12660.

    Article  Google Scholar 

  53. Castillo C, Muñoz L, Carrillo I, Liempi A, Medina L, Galanti N, et al. Toll-like receptor-2 mediates local innate immune response against Trypanosoma cruzi in ex vivo infected human placental chorionic villi explants. Placenta. 2017;60:40–6.

    Article  CAS  Google Scholar 

  54. Regidor-Cerrillo J, Gomez-Bautista M, Sodupe I, Aduriz G, Alvarez-Garcia G, Del Pozo I, et al. In vitro invasion efficiency and intracellular proliferation rate comprise virulence-related phenotypic traits of Neospora caninum. Vet Res. 2011;42:41.

    Article  Google Scholar 

  55. Almería S, Serrano-Perez B, López-Gatius F. Immune response in bovine neosporosis: protection or contribution to the pathogenesis of abortion. Microb Pathog. 2017;109:177–82.

    Article  Google Scholar 

  56. Regidor-Cerrillo J, Gómez-Bautista M, Pereira-Bueno J, Adúriz G, Navarro-Lozano V, Risco-Castillo V, et al. Isolation and genetic characterization of Neospora caninum from asymptomatic calves in Spain. Parasitology. 2008;135:1651–9.

    Article  CAS  Google Scholar 

  57. Rojo-Montejo S, Collantes-Fernández E, Regidor-Cerrillo J, Álvarez-García G, Marugan-Hernández V, Pedraza-Díaz S, et al. Isolation and characterization of a bovine isolate of Neospora caninum with low virulence. Vet Parasitol. 2009;159:7–16.

    Article  Google Scholar 

  58. Pérez-Zaballos FJ, Ortega-Mora LM, Álvarez-García G, Collantes-Fernández E, Navarro-Lozano V, García-Villada L, et al. Adaptation of Neospora caninum isolates to cell-culture changes: an argument in favor of its clonal population structure. J Parasitol. 2005;91:507–10.

    Article  Google Scholar 

  59. Bridger PS, Menge C, Leiser R, Tinneberg HR, Pfarrer CD. Bovine caruncular epithelial cell line (BCEC-1) isolated from the placenta forms a functional epithelial barrier in a polarised cell culture model. Placenta. 2007;28:1110–7.

    Article  CAS  Google Scholar 

  60. Hambruch N, Haeger JD, Dilly M, Pfarrer C. EGF stimulates proliferation in the bovine placental trophoblast cell line F3 via Ras and MAPK. Placenta. 2010;31:67–74.

    Article  CAS  Google Scholar 

  61. Menzies M, Ingham A. Identification and expression of Toll-like receptors 1–10 in selected bovine and ovine tissues. Vet Immunol Immunopathol. 2006;109:23–30.

    Article  CAS  Google Scholar 

  62. Puech C, Dedieu L, Chantal I, Rodrigues V. Design and evaluation of a unique SYBR Green real-time RT-PCR assay for quantification of five major cytokines in cattle, sheep and goats. BMC Vet Res. 2015;11:65.

    Article  Google Scholar 

  63. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3:1101–8.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the Spanish Ministry of Economy and Competitiveness (AGL2013-44694-R) and the Community of Madrid (PLATESA2-CM P2018/BAA-4370). Laura Jiménez-Pelayo was financially supported by a fellowship from the Complutense University of Madrid and Marta García-Sánchez was financially supported through a grant from the Spanish Ministry of Economy and Competitiveness (BES-2014-070723). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Availability of data and materials

Not applicable.

Authors’ contributions

JRC, PH, ECF, MGB and LMO conceived the study and participated in its design. LJP and MGS wrote the manuscript, with interpretation of results and discussion input from JRC, ECF, MGB, LMO, NH and CP. LJP and MGS performed in vitro infection of the cultures, collection of the samples and ELISA assays. JRC, PH, LJP and MGS designed and performed RT-qPCR analyses. NH and CP isolated bovine trophoblast and caruncular cell lines used in the assays. LJP and MGS carried out statistical analyses and interpreted the results. All authors read and approved the final manuscript.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Luis Miguel Ortega-Mora.

Additional file

Additional file 1:Table S1.

Statistical test results for mRNA expression levels. Table S2. Statistical test results for protein secretion.

Rights and permissions

Open Access This 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jiménez-Pelayo, L., García-Sánchez, M., Regidor-Cerrillo, J. et al. Immune response profile of caruncular and trophoblast cell lines infected by high- (Nc-Spain7) and low-virulence (Nc-Spain1H) isolates of Neospora caninum. Parasites Vectors 12, 218 (2019). https://doi.org/10.1186/s13071-019-3466-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-019-3466-z

Keywords