- Research
- Open access
- Published:
Identification of a novel methyltransferase-type 12 protein from Haemonchus contortus and its effects on functions of goat PBMCs
Parasites & Vectors volume 13, Article number: 154 (2020)
Abstract
Background
Methyltransferases (MTFs) are broad range of enzymes, which are ubiquitously expressed in diverse organisms ranging from bacteria to animals. MTFs proteins have been associated with various biological/cellular processes including transcriptional regulation, subcellular protein and RNA localization, signal transduction and DNA-damage repair. However, the role of MTFs in immune mechanism during host–parasite interaction has not been addressed yet.
Results
An open reading frame (764 bp) of methyltransferase-type 12 gene of H. contortus denoted as HcMTF-12, was successfully cloned using reverse transcriptase-polymerase chain reaction (RT-PCR) followed by prokaryotic expression in Escherichia coli BL21 (DE3 strain). The recombinant HcMTF-12 protein (rHcMTF-12) was about 47 kDa along with a fusion vector protein of 18 kDa. Immunoblot results identified the native protein MTF-12 with antibodies produced in rats against rHcMT-12, whereas rHcMTF-12 protein was recognized with sera of goat experimentally infected with H. contortus. Immunohistochemical analysis revealed that the native MTF-12 protein was mainly located in the periphery (cuticle) of parasite sections as well as within the pharynx and intestinal region. An immunofluorescence assay validated that rHcMTF-12 attached to the surface of goat PBMCs. Furthermore, the cytokines transcription of IL-2, IFN-γ and IL-4 transcripts of PBMCs incubated with rHcMTF-12 were enhanced in a dose-dependent manner. The secretion of TGF-β1 and IL-10 was significantly decreased. However, IL-6 production was not significantly different as compared to the control groups. Moreover, the migration activity and nitric oxide (NO) production by PBMCs were induced considerably, whereas the proliferation of PBMCs cells was negatively affected when incubated with the rHcMTF-12 protein.
Conclusions
Our findings suggest that HcMTF-12 significantly mediated the functions of PBMCs, and it might be a potential candidate for therapeutic interventions against haemonchosis.
Background
Methyltransferases (MTFs) are a broad range of enzymes, which are ubiquitously expressed in varied organisms ranging from bacteria to animals [1]. MTFs are membrane-bound enzymes that transfer a methyl group from S-adenosyl-l-methionine to thiol, amino or catechol acceptor molecules and this reaction modifies RNA, DNA and protein molecules for regulatory activities [2, 3]. Previously, more than 150 methyltransferase enzymes have been documented with structural differences but similar characteristics within the local active sites [4]. MTFs have been implicated in transcriptional regulation, subcellular protein and RNA localization, DNA-damage repair and signal transduction [5]. The multiple roles of methylated DNA in restriction-modification systems of prokaryotes and various cellular processes including gene differentiation and regulation in eukaryotes are well documented. However, apart from biochemical/structural properties, the role of methyltransferases participating in immune mechanism during host–parasite interaction has not been addressed yet. S-adenosyl-l-methionine (Ado-Met), a methyl donor in methyl transfer system has been suggested as an important pathway for metabolism of endogenous as well as xenobiotic compounds [6]. A previous study highlighted that methylation is an essential process to control gene transcription through which many drugs and small molecules undergo methylation, including nucleic acids, polysaccharides, lipids, proteins, hormones and neurotransmitters [7]. A recent research in the context of drug discovery and drug development revealed that MTFs are involved in using of S-Adenosyl methionine (SAM) analogs to methylate naturally produced anticancer agents which keeps alternate alkyl group as a substitute for methyl [8].
Haemonchus contortus is a parasite of veterinary importance globally by inducing huge economic and production losses due to severe anemia, hemorrhagic gastritis, oedema and even death of infected small ruminants [9]. After three successive developmental stages that occur in the environment, the third-stage (L3) larvae are taken into hostʼs digestive tract during grazing where fourth-stage larvae (L4) develop into mature adult worms feeding on blood from capillaries in the abomasum [10]. The widespread use of anthelmintic drugs as an effective method for control of H. contortus infections and generation of parasitic resistance to these commercially available drugs has become a serious problem worldwide [11,12,13]. During host–pathogen interactions, many proteins and other molecules of pathogens called excretory/secretory (ES) proteins target the host proteins to facilitate the infection process and modulate host immune response [14]. ES antigens of H. contortus possibly perform significant roles in infection process at the parasite-host interface [15].
Our previous proteomic analysis of H. contortus ES proteins bound to peripheral blood mononuclear cells (PBMCs) of goats revealed that methyltransferase type-12 is an important ES interacting antigen bound to goat immune cells at L4 and L5 developmental stages of H. contortus in vivo [16], indicating that this binding may have vital role in modulation of host immune responses during infection process. Here, we cloned and characterized the gene MTF-12 from H. contortus and describe the effects of the recombinant protein rHcMTF-12 on host PBMCs.
Methods
Experimental animals
Haemonchus contortus-free goats (local crossbreed) 3–6 months of age were housed indoors at the animal experimental station of Nanjing Agricultural University. Animals were given whole-shelled corn and hay, and water ad libitum. The goats have been dewormed with levamisole orally (8 mg/kg body weight) at a 2 week dose gap. In accordance with standard parasitological techniques, faecal samples were examined microscopically twice a week for helminth infection and goats having no sign of parasitic infection were used for the functional study. Three biological repeats (three goats), each comprising of three technical repeats (three replicates for each goat), were used in immunological functional experiments such as cytokine transcriptional analysis, immunofluorescence assays, nitric oxide (NO) production, cell proliferation and migration assays.
Specific anti-protein antibodies were produced in female Sprague-Dawley (SD) (~ 150 g body weight) rats bought from the Experimental Animal Center of Jiangsu Province, PR China (certified: SCXK 2008-0004). The rats were kept in a microbe-free state.
Parasite collection and mononuclear cell separation/culture
Helminth-free goats (6 months of age) from the research and teaching flock at Nanjing Agricultural University were administered orally with approximately 10,000 L3 larvae of H. contortus (Nanjing strain 2005). The microscopical examination of fecal samples was performed on a weekly basis for the existence of parasite infection. The collection and preservation of parasite eggs, larvae, and adult worms were performed in accordance with the previously described methods [17].
Blood samples were collected from dewormed healthy goats and cultured as previously described [18]. Briefly, whole blood was taken by venipuncture of the jugular vein into a 10 ml vacutainer coated with ethylenediaminetetraacetic acid (EDTA) and brought to laboratory for PBMCs isolation and culture. PBMCs were separated by Ficoll-hypaque (GE Healthcare, Munich, USA) using a density gradient centrifugation protocol [19]. After washing twice with phosphate-buffered saline (PBS: Ca2+/Mg2+-free, pH 7.4), cell viability assessed by the trypan blue exclusion test was > 95% in all experiments. The isolated PBMCs were then cultured in complete medium RPMI 1640 (Roswell Park Memorial Institute 1640) or DMEM (Dulbeccoʼs Modified Eagleʼs Medium) containing 100 U/ml penicillin and 100 mg/ml streptomycin together with 10% heat-inactivated fetal bovine serum (FBS), (Gibco, Paisley, UK). Before performing cell isolation, strict precautions were taken and all solutions were passed through a syringe-supported Millex®GP 0.22 μm filter (Merck Millipore Ltd., Cork, Ireland) to avoid contamination. For functional analysis, PBMCs were cultured in 24-well flat bottomed culture plates (Costar, Cambridge, MA, USA) at 37 °C in 5% CO2. To investigate the effects of rHcMTF-12 protein on certain functions of goat PBMCs, cultured cells were subjected to rHcMTF-12 protein with discrete concentrations of 5 μg/ml, 10 μg/ml, 20 μg/ml and 40 μg/ml. PBMCs treated with empty-pET32a protein (10 μg/ml) and treated with PBS (pH 7.4) were used as controls in all triplicate experiments.
RNA extraction and synthesis of 1st strand cDNA
Approximately 500 adult H. contortus parasites were crushed and homogenized at 4 °C in a pestle containing 1 ml Trizol (Invitrogen, Carlsbad, CA, USA), and total RNA was isolated under RNase-free conditions according to a previously described procedure [20]. Briefly, the homogenized mixture was transferred to sterile 2 ml Eppendrof tube and RNA was precipitated with the addition of 0.25 ml isopropyl alcohol per 1 ml of Trizol and centrifuged at 10,000×g for 10 min. After treatment with 70% ethanol in diethyl pyrocarbonate (DEPC), the pellet was finally dispensed in RNase/DNase free water and immediately used for cDNA synthesis. The RNA quantity was assessed using a spectrophotometry at optical density of 260 nm (OD260) and the quality was determined by the OD260/OD280 ratio. Only samples with an OD260/OD280 ratio between 1.8–2.0 were used. cDNA was synthesized using the PrimeScriptTM II 1st strand cDNA synthesis kit (TaKaRa Biotechnology, Shiga, Japan) as per the manufacturer’s instructions and stored at – 20 °C until further use [17].
Cloning and characterization of HcMTF-12 gene
The open reading frame (ORF) of the H. contortus MTF-12 gene (GenBank/Uniprot: CDJ87424.1/ U6P1M9) was used to design the sense primer (5ʹ-GAA TTC ATG CGC AAG GAC GCT G-3ʹ) and anti-sense primer (5ʹ-CGA GCT CTT AAC CTG CAA CCA CTG TC-3ʹ) by using Primer Premier v5 software. The EcoRI and SacI sites (underlined) were contained in the primers, respectively. The PCR reaction mixture (total volume of 50 μl) contained 50 μM 10× LA PCR buffer (Mg2+ Free) (TaKaRa Biotech, Dalian, China), 400 μM dNTP mixture, 3 mM MgCl2, 1 U Taq DNA polymerase (TaKaRa Biotech), 2 μl cDNA template and 400 nM forward and reverse primers. The PCR amplification cycles were as follows: 95 °C for 3 min followed by 40 cycles (at 95 °C for 35 s, 55 °C for 35 s, 72 °C for 60 s), followed by a final extension step at 72 °C for 8 min. The amplified PCR product was run on a 1% agarose G-10 (Biowest, Gene Co., Ltd., Chai Wan, Hong Kong) gel electrophoresis and purified with a gel purification kit (Omega, Norcross, GA, USA). The target fragment was ligated in the frame of the pMD-19T cloning vector (TaKaRa Biotech) followed by transformation into Trans5α competent cells (TransGen Biotech, Beijing, China). The amplified MTF-12 gene was subjected to restriction endonuclease digestion with EcoRI and SacI, and sequenced by GenScript (Nanjing, Jiangsu, China). Subsequently, the identified target gene (HcMTF-12) was cloned into pET-32a (+) prokaryotic expression vector (Novagen, Madison, USA) to produce pET-32a (+)-HcMTF-12 and ligated into E. coli BL21. Positive clones were digested with EcoRI and SacI and confirmed again by sequencing (GenScript, Nanjing, Jiangsu, China). The target gene was further subjected to bioinformatics analysis with available online sequences on GenBank (http://www.ncbi.nlm.nih.gov/BLAST/).
The comparison of the HcMTF-12 protein fragment with that of other parasites was carried out using BLASTx and BLASTp (https://blast.ncbi.nlm.nih.gov/Blast.cgi), and all sequences were aligned using ClustalW v1.8 software (http://www.clustal.org/). A phylogenetic tree was constructed using the neighbor-joining method in MEGA v6.0 program [21].
Purification of H. contortus MTF-12 recombinant protein
The bacterial culture containing pET-32a (+)-HcMTF-12 was re-cultured in 1 l Luria Bertini (LB) medium containing ampicillin (100 μg/ml) for 2–3 h at 37 °C and 120× g. After the OD600 of the culture reached 0.6, 1 mM isopropyl-b-d-thiogalactopyranoside (IPTG; Solarbio, Beijing, China) was added to induce protein expression for 7 h and samples were taken at different time points at 1-h intervals to check the ascending expression level. The cells were collected by centrifugation at 4 °C and 8000×g for 20 min and washed twice with PBS. The pellet was lysed with 10 μg/ml lysozyme (Sigma-Aldrich, Shanghai, China) followed by sonication at 600 W for 30 min and centrifuged at 8000×g for 10 min at 4 °C. The rHcMTF-12 protein was eluted using a Ni2+-nitrilotriacetic acid (Ni-NTA) column (GE Healthcare). Thereafter, purified protein was assessed on 12% (w/v) sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) followed by staining with R-250 Coomassie brilliant blue (Sigma-Aldrich). The fused protein of pET-32a with the 109aa Trx·Tag™ thioredoxin protein and six histidines was extracted using the same procedure described for the recombinant protein. The endotoxin level in samples was quantified using a Limulus Amebocyte Lysate (LAL) gel clot assay using a Pyrosate® Kit (Cape Cod Inc., East Falmouth, MA, USA). Finally, the proteins were quantified using the Bradford method [22], according to the described procedure, and stored at − 20 °C for functional analysis.
Antibody generation and confirmation by western blot
Polyclonal antibodies against rHcMTF-12 protein (anti-rHcMTF-12) were produced using 300 μg of rHcMTF-12 protein homogenized with an equivalent volume of Freund’s complete adjuvant injected intradermally into 2 female SD rats as described previously [23]. Two weeks after the primary dose, the rats were given three booster doses, seven days apart, using the same quantity and administration route as previously of rHcMTF-12 mixed with Freund’s incomplete adjuvant. Ten days after the last booster immunization, the rats were anesthetized with 25% isoflurane and sera with anti-rHcMTF-12 antibodies were separated from the blood and stored at − 20 °C for further applications. The sera against H. contortus infection were separated from blood of experimentally infected goats [17]. Sera separated from blood of normal goats/rats (un-immunized) were used as a control. The antibodies were quantified by ELISA.
For assessment of the specific reactivity with rHcMTF-12, an immuno-blot analysis was performed. Briefly, the rHcMTF-12 protein was separated on 12% SDS-PAGE and shifted onto Hybond® ECL™ nitrocellulose membranes (Sigma-Aldrich). Unnecessary binding sites were prevented by 5% skimmed milk mixed with Tris-buffer saline containing 0.1% Tween-20 (TBST). The membrane was washed thrice (5 min per wash) in TBST, and further incubated with the primary antibody (1: 300 dilution in TBST) anti-rHcMTF-12 at 4 °C overnight followed by the secondary antibody (1: 3000 dilution in TBST), horseradish peroxidase (HRP)-conjugated goat anti-rat IgG (Sigma-Aldrich), for 2 h 37 °C. The reactions were visualized by exposing to 3, 3-diaminobenzidine tetrahydrochloride (DAB) (Boster Bio-technology).
Immunofluorescence assay (IFA)
PBMCs from freshly collected goat blood were incubated with rHcMTF-12 protein (treatment group) or pET32a empty protein (positive control group) and PBS (negative control group) for 2 h at 37 °C with 5% CO2. PBMCs were fixed on poly-l-lysine-treated slides with 4% paraformaldehyde and permeabilized using 1% TritonX-100 in PBS. Thereafter, cells were treated with the primary antibody, rat anti-rHcMTF-12 IgG (1:300 dilution) or negative rat sera, followed by the secondary antibody (1:3000 dilution), Cy3-coupled goat anti-rat IgG (Beyotime Biotechnology, Shanghai, China) for 1 h at 37 °C in the dark. Cells were stained for 5 min in the dark with 1.5 μM 2-(4-amidinophenyl)-6-indolecarbamidine dihydrochloride (DAPI: Sigma-Aldrich) and visualized under a laser microscope (LSM710; Zeiss, Jena, Germany). Each step was performed after washing slides three times with PBS.
Distribution of MTF-12 protein in worms
Adult male/female H. contortus, suspended in Tissue-Tek® OCT (Sakura, Torrance, CA, USA), were solidified in liquid nitrogen. The samples were sliced into 10 μm thickness using Cryotome (CM1950, Wetzlar, Germany), and allowed to settle on poly-l-lysine hydrobromide slides (Merck Millipore, Darmstadt, Germany). The slides were incubated with first and secondary antibodies, rat-anti-HcMTF-12 (1:300 dilution) or negative rat sera and Cy3-labeled goat anti-rat IgG (1:3000 dilution) (Beyotime Biotechnology), respectively, for 2 h at 37 °C. Cell nuclei were stained with DAPI (Sigma-Aldrich) and examined under a confocal microscope (LSM710; Zeiss).
Cytokine transcriptions evaluated by real-time PCR
PBMCs were separated at a density of 1 × 106 cells/ml using Ficoll® Hypaque Lymphocyte Separation Medium (GE Healthcare). Then cells were added into a 24-well plate (1 ml/well) containing concanavalin A (ConA: 10 μg/ml) and cultured with control buffer (PBS) and pET-32a control protein (10 μg/ml) or rHcMTF-12 protein (5, 10, 20 and 40 μg/ml) at 37 °C and 5% humidity for 24 h. After washing with PBS, the cells were subjected to total RNA extraction with PrimeScriptTM II 1st strand cDNA synthesis kit (TaKaRa Biotech) following the manufacturer’s guidelines. The quantity and quality of extracted RNA was analyzed using a spectrophotometer (Eppendorf, Hamburg, Germany) at absorption values between 1.8 and 2.0 (OD260/OD280). The primers used for cytokine detection (IL-2, IL-4, IL-6, IL-10, TGF-β1 and IFN-γ) including an endogenous reference gene (β-actin), their amplification efficiencies as well as standard reaction conditions, are provided in Additional file 1: Tables S1–S3. Transcriptional analysis was carried out on an ABI 7500 real-time PCR system (Applied Biosystems, Foster City, USA). The data were analyzed based on raw cycle thresholds (Cq), obtained from the ABI Prism 7500 software [24] using the comparative Cq method (2−ΔΔ Cq method).
Cell proliferation assay
Anti-proliferative efficiency of rHcMTF-12 was assessed by culturing goat PBMCs in a 96-well tissue culture plate (Corning Life Sciences, New York, USA) as described previously [25]. Briefly, 100 µl of cells (1 × 106 cells/well) was activated with ConA and cultured with different rHcMTF-12 protein dilutions (5, 10, 20 and 40 µg/ml) or with pET-32a protein and PBS control, in triplicate. The mixture was cultured at 37 °C and 5% CO2 for 72 h, and 10 µl of CCK-8 solution (Beyotime Biotechnology) was mixed by pipetting in each well 4 h before harvesting. The optical density of the mixture was checked at 450 nm (OD450) using a spectrophotometer.
Cell migration activity of PBMCs
Migration assessment of the cells was evaluated using sterilized 8.0 μm pore polycarbonate membrane inserts [26]. Briefly, 200 μl of 1 × 106 cells in RPMI medium either rHcMTF-12-tretaed or untreated (control) were added into the Millicell® insert (Merck Millipore, Hayward, CA, USA) and placed on a Transwell tissue culture plate (Corning). Subsequently, the lower section was filled with 1300 μl of RPMI medium. The plate was incubated for 4 h in 5% CO2 at 37 °C. After incubation, PBMCs that moved through the membrane into the bottom of the well were determined using an automated cell counter (Countstar, Shanghai, China). The results are presented as percentages of the seeded PBMCs.
Nitric oxide assay
One hundred µl of 1 × 106 goat PBMCs per ml in DMEM medium were added into 96-well plate with varying concentrations of rHcMTF-12 (5, 10, 20 and 40 µg/ml) or pET32a and control buffer (PBS). After 48 h of cell culture at 37 °C and 5% CO2, the intracellular nitrite secretion of PBMCs was measured using the Griess assay [27] in accordance with total nitric oxide assay kit (Beyotime Biotechnology). Absorbance values were quantified on a spectrophotometer at 540 nm (OD540). Data were obtained from three independent experiments.
Data analysis
Data were presented as the mean ± standard error of the mean (SE). Statistical analysis for significant differences were carried out using a one-way analysis of variance (ANOVA) using Graphpad PrismTM v6 software (San Diego, California, USA), followed by a Tukeyʼs post-hoc test. The Student’s t-test was used for parametric samples and considered statistically significant at P < 0.05.
Results
Amplification and cloning of the HcMTF-12 gene
Amplification of the HcMTF-12 gene was successfully performed by PCR from H. contortus cDNA with specific primer pairs with the fragment detected at the desired size of 764 bp (Fig. 1). The product was subsequently cloned and confirmed as recombinant plasmid pMD-19T-HcMTF-12 by endonuclease enzymatic digestion with EcoRI and SacI (Fig. 1). Further validation was carried out by successful generation of recombinant expression plasmid pET-32a (+)-HcMTF-12 on 1% agarose gel electrophoresis, and sequencing confirmed its correct insertion into accurate reading frame (Fig. 1).
Sequence analysis using bioinformatics tools
The nucleotide and translated amino acid sequences of HcMTF-12 were analyzed using BLASTx and BLASTp, and compared with known protein sequences from other parasite species available on GenBank. Analysis of the amino acids sequences showed considerable identity (82%) with Teladorsagia circumcincta methyltransferase domain protein (GenBank: PIO70527.1), 79% with Ancylostoma ceylanicum SCP-like protein (GenBank: EPB79002.1), 78% with Necator americanus methyltransferase domain protein (GenBank: XP_013304755.1), 77% with Oesophagostomum dentatum methyltransferase domain protein (GenBank: KHJ75858.1), 74% with Dictyocaulus viviparus methyltransferase domain protein (GenBank: KJH40840.1), 60% with Caenorhabditis elegans methyltransferase-like protein (GenBank: NP_001040827.1) and 59% with Toxocara canis methyltransferase-like protein 2-A (GenBank: KHN87041.1) (Fig. 2). The distinctive characteristics of HcMTF-12 were verified as S-adenosyl-l-methionine binding sites in its structure (Fig. 2). The target sequence was further analyzed by multiple alignment with a number of published sequences and determined that the cloned ORF of HcMTF-12 belongs to the H. contortus Adenosyl-methionine superfamily. The phylogenic tree analysis using MEGA 5.10 software indicated that HcMTF-12 was accurately allied with MTF of homologous protein sequences obtained from NCBI and represented as bootstrap values for 1000 replicates (Fig. 2). The bootstrap consensus calculation showed high amino acid sequence similarity between HcMTF-12 and the available data on GenBank for T. circumcincta, A. ceylanicum, N. americanus, O. dentatum, D. viviparus, C. elegans and T. canis.
rHcMTF-12 protein expression and purification
The recombinant plasmid rHcMTF-12 inserted into E. coli strain BL21was expressed as a fusion protein after 0.8 mM IPTG induction and purified by chromatography on a Ni-NTA super column. The purified product was electrophoresed on a 12% SDS-PAGE and stained with Coomassie brilliant blue (Fig. 3). The rHcMTF-12 protein was expressed at a molecular weight of 47 kDa rather than calculated mass of 29 kDa due to the extra-size of the vector protein (pET32a), which was approximately 18 kDa fused with the target protein (Fig. 3).
Antibodies detected by western blot analysis
Immunoblotting was utilized to detect polyclonal antibodies against rHcMTF-12, which showed that rHcMTF-12 could be detected with antibodies in sera of goats infected with H. contortus (Fig. 3). The native MTF-12 antigen in whole soluble extracts of H. contortus was distinguished with antibodies raised in SD rats immunized with rHcMTF-12 and the presence of a single band indicated that these antibodies had only specificity against the MTF-12 protein of H. contortus (Fig. 3). However, sera taken from un-immunized goats/rats did not show any antibodies against the protein.
Detection of rHcMTF-12 on the surface of goat PBMCs
The interactions of rHcMTF-12 with host immune cells were detected by immunofluorescence assay (IFA) according to previously described procedures [17]. Results showed red fluorescence at the periphery of incubated PBMCs thus confirming Cy3-labeled rHcMTF-12 binding in the treatment group along with the DAPI labelled cell nuclei displaying in blue color. However, in the control groups, no red visualization was detected either in cells incubated with pET32a control protein or the PBS control. These results indicate that rHcMTF-12 could attach on PBMCs surface in shape of surface ligand complex and suppress or modulate the immune functions of host PBMCs (Fig. 4).
Distribution of MTF-12 protein within the worm
As shown in Fig. 5, adult male and female worms were used to detect the target protein within partial body sections (not whole body length) by immunohistochemical assay. The distribution of the native MTF-12 protein within the worm was indicated as red fluorescence in the cuticle as well as within the pharynx and intestinal region in treatment groups. However, in control sections, no protein localization was detected (Fig. 5).
rHcMTF-12 mediated cytokine expression of goat PBMCs
The mRNA levels of cytokines in PBMCs cultured with rHcMTF-12 protein were evaluated, and the results showed that transcriptions of interleukin-2 (IL-2) (ANOVA: F(5, 12) = 243.1, P < 0.0001), IL-4 (ANOVA: F(5, 12) = 28.31, P < 0.0001) and IFN-γ (ANOVA: F(5, 12) = 266.5, P < 0.0001) were augmented at a considerable level in a dose-dependent manner. The secretions of TGF-β1 (ANOVA: F(5, 12) = 18.98, P < 0.0001) and IL-10 (ANOVA: F(5, 12) = 117.2, P < 0.0001) were significantly decreased in PBMCs incubated with rHcMTF-12 in comparison with the pET32a control group and PBS control group. However, IL-6 was not significantly different (ANOVA: F(5, 12) = 0.410, P < 0.833) when compared to the control group (Fig. 6).
rHcMTF-12 downregulated multiplication of goat PBMCs
The anti-proliferative effect of rHcMTF-12, in comparison with the PBS control and pET32a protein groups, was estimated using a cell counting kit (CCK8). Our results highlighted that PBMCs activated with ConA and treated with rHcMTF-12 protein downregulated PBMCs multiplication efficiency at 5 µg/ml (ANOVA: F(5, 12) = 30.44, P = 0.0014), 10 µg/ml (ANOVA: F(5, 12) = 30.44, P = 0.0012), 20 µg/ml (ANOVA: F(5, 12) = 30.44, P = 0.0013) and 40 µg/ml (ANOVA: F(5, 12) = 30.44, P = 0.0002) multiple protein concentrations dose-dependently (Fig. 7). However, cell proliferation was significantly unchanged between the pET32a group and the control buffer group (ANOVA: F(5, 12) = 30.44, P = 0.692).
rHcMTF-12 promoted migration efficiency of PBMCs
The change in migratory efficiency of host PBMCs in response to rHcMTF-12 was assessed by the percentile of shifted PBMCs through the Millipore polycarbonate membrane into the lower chamber. Our results revealed that no significant cell movement was achieved in the control buffer and protein (pET32a) group (ANOVA: F(5, 12) = 12.51, P = 0.783). However, in treatment groups of 5 µg/ml (29.33 ± 1.453%) (P = 0.003), 10 µg/ml (28.00 ± 2.082%) (P = 0.008), 20 µg/ml (29.00 ± 2.887%) (P = 0.0014) and 40 µg/ml (25.67 ± 1.202%) (P = 0.014) of rHcMTF-12, a significant percentage of cells migrated through the membrane (Fig. 8).
rHcMTF-12 mediated NO production in goat PBMCs
Using the Griess assay, no significant NO production occurred between control group and 5 µg/ml protein concentration (ANOVA: F(5, 12) = 98.41, P = 0.716). The rHcMTF-12 treatment groups showed a significant difference in NO production by PBMCs at 10 µg/ml (ANOVA: F(5, 12) = 98.41, P < 0.001), 20 µg/ml (ANOVA: F(5, 12) = 98.41 P < 0.001) and 40 µg/ml (ANOVA: F(5, 12) = 98.41, P < 0.001) when compared to the control group. Meanwhile, there was no significant difference between the PBS and pET32a groups (ANOVA: F(5, 12) = 98.41, P = 0.682) (Fig. 9).
Discussion
S-Adenosyl methionine (SAM)-dependent methyltransferases genes have been implicated in numerous biological processes through various cellular pathways, which regulate gene expression and/or functions [28]. However, nothing is known about methyltransferases regulation of immunity during host–parasite interactions at the cellular level. Here, we report the host cell binding protein of H. contortus parasite (HcMTF-12) and its immuno-regulatory characteristics with host PBMCs in vitro for the first time.
Previously, it was reported that a large class of SAM-dependent methyltransferases shared a conserved catalytic domain structurally in interaction with S-adenosylmethionine as a common cofactor [29]. In addition, nine AdoMet-dependent methyltransferases had been found to have same folding arrangement in a central parallel β-sheet encircled by α-helicals [30,31,32], and a large number of AdoMet-dependent methyltransferases shared conserved sequence motifs in their three-dimensional structure [33,34,35]. In this study, we cloned a SAM-dependent methyltransferase gene from H. contortus (HcMTF-12), which showed highly similar sequence identity and structural characteristics among various sequences of methyltransferases enzymes from other parasitic species. The size of the native HcMTF-12 was c.29 kDa in somatic extracts of H. contortus detected by the antibodies of rat against rHcMTF-12 and the molecular weight of rHcMTF-12 was c.47 kDa. Apart from the Tag peptides of 18 kDa of the vector, the size of the rHcMTF-12 was identical to that of the native one. Immunoblotting showed that rHcMTF-12 was identified with sera from experimentally infected with H. contortus goats. The anti-rHcMTF-12 antibodies used to probe the whole soluble proteins of H. contortus produced only a single band, which indicated that these antibodies did not show cross-reactivity with other proteins in this parasite. Our results suggest that these antibodies had specificity against the MTF-12 protein of H. contortus. These results indicate that the native MTF-12 could enter into the circulatory system of the goat and be recognized by the host immune system.
The ES proteins are released at the host–parasite interface and hence the protein-cell interaction is the primary requirement for immune functions [25]. A previous study demonstrated that methyltransferase from Plasmodium falciparum localized in membrane structures where it proposed to play an important role in parasite membrane biogenesis [36]. In the present study, we detected HcMTF-12 protein interaction on the surface of host PBMCs using IFA. Its localization was detected within the parasite structure by immunohistochemical analysis, which confirmed that HcMTF-12 protein was mainly located in the cuticle as well as within the pharynx and intestinal region.
T helper cells 1 (Th1) and Th2 along with inflammatory responses are linked with production of numerous cytokines that play vital anti-parasitic roles against helminths particularly H. contortus infection [37]. Previously, it was noted that a type-2 immune response was responsible to helminth infections, which were usually associated with IL-10, IL-5 and IL-4 secretions [38]. In addition, interleukin 10 (IL-10), a key element generated through inducible T regulatory cells (Treg), was responsible to prevent production of allergic Th2 responses [39]. In our study, IL-4 secretion by PBMCs in response to rHcMTF-12 protein was increased, while the secretion of the immunosuppressive IL-10 was decreased. This indicates that rHcMTF-12 might be a protective antigen which diminished the effects of Treg cells and induced Th2 humoral responses against H. contortus infection.
Apart from humoral responses, cellular immunity is another important mechanism against parasitic infection regulated by Th1 cells which produce IFN-γ [40, 41]. Moreover, progression of Th2 cells in immune responses was found to be negatively regulated by IFN-γ [42]. According to Schallig et al. [43], sheep experimentally infected with H. contortus retorted to a cell mediated immune response, which was characterized with high levels of IL-2 and IFN-γ. In the present study, higher expressions of IFN-γ and IL-2 cytokines were detected in goat PBMCs after incubation with rHcMTF-12 protein. In this regard, a higher level of IL-4 and IFN-γ suggested a balanced Th1 and Th2 environment which might either be a collective effect of immune cells (PBMCs) or an individual cell playing a distinctive role; this needs to be investigated further.
Transforming growth factor beta 1 (TGF-β1) is a multifunctional interleukin which differentially operates several immune system events by various cell types and potentially regulates many biological functions [44]. Moreover, it was suggested that TGF-β1 also played a dynamic role in inhibition of cell multiplication and stimulation of programmed cell death of numerous immune cell subsets [45]. The production of TGF-β1 after incubation of goat PBMCs with rHcMTF-12 increased at a significant level, which might be another mechanism of rHcMTF-12 antigen to promote immune evasion at the parasite-host interface.
Interleukin 6 (IL-6), a pro-inflammatory cytokine, was released in reaction to infection and contributed in the host defense mechanism by generation of an acute stage immune response. In this study, the IL-6 did not show any changes indicating that rHcMTF-12 protein could not stimulate IL-6 production.
Modulation of immune responses during the infection process is also linked with active stimulation and recruitment of eosinophils and other lymphocytes by helminths, resulting in stability of worms within the host [46, 47]. Previously, we determined stimulatory roles of HcESPs collectively and individually in the case of rHcftt-2 and rHcES-24 to induce the migration of goat PBMCs [20, 48, 49]. Here, we also found that rHcMTF-12, as an interacting protein of HcESPs, showed significant enhancement ability to the migration of host PBMCs. These results indicate that rHcMTF-12 protein might be an active contributor of HcESPs in the modulation of PBMCs migration. However, the underlying mechanism involved behind this strategy needs to be further investigated.
Proliferation of immune cells plays an important role in regulation of immune responses during helminth infection [50, 51]. Previously, we found a suppressive effect of various H. contortus ES proteins on PBMC multiplication [20, 48, 49, 52]. In accordance with earlier studies, the present study revealed that rHcMTF-12 could also significantly inhibit the proliferation of goat PBMCs.
Previous molecular approaches suggested that protein methyltransferase is actively involved in regulation of NO release by endothelial cells [53]. NO was involved in non-specific host resistance and cytotoxic activities against a variety of parasitic infections, and its release could be enhanced by IFN-γ [54, 55]. The results from the present study implicate that rHcMTF-12 could stimulate the NO production. This result together with the increased level of IFN-γ, suggest that rHcMTF-12 might be an important contributor in the host defense mechanism against H. contortus infection.
Conclusions
We first cloned and characterized HcMTF-12, an important ES antigen of H. contortus. Our results show that rHcMTF-12 protein could bind to the host immune cells and modulate a variety of immune functions, including upregulation of IFN-γ, IL-4 and IL-2, and enhancing PBMCs migratory activity and NO production. However, IL-10 and TGF-β1 production and cell multiplication were suppressed at different rHcMTF-12 protein doses. These results indicate that the HcMTF-12 protein plays a modulatory role on multiple functions of host PBMCs. However, the underlying regulatory mechanisms of HcMTF-12 in vivo during host–parasite interactions need to be explored in future studies.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and in Additional file 1: Tables S1–S3.
Abbreviations
- MTF-12:
-
methyltransferase-type 12
- LB:
-
Luria Bertani medium
- IPTG:
-
isopropyl-Bd-thiogalactopyranoside
- Ni-NTA:
-
Ni2+-nitrilotriacetic acid
- SDS-PAGE:
-
sodium dodecyl sulfate polyacrylamide gel electrophoresis
- PVDF:
-
polyvinylidene difluoride
- IgG-HRP:
-
horseradish peroxidase labeled immunoglobulin G
- DAB:
-
3, 3-diaminobenzidine tetrahydrochloride
- DEPC:
-
diethyl pyrocarbonate
- OCT:
-
optimal cutting temperature
- BSA:
-
bovine serum albumin
- PBMCs:
-
peripheral blood mononuclear cell
- PBS:
-
phosphate buffered saline
- RPMI 1640:
-
Roswell park memorial institute 1640 culture media
- FBS:
-
fetal bovine serum
- con A:
-
concanavalin A
- Cy3:
-
cyanine dyes 3
- DAPI:
-
2-(4-amidinophenyl)-6-indolecarbamidine dihydrochloride
- Th1:
-
helper T cell 1
- Th2:
-
helper T cell 2
- IL-10:
-
interleukin-10
- IL-6:
-
interleukin-6
- IL-4:
-
interleukin-4
- IL-2:
-
interleukin-2
- IFN-γ:
-
interferon-γ
- TGF-β1:
-
transforming growth factor-β1
- IL-17:
-
interleukin-17
References
Posfai J, Bhagwat AS, Posfai G, Roberts RJ. Predictive motifs derived from cytosine methyltransferases. Nucleic Acids Res. 1989;17:2421–35.
Weisiger RA, Jakoby WB. Thiol S-methyltransferase. Methods Enzymol. 1981;77:257–62.
Lyon ES, Jakoby WB. Arylamine N-methyltransferase. Methods Enzymol. 1981;77:263–6.
Petrossian T, Clarke S. Bioinformatic identification of novel methyltransferases. Epigenomics. 2009;1:163–75.
Bedford MT. Arginine methylation at a glance. J Cell Sci. 2007;120:4243–6.
Sharma A, Gerbarg P, Bottiglieri T, Massoumi L, Carpenter LL, Lavretsky H, et al. S-Adenosylmethionine (SAMe) for neuropsychiatric disorders: a clinician-oriented review of research. J Clin Psychiatry. 2017;78:656–67.
Bottiglieri T, Laundy M, Martin R, Carney MWP, Nissenbaum H, Toone BK, et al. S-adenosylmethionine influences monoamine metabolism. Lancet. 1984;324:224.
Singh S, Zhang J, Huber TD, Sunkara M, Hurley K, Goff RD, et al. Facile chemoenzymatic strategies for the synthesis and utilization of S-adenosyl-(l)-methionine analogues. Angew Chem Int Ed Engl. 2014;53:3965–9.
Padula A. Gastrointestinal nematodes of sheep and cattle: biology and control. Aust Vet J. 2011;89:54.
Sutherland I, Scott I. Gastrointestinal nematodes of sheep and cattle: biology and control. Chicester: Wiley; 2009.
Mortensen LL, Williamson LH, Terrill TH, Kircher RA, Larsen M, Kaplan RM. Evaluation of prevalence and clinical implications of anthelmintic resistance in gastrointestinal nematodes in goats. J Am Vet Med Assoc. 2003;223:495–500.
Kaplan RM, Vidyashankar AN. An inconvenient truth: global worming and anthelmintic resistance. Vet Parasitol. 2012;186:70–8.
Gilleard JS. Understanding anthelmintic resistance: the need for genomics and genetics. Int J Parasitol. 2006;36:1227–39.
Schmid-Hempel P. Parasite immune evasion: a momentous molecular war. Trends Ecol Evol. 2008;23:318–26.
Mohandas N, Young ND, Jabbar A, Korhonen PK, Koehler AV, Amani P, et al. The barber’s pole worm CAP protein superfamily - a basis for fundamental discovery and biotechnology advances. Biotechnol Adv. 2015;33:1744–54.
Gadahi JA, Wang S, Bo G, Ehsan M, Yan R, Song X, et al. Proteomic analysis of the excretory and secretory proteins of Haemonchus contortus (HcESP) binding to goat PBMCs in vivo revealed stage-specific binding profiles. PLoS ONE. 2016;11:e0159796.
Ehsan M, Gao W, Gadahi JA, Lu M, Liu X, Wang Y, et al. Arginine kinase from Haemonchus contortus decreased the proliferation and increased the apoptosis of goat PBMCs in vitro. Parasites Vectors. 2017;10:311.
Yanming S, Ruofeng Y, Muleke CI, Guangwei Z, Lixin X, Xiangrui L. Vaccination of goats with recombinant galectin antigen induces partial protection against Haemonchus contortus infection. Parasite Immunol. 2007;29:319–26.
Nicholson IC, Mavrangelos C, Fung K, Ayhan M, Levichkin I, Johnston A, et al. Characterisation of the protein composition of peripheral blood mononuclear cell microsomes by SDS-PAGE and mass spectrometry. J Immunol Methods. 2005;305:84–93.
Gadahi JA, Ehsan M, Wang S, Zhang Z, Wang Y, Yan R, et al. Recombinant protein of Haemonchus contortus 14-3-3 isoform 2 (rHcftt-2) decreased the production of IL-4 and suppressed the proliferation of goat PBMCs in vitro. Exp Parasitol. 2016;171:57–66.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6. Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72:248–54.
Han K, Xu L, Yan R, Song X, Li X. Molecular cloning, expression and characterization of enolase from adult Haemonchus contortus. Res Vet Sci. 2012;92:259–65.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25:402–8.
Wang W, Wang S, Zhang H, Yuan C, Yan R, Song X, et al. Galectin Hco-gal-m from Haemonchus contortus modulates goat monocytes and T cell function in different patterns. Parasites Vectors. 2014;7:342.
Ehsan M, Wang W, Gadahi JA, Hasan MW, Lu M, Wang Y, et al. The serine/threonine-protein phosphatase 1 from Haemonchus contortus is actively involved in suppressive regulatory roles on immune functions of goat peripheral blood mononuclear cells. Front Immunol. 2018;9:1627.
Sun J, Zhang X, Broderick M, Fein H. Measurement of nitric oxide production in biological systems by using Griess reaction assay. Sensors. 2003;3:276.
Siedlecki P, Zielenkiewicz P. Mammalian DNA methyltransferases. Acta Biochimica Polonica. 2006;53:245–56.
Schluckebier G, O’Gara M, Saenger W, Cheng X. Universal catalytic domain structure of AdoMet-dependent methyltransferases. J Mol Biol. 1995;247:16–20.
Cheng X, Kumar S, Posfai J, Pflugrath JW, Roberts RJ. Crystal structure of the Hhal DNA methyltransferase complexed with S-adenosyl-l-methionine. Cell. 1993;74:299–307.
Vidgren J, Svensson LA, Liljas A. Crystal structure of catechol O-methyltransferase. Nature. 1994;368:354.
Reinisch KM, Chen L, Verdine GL, Lipscomb WN. The crystal structure of Haelll methyltransferase covalently complexed to DNA: an extrahelical cytosine and rearranged base pairing. Cell. 1995;82:143–53.
Ingrosso D, Fowler AV, Bleibaum J, Clarke S. Sequence of the d-aspartyl/l-isoaspartyl protein methyltransferase from human erythrocytes. Common sequence motifs for protein, DNA, RNA, and small molecule S-adenosylmethionine-dependent methyltransferases. J Biol Chem. 1989;264:20131–9.
Kagan RM, Clarke S. Widespread occurrence of three sequence motifs in diverse S-Adenosylmethionine-dependent methyltransferases suggests a common structure for these enzymes. Arch Biochem Biophys. 1994;310:417–27.
Malone T, Blumenthal RM, Cheng X. Structure-guided analysis reveals nine sequence motifs conserved among DNA amino-methyl-transferases, and suggests a catalytic mechanism for these enzymes. J Mol Biol. 1995;253:618–32.
Witola WH, Pessi G, El Bissati K, Reynolds JM, Mamoun CB. Localization of the phosphoethanolamine methyltransferase of the human malaria parasite Plasmodium falciparum to the Golgi apparatus. J Biol Chem. 2006;281:21305–11.
Shakya KP, Miller JE, Lomax LG, Burnett DD. Evaluation of immune response to artificial infections of Haemonchus contortus in gulf coast native compared with Suffolk lambs. Vet Parasitol. 2011;181:239–47.
Munoz-Guzman MA, Cuenca-Verde C, Valdivia-Anda G, Cuellar-Ordaz JA, Alba-Hurtado F. Differential immune response between fundic and pyloric abomasal regions upon experimental ovine infection with Haemonchus contortus. Vet Parasitol. 2012;185:175–80.
Grencis RK, Humphreys NE, Bancroft AJ. Immunity to gastrointestinal nematodes: mechanisms and myths. Immunol Rev. 2014;260:183–205.
Murata Y, Ohteki T, Koyasu S, Hamuro J. IFN-gamma and pro-inflammatory cytokine production by antigen-presenting cells is dictated by intracellular thiol redox status regulated by oxygen tension. Eur J Immunol. 2002;32:2866–73.
Cope A, Le Friec G, Cardone J, Kemper C. The Th1 life cycle: molecular control of IFN-γ to IL-10 switching. Trends Immunol. 2011;32:278–86.
Coomes SM, Pelly VS, Kannan Y, Okoye IS, Czieso S, Entwistle LJ, et al. IFNγ and IL-12 restrict Th2 responses during helminth/plasmodium co-infection and promote IFNγ from Th2 cells. PLoS Pathog. 2015;11:e1004994.
Schallig HD. Immunological responses of sheep to Haemonchus contortus. Parasitology. 2000;120:S63–72.
Li MO, Wan YY, Sanjabi S, Robertson AK, Flavell RA. Transforming growth factor-beta regulation of immune responses. Ann Rev Immunol. 2006;24:99–146.
Lebman DA, Edmiston JS. The role of TGF-beta in growth, differentiation, and maturation of B lymphocytes. Microbes Infect. 1999;1:1297–304.
Turner DG, Wildblood LA, Inglis NF, Jones DG. Characterization of a galectin-like activity from the parasitic nematode, Haemonchus contortus, which modulates ovine eosinophil migration in vitro. Vet Immunol Immunopathol. 2008;122:138–45.
McGovern KE, Wilson EH. Role of chemokines and ttrafficking of immune cells in parasitic infections. Curr Immunol Rev. 2013;9:157–68.
Gadahi JA, Yongqian B, Ehsan M, Zhang ZC, Wang S, Yan RF, et al. Haemonchus contortus excretory and secretory proteins (HcESPs) suppress functions of goat PBMCs in vitro. Oncotarget. 2016;7:35670–9.
Gadahi JA, Li B, Ehsan M, Wang S, Zhang Z, Wang Y, et al. Recombinant Haemonchus contortus 24 kDa excretory/secretory protein (rHcES-24) modulate the immune functions of goat PBMCs in vitro. Oncotarget. 2016;7:83926–37.
Allen JE, MacDonald AS. Profound suppression of cellular proliferation mediated by the secretions of nematodes. Parasite Immunol. 1998;20:241–7.
Loke PN, MacDonald AS, Robb A, Maizels RM, Allen JE. Alternatively activated macrophages induced by nematode infection inhibit proliferation via cell-to-cell contact. Eur J Immunol. 2000;30:2669–78.
Wen Y, Wang Y, Wang W, Lu M, Ehsan M, Tian X, et al. Recombinant Miro domain-containing protein of Haemonchus contortus (rMiro-1) activates goat peripheral blood mononuclear cells in vitro. Vet Parasitol. 2017;243:100–4.
Han Z, Druhan LJ, Cardounel AJ, Chicoine LG. Regulation of endothelial nitric oxide production by protein methyltransferase’s. FASEB J. 2011;25:lb490.
Kotze AC, McClure SJ. Haemonchus contortus utilises catalase in defence against exogenous hydrogen peroxide in vitro. Int J Parasitol. 2001;31:1563–71.
Wandurska-Nowak E. The role of nitric oxide (NO) in parasitic infections. Wiad Parazytol. 2004;50:665–78.
Acknowledgements
We are thankful to all team members at the MOE Joint International Research Laboratory of Animal Health and Food Safety, College of Veterinary Medicine, Nanjing Agricultural University, for their valuable suggestions and help.
Funding
This work was funded by the “National Key Basic Research Program (973 Program) of P.R. China” (Grant No. 2015CB150300) and by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).
Author information
Authors and Affiliations
Contributions
XL was engaged in project administration, funding acquisition and supervision. XQZ participated in management and coordination throughout the study. ME prepared original draft of manuscript and conducted experiments. JAG and ML performed data curation and analysis. YW and MWH collected blood samples and helped in some immunological experiments. TL and MH provided input in computational analysis. RY, XS and LX provided all research resources during the study period. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All animal and laboratory experiments were strictly executed with guidelines and instructions of Animal Ethics Committee, Nanjing Agricultural University, China, and abided by the guidelines of the Animal Welfare Council of China. All experimental procedures were approved by the Science and Technology Agency of Jiangsu Province (ID: SYXK (SU) 2010-0005).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1: Table S1.
Primers sequences for real-time PCR. Table S2. Reagents used for real-time PCR. Table S3. Thermal cycling conditions for real-time PCR.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Ehsan, M., Gadahi, J.A., Liu, T. et al. Identification of a novel methyltransferase-type 12 protein from Haemonchus contortus and its effects on functions of goat PBMCs. Parasites Vectors 13, 154 (2020). https://doi.org/10.1186/s13071-020-04028-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13071-020-04028-y