Identification and molecular characterization of a novel Babesia orientalis thrombospondin-related anonymous protein (BoTRAP1)

The thrombospondin-related anonymous protein (TRAP) family, a kind of transmembrane protein, is widely distributed with a conserved feature of structure in all apicomplexan parasites and plays a crucial role in the gliding motility and survival of parasites. The Babesia orientalis TRAP1 gene (BoTRAP1) was truncated and cloned into a pET-42b expression vector and expressed as a GST-tag fusion protein with a TEV protease site. Rabbit anti-rBoTRAP1 antibody was produced and purified using a protein A chromatography column. Western blot analysis was performed to identify the native protein of BoTRAP1 and differentiate B. orientalis-infected positive from negative serum samples. The localization of BoTRAP1 on merozoites was identified by the indirect florescent antibody test (IFAT). The partial sequence of the TRAP1 gene was cloned from B. orientalis cDNA and identified to contain a von Willebrand factor A (vWFA) region and a thrombospondin type-1 (TSP-1) domain; it had a length of 762 bp, encoding a polypeptide of 254 amino acid residues with a predicted size of 28.2 kDa. The partial sequence was cloned into a pET-42b expression vector and expressed in E. coli as a GST fusion protein. Western blot indicated that rBoTRAP1 has a high immunogenicity and can differentiate B. orientalis-infected positive and negative serum samples collected from water buffaloes. IFAT showed that BoTRAP1 is mainly localized on the apical end of intracellular parasites by using polyclonal antibodies (PcAb) against rBoTRAP1. Meanwhile, the PcAb test also identified the native BoTRAP1 as a ~65 kDa band from B. orientalis lysates. The predicted 3D structure of BoTRAP1 contains a metalion-dependent adhesion site (MIDAS), which could be important for interaction with ligand on the surface of the host cells. Like all known protozoa, B. orientalis has a TRAP family, comprising TRAP1, TRAP2, TRAP3 and TRAP4. The newly identified and characterized BoTRAP1 may play a key role in the invasion of B. orientalis into water buffalo erythrocytes.


Introduction
Babesia orientalis is an apicomplexan parasite that is widespread in southern China and causes babesiosis in water buffaloes, leading to an enormous economic loss [1,2]. The clinical symptoms in water buffalo include anemia, fever, icterus, hemoglobinuria and even death [2,3]. Currently, no vaccine is available to control B. orientalis infection, and drugs for treating B. orientalis are also scarce, suggesting the importance and necessity to explore potential vaccines based on related antigen molecules.
All the thrombospondin-related anonymous protein (TRAP) family members are secreted by micronemes as a membrane protein, and TRAPs with conserved structures are present in all protozoans, with one or more von Willebrand factor A (vWFA) and thrombospondin type-1 repeat (TSR) domain in their extracellular region, as well as a cytoplasmic tail domain (CTD) with a tryptophan residue [4]. In malaria parasites, the TRAPs were first identified in Plasmodium falciparum, and their homologues were found in all other Plasmodium species [5,6]. Subsequent studies have shown that the TRAPs are expressed in different plasmodial stages, such as sporozoite, merozoite and ookinete, and their orthologues are also present in other protozoa, including Toxoplasma gondii, Neospora caninum, Babesia spp., Cryptosporidium spp. and Eimeria spp. [7,8]. In B. bovis and B. gibsoni, the TRAPs were only expressed during the asexual stage, and the antibodies of their recombinant proteins have an obvious influence on Babesia invasion into the host red blood cells (RBCs) [9,10].
In the life-cycle of apicomplexan parasites, host cell invasion is a crucial step for survival, and the process is highly dependent on the interaction between the parasite-and host-surface molecules [11]. In Plasmodium spp., the first step in the invasion of the extracellular merozoites is the attachment to the host cells. In this process, the initial adhesion with host cells based on glycosyl phosphatidylinositol anchor (GPI) of merozoite surface proteins (MSPs) is invertible, followed by reorientation to link the anterior tip of merozoites with the plasma membrane of host cells, leading to the formation of tight junctions from higher-affinity transmembrane proteins secreted by micronemes and rhoptries of parasites; this attachment to the surface of host cells is irreversible. Finally, the parasites invade host cells via a moving complex that involves both apical membrane antigen 1 (AMA1) and rhoptry neck proteins (RONs); this motor process is driven by an actomyosin motor [12]. During the invasion, TRAPs play an important role in the formation of actomyosin motor by linking to actin through their cytoplasmic tail domains (CTD) while binding to host cells via their vWFA domains [7,13]. Subsequent studies have demonstrated that the interaction between TRAP CTD and actin-myosin is connected by aldolase and depends on the sub-terminal tryptophan residue of cytoplasmic tail [14].
Currently, vaccine development efforts have shifted toward the use of antigenically defined immunogens, particularly the molecules interacting or disrupting the process of parasite invasion into host RBCs [10,[15][16][17]. Therefore, identification and characterization of these genes encoding TRAPs in Babesia spp. would facilitate the discovery of novel Babesia vaccine candidate antigens.

Parasites
Babesia orientalis (Wuhan strain) was isolated from Wuhan city, Hubei Province, China, and preserved in liquid nitrogen with the additive of dimethyl sulfoxide (DMSO) in the State Key Laboratory of Agricultural Microbiology, Huazhong Agricultural University, China.
Two water buffaloes (1.5 years-old) were purchased from a Babesia-free area, and were confirmed to be free of Babesia and Theileria by microscope examination and real-time PCR [18]. The water buffalos were splenectomized two weeks before injection of 4 ml of B. orientalis infected blood with the percentage of parasitized erythrocytes (PPE) being 1%. Blood samples were collected every day to monitor the parasitemia until reaching 3%.

Preparation of RNA and cDNA
Blood from the jugular vein of experimentally infected water buffaloes was collected in 50 ml centrifuge tubes containing EDTA-K 2 (Sigma, Shanghai, China). Total RNA was extracted from purified B. orientalis merozoites by using the TransZol Up (TransGen Biotech, Beijing, China) and treated with RNase-free DNaseI (TAKARA, Dalian, China). The cDNA was prepared from 1 μg of the total RNA using PrimeScript TM RT reagent Kit with gDNA eraser (TAKARA, Dalian, China) according to the manufacturer's instructions.

Preparation of recombinant plasmid
Primer pairs for the full-length and partial BoTRAP1 sequences including a vWFA region and a TSP-1 domain were designed based on the fragment of BoTRAP1 screened from B. orientalis genome database ( Table 1). The PCR reaction was performed using the following cycling parameters: 94°C for 5 min, followed by 35 cycles (94°C for 30 s, 58°C for 30 s, 68°C for 1 min), and a final extension at 68°C for 10 min. The resulting PCR product was purified by using EasyPure® PCR Purification Kit (TransGEN, Beijing, China), and then cloned into a pET-42b expression vector. All the constructs were confirmed by DNA sequencing.

Protein expression and purification
The recombinant plasmid was transformed into E. coli BL21 expression host cells, followed by incubation at 37°C in LB medium containing 100 mg/ml kanamycin for 3 h. At an optical density of 0.6 to 0.8 at 600 nm (OD 600 ), the cells were induced with 0.8 mM IPTG (Biosharp, Anhui, China), further cultured at 25°C for another 12 h and then harvested.
For protein purification, cells were harvested by centrifugation at 7000 rpm for 10 min in a high-speed refrigerated centrifuge (Hitachi, Tokyo, Japan), resuspended in PBS (pH 7.5) and lysed by passing through a high-pressure homogenizer at 1000 Bar. After centrifugation at 10,000× rpm and 4°C for 10 min, the supernatant was filtered through a 0.45 μm pore size filter and loaded onto glutathione sepharose beads (GE Healthcare, Uppsala, Sweden). The proteins were eluted with elution buffer (9.5 ml ddH 2 O, 50 mM Tris-HCl, 10 mM L-glutathione reduced glutathione). The GST-tag on the N-terminus was cleaved using TEV (Tobacco Etch Virus) protease (Solarbio, Shanghai, China) at 4°C for 8 h.

Preparation of anti-rBoTRAP1 immune serum
Polyclonal antibodies against BoTRAP1 protein were prepared in Japanese white rabbits (n = 2) according to the established immune procedure. Briefly, seven days after purchase, the rabbits were subcutaneously immunized with the purified 500 μg recombinant protein emulsified in equal amounts of Freund's complete adjuvant (Sigma, Shanghai, China), followed by a second immunization with one half of the same recombinant protein emulsified in equal amounts of Freund's incomplete adjuvant (Sigma, Shanghai, China) at a 14-day interval. This was followed by another three immunizations each at 7-day intervals. After the titer of the antisera was assayed by ELISA (the recombinant protein was coated), blood samples were collected from the carotid artery of the rabbits to prepare the polyclonal antibodies. Total immunoglobulin Gs (IgGs) were purified from rabbit sera through a Protein A chromatography column according to the manufacturer's instructions (Beyotime Biotechnology, Shanghai, China).

Identification of immunogenicity and native BoTRAP1 by Western blot
The B. orientalis-infected RBCs (1 ml) were suspended in an equal amount of phosphate-buffered saline (PBS) buffer, and supplemented with red blood cell lysis buffer preheated at 37°C (18 ml), followed by heating in a water bath for 5 min and centrifugation at 2450× rpm for 5 min. The resulting supernatant was collected and centrifuged at 12,000× rpm for 20 min. After a high-speed centrifugation, pellets (parasite fraction) from previous steps were collected and washed three times in cold phosphate-buffered saline (PBS) buffer, followed by the addition of an equal amount of PBS and centrifugation at 15,000× rpm for 20 min. After discarding the supernatant, the pellets were suspended in PBS and stored at -20°C.
The rBoTRAP1 and B. orientalis lysates were separately subjected to 12% SDS-PAGE using the standard method, followed by electroblotting onto a nitrocellulose membrane. The membranes for the Western blot were blocked with 0.05% Tween-20 in TBS (TBST) plus 1% BSA overnight at 4°C, and then separately probed with the anti-BoTRAP1 PcAb or positive serum of B. orientalis diluted with TBST (1:200) at 37°C for 1 h. The membranes were washed 5 times in TBST and incubated with secondary antibodies diluted with TBST (1:1000, HRP labeled goat anti-rabbit IgG) at 37°C for 1 h. After rinsing 5 times in TBST, the positive bands on membranes were visualized using the ECL method.

Localization of BoTRAP1 by indirect florescent antibody test (IFAT)
The B. orientalis-infected water buffalo RBCs (3% parasitemia) were smeared on glass slides and fixed in cold 100% methanol (-20°C). After three washes, blood smears were permeabilized by using 0.1% Triton-100 for 30 min, followed by incubation separately with anti-BoTRAP1 PcAb and pre-immune rabbit sera diluted 200× with 1× PBS/1% BSA for 30 min. The secondary antibody was anti-rabbit Alexa-594 diluted 1000× with 1× PBS/1% BSA and the parasite nucleus was stained with a Hoechst stain. Finally, the coverslips with cells were mounted on a slide in 10 μl anti-fluorescence quenching agent and observed with a confocal laser scanning microscope.

Sequence analysis and crystal structure modeling of BoTRAP1
The BoTRAP1 gene has a significant degree of similarity to the previously reported BbTRAP1 (XM_001609738.1), and was analyzed by Conserved Domain Search Service (CD Search) of NCBI (https://www.ncbi.nlm.nih.gov/Structure/

Expression of the recombinant BoTRAPl
The full and partial BoTRAP1 sequences were obtained separately from B. orientalis cDNA by PCR (Fig. 2). For molecular characterization of BoTRAP1, the truncated fragment was expressed in E. coli with a TEV protease site, and the size of rBoTRAP1 was~55 kDa including a 26 kDa GST-tag (Fig. 3, Lane 1). The rBoTRAP1 was purified by glutathione sepharose beads (Fig. 3, Lane 4), and the GST tag on the N-terminus was removed by using the TEV protease (Fig. 3, Lane 5).

Immune reactivity of rBoTRAP1 and identification of native BoTRAP1
The potential of BoTRAP1 as vaccine candidate antigen for preventing and controlling B. orientalis infection in water buffalo was tested by Western blot, and two bands corresponding to 55 kDa (GST fusion) and 28 kDa (GST removed) rBoTRAP1 were detected by using positive serum in contrast to no signal observed with negative sera used as control (Fig. 4a). The results indicated that rBoTRAP1 could differentiate B. orientalis positive and negative serum samples from water buffalo. To identify native BoTRAP1 in B. orientalis, anti-rBoTRAP1 immune sera were produced in rabbits, purified with a Protein A chromatography column, and used to detect B. orientalis lysates by Western blot. One band corresponding to~65 kDa native BoTRAP1 was detected (Fig. 4b), and no signal was found in uninfected water buffalo RBCs and pre-immune sera (data not shown).
For further characterization of BoTRAP1, we performed an IFAT separately using anti-BoTRAP1 rabbit immune serum and pre-immune sera. The results showed that BoTRAP1 is localized on the apical end of intracellular or extracellular parasite under a confocal laser scanning microscope (Fig. 5a). No signals were observed in parasites when tested using pre-immune sera (Fig. 5b).

Homology modeling of BoTRAP1
To further confirm the potential role of BoTRAP1 in the invasion of parasites into host cells, the 3D structures of BoTRAP1 were homology-modeled using the crystal structure of P. vivax TRAP protein as a template. As shown in Fig. 6b, the predicted BoTRAP1 3D model contain vWFA and TSP-1 domains, with an inserted domain in the integrin of the vWFA domain, containing a metalion-dependent adhesion site (MIDAS) for interaction with ligand. While the amino acid sequences are not well conserved between BoTRAP1 and PvTRAP (28% similarity), their overall structures share a conserved feature (Fig. 6a, b).

Discussion
Four TRAP genes have been identified in B. bovis and named BbTRAP1-4 (GenBank accession numbers XM_001609738, XM_001609762, XM_001609736 and XM_001609760, respectively). The BbTRAPs share a typical structure with other apicomplexan TRAPs, including vWFA and TSP domains as well as a transmembrane region followed by a cytoplasmic Cterminal tail domain, except for the absence of the TSP domain from BbTRAP4 [19]. The TRAP family has also been identified in B. orientalis and named as TRAP1-4 (Data for TRAP2-4 not shown). Sequence analysis indicated that the four TRAP genes are very similar in both B. bovis and B. orientalis. The TRAPs are secreted from micronemes as a transmembrane protein and are present in all apicomplexan parasites, with a high similarity for their cellular localization in Plasmodium spp., Neospora caninum, B. bovis and B. Pre-immune sera were used as negative control for the validation test (c). PcAb-BoTRAP1 (Red) reacted with native TRAP1 on merozoites. Nuclei were counterstained by Hochest (blue). Scale-bars: 1 μm gibsoni [7,10,19,20]. Confocal laser microscopy revealed that the native BoTRAP1 is located in the apical end of parasites, which is consistent with the TRAP localization of other parasites. On Western blot, the anti-rBoTRAP1 antibody detected the native protein as a~65 kDa band from B. orientalis lysates. This protein size is in agreement with the theoretical mass of native BoTRAP1 (65 kDa) and is approximately the size of BbTRAP1 (73 kDa). Together, these results further verified that the sequence we obtained from B. orientalis cDNA precisely is the BoTRAP1 gene.
Furthermore, we tried to determine the ectodomain crystal structure of BoTRAP1 (vWFA and TSP-1 domains). Despite a high purity and high concentration, the rBoTRAP1 failed to grow crystals in crystal buffer (576-well) (Hampton, California, USA). This result is probably attributed to the use of the prokaryotic, rather than the eukaryotic, expression system, leading to the absence of glycosylation on the surface of the recombinant protein. Fortunately, the overall structure of TRAP family is quite conservative across all apicomplexan parasites [21,22].
For further characterization, the BoTRAP1 ectodomain structure was simulated by homologous modeling.
The interaction between the parasite ligands and the corresponding host-cell receptors is essential to the invasion of apicomplexan parasites [4,7]. Currently, the invasion mechanism of apicomplexan parasites is poorly understood. A related study has shown that TRAPs participate in the formation of actin-myosin motor and drive the parasites into the host cells, enabling the translocation of TRAPs to the posterior pole of the parasites along the cytomembrane surface [23]. However, the actin-myosin motor was mediated by a sub-terminal tryptophan residue within the cytoplasmic tail of TRAPs, and the tryptophan residue connected to aldolase is rather conserved across apicomplexan parasites [4]. Interestingly, in order to release the parasites into host cells completely, the TRAPs translocated to the posterior end from the apical end of parasites need to be cleaved at their transmembrane domains by the rhomboid proteases (ROMs) [24,25]. In this study, we predicted the  . The structure of PvTRAP is shown as a cartoon (left) and surface structure (right). The Mg 2+ ion at the ligand binding site was colored in yellow. The α-helices were shown as red, and yellow indicates β-strands. b Predicted BoTRAP1 structure (residues 44-297). The 3D structure in identical orientations contains six α-helices (α1 and α6) and seven β-strands (β1 to β7). The BoTRAP1 vWFA domain also includes a metalion-dependent adhesion site (MIDAS) for binding ligand. The figure on the left shows the cartoon structure of BoTRAP1, and the figure on the right shows BoTRAP1 surface structure Fig. 7 The predicted transmembrane regions and the cytoplasmic tails of apicomplexan TRAP-family proteins. The putative rhomboid cleavage sites are marked with asterisks, the predicted transmembrane domains by rectangles, and the sub-terminal tryptophan residues with an arrow transmembrane region and the putative rhomboid cleavage site of BoTRAP1, as well as detected the conserved tryptophan residue present, in the cytoplasmic tail domain (Fig. 7).
In malarial parasites, the GPI-linked protein semaphorin-7A (CD108) had been identified as an erythrocyte surface receptor for the P. falciparum MTRAP [26]. However, the erythrocyte receptor is still unknown for the BoTRAPs. Therefore, exploring the interaction of BoTRAPs with erythrocyte surface would facilitate the understanding of the invasion mechanism of B. orientalis and the discovery of new blocking agents or methods for controlling B. orientalis.

Conclusions
Overall, BoTRAP1 was cloned, sequenced and expressed, with the native BoTRAP1 being characterized as~65 kDa by Western blot and located in the apical end of the parasites. The simulation structure and functional analysis revealed that BoTRAP1 might perform critical functions in aiding parasite invasion and interaction of parasite-host RBCs. Blocking TRAP1 function might prevent the parasites from invading RBCs. Therefore, the newly identified and characterized BoTRAP1 may have potential as drug target or vaccine candidate antigen for preventing and controlling B. orientalis infection in water buffalo.