Functional characterization of acyl-CoA binding protein in Neospora caninum
Parasites & Vectors volume 13, Article number: 85 (2020)
Lipid metabolism is pivotal for the growth of apicomplexan parasites. Lipid synthesis requires bulk carbon skeleton acyl-CoAs, the transport of which depends on the acyl-CoA binding protein (ACBP). In Neospora caninum, the causative agent of neosporosis, the FASII pathway is required for growth and pathogenicity. However, little is known about the fatty acid transport mechanism in N. caninum.
We have identified a cytosolic acyl-CoA binding protein, with highly conserved amino acid residues and a typical acyl-CoA binding domain in N. caninum. The recombinant NcACBP protein was expressed to verify the binding activities of NcACBP in vitro, and the heterologous expression of NcACBP in Δacbp yeast in vivo. Lipid extraction from ΔNcACBP or the wild-type of N. caninum was analyzed by GC-MS or TLC. Furthermore, transcriptome analysis was performed to compare the gene expression in different strains.
The NcACBP recombinant protein was able to specifically bind acyl-CoA esters in vitro. A yeast complementation assay showed that heterologous expression of NcACBP rescued the phenotypic defects in Δacbp yeast, indicating of the binding activity of NcACBP in vivo. The disruption of NcACBP did not perturb the parasite’s growth but enhanced its pathogenicity in mice. The lipidomic analysis showed that disruption of NcACBP caused no obvious changes in the overall abundance and turnover of fatty acids while knockout resulted in the accumulation of triacylglycerol. Transcriptional analysis of ACBP-deficient parasites revealed differentially expressed genes involved in a wide range of biological processes such as lipid metabolism, posttranslational modification, and membrane biogenesis.
Our study demonstrated that genetic ablation of NcACBP did not impair the survival and growth phenotype of N. caninum but enhanced its pathogenicity in mice. This deletion did not affect the overall fatty acid composition but modified the abundance of TAG. The loss of NcACBP resulted in global changes in the expression of multiple genes. This study provides a foundation for elucidating the molecular mechanism of lipid metabolism in N. caninum.
Neospora caninum is an apicomplexan parasite responsible for neosporosis, a disease characterized by abortions and stillbirths (mainly among cattle) and by neuromuscular diseases in dogs [1, 2], that results in substantial economic losses to the beef and dairy industries worldwide .
Lipids are crucial to the biology of all cells and organisms, acting not only as primary sources of energy but also as regulators of metabolism and growth, participating in various signaling networks . Lipids are the major structural elements of all biological membranes, serving as signaling molecules within and between cells. As a highly efficient store, the lipids could reduce the energy of all cells and organisms [5, 6]. Lipids can also act as pivotal pathogenetic factors that allow pathogens to escape immune responses, manipulate host processes, and develop disease .
Fatty acids provide the acyl skeleton for lipid synthesis. In apicomplexan parasites, the FASII de novo synthesis pathway is localized to a special metabolic organelle, the apicoplast, and is responsible for the production of long-chain fatty acids, processed for fatty acids elongation and desaturation in the ER, such as myristic acid and palmitic acid [8, 9]. In addition, parasites can scavenge lipids from the host cells and the surrounding environment . Before entering lipid synthesis pathways, different kinds of acyl-CoA transporters which acts as the major carrier of acyl-CoAs, such as fatty acid binding protein (FABP), sterol carrier protein 2 (SCP2) and acyl-CoA-binding protein (ACBP), would activate and convert fatty acids to fatty acyl-CoA esters via a reaction catalyzed by fatty acyl-CoA synthetase and transported to various metabolic locations . ACBP is a protein of approximately 10 kDa that is ubiquitously expressed and highly conserved in tissues with active lipid metabolism such as liver and adipose tissues  from humans to protozoans  and binds long-chain fatty acid (LCFA) CoA esters (C14-C22) through its acyl-CoA-binding domains (ACBDs) with high affinity and specificity . Several studies have demonstrated the precise biochemical functions of ACBP, such as binding and transporting acyl-CoAs, maintaining intracellular acyl-CoA pools , and participating in membrane biosynthesis , fatty acid elongation and sphingolipid synthesis, in eukaryotes . Consistent with its biochemical properties, ACBP deficiency decreases the intracellular LCFA-CoA pool while ACBP overexpression has the opposite effect in yeasts  and mice . In addition, it has been shown that ACBP is required for LCFA esterification into triglycerides and phospholipids (PLs) [19, 20] and for oxidation . Furthermore, functional loss of ACBP decreases the very-long-chain fatty acid (VLCFA) content, suggesting that ACBP regulates LCFA elongation and hence VLCFA levels .
Functional investigations of ACBPs in apicomplexan parasites have been reported. In Cryptosporidium parvum, ACBP is presumed to play roles in lipid metabolism and trafficking in parasitophorous vacuoles (PV) . Recently, our laboratory showed that acyl-CoA binding protein and sterol carrier protein 2 in Toxoplasma cooperate in lipid metabolism . As Neospora structurally and morphologically resembles Toxoplasma, we hypothesized that Neospora may depend on the same fatty acid metabolic pathways in which acyl-CoA binding protein determines the metabolic fate of fatty acids. However, little is known about the role of ACBP in N. caninum. We identified a candidate gene (NCLIV_066640) encoding a Neospora acyl-CoA binding protein via bioinformatics analysis. Then, we investigated the expression, localization and activity of this protein in vitro and in vivo. The combined results of genetic disruption, lipidomic analysis and transcriptome profiling revealed roles of ACBP in lipid metabolism. These data may provide a reference for the further investigation of lipid metabolism in N. caninum.
Parasites and cell culture
Human foreskin fibroblasts (HFF, SCSP-106) cells and Vero cells were obtained from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). The cells were cultured in DMEM (Dulbecco’s modified Eagle’s medium) supplemented with 20% (HFF cells) and 8% (Vero cells) heat-inactivated fetal bovine serum (FBS; Gibco, New York, USA), respectively in a humidified incubator containing 5% CO2 at 37 °C. The N. caninum wild-type strain (Nc-1) and the constructed NcACBP-deficient strain were cultured as tachyzoites by serial passages in HFF or Vero cells supplemented with 2% FBS, 10 units/ml penicillin and 100 mg/ml streptomycin. The parasites were harvested from freshly lysed Vero cells and washed twice with cold phosphate-buffered saline (PBS), and centrifuged at 1800× rpm for 10 min for collection as previously described .
Mice and virulence assay
A virulence assay was performed on six-week-old female BALB/c mice (Peking University Health Science Center, China). The animals were housed under specific-pathogen-free conditions with ad libitum access to feed and water. Groups of BALB/c mice (n = 5) were infected with 5 × 106 freshly harvested tachyzoites of different strains intraperitoneally.
Sequence analysis and recombinant protein expression
To obtain detailed information on the acyl-coenzyme A binding protein in N. caninum, the Toxoplasma genomic resource database (ToxoDB ver.9.0) was used to search for ACBP-related genes. The NcACBP (NCLIV_066640) gene sequence was obtained from the ToxoDB (http://toxodb.org/toxo/) website. Alignment with the corresponding sequence in Saccharomyces cerevisiae and in other apicomplexan parasites was performed by Clustal W using DNAMAN (Lynnon Biosoft, San Ramon, USA) and bootstrap consensus trees were then generated. The physical and chemical properties of ACBPs were obtained using the online tool ExPAS-ProtParam (https://web.expasy.org/protparam/); the conserved domains of ACBPs were analyzed by SMART (http://smart.embl-heidelberg.de/) and mapped using IBS 1.0 software.
DNA extraction kit (Aidlab Biotechnologies Co., Ltd., Beijing, China) was used to extract the genomic DNA of the wild-type strain. The open reading frame was amplified by PCR using the primer pairs F1/R1, the 5′UTR, F2/R2, and the 3′UTR, F3/R3. All primers used in this study are listed in Additional file 1: Table S1. Then, we used the pET-28a vector (Novagen, Madison, Germany) to clone the complete coding sequence and transformed into Escherichia coli (Transetta, TansGenBiotech Co., Ltd., Beijing, China) for recombinant protein expression. The rNcACBP-His and rNcACBP-GST proteins were purified by affinity chromatography using Ni-IDA agarose in accordance with the manufacturer’s standard protocol. Purified recombinant proteins were assessed by SDS-PAGE analysis. The gel was consisted of running gel (12%, w/v, acrylamide) and stacking gel. The electrophoresis buffer was 25 mM Tris, 192 mM glycine, 0.1% SDS, pH 8.3. After electrophoresis, the gels were stained using Coomassie brilliant blue (Thermo Fisher Scientific Inc., Waltham, MA, USA) for protein detection. They were treated in protein destaining solution (10% acetic acid, 30% methyl alcohol, 60% distilled water, v:v:v) overnight. The relative molecular mass of each of the recognized bands was determined by comparison with standard markers.
Quantitative real-time PCR (qRT-PCR)
Total RNA was extracted from 1 × 108 tachyzoites of the wild-type strain and ΔNcACBP strain with TRIzol reagent and converted to cDNA using an EasyScript First-Strand cDNA Synthesis SuperMix kit (TransGen, Beijing, China) in accordance with the manufacturer’s instructions. The NcActin (NcLiv_061190) gene was selected as the endogenous reference gene to normalize the P-value in each sample . As per the manufacturer’s instructions, RT-PCR was performed in triplicate with three independent samples for each experimental group in an ABI Prism 7500 System (Applied Biosystems Inc., Foster, USA) with SYBR Green II (Takara Biotechnology, Dalian, Co., Ltd, China). The RT-PCR conditions were as follows: 94 °C for 5 s, followed by 40 cycles at 94 °C for 5 s and 60 °C for 30 s. The relative expression levels of genes were calculated from the quantification cycle (Cq) value and standardized by the 2−ΔΔCq method . All primers used in this study are listed in Additional file 1: Table S1.
Enzyme kinetics and substrate preference assays
The in vitro binding kinetics and substrate preference of the rNcACBP-GST protein were determined by an NBD-C16:0-CoA-based assay as described previously [23, 28]. The reaction components consisted of 0.25 μM fluorescently labeled NBD-C16:0-CoA substrate, 1 μM rNcACBP-GST or GST protein (as the negative control group), and PBS in a final volume of 100 μl. The enzyme kinetic assays were performed using 1 μM rNcACBP-GST protein and NBD-C16:0-CoA (0–0.3 μM) and PBS in a final volume of 100 μl. The reactions were performed in a 96-well plate and incubated for 5 min at room temperature to ensure maximum binding before proceeding with fluorescence measurements. The fluorescence intensity of the conjugates was measured at 460 nm and 538 nm by a SpectraMax M5.
Heterologous complementation of the yeast ACBP gene by the NcACBP gene
To determine whether the NcACBP gene could complement the yeast ACBP gene, the NcACBP gene was inserted into the yeast expression vector p405ADH1 under the control of the ADH1 promoter and CYC1 terminator. The yeast wild-type and ACBP null mutant (ΔScACBP) strains (Dharmacon Inc., Lafayette, CO, USA) were obtained and grown on YPD (1% yeast extract, 2% peptone, 2% dextrose) medium. 2% agar was added to obtain a solid medium if necessary. To express N. caninum proteins in a yeast strain defective in ScACBP formation, the putative NcACBP sequence was released from the p405ADH vector by restriction digestion with BamHI and XhoI, purified and ligated into the p405ADH-NcACBP plasmid digested with SacI. The linearized p405ADH-NcACBP plasmid was transformed into the ΔScACBP strain using a commercial kit (FunGenome Company, Beijing, China) according to the manufacturer’s instructions. Candidate transformants were picked and streaked on plates lacking leucine (SCGal-Glu-Leu) to select against the presence of the p405ADH-NcACBP plasmid. Analysis of the yeast vacuolar structure was performed using the fluorescent vital dye FM4-64 [29, 30], and visualized by a Leica confocal microscope system (Leica TCS SP52, Wetzlar, Germany).
The subcellular localization of NcACBP and apicoplast were detected by IFA. Tachyzoites that freshly released or infected HFF cells were fixed with 4% paraformaldehyde for 30 min, as previously described . Samples were blocked with 3% BSA-PBS after permeabilized with 0.1% Triton X-100 and incubated with primary antibodies for 1 h. Rabbit anti-NcSRS2 (1:500), mouse anti-HA (1:500), mouse anti-NcENR (1:500) were used as primary antibodies in this study. Then, FITC-conjugated goat-anti mouse IgG (Sigma-Aldrich, Louis, MO, USA) and Cy3-conjugated goat-anti-rabbit IgG (Sigma-Adrich) were used as secondary antibodies at 1:1000 dilution for labeling. The nuclear was stained with Hoechst (1:100) (Sigma-Aldrich), and the lipid bodies were stained with Nile red (1:50). Mouse anti-HA monoclonal antibody was purchased from Sigma-Aldrich. Mouse anti-NcENR and rabbit anti-NcSRS2 were all polyclonal antibodies stored in our laboratory.
Construction of the NcACBP knockout N. caninum strain
We used the homologous recombination strategy and a CRISPR/Cas9 plasmid to construct the NcACBP deletion strain. The pTCR-CD plasmid contains the chloramphenicol resistance gene (CmR), red fluorescence protein gene (RFP), bacterial cytosine deaminase gene (CD), and ampicillin resistance gene (Amp), and was modified as previously described . The expression of the CmR-RFP fusion gene and CD gene were under the control of the NcTublin promoter. The NcACBP 5′ and 3′ untranslated region (UTR) fragments were amplified and inserted into the pTCR-CD plasmid to delete the NcACBP gene. To produce the NcACBP gene deletion plasmid, correctly sequenced plasmids were double digested with HindIII and XhoI for the 5′ UTR and with XmaI and SpeI (NEB, Ipswich, USA) for the 3′ UTR, named pTCR-NcACBP-CD KO. The linearized pTCR-NcACBP-CD KO plasmid, digested with NotI (NEB), was purified using ethanol precipitation and then resuspended by cytomix . The single guide RNA (5′-GCT TAC AAA AGG CTC ATT CA-3′) was designed using the E-CRISPR website to recruit Cas9 and cut the NcACBP gene. In addition, 19-bp sequences, up- and downstream in the CRISPR/Cas9 plasmid were selected as overlapping regions. Correctly sequenced CRISPR/Cas9 plasmids were double digested with KpnI and AvrII to produce the NcACBP gene deletion plasmid, named CRISPR/Cas9-NcACBP. All constructs were verified by sequencing (Beijing Ruibiotech Co. Ltd, Beijing, China). Fifty μg of the above plasmids were electroporated into wild-type tachyzoites (1 × 107). The transgenic parasites were grown under chloramphenicol (20 mM) and 5-fluorocytosine (40 mM) selection pressure and then screened to confirm the purity of the selected strains until cultured consecutively to the 10th generation.
The size and number of plaques represent parasites’ successive rounds of lytic cycles, including invasion, replication, and egress, and thus can be used to evaluate the exhaustive fitness of tachyzoites. HFF cells were seeded in 6-well plates previously, then infected with 800 parasites per well and incubated for 7 days in a 37 °C incubator with 5% CO2. Subsequently, cells were stained with 2% crystal violet for 15 min after fixed in PBS containing 4% paraformaldehyde for 30 min. The stained wells were washed with deionized water, airdried and visualized by microscopy (Olympus Co., Tokyo, Japan) using image acquisition software. Plaque area was measured using photoshop by statistical pixels .
HFF cells were prepared into 24-well plates with coverslips before infected with 800 tachyzoites per well and cultured at 37 °C in DMEM with 20% FBS for 30 min, then washed with PBS for three times. Then cells were washed with PBS and fixed with 4% paraformaldehyde after adhered overnight in a 37 °C incubator with 5% CO2, followed by IFA assay to counted the number of parasites per vacuole with a fluorescence microscope.
Intracellular tachyzoites (1 × 108) were harvested and total lipids were extracted in chloroform/methanol (1:2, v/v) for 30 min at 60 °C, as previously described [33, 34]. Then, polar and nonpolar metabolites were separated by phase partitioning. The organic phase was dried under N2 gas and dissolved in chloroform/methanol (2:1, v/v) for lipid analysis. Then, the lipids were mixed with 1 nM lauric acid (C12:0) as the internal standard and derivatized using MethPrep II (Alltech, Chicago, Illinois, USA). The resulting fatty acid methyl esters were analyzed using GC-MS [33, 35]. The comparison of retention times and mass spectra from GC-MS could identify all the fatty acid methyl esters profile with authentic chemical standards. The data were analyzed with the Agilent workstation software MSD Chemstation D.01.02.16 and compared with the spectral library. The peak areas of fatty acids from N. caninum were analyzed and standardized with respect to the internal standard. Finally, the peak area ratio of total fatty acids from different strains was calculated.
Stable isotope metabolic labeling of N. caninum fatty acids
Stable isotope metabolic labeling experiments followed by lipid extraction and GC-MS analysis were performed as previously described [33,34,35]. Infected HFF cells were cultured in medium in which the unlabeled glucose was replaced with 8 mM U-13C-glucose (Cambridge Isotope Laboratories, Andover, MA, USA) for polar metabolites 24 h before the tachyzoites egress. All lipids were analyzed by GC-MS after derivatization using MethPrep II (Alltech). The shift in the mass spectrum of each fatty acid was analyzed to assess the incorporation of 13C into fatty acids.
Thin layer chromatography
We adopted the Folch method to extract the total lipids as previously described . In short, fatty acids were released as their methyl esters after treated with chloroform/methanol (2:1 v/v) and dried under N2. Each sample was suspended in 50 μl of chloroform, and 10 μl was loaded onto TLC plates (Whatman, Maidstone, Kent, UK). The samples were delivered by small drops, and the esterification mix was loaded as the standard marker to quantify total palmitate. On TLC Silica Gel 60 plates (Merck, Gibbsboro, NJ, USA), total lipids were separated by a solvent for neutral lipids (hexane: diethyl ether: acetic acid (90:10:1, v/v/v)) and run with lipid standards for the separation of total phospholipids from mono-, di-, and triacylglycerols (TAG). Canon digital scanner (model F917500; Tokyo, Japan) was used to image the plates, and the intensities of TAG bands were measured by densitometry.
Lipid body staining
We used fluorescence microscopy to detected the lipid bodies in N. caninum, intravacuolar parasites were fixed in 4% paraformaldehyde and following the IFA protocol. Lipid bodies were stained with Nile red , a fluorescent dye that preferentially binds neutral lipids such as triglycerides to assess the effect of NcACBP on lipid accumulation. Infected cells and extracellular tachyzoites were fixed with 4% paraformaldehyde, and permeabilize with 0.1% Triton X-100 for 15 min, then incubated with the Nile red in the dark at RT for 20 min. Cells were washed in PBS and imaged using a Leica confocal microscope system (TCS SP52; Leica, Wetzlar, Germany).
Transcriptome sequencing was performed according to the manufacturer’s recommendations. Briefly, a total amount of 3 μg RNA per sample was used as input material for the RNA sample preparation and subjected to poly-T oligo-attached magnetic bead enrichment. Sequencing libraries were generated using the NEBNext® UltraTM Directional RNA Library Prep Kit (Illumina, NEB, Ipswich, USA) following the manufacturer’s instructions and index codes were added to attribute sequences to each sample. Library quality was assessed on the Agilent Bioanalyzer 2100 system. The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumina), according to the manufacturer’s recommendations. After cluster generation, the library preparations were sequenced on an Illumina Hiseq platform and paired-end reads were generated. All treatments and subsequent analyses were performed on individual transcripts.
Differential expression analysis
Differential expression analysis was performed using the DESeq R package (1.18.0). DESeq provides statistical routines for determining differential expression in digital gene expression data using a model based on the negative binomial distribution. Raw counts were normalized using the Benjamini and Hochberg’s approach for controlling the false discovery rate (FDR) . A corrected P-value of 0.005 and log2(Fold change) of 1.2 were set as the threshold for significantly differential expression [39, 40].
GO and KEGG enrichment analysis of differentially expressed genes
Gene Ontology (GO) enrichment analysis of differentially expressed genes was implemented by the GOseq package in R software , in which gene length bias was corrected. GO terms with a corrected P-value less than 0.05 were considered significantly functional enrichment in differential expressed genes.
KEGG is a database resource for understanding high-level functions and utilities of the biological system, such as the cell, the organism and the ecosystem, from molecular-level information, especially large-scale molecular datasets generated by genome sequencing and other high-throughput experimental technologies (http://www.genome.jp/kegg/) . We used KOBAS software to test the statistical enrichment of differential expression genes in KEGG pathways.
Graphs and statistical analyses were made using GraphPad Prism (GraphPad, San Diego, CA, USA). All data were analyzed using Student’s t-test and univariate analysis of survival using Log-rank (Mantel-Cox) test. P-values are represented in the figures as follows: *P < 0.05; **P < 0.01; ***P < 0.001; ns, not significant.
NcACBP is a conserved protein in apicomplexan parasites
Only one protein containing the acyl-coenzyme A binding domain was found on chromosome XII in N. caninum, also named the diazepam binding inhibitor based on the annotations. Multiple alignments performed by Clustal W showed there were many highly conserved amino acid residues between S. cerevisiae and other apicomplexan species which were reported to play important roles in binding acyl-CoA ligand . The conserved Lys residues in H2 and H3 α-helices clusters and Tyr residue in H2 are presumed to interact with 3′-phosphate group of the CoA part in the acyl-CoA ligand, while the conserved Tyr in H4 may contribute to the stacking of its own aromatic ring and the adenine ring of the CoA part of ligand. The triangles indicate the potential binding sites for acyl-CoA esters, and H1-H4 indicates the positions of four putative alpha-helices (Fig. 1a). TgACBP1 shares the highest amino acid sequence identity (85.26%) with NcACBP by phylogenetic analysis (Fig. 1b). NcACBP consists of 95 amino acids and its predicted molecular weight is ~ 10.7 kDa. The 3D structure of NcACBP showed four α-helix bundles that constitute the binding pocket for acyl-CoA esters (Fig. 1c). In the large molecular weight ACBPs (> 12 kDa), the conserved acyl-CoA-binding domain (ACBD) lies at the N-terminus while other structural domains, such as the ankyrin repeats (ANK) domain, reside at the C-terminus, and some of these proteins contain a transmembrane region (TMR) at the N-terminus (i.e. EtACBP2). In contrast, the small molecular weight ACBPs (10-kDa) contain only one conserved acyl-CoA-binding domain (Fig. 1d). These data indicate that NcACBP is a conserved protein in apicomplexan parasites.
Identification and cellular localization of NcACBP
To assess the expression and localization of NcACBP, recombinant proteins fused with a GST tag and/or histidine tag were expressed in E. coli successfully. The recombinant protein rNcACBP-His was identified by SDS-PAGE (Fig. 2a-i). The native NcACBP expression in N. caninum was identified by western blot. The expected bands were elicited by anti-NcACBP polyclonal antibody, recognized an ~ 11 kDa protein in the lysate of tachyzoites and NcActin served as the loading control (Fig. 2a-ii). To localize the NcACBP, the C-terminal of NcACBP was fused with a triple hemagglutinin (3× HA) epitope tag by single homologous recombination (Fig. 2b). The IFA showed that NcACBP was distributed in the cytosol of intracellular and extracellular parasites (Fig. 2d), suggesting that NcACBP is a cytoplasmic protein in N. caninum.
Determination of NcACBP binding activity in vitro and in vivo
To verify the binding activities of NcACBP in vitro, we expressed recombinant NcACBP-GST in E. coli and performed fluorescent substrate binding assay using NBD-C16:0-CoA. We observed an increased fluorescence upon the binding of NBD-C16:0-CoA to rNcACBP, while there is no signal in the GST control group (Fig. 3a). Through this fluorometric assay, we determined the dissociation constant (KD) of NcACBP 11.65 nM for NBD-C16:0-CoA (Fig. 3b). These results indicate that rNcACBP has an acyl-CoA binding activity in vitro. Additionally, the function of NcACBP in vivo was analyzed. The NcACBP gene was inserted into the yeast expression vector p405ADH1 under the control of the ADH1 promoter and CYC1 terminator and transforming into ΔScACBP mutant yeast (Fig. 3c–d). Transformed yeast cells were stained with the vacuolar marker FM4-64 and observed with a fluorescence microscope. Lack of ScACBP caused disintegration of yeast vacuoles, leading to an increased number of cells with the multilobed vacuole phenotype. Cells transformed with an empty vector did not alter the vacuolar disintegration phenotype; however, yeast cells expressed NcACBP protein complemented the phenotype (Fig. 3e–f). The significance was determined by Student’s t-test: t(5) = 7.481, P < 0.001. These results demonstrated that NcACBP can bind acyl-CoA in vitro and complement the function of yeast ACBP in vivo.
NcACBP is not essential during the tachyzoite stage
To characterize the role of NcACBP in N. caninum, we generated a complete knockout mutant of NcACBP (ΔNcACBP) via the CRISPR/Cas9 system. By targeting the native NcACBP locus in the wild-type strain, the gene was replaced by the CAT-RFP cassette through homologous recombination (Fig. 4a). The deletion of ACBP gene was confirmed by PCR (Fig. 4b). The RT-PCR results also showed that the transcriptional level of NcACBP was significantly reduced after deletion (Fig. 4c). The phenotype assays showed that knockout of NcACBP did not affect the plaque number and size (Fig. 4d) and did not affect the replication of parasites either (Fig. 4e). These data showed that NcACBP is not an essential gene in N. caninum during the tachyzoite stage.
Knockout of NcACBP enhanced the pathogenicity of N. caninum in mice
To evaluate the effects of NcACBP on N. caninum pathogenicity, BALB/c mice were intraperitoneally infected with 5 × 106 freshly released tachyzoites of the wild-type or ΔNcACBP. Signs of illness, such as coat ruffling, inactivity, and mental depression, were observed after 4 days post-infection. The mice infected with the ΔNcACBP strain began to die at 8 days post-infection, while mice infected with the wild-type strain showed a significant delay in the time of death (dpi = 18). The survival rate of ΔNcACBP infected mice was lower than the wild-type infected mice significantly. Significance was determined by Log-rank (Mantel-Cox) test: χ2 = 9.67, df = 15, P = 0.0146 (Fig. 4f). These results demonstrate that the pathogenicity of N. caninum in mice was enhanced after the deletion of NcACBP.
NcACBP disruption does not impair the biogenesis of the apicoplast
Since ACBPs were indicated to participate in acyl-CoA binding and transport, the maintenance of intracellular acyl-CoA pools and membrane biosynthesis in eukaryotes, it is possible that loss of NcACBP may affect the synthesis of fatty acids. To determine the potential role of NcACBP in fatty acid synthesis, especially de novo fatty acid synthesis in the apicoplast, the apicoplast was stained with mouse anti-NcENR antibodies, which showed that disruption of NcACBP did not affect the morphology of the apicoplast (Fig. 4g), demonstrating that loss of NcACBP does not impair the biogenesis of the apicoplast during the tachyzoite stage.
NcACBP disruption does not affect the total fatty acid composition
To validate our hypothesis that loss of the NcACBP may alter the lipid profile in parasites, the fatty acids were extracted from both the wild-type and ΔNcACBP strains and its species and abundance were analyzed by GC-MS. The results showed that there was no obvious difference in the total fatty acid content between the parental and ΔNcACBP strains (Fig. 5a), indicating that disruption of NcACBP did not significantly affect the abundance of fatty acids. To investigate the effect of disruption of NcACBP on the fatty acid synthesis, we labeled parasites with U-13C-glucose as previously reported . The GC-MS results showed that this disruption did not cause any changes in the incorporation of 13C into fatty acids (Fig. 5b). These data demonstrate that NcACBP is dispensable for fatty acid production.
NcACBP disruption leads to the accumulation of neutral lipids in parasites
To verify whether disruption of NcACBP affects lipid accumulation, we extracted lipids from wild-type and ΔNcACBP parasites and performed TLC analysis, which showed an increased abundance of TAG in ΔNcACBP parasites compared to that of parental parasites (Fig. 5c). Then, we focused on the delivery of neutral lipids stored in host lipid bodies to the PV and the parasite, as several pathogens are able to recruit lipid bodies around their vacuoles [43, 44]. The infected cells were incubated with Nile red, a fluorescent dye that preferentially binds neutral lipids such as triglycerides. Compared with cells infected with the wild-type strain, numbers and size of the lipid bodies were increased in HFF cells infected with ΔNcACBP strains; however, there was no statistically significant difference (Fig. 5d).
Loss of NcACBP results in global changes in the expression of multiple genes
Transcriptome analysis was performed to compare the gene expression levels in wild-type and ΔNcACBP, so as to investigate the mechanism of the phenotypic changes in ΔNcACBP. Statistical analysis showed that 1474 genes were expressed differentially, including 623 upregulated genes and 851 downregulated genes with a log2 fold change of ≥ 1.2 (Fig. 6a) . The regulated genes with a log2 fold change of ≥ 4 between Nc-1 and ΔNcACBP strains are provided in Additional file 1: Table S2. To characterize 1474 differentially expressed genes, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was executed. The most highly enriched canonical pathways were the ubiquitin-mediated proteolysis, aminoacyl-tRNA biosynthesis, glycolysis and gluconeogenesis, proteasome, ribosome and regulation of autophagy pathways (Fig. 6b). We also found several differentially expressed genes involved in multiple metabolic pathways, including oxidative phosphorylation, peroxisome, protein processing in the endoplasmic reticulum and glycerophospholipid metabolism pathways. Several genes involved in fatty acid biosynthesis and metabolism were upregulated, while genes involved in ER protein processing and glycerophospholipid metabolism were downregulated (Fig. 6c). These results suggest that loss of NcACBP affects the metabolism of fatty acids and glycerol phospholipids in N. caninum and results in global changes in the expression of multiple genes involved in various biological functions and cell components. Gene Ontology enrichment analysis was executed to identify the functions of the differentially expressed genes. These genes were classified into three categories: biological process, cellular component and molecular function. We observed that 1344 genes were involved in biological processes, such as glycoprotein biosynthetic and metabolic processes, cellular lipid biosynthesis and metabolic processes, phosphorylation and glycosylation; 621 genes participated in molecular function (i.e. transferase activity, transfer pentosyl groups, kinase activity and transporter activity); and 331 genes took part in cellular components, particularly in the components of extracellular region, membrane region and cytoplasm (Fig. 6d). These data demonstrate that the deletion of NcACBP damages the profile of the cell membrane, the metabolism of glycoproteins and lipids, phosphorylation and glycosylation.
ACBP, originally identified as a mammalian diazepam-binding inhibitor , mainly functions as an intracellular acyl-CoA transporter and pool former  and is critical to lipid metabolism in cells . Among apicomplexans, NcACBP functions as a short protein with a unique acyl-CoA binding domain and is dispersed in the cytosol of either intracellular or extracellular parasites. However, other apicomplexans may have multiple ACBP proteins of various types (i.e. short proteins or long proteins fused with ankyrin repeats or with containing transmembrane region domains), revealing that ACBP-mediated metabolic pathways may be highly divergent in the phylum Apicomplexa. Our biochemical data indicate that rNcACBP is capable of binding to palmitoyl-CoA (25 μM) with the highest binding affinity (KD = 11.65 nM) in vitro. In addition, the abnormal phenotype of slow growth and an increased number of multilobed vacuoles in ΔScACBP mutant yeast  can be rescued by NcACBP, similar to the effect of RpACBP-1 , demonstrating that NcACBP possibly plays a similar role to that of ScACBP as a homologous protein.
Although the deletion of NcACBP by CRISPR/Cas9 using homologous recombination did not perturb the survival and morphology of the parasites in vitro, its pathogenicity in mice was enhanced. We have compared the transcriptomes between ΔNcACBP and wild-type parasites and tried to reveal overexpressed genes that may result in enhanced virulence of ΔNcACBP. Among 8 genes with a 6-fold higher increase in ΔNcACBP parasites, we observed that NCLIV_033380, NCLIV_069820, NCLIV_019580, and NCLIV_020100 are presumed to be involved in cellular membrane biogenesis, while there are still several genes with unknown protein domains and functions such as NCLIV_052460, NCLIV_020100, NCLIV_023980, and NCLIV_002870. These four proteins have no signal peptides, indicative of the low possibility of secretion into PV and host cells to regulate host-parasite interaction. In addition, we also analyzed 25 genes with a 4–6-fold transcriptional increase and also found that most of these genes may be involved in membrane biogenesis. However, transcriptome sequencing did not identify any genes related to known Neospora virulence factors, such as ROP5, ROP16 and GRA17.
Reportedly, lipids as an important pathogenic factor can help parasites escape the host immune response and cause diseases [5, 49]. In our study, the TAG abundance was increased in the ΔNcACBP strains compared with that in the wild-type strain. Moreover, an increased amount of neutral lipids derived from host lipid bodies was observed surrounding the ΔNcACBP strains compared with that in the wild-type strain. Although these lipid bodies were mainly detected under experimental conditions, they may be relevant physiologically and involved in nutrient uptake and delivery to the parasite’s cytoplasm potentially. Lipid body formation is closely linked to the biosynthesis of neutral lipids, such as sterols and TAGs that are packaged in the lipid droplet core [50, 51]. TAG is formed by redundant non-esterified free fatty acids, which are cytotoxic at high concentration, in a self-protective manner . In addition, our laboratory has previously shown that disruption of ACBP alone in Toxoplasma did not affect the growth ability and intracellular replication, which is the same as NcACBP knockout phenotypes in vitro . However, TgACBP disruption did not affect the virulence to mice while NcACBP knockout enhanced the pathogenicity to mice. Notably, TgACBP and TgSCP2 double disruption reduced the TAG overall abundance while NcACBP knockout resulted in the accumulation of TAG. Therefore, we hypothesize that one possibility of the virulence difference lies in the metabolic regulation mechanism of TAG between Toxoplasma and Neospora. TAG could provide acyl skeleton for phospholipids and other many kinds of lipid derivatives, which can be involved in membrane synthesis. Transcriptional levels of many genes related to membrane biogenesis have been identified, indicating that NcACBP disruption may lead to membrane biogenesis disorder and thus affect the immunological recognition of specific antigens in parasite membranes by the host. In this perspective, we should apply lipidomics to reveal lipid components in cellular membranes of ΔNcACBP parasites and pay more attention to the immunological response of host cells to ΔNcACBP infection in future work.
Stable isotope labeling and metabolomic analysis showed that loss of NcACBP did not significantly alter the total abundance of fatty acids or affect the synthesis of fatty acids, nor did it impair the synthesis of unsaturated long-chain fatty acids or the ability to salvage short-chain fatty acids from the host. Furthermore, because U-13C-glucose labeled fatty acids were synthesized by the FASII pathway, there may be alternate substrates for fatty acid synthesis, and the sources of fatty acids are likely complicated and diverse. As intracellular pathogens acquire essential non-diffusible host metabolites , perhaps the fatty acid uptake of NcACBP-deficient parasites is increased in host cells. However, we hypothesize that proteins functionally compensatory to NcACBP are expressed to compensate for fatty acid metabolism in NcACBP-deficient parasites. The results of the present studies strongly suggest that in addition to the essential fatty acid synthesis pathways (FASI, FASII, and FAE), intracellular tachyzoites can salvage fatty acids from the host and the surrounding environment in a fine-tuned manner to meet their own needs [4, 35, 54]. Since there are few studies on fatty acid biosynthesis and metabolism in N. caninum, we propose that the fatty acid biosynthesis and metabolic pathways in N. caninum may be similar to those in T. gondii and P. falciparum, based on our results and previous studies [35, 54]. Acetyl-CoA is a key metabolite involved in the TCA cycle, fatty acid synthesis, fatty acid elongation, etc. The apicoplast FASII pathway generates its own pool of acetyl-CoA, mitochondria generate acetyl-CoA for the FASII pathway in the cytosol, and the elongation pathway in the ER relies on the TCA cycle. In addition, fatty acids and lipids scavenged from the host may be a possible redundant source for parasites.
We demonstrated that NcACBP is not an essential gene during the tachyzoites stage and is present as a short ACBP dispersed in the cytoplasm. Genetic ablation of NcACBP did not impair the survival and growth phenotype of N. caninum but enhanced its pathogenicity in mice. This deletion did not affect the overall fatty acid composition but did modify the abundance of TAG and loss of NcACBP resulted in global changes in the expression of multiple genes.
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional file.
clustered regularly interspaced short palindromic repeats/CRISPR associated protein 9
fatty acid synthesis
fatty acid elongation
fatty acid binding protein
sterol carrier protein 2
long-chain fatty acid
very long-chain fatty acid
fetal bovine serum
bovine serum albumin
open reading frame
green fluorescence protein
gas chromatography-mass spectrometry
thin layer chromatography
false discovery rate
Kyoto Encyclopedia of Genes and Genomes
tricarboxylic acid cycle
nuclear localization signal
Almeria S, Serrano-Perez B, Lopez-Gatius F. Immune response in bovine neosporosis: protection or contribution to the pathogenesis of abortion. Microb Pathog. 2017;109:177–82.
Dubey JP, Carpenter JL, Speer CA, Topper MJ, Uggla A. Newly recognized fatal protozoan disease of dogs. J Am Vet Med Assoc. 1988;192:1269–85.
Calarco L, Barratt J, Ellis J. Genome-wide identification of mutational hotspots in the apicomplexan parasite Neospora caninum and the implications for virulence. Genome Biol Evol. 2018;10:2417–31.
Mancio-Silva L, Slavic K, Grilo Ruivo MT, Grosso AR, Modrzynska KK, Vera IM, et al. Nutrient sensing modulates malaria parasite virulence. Nature. 2017;547:213–6.
Ramakrishnan S, Serricchio M, Striepen B, Butikofer P. Lipid synthesis in protozoan parasites: a comparison between kinetoplastids and apicomplexans. Prog Lipid Res. 2013;52:488–512.
Fujimoto T, Parton RG. Not just fat: the structure and function of the lipid droplet. Cold Spring Harb Perspect Biol. 2011;3:a004838.
Ehrt S, Schnappinger D. Mycobacterium tuberculosis virulence: lipids inside and out. Nat Med. 2007;13:284–5.
McFadden GI, Reith ME, Munholland J, Lang-Unnasch N. Plastid in human parasites. Nature. 1996;381:482.
Ramakrishnan S, Docampo MD, MacRae JI, Ralton JE, Rupasinghe T, McConville MJ, et al. The intracellular parasite Toxoplasma gondii depends on the synthesis of long-chain and very long-chain unsaturated fatty acids not supplied by the host cell. Mol Microbiol. 2015;97:64–76.
Amiar S, MacRae JI, Callahan DL, Dubois D, van Dooren GG, Shears MJ, et al. Apicoplast-localized lysophosphatidic acid precursor assembly is required for bulk phospholipid synthesis in Toxoplasma gondii and relies on an algal/plant-like glycerol 3-phosphate acyltransferase. PLoS Pathog. 2016;12:e1005765.
Neess D, Bek S, Engelsby H, Gallego SF, Faergeman NJ. Long-chain acyl-CoA esters in metabolism and signaling: role of acyl-CoA binding proteins. Prog Lipid Res. 2015;59:1–25.
Faergeman NJ, Wadum M, Feddersen S, Burton M, Kragelund BB, Knudsen J. Acyl-CoA binding proteins; structural and functional conservation over 2000 MYA. Mol Cell Biochem. 2007;299:55–65.
Burton M, Rose TM, Faergeman NJ, Knudsen J. Evolution of the acyl-CoA binding protein (ACBP). Biochem J. 2005;392:299–307.
Xiao S, Chye ML. An Arabidopsis family of six acyl-CoA-binding proteins has three cytosolic members. Plant Physiol Biochem. 2009;47:479–84.
Faergeman NJ, Feddersen S, Christiansen JK, Larsen MK, Schneiter R, Ungermann C, et al. Acyl-CoA-binding protein, Acb1p, is required for normal vacuole function and ceramide synthesis in Saccharomyces cerevisiae. Biochem J. 2004;380:907–18.
Bloksgaard M, Neess D, Faergeman NJ, Mandrup S. Acyl-CoA binding protein and epidermal barrier function. Biochim Biophys Acta. 2014;1841:369–76.
Knudsen J, Faergeman NJ, Skott H, Hummel R, Borsting C, Rose TM, et al. Yeast acyl-CoA-binding protein: acyl-CoA-binding affinity and effect on intracellular acyl-CoA pool size. Biochem J. 1994;302:479–85.
Huang H, Atshaves BP, Frolov A, Kier AB, Schroeder F. Acyl-coenzyme A binding protein expression alters liver fatty acyl-coenzyme A metabolism. Biochemistry. 2005;44:10282–97.
Gaigg B, Neergaard TB, Schneiter R, Hansen JK, Faergeman NJ, Jensen NA, et al. Depletion of acyl-coenzyme A-binding protein affects sphingolipid synthesis and causes vesicle accumulation and membrane defects in Saccharomyces cerevisiae. Mol Biol Cell. 2001;12:1147–60.
Elle IC, Simonsen KT, Olsen LC, Birck PK, Ehmsen S, Tuck S, et al. Tissue- and paralogue-specific functions of acyl-CoA-binding proteins in lipid metabolism in Caenorhabditis elegans. Biochem J. 2011;437:231–41.
Harris FT, Rahman SM, Hassanein M, Qian J, Hoeksema MD, Chen H, et al. Acyl-coenzyme A-binding protein regulates Beta-oxidation required for growth and survival of non-small cell lung cancer. Cancer Prev Res. 2014;7:748–57.
Neess D, Bek S, Bloksgaard M, Marcher AB, Faergeman NJ, Mandrup S. Delayed hepatic adaptation to weaning in ACBP−/− mice is caused by disruption of the epidermal barrier. Cell Rep. 2013;5:1403–12.
Zeng B, Cai X, Zhu G. Functional characterization of a fatty acyl-CoA-binding protein (ACBP) from the apicomplexan Cryptosporidium parvum. Microbiology. 2006;152:2355–63.
Fu Y, Cui X, Liu J, Zhang X, Zhang H, Yang C, et al. Synergistic roles of acyl-CoA binding protein (ACBP1) and sterol carrier protein 2 (SCP2) in Toxoplasma lipid metabolism. Cell Microbiol. 2019;21:e12970.
Li WS, Liu J, Wang J, Fu Y, Nan HZ, Liu Q. Identification and characterization of a microneme protein (NcMIC6) in Neospora caninum. Parasitol Res. 2015;114:2893–902.
Ma L, Liu G, Liu J, Li M, Zhang H, Tang D, et al. Neospora caninum ROP16 play an important role in the pathogenicity by phosphorylating host cell STAT3. Vet Parasitol. 2017;243:135–47.
Azadi Y, Ahmadpour E, Hamishehkar H, Daryani A, Spotin A, Mahami-Oskouei M, et al. Quantification of Toxoplasma gondii in the tissues of BALB/c mice after immunization with nanoliposomal excretory-secretory antigens using real-time PCR. Comp Immunol Microbiol Infect Dis. 2018;59:52–6.
Fritzler JM, Zhu G. Novel anti-Cryptosporidium activity of known drugs identified by high-throughput screening against parasite fatty acyl-CoA binding protein (ACBP). J Antimicrob Chemother. 2012;67:609–17.
Zal T, Zal MA, Lotz C, Goergen CJ, Gascoigne NR. Spectral shift of fluorescent dye FM4-64 reveals distinct microenvironment of nuclear envelope in living cells. Traffic. 2006;7:1607–13.
Fischer-Parton S, Parton RM, Hickey PC, Dijksterhuis J, Atkinson HA, Read ND. Confocal microscopy of FM4-64 as a tool for analysing endocytosis and vesicle trafficking in living fungal hyphae. J Microsc. 2000;198:246–59.
Li M, Wang H, Liu J, Hao P, Ma L, Liu Q. The apoptotic role of metacaspase in Toxoplasma gondii. Front Microbiol. 2015;6:1560.
Wang H, Lei T, Liu J, Li M, Nan H, Liu Q. A nuclear factor of high mobility group box protein in Toxoplasma gondii. PLoS ONE. 2014;9:e111993.
MacRae JI, Sheiner L, Nahid A, Tonkin C, Striepen B, McConville MJ. Mitochondrial metabolism of glucose and glutamine is required for intracellular growth of Toxoplasma gondii. Cell Host Microbe. 2012;12:682–92.
Ramakrishnan S, Docampo MD, Macrae JI, Pujol FM, Brooks CF, van Dooren GG, et al. Apicoplast and endoplasmic reticulum cooperate in fatty acid biosynthesis in apicomplexan parasite Toxoplasma gondii. J Biol Chem. 2012;287:4957–71.
Dubois D, Fernandes S, Amiar S, Dass S, Katris NJ, Botte CY, et al. Toxoplasma gondii acetyl-CoA synthetase is involved in fatty acid elongation (of long fatty acid chains) during tachyzoite life stages. J Lipid Res. 2018;59:994–1004.
Taib B, Bouyakdan K, Hryhorczuk C, Rodaros D, Fulton S, Alquier T. Glucose regulates hypothalamic long-chain fatty acid metabolism via AMP-activated kinase (AMPK) in neurons and astrocytes. J Biol Chem. 2013;288:37216–29.
Greenspan P, Mayer EP, Fowler SD. Nile red: a selective fluorescent stain for intracellular lipid droplets. J Cell Biol. 1985;100:965–73.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.
Singhania A, Graham CM, Gabrysova L, Moreira-Teixeira L, Stavropoulos E, Pitt JM, et al. Transcriptional profiling unveils type I and II interferon networks in blood and tissues across diseases. Nat Commun. 2019;10:2887.
Kobayashi K, Umeda K, Ihara F, Tanaka S, Yamagishi J, Suzuki Y, et al. Transcriptome analysis of the effect of C-C chemokine receptor 5 deficiency on cell response to Toxoplasma gondii in brain cells. BMC Genom. 2019;20:705.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11:R14.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–4.
Toledo DA, D’Avila H, Melo RC. Host lipid bodies as platforms for intracellular survival of protozoan parasites. Front Immunol. 2016;7:174.
Saka HA, Thompson JW, Chen YS, Dubois LG, Haas JT, Moseley A, et al. Chlamydia trachomatis infection leads to defined alterations to the lipid droplet proteome in epithelial cells. PLoS ONE. 2015;10:e0124630.
Zheng M, Li D, Zhao Z, Shytikov D, Xu Q, Jin X, et al. Protein phosphatase 2A has an essential role in promoting thymocyte survival during selection. Proc Natl Acad Sci USA. 2019;116:12422–7.
Guidotti A, Forchetti CM, Corda MG, Konkel D, Bennett CD, Costa E. Isolation, characterization, and purification to homogeneity of an endogenous polypeptide with agonistic action on benzodiazepine receptors. Proc Natl Acad Sci USA. 1983;80:3531–5.
Knudsen J, Neergaard TB, Gaigg B, Jensen MV, Hansen JK. Role of acyl-CoA binding protein in acyl-CoA metabolism and acyl-CoA-mediated cell signaling. J Nutr. 2000;130:294s–8s.
Majerowicz D, Hannibal-Bach HK, Castro RSC, Bozaquel-Morais BL, Alves-Bezerra M, Grillo LAM, et al. The ACBP gene family in Rhodnius prolixus: expression, characterization and function of RpACBP-1. Insect Biochem Mol Biol. 2016;72:41–52.
Melo RC, Dvorak AM. Lipid body-phagosome interaction in macrophages during infectious diseases: host defense or pathogen survival strategy? PLoS Pathog. 2012;8:e1002729.
Brown DA. Lipid droplets: proteins floating on a pool of fat. Curr Biol. 2001;11:R446–9.
Guo Y, Cordes KR, Farese RV Jr, Walther TC. Lipid droplets at a glance. J Cell Sci. 2009;122:749–52.
Bartz R, Li WH, Venables B, Zehmer JK, Roth MR, Welti R, et al. Lipidomics reveals that adiposomes store ether lipids and mediate phospholipid traffic. J Lipid Res. 2007;48:837–47.
Pernas L, Bean C, Boothroyd JC, Scorrano L. Mitochondria restrict growth of the intracellular parasite Toxoplasma gondii by limiting its uptake of fatty acids. Cell Metab. 2018;27:886–97.
Botte CY, Yamaryo-Botte Y, Rupasinghe TW, Mullin KA, MacRae JI, Spurck TP, et al. Atypical lipid composition in the purified relict plastid (apicoplast) of malaria parasites. Proc Natl Acad Sci USA. 2013;110:7506–11.
We thank Sai Fan and Xia Cui (Beijing Research Center for Preventive Medicine, China) for support in GC-MS analysis. We are grateful to Professor Bang Shen (Huazhong Agricultural University) for kindly providing the CRISPR/CAS9-UPRT vector.
This study was supported by the National Key Research and Development Program of China (2017YFD0501304), the National Natural Science Foundation of China (31772730, 31972700), the National Key Basic Research Program (973 program) of China (Grant 2015CB150300), and the Beijing Municipal Natural Science Foundation (6172023).
Ethics approval and consent to participate
The experiments were performed in strict according to the recommendations of the Guide for the Care and Use of Laboratory Animals of the Ministry of Science and Technology of China. All experimental animal procedures and protocols were approved by the Institutional Animal Care and Use Committee of China Agricultural University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhou, B., Fu, Y., Zhang, H. et al. Functional characterization of acyl-CoA binding protein in Neospora caninum. Parasites Vectors 13, 85 (2020). https://doi.org/10.1186/s13071-020-3967-9