Rhoptry protein 5 (ROP5) Is a Key Virulence Factor in Neospora caninum

Neospora caninum, of the Apicomplexa phylum, is a common cause of abortions in cattle and nervous system dysfunction in dogs. Rhoptry proteins of Apicomplexa play an important role in virulence. The objectives of this study were to study functions of NcROP5 in N. caninum by deleting the NcROP5 gene from the wild Nc-1 strain. We selected NcROP5 in ToxoDB and successfully constructed an NcROP5 gene-deleted vector, pTCR-NcROP5-CD KO. Then we screened the NcROP5 knockout strains (ΔNcROP5) at the gene, protein and transcription levels. Plaque assay, host cell invasion assay and intracellular proliferation test showed that the ΔNcROP5 strain had less plaque space, weakened invasion capacity and slower intracellular growth. Animal testing showed significantly lower cerebral load of ΔNcROP5 than the load of the Nc-1 strain, as well as a loss of virulence for the ΔNcROP5 strains. Phenotypic analyses using the label-free LC-MS/MS assay-based proteomic method and KEGG pathway enrichment analysis showed a reduction of NcGRA7 transcription and altered expression of multiple proteins including the apicomplexan family of binding proteins. The present study indicated that ROP5 is a key virulence factor in N. caninum in mice. The proteomic profiling of Nc-1 and ΔNcROP5 provided some data on differential proteins. These data provide a foundation for future research of protein functions in N. caninum.


INTRODUCTION
Neospora caninum is a common cause of abortions in cattle and nervous system dysfunction in dogs (Hall et al., 2005;Lyon, 2010). The parasite has become an international concern due to the connection of the parasite infection to abortions in dairy and beef cattle worldwide (Dubey, 1999). Antibodies to N. caninum have been detected in humans in Brazil, Korea, Northern Ireland, and the United States, although viable parasites have not been isolated (Dubey et al., 2007). Investigations are ongoing to determine whether the pathogen is a potential threat to human health.
Neospora caninum and Toxoplasmosis gondii belong to the same Apicomplexa family and share similar morphology, host range and clinical symptoms, although they appear to be divergent in their pathogenicity in mice (Lyon, 2010). A tachyzoite of the T. gondii RH strain can lead to death in a mouse (Saeij et al., 2006;Taylor et al., 2006), whereas the same outcome requires more than 10 6 tachyzoites of the N. caninum Arranz-Solís et al., 2015). Rhoptry proteins (ROPs), including TgROP5, TgROP16, TgROP17, and TgROP18 are impotant virulence factors of T. gondii (TgROPs) (Yamamoto et al., 2009;Etheridge et al., 2014;Behnke et al., 2015;Shwab et al., 2015). Deletion of the TgROP5 gene, which forms a complex with TgROP18, results in a complete loss of virulence (Reese and Boothroyd, 2011). We hypothesized that the above genes had similar functions in N. caninum. Unexpectedly, we found that ROP18 in N. caninum was a pseudogene . We suspected that another rhoptry gene of N. caninum might play a crucial role in virulence. Transgenic strains of T. gondii and Plasmodium falciparum of Apicomplexa parasites have been widely used to study the functions of parasite proteins (Ghorbal et al., 2014;Hui et al., 2014;Shen et al., 2014;Hammoudi et al., 2015;Williams et al., 2015;Li et al., 2016). However, no laboratory strains and genomeediting techniques have been developed and no transgenic strains constructed for the study of protein functions of N. caninum. We sought to construct a ROP5-knockout strain of N. caninum and determine the function of ROP5 in this parasite.
Firstly, we constructed a recombinant plasmid containing the untranslated region of the NcROP5 gene and obtained a NcROP5-deficienct strain ( NcROP5). The NcROP5-deficienct parasite exhibited, weakened host invasion, proliferation and virulence, compared with the wild parent strain. The technique we developed and the new transgenic strain laid ground work for future research of this highly pathogenic parasite.

Ethics Statement
The experiments were performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals, Ministry of Science and Technology, China. All experimental proceduresw ere approved by the Institutional Animal Care and Use Committee of China Agricultural University (under the certificate of Beijing Laboratory Animal employee ID: 18049). The mice were humanely euthanized by cervical dislocation after anesthetization by subcutaneous injection of atropine (0.02 mg/kg) when they were unable to reach food or water for more than 24 h and lost 20% body weight.

Cell Culture
Vero cells and human foreskin fibroblast (HFF) cells were purchased from the ATCC (Manassas, VA, USA) and serially passaged in our laboratories as previously described (Gaskell et al., 2009). The cell lines were cultured in Dulbecco's modified Eagle's medium (DMEM, MACGENE, Beijing, China) supplemented with 10% heat-inactivated fetal bovine serum (FBS, Gibco, USA) in a humidified incubator containing 5% CO 2 at 37 • C. The cells were passaged every 3-6 days by trypsinization.

Construction of the NcROP5-Deficient N. caninum Strain
The pTCR-CD plasmid was used at the Key Laboratory of Animal Parasitology (Beijing City, China) as previously described (Li et al., 2016). The plasmid contains the chloramphenicol resistance gene (CmR), red fluorescence protein gene (RFP), 5-fluorine cytosine gene (CodA), and ampicillin resistance gene (Amp). Nctublin promoter was used to drive the CmR-RFP fusion gene and CodA gene expression. The fragments of the NcROP5 5 (2826 bp) and 3 UTR regions (2297 bp) were amplified using the primer pairs F1/R1 and F2/R2 (Supplementary Table S1), which were designed based on the N. caninum Liverpool protein coding gene sequence 1 (FR823393, Chromosome: XII). Primers F1, R1, F2, and R2 introduced KpnI, XhoI, XmaI, and SpeI sites, respectively (Supplementary  Table S1). After amplification, the fragments were cloned into a pEASY-T-Blunt vector (TransGen Biotech Co., Ltd, Beijing). Then the DNA insert was sequenced (Beijing Sunbiotech Co. Ltd., China) and analyzed by blast in ToxoDB 1 to confirm the authenticity of the cloned sequence. All alignments were performed using DNAMAN version v5.2.2 and the correct plasmids were double digested with KpnI and XhoI for the 5 UTR and XmaI and SpeI for the 3 UTR (NEB, USA). Then, the fragments were ligated to the pTCR-CD vector to produce the NcROP5 gene deleted vector, pTCR-NcROP5-CD KO.
All constructs were verified by sequencing (Beijing Sunbiotech Co. Ltd, China). After digestion with NotI (NEB, Beijing), the linearized pTCR-NcROP5-CD KO plasmids were purified using ethanol precipitation and then resuspended using cytomix . Nc-1 tachyzoites (1 × 10 7 ), mixed with 50 µg of the above plasmids, were transferred to a 0.4 mm gap cuvette and electroporated with 2 kV at 25 µFd and 50 with the Gene Pulser Xcell electroporation system (Bio-Rad, USA). Then the parasites cultured in an atmosphere containing 5% CO 2 at 37 • C for 24 h prior to the addition of chloramphenicol (20 µM) and 5-fluorine cytosine (40 µM). The parasites were cultured consecutively to the 8th generation and then screened to confirm the purity of the selected strains.
For complementation analysis, pDMG plasmid Tao et al., 2014;Li et al., 2016) was used to express NcROP5 in the NcROP5 strain. The complete coding sequence of NcROP5 was amplified using primers F3/R3 (Supplementary Table S1) that appended the flanking EcoRV and AvrII restriction endonuclease sites. The amplification products were introduced into the same sites of a modified pDMG plasmid with GFP replaced by HA (Li et al., 2016). The pDMG-NcROP5-HA vector was electroporated into the NcROP5 strain.
The transgenic parasites were grown under pyrimethamine selection pressure. The selected NcROP5 complementary strain, named as i NcROP5, was identified by western blotting and immunofluorescence assays (IFAs).

PCR
To screen the NcROP5 knockout parasites ( NcROP5), we designed primers targeting the NcROP5 fragment (ToxoDB: NCLIV_060730, 800 bp) to detect if the gene could be amplified from different clones. Polymerase chain reaction (PCR) was conducted using the T100 TM Thermal Cycler (Bio-Rad, USA) with PCR SuperMix (TransGen, China) following manufacturer's instructions. The primer sequences were F4/R4 (Supplementary Table S1). The Nc5 gene (336 bp) served as the internal reference with the primer pair F5/R5 (Supplementary Table S1). The PCR conditions were as follows: 95 • C for 5 min, 30 cycles of 95 • C for 30 s, 56 • C for 30 s and 72 • C for 1 min, and 72 • C for 10 min. The PCR product was identified using electrophoresis.

Western Blot
Purified parasites (Nc-1, three NcROP5 clones and one i NcROP5) were lysed using RIPA buffer (Beyotime, Beijing) supplemented with a cocktail of protease inhibitors (Sigma, USA). The lysates were resolved on a 10% (w/v) SDS-PAGE gel. After electrophoresis, separated strips were transferred onto polyvinylidene fluoride (PVDF) membranes (Millipore, MA, USA) together with a visible pre-stained protein marker (TransGen Biotech Co., Ltd, Beijing). The membranes were blocked with 5% (w/v) skim milk (BD Difco, USA) in PBS for 1 h at 37 • C and then incubated with the mouse anti-rNcROP5 antibody  (National Animal Protozoa Laboratory in China Agricultural University, 1:500) for NcROP5 and the mouse anti-HA antibody (Sigma-Aldrich, 1:500) for i NcROP5. The N. caninum F-actin subunit beta (NcActin, National Animal Protozoa Laboratory in China Agricultural University) was used as a loading control and was incubated with the rabbit anti-rNcActin antibody (1:2000). After five washes with PBST (1% Tween-20), the membranes were incubated for 1 h with a horseradish peroxidase (HRP)labeled goat anti-mouse IgG (H+L) secondary antibody (Sigma, USA) diluted 1:5000 and an HRP-labeled goat anti-mouse IgG secondary antibody (Sigma, USA) diluted 1:10,000 in PBS with 5% BSA. Finally, the reactive bands were visualized using enhanced chemiluminescence reagents (CoWin Biotech Co., Ltd, Beijing).

Immunofluorescence Assay (IFA)
Immunofluorescence Assay was used to detect the NcROP5 subcellular localization as previously described (Li et al., 2016). Appropriate numbers of parasites were coated on glass coverslips in 12-well plates and then fixed in 4% paraformaldehyde, permeabilized with 0.25% Triton X-100 and blocked with 3% BSA. Subsequently, the cells were incubated with a mouse anti-HA monoclonal antibody, mouse anti-NcROP5 and rabbit anti-NcSRS2 antibody (National Animal Protozoa Laboratory in China Agricultural University) followed by FITCconjugated goat-anti mouse IgG (H+L), Cy3-conjugated goatanti mouse IgG (H+L), and Cy3-conjugated goat-anti rabbit IgG (H+L) (Sigma, USA). The nuclear DNA was stained with Hoechst33258 (Sigma, USA). Extracellular i NcROP5 parasites were identified by IFA. The images were obtained using a Leica confocal microscope system (Leica, TCS SP52, Germany).

Invasion Assays
Cell invasion assays were performed using the following protocol . Freshly egressed parasites (Nc-1, NcROP5 and i NcROP5 strains) were allowed to settle onto the HFF monolayers on 12-well plates for 30 min at 37 • C prior to incubation. Extracellular parasites were removed by washing three times with PBS. The cells were incubated for 20 h prior to IFA using anti-NcSRS2 primary and FITC-conjugated secondary antibodies (1:100 dilution). The number of vacuoles representing successful invasion parasites was counted under a fluorescence microscope. A total of three replicates (each clone in duplicate) were performed.

Plaque Assay
Plaque assays were set up immediately following the initiation of the invasion assays using the same parasite suspensions. Nc-GFP strain (kindly provided by Professor Xuenan Xuan, Obihiro University of Agriculture and Veterinary Medicine, Japan), as a null vector constructing group, which is a transgenic parasite by transfecting Nc-1 wild-type strain with the pDMG plasmid, was propagated as tachyzoites by serial passages in HFF cells as previously described (Gaskell et al., 2009). The plaque assay was performed as previously described . Five hundred parasites were added to each well of a 6-well plate containing confluent host cell HFF monolayers. The plate was incubated for 7 days in a 37 • C incubator with 5% CO 2 , fixed with 4% paraformaldehyde for 15 min and stained with 2% crystal violet for 15 min at room temperature. The stained wells were washed with deionized water, air dried and visualized using a fluorescence microscope (IX71, Olympus, Japan). The plaque area measurement was performed as previously described (Roos et al., 1994).

Replication Assay
The Nc-1, NcROP5 and i NcROP5 parasites were harvested and counted. Five-hundred parasites were added to confiuent vero and HFF cells on 25 mm coverslips. The coverslips were allowed to sit at 37 • C for 30 min in DMEM with 20% FBS and then washed three times with PBS. The coverslips were incubated at 37 • C for 24 h in DMEM with 2% FBS and then fixed with 4% paraformaldehyde for 15 min. IFAs were performed using anti-NcSRS2 primary  and FITC-conjugated secondary antibodies (Proteintech, USA). The number of parasites per vacuole was counted under a fluorescence microscope. The experiments were performed in triplicates and repeated three times, and 200 vacuoles were counted per coverslip.

Egress Assay
Neospora caninum tachyzoites, as well as T. gondii, must form parasitophorous vacuole (PV) after invading the host cell (Tilley et al., 1997). They proliferate in the PV and then egress out of PV. Then, the tachyzoites invades other cells again (Soldati and Meissner, 2004). In our experiments, parasites were inoculated in 12-well plates containing prepared HFF cells for 30 h. Egress was stimulated with 10 µM of the Ca 2+ ionophore A23187 from Streptomyces chartreusensis (Sigma, USA) for 2 min at 37 • C prior to fixation with paraformaldehyde (Williams et al., 2015). The IFA was performed as described above using the mouse anti-NcSRS2 primary and FITC-conjugated secondary molecules. The average number of ruptured vacuoles was determined by counting a minimum of 100 vacuoles per slide and three slides for each experiment. Three independent experiments were performed.

Gene Expression Assay
Total RNA was extracted from 1 × 10 7 tachyzoites of the Nc-1 wild-type strain and three NcROP5 clones with TRIzol reagent (Invitrogen, USA). cDNA was synthesized using the EasyScript First-Strand cDNA Synthesis SuperMix kit (TransGen, China). Specific primers were designed using Primer Premier 5.0 , including primers for rhoptry necks (RON2 and RON4), rhoptrys (ROP4, ROP5, ROP7, ROP16, and ROP17), dense granules (GRA2, GRA6, and GRA7) and the endogenous reference gene NcActin (Supplementary Table S1). The specificity of these primers was evaluated using conventional quantitative real-time PCR (qRT-PCR). The qRT-PCR was conducted using the ABI Prism 7500 System (Biosystems Inc., USA) with SYBR Green II (Takara Biotechnology, Dalian, Co., Ltd) following manufacturer's instructions. The resulting RNA concentrations were normalized using Ncactin (Wang et al., 2017), and the relative expression levels of the target genes were analyzed using the ABI Prism 7500 software v2.0.5 (Biosystems Inc., USA). The RT-PCR conditions were as follows: 94 • C for 5 s, 40 cycles of 94 • C for 5 s and 60 • C for 30 s. The relative expression of genes was calculated using the 2 − Ct method and standard deviation was calculated from three replicates (Cárdenas-Mondragón et al., 2016;Gagnaire et al., 2016).

Label-Free LC-MS/MS Assay
The N. caninum Nc-1 and NcROP5 strains were cultured as tachyzoites in HFF cells as previously described (Gaskell et al., 2009) for 96 h. For each parasite, a total of 5 × 10 7 tachyzoites were harvested and isolated by washing in cold PBS, centrifugation, resuspension in cold PBS, syringing three times through a 27-gaure needle, filtering through a 5.0 µm pore filter (Millipore, MA, USA), washing three times with PBS, and finally centrifugation at 1800 rpm for 10 min . Then the parasites were ultrasonically lysed in protein lysis buffer (8.0 M urea, 100 mM pH 8.0 Tris-HCl, 1×cocktail) using an ultrasonic liquid processor (Sonic & Materials INC., USA). The lysates were sonicated five times for 2 s each at 30% amplitude. Then the sonicated lysates were centrifuged for 30 min at 4 • C at 12000 g and obtained the supernatants (Deutsch et al., 2014). After quantitative assay by the Bradford method (Bradford Protein Assay kit) (Deutsch et al., 2014), 60 µg protein was reduced, alkylated and trypsin digested. Each sample was separated by a high performance liquid chromatography (HPLC) system (Ultimate 3000, Thermo Scientific), and analyzed by mass spectrometry (Q-Exactive HF, Thermo Scientific). The spectra were identified and analyzed by the Mascot software (Thermo Scientific, USA) and Proteome discoverer 2.1 (Thermo Scientific, USA) on the basis of database listed in Supplementary Table S2 (Sandberg et al., 2012).

Mouse Virulence Assay
To compare the virulence of the NcROP5 strain with the parental Nc-1 strain and complemented i NcROP5 strain, tachyzoites (3 × 10 6 ) were injected i.p. into 6-week-old female BALB/c mice (five mice per strain), which were purchased from the Laboratory Animal Center of Academy of Military Medical Sciences (Beijing, China). Survival was monitored for 30 days. To examine the cerebral load of parasites, tachyzoites were injected i.p. into 6-week-old female BALB/c mice at a dose of 1 × 10 6 (five mice). The mouse brains were examined for parasites on the 30th day (Pan et al., 2014).

Statistical Analysis
The statistical analysis of all of the data was performed using GraphPad Prism 5 v. 5.01 (San Diego, CA, USA). The results were expressed as mean ± SD and evaluated by non-parametric tests. Values of P < 0.05 and P < 0.01 were considered statistically significant and very significant, respectively.

Successful Construction of the NcROP5 Knockout Strain
Using the knockout plasmid, pTCR-CD, we successfully constructed the NcROP5 knockout strain. To use the plasmid for deleting a gene, we ligated the 5 UTR and 3 UTR regions of the target gene. Then, we analyzed the genome of N. caninum Liverpool strain in ToxoDB 2 and identified the sequences with high identity to TgROP5, which were NcLIV_060730, NcLIV_060740, and NcLIV_060741. Compared to the tandem cluster of polymorphic alleles of TgROP5 (TgROP5A, TgROP5B, and TgROP5C) (Behnke et al., 2011), we believed that NcROP5 would contain NcROP5A, NcROP5B, and NcROP5C, coded by the above three genes, respectively. Homology of nucleotide and protein sequence between NcLIV_060730 and TgROP5A were 52.72 and 76.76%, respectively. Since TgROP5 is an important virulence factor of T. gondii (Behnke et al., 2012), we suspected that NcROP5 might play an active role in N. caninum virulence. To examine the physiological role of ROP5 in the Nc-1 strain, NcROP5 KO parasites were generated by targeted gene disruption. The non-coding region (UTR) sequences of NcROP5, NcROP5 5 UTR and 3 UTR, were cloned (Supplementary Figure S1A) and ligated to the pTCR-CD vector, forming pTCR-NcROP5-CD KO (Supplementary Figure S1B). The plasmid was identified by restriction enzyme digestion (Supplementary Figure S1A) and sequenced and then transfected into Nc-1. Theoretically, the entire ROP5 coding region, containing NcROP5A, NcROP5B, and NcROP5C, could be replaced with the CmR and RFP genes (Supplementary Figure S1C). RFP was expressed at normal levels in the suspected ROP5 mutants ( Figure 1A). Following CmR selection of the RFP and 5-fluorine cytosine-negative parasites, we isolated three mutants out of 218 clones.
To generate the ROP5 gene-deficient Nc-1 strain, the NcROP5 gene was targeted and replaced by CmR-RFP. PCR was used to confirm the deletion of the NcROP5 gene. NcROP5 gene amplification was negative in the parasite genomes ( Figure 1B). NcROP5 deficiency was confirmed using western blotting ( Figure 1C) and RT-PCR ( Figure 1D). We also successfully generated a strain with complete knockout of the NcROP5 gene. These results suggested that NcROP5 was not an essential, surviving gene in the parasites.

Successful Construction of the NcROP5 Complementary Strain
To verify the biological role of NcROP5, we obtained the complete coding sequence of NcROP5A without the termination codon TAA (NcLIV_060730, 1647 bp, Supplementary Figure S2A) and generated the complementary plasmid pDMG-NcROP5-HA. The expression of NcROP5 with an HA-tag and stop codon TAA at the C-terminus was driven under the control of the NcGRA1 promoter (Supplementary Figure S2B). We correctly identified the plasmid by restriction enzyme digestion and DNA sequencing (Supplementary Figure S2A), and then electroporated it into one of the NcROP5 strains. After repeated passage, we selected strains that stably expressed NcROP5 (named NcActin was used as the control. (C) IFA analysis of i NcROP5 tachyzoites. The HA-tagged ROP5 and ROP5 proteins (green) were localized at the anterior end of the parasites, and HA co-localized with ROP5. The shapes of parasites were stained with anti-NcSRS2 (red), and the nuclear DNA was stained with Hoechst (blue). Scale bar, 1.0 µm. as i NcROP5). NcROP5 expression in the complementary strain based on WB and RT-PCR was comparable to the expression of the wild strain (Figures 2A,B). IFA was used to analyze the subcellular localization of NcROP5 in the i NcROP5 strains using a mouse anti-HA monoclonal antibody. NcROP5 was colocalized with HA at the anterior end of the extracellular parasites ( Figure 2C).

The NcROP5 Gene Is Critical in the Invasion and Proliferation of N. caninum
Plaque formation measures the viability of parasites in cell culture and reflects parasite motility on the surface of the host cell layer, invasion, intracellular growth and egress (Plattner et al., 2008). We assessed the NcROP5 strain using the plaque assay. The NcROP5 strain formed noticeably smaller and considerably fewer plaques than the parental strain and the complementary strain (Figures 3A,B). The NcROP5 strain showed significant weakening of host cell invasion ( Figure 3C). Intracellular replication assays were performed to investigate if deleting ROP5 in the Nc-1 strain weakened the ability of the parasite to proliferate in host cells. Parasites that invaded host cells 24 h after inoculation were analyzed for intracellular growth by counting the number of parasites per vacuole. Due to the lack of synchronization of host cell invasion, the intracellular vacuoles contained 1, 2, 4, or 8 FIGURE 3 | Deletion of the NcROP5 gene affects the Nc-1 strain phenotype. (A,B) Plaque assays showed very low plaque formation of the NcROP5 strain compared with the strains harboring NcROP5. The assay was performed three times independently, and the plaque area was quantified from three independent experiments. At least 10 fields were quantified per strain in every experiment. Each bar indicates mean ± SD; "Pi" indicates pixels of each plaque. (C) The NcROP5 strain possess significantly reduced cell invasion compared with the Nc-1, Nc-GFP, and i NcROP5 strains. The number of parasite vacuoles was counted after immunofluorescent staining with anti-NcSRS2 antibodies in the cell invasion assay. (D) Intracellular replication assay showed that NcROP5 strain proliferated more slowly than the other three strains. The IFA was performed with anti-NcSRS2 antibodies. At least 100 vacuoles were examined for each per strain. (E) Egress assay showed no difference between strains indicating that NcROP5 is not involved in parasite egress. * P < 0.01 and * * P < 0.001 (Non-parametric tests). Scale bar, 10 µm.
parasites. However, intracellular replication of NcROP5 in Vero and HFF cells was lower than that of the wild type parasite (Nc-1). As shown in Figure 3D, the distribution of the eight parasites per vacuole was significantly lower for the ROP5 knockout strain than that for Nc-1, Nc-GFP, and i NcROP5, indicating that the intracellular replication of the NcROP5 strain was slower than that of the other three strains. However, there was no significant difference in egress between the four strains ( Figure 3E). These findings indicated that NcROP5 expression in N. caninum is a critical determinant of parasite phenotype, including host cell invasion and proliferation.

Identification of Peptides and Proteins in Nc-1 and NcROP5
To analyze the cause of the NcROP5 phenotype changes, we compared protein expression levels in Nc-1 and NcROP5. Nc-1 and NcROP5 tachyzoites were analyzed by label-free LC-MS/MS for quantification of peptides. A total of 18,809 peptides with FDR ≤ 0.01 were identified, and a total of 1149 proteins with expression difference ≥ 1.2-fold between the two strains were functionally annotated according to the N. caninum genome annotation in the ToxoDB database (Supplementary Figure S3A and Table S3). The three main annotation types were biological process, cellular component and molecular function (Supplementary Figures S3B-D). The differentially expressed proteins were classified into 31 functional groups (Supplementary Figure S3E), with biological process possessing 12 GO (gene ontology annotation) terms, cellular component 10 GO terms, and molecular function 9 GO terms. Cellular process and metabolic process were major biological processes (32 and 31%, respectively), cell part and cell are major cellular components (25% each), and binding and catalytic activities were major molecular functions (53 and 37%, respectively).

Molecularly Enriched Pathways Associated with Differentially Expressed Proteins in Nc-1 and NcROP5
To further identify biological pathways of the differentially expressed proteins, the number of proteins expressed in Nc-1 and NcROP5 and differences between the two strains are shown in Supplementary Figure S4A. Demographic analysis demonstrated that 76 proteins were increased and 89 proteins were decreased with a 2.0-fold change; 198 proteins were increased and 252 proteins were decreased with a 1.5-fold change; and 477 proteins were increased and 671 proteins were decreased with a 1.2-fold change. We predicted that 477 increased proteins participated in 92 biological pathways and 671 decreased proteins participated in 88 biological pathways by KEGG pathway analysis 3 . The differentially expressed proteins were matched with the proteins annotated in the KEGG pathway database, and they were compared to determine their involvement in the KEGG pathway. We found that five KEGG pathways were clearly affected (P ≤ 0.05) (Supplementary Figure S4B and Table S4), followed by mismatch repair, SNARE interactions in vesicular transport, glycosylphosphatidylinositol (GPI)-anchor biosynthesis, ABC transporters and RNA degradation.

Knockout of NcROP5 Reduced Parasite Virulence
To investigate if NcROP5 is necessary for Nc-1 virulence, mice were infected intraperitoneally with 1 × 10 6 Nc-1, Nc-GFP, NcROP5, or i NcROP5 tachyzoites, followed by observation of parasite load in the brain and survival rate. Significant differences were observed between the NcROP5 group and the other groups (Nc-1, Nc-GFP and i NcROP5). Signs of illness such as rough coats, inactivity, or nervous system signs (hind limb weakness, head tilting, or walking in circles) were observed 2 or Frontiers in Microbiology | www.frontiersin.org FIGURE 4 | The transcription of differential genes in Nc-1, NcROP5 and i NcROP5. Quantitative RT-PCR analyses of the transcription levels of the indicated genes in NcROP5 parasites compared with the parental strain (Nc-1) and the complementary strain (i NcROP5). Each bar indicates the relative quantity (RQ) ±SD. The RQs of the NcROP5 and i NcROP5 parasites were expressed as fold of the expression by Nc-1, which was artificially set as 1. The data presented are representative of three independent experiments each performed in duplicate. * P < 0.05.
FIGURE 5 | The ROP5 gene is critical in N. caninum virulence. (A) Brain parasite load of NcROP5 was significantly reduced compared to Nc-1, Nc-GFP, and i NcROP5. BALB/c mice (n = 5) were injected i.p. with 1 × 10 6 tachyzoites and the number of parasites in the brain was determined 30 days after parasite inoculation. (B) All NcROP5-inoculated mice survived compared with 40-60% survival of Nc-1, Nc-GFP and i NcROP5-inoculated animals. BALB/c mice (n = 5) were infected with 3 × 10 6 tachyzoites and monitored for 30 days. Two independent experiments were performed.
5 days before death. In addition, rough hair coats and reduced activity were observed in 2 Nc-1-infected mice from day 28 postinfection, 1 Nc-GFP-infected mouse from day 29 post-infection and 1 i NcROP5-infected animal from day 26 post-infection. There were no clinical signs in NcROP5-infected mice during the 30 days observation period. The cerebral parasite load of the NcROP5 group was approximately one-fourth of the load in the Nc-1 and Nc-GFP groups. The i NcROP5 group had the most parasites in the brain (Figure 5A). Consistent with the low parasite load in the brain, the survival rate of the NcROP5infected animals was higher than that of the other groups. Actually, all NcROP5-infected mice survived, compared with the survival rate of 40-60% for the other groups ( Figure 5B). These results demonstrate a requirement of NcROP5 for Nc-1 virulence.

DISCUSSION
Neospora caninum and T. gondii are obligate intracellular parasites that belong to the phylum Apicomplexa, which includes Cryptosporidium, Eimeria, Plasmodium, and Theileria (Hajj et al., 2007;Plattner et al., 2008). Toxoplasma is a widely used model organism for studying protein functions in Apicomplexa (Shen et al., 2014). In T. gondii, the methods of gene deletion and protein expression regulation have been developed to study protein functions and provided valuable information (Payne et al., 2011). However, related research reports about functions of proteins in N. caninum are limited to epidemiological studies, the expression and localization of cloned genes, and immunogenicity (Liu et al., 2013;Li et al., 2014;Hecker et al., 2015;Pastor-Fernández et al., 2016). The overexpression of some genes and addition of labeling with biotin in N. caninum have been reported Mota et al., 2016). Few reports on N. caninum have investigated the loss of one or more genes and their effects on gene functions (Mota et al., 2016).
The genome-editing technique is very important and has advanced researches on the roles of gene functions in some organisms. We constructed a pTCR-CD knockout plasmid and used the plasmid to edit genes in N. caninum successfully with the example of the NcROP5 gene. This step was a very positive development for research on N. caninum protein functions. The roles of NcROP5 were partially deduced using the strain lacking the NcROP5 gene based on the PCR assay results, gene transcription levels and protein expression levels. The study provides a reference method and lays the foundation for further studies on NcROP5 functions and other tachyzoite proteins. We selected three NcROP5 monoclonal strains and observed their phenotypes. The phenotypes of the NcROP5 strains were similar, with a weakened invasion capacity, slower intracellular proliferation rate and smaller sized plaques. However, there were significant differences compared to the Nc-1, Nc-GFP, and NcROP5 complement strains. The slight phenotypic differences between i NcROP5 and the reference strains might be due to the efficient promoter NcGRA1 of pDMG-NcROP5-HA, while ROP5 is encoded in a multicopy gene in the reference strain Nc-1. The results indicated that the NcROP5 gene played an important role in N. caninum virulence. The phenotypic changes in the NcROP5 strains were probably caused by the abnormal expression of other genes associated with NcROP5.
Proteomics analysis was performed to explore the proteins expression characteristics of Nc-1 and NcROP5, but regrettably there were few interesting proteins secreted by microneme, rhoptry, or dense granule. GRA7 transcription level was reduced in NcROP5 parasites, and for the reason that qRT-PCR was more sensitive than label-free analysis. In this part, we noticed some mysterious proteins, AP2 protein family. The apicomplexan AP2 (ApiAP2) family of binding proteins might be a major class of transcriptional regulators in Apicomplexan (Radke et al., 2013). There are 68 ApiAP2-encoding genes predicted in the Toxoplasma genome (Altschul et al., 2010). We crudely inferred the functions of the genes according to other Apicomplexan parasites, because they have not been fully studied in N. caninum. In T. gondii, AP2XI-5 targets more than 300 gene promoters and actively controls the transcription of these genes. AP2XI-5 was reported to be enriched at numerous active promoters, especially crucial virulence factors such as the rhoptry and microneme proteins (Walker et al., 2013b). In addition, TgAP2XI-4 regulates bradyzoite gene expression during parasite differentiation and cyst formation (Walker et al., 2013a), and TgAP2XI-3 dampens proliferation of parasites (Lai et al., 2012). In our study, the expression levels of AP2VIII-2, AP2XII-4, and AP2XII-5 were increased and AP2VIII-3 was not detected in the NcROP5 group. These findings might be related to the loss of NcROP5 or abnormal expression of virulence factors NcROPs, MICs, GRAs and some metabolic enzymes as result of NcROP5 deletion. Meanwhile expression of some proteins, such as actin, myosin, and cell division protein, was decreased in the NcROP5 strain, consistent with the slow proliferation of NcROP5. The lower expression of microneme protein 5 (MIC5) was one possible reason for the weaker invasion capacity of NcROP5, because MIC5 as a secreted protein modulates parasite adhesion and invasion (Brydges et al., 2008;Gaji et al., 2011).
Toxoplasmosis gondii strains are primarily divided into three types based on genomic analysis (type I, type II, and type III) (Grigg et al., 2001;Su et al., 2003). Previous studies identified a closely related ROP5 gene of T. gondii that was absolutely critical for pathogenesis in mice; the deletion of this gene resulted in a complete loss of virulence in mice even following infection with 10 6 parasites (Reese and Boothroyd, 2011). For N. caninum, various authors have described evidence of different toxicity among isolates. They demonstrated Nc-liv and Nc-Spain 7 had a higher virulence than Nc-1 Caspe et al., 2012), and Nc-Spain 2H, Nc-Spain 3H, Nc-Spain 5H, Nc-Spain 6, Nc-Spain 7, and Nc-Spain 9 isolates showed different parasite burdens and brain lesions during the late phase of infection (Pereira et al., 2010). Although among N. caninum isolates, Nc-1 is less virulent in cattle than other isolates, it still leads to nervous system deficiency and deaths (Pan et al., 2014). In our study, the pathogenicity of the NcROP5 strain from Nc-1 was greatly decreased, similar to TgROP5. The phylogenetic tree was showed that NcROP5 is closely related to TgROP5 (Ghorbal et al., 2014). ROP5 likely has similar functions in Apicomplexa and mediates parasite virulence (Reese and Boothroyd, 2011;Walzer et al., 2013). In Toxoplasma, ROP5, ROP18, and GRA7 form a complex that plays critical roles in parasite virulence (Hermanns et al., 2016). ROP5 and ROP17 were also reported to control T. gondii RH strain virulence through the complex (Etheridge et al., 2014). However, ROP18 has been proven to be a pseudogene and not a virulence factor in N. caninum (Reid et al., 2012). Surprisingly, the loss of NcROP5 led to a reduction in NcGRA7 transcription in our study. In contrast, TgGRA7 bound the host cell immunity-related GTPase A6 and enhanced its polymerization, rapid turnover, and eventual disassembly . TgGRA7 was considered a virulence factor, which was in contrast to the NcGRA7 immunogenicity analysis (Vemulapalli et al., 2007;Alaganan et al., 2014). Therefore, further research is needed to clarify the association between NcROP5, NcROP17, and NcGRA7, although we have defined NcROP5 as a factor that affects N. caninum virulence in mice.

CONCLUSION
Our findings here highlight the role of ROP5 as an important virulence factor in N. caninum and may contribute to a better understanding of the possible mechanisms underlying the interactions between parasites and their hosts.

ACKNOWLEDGMENTS
We are grateful to Professor Silvia N. J. Moreno (University of Georgia, USA) for kindly providing the pTCY vector. We thank Professor Xuenan Xuan (Obihiro University of Agriculture and Veterinary Medicine, Japan) for kindly providing the Nc-GFP strain. We are grateful to Dr. Jin Zhu (Therapeutic Goods Administration, Australia) for his assistance in the preparation of the manuscript.