Toxoplasma gondii Dense Granule Proteins 7, 14, and 15 Are Involved in Modification and Control of the Immune Response Mediated via NF-κB Pathway

Toxoplasma gondii infects almost all warm-blooded animals, including humans, leading to both cellular and humoral immune responses in the host. The virulence of T. gondii is strain specific and is defined by secreted effector proteins that disturb host immunity. Here, we focus on nuclear factor-kappa B (NFκB) signaling, which regulates the induction of T-helper type 1 immunity. A luciferase assay for screening effector proteins, including ROPs and GRAs that have biological activity against an NFκB-dependent reporter plasmid, found that overexpression of GRA7, 14, and 15 of a type II strain resulted in a strong activity. Thus, our study was aimed at understanding the involvement of NFκB in the pathogenesis of toxoplasmosis through a comparative analysis of these three molecules. We found that GRA7 and GRA14 were partially involved in the activation of NFκB, whereas GRA15 was essential for NFκB activation. The deletion of GRA7, GRA14, and GRA15 in the type II Prugniaud (Pru) strain resulted in a defect in the nuclear translocation of RelA. Cells infected with the PruΔgra15 parasite showed reduced phosphorylation of inhibitor-κBα. GRA7, GRA14, and GRA15 deficiency decreased the levels of interleukin-6 in RAW246.7 cells, and RNA-seq analysis revealed that GRA7, GRA14, and GRA15 deficiency predominantly resulted in downregulation of gene expression mediated by NFκB. The virulence of all mutant strains increased, but PruΔgra14 only showed a slight increase in virulence. However, the intra-footpad injection of the highly-virulent type I RHΔgra14 parasites in mice resulted in increased virulence. This study shows that GRA7, 14, and 15-induced host immunity via NFκB limits parasite expansion.


INTRODUCTION
The obligate intracellular protozoan parasite Toxoplasma gondii can cause congenital toxoplasmosis, opportunistic infections in immunocompromised patients, and ocular disease (1)(2)(3). Epidemiological investigation of toxoplasmosis revealed that the majority of European and North American strains of the parasite belong to three distinct clonal lineages: type I, II, and III (4). These strains differ in virulence in mice: type I strains are the most virulent with a lethal dose (LD 100 ) of one parasite, whereas the LD 50 of type II and III strains are ∼10 3 and 10 5 , respectively (5). Previous studies demonstrated that virulence is largely mediated by several families of secretory pathogenesis determinants (6). These secreted effector proteins originate from different organelles, namely the rhoptries, known as rhoptry proteins (ROPs), and dense granules, known as dense granule proteins (GRAs) (7). Recently, it has become clear that T. gondii manipulates and modulates host resistance mechanisms at multiple points along pro-inflammatory pathways, which in turn dictates parasite burden and disease (8).
Nuclear factor-kappa B (NFκB), the central mediator of inflammatory responses and immune function, comprises homoand heterodimers of five members: NFκB1 (p50), NFκB2 (p52), RelA (p65), RelB, and c-Rel (9,10). The NFκB complex structure resides in the cytoplasm of unstimulated cells, where it is complexed with the inhibitor-κB (IκB) family of proteins, such as IκBα, IκBβ, and IκBε, which bind to the NFκB DNA binding domain and dimerization domain, the Rel homology domain, and thereby interfere with the function of the nuclear localization signal (11). Upon exposure to various infectious and inflammatory stimuli, the inhibitor proteins are phosphorylated, resulting in their ubiquitination and degradation, allowing the nuclear translocation of NFκB dimers to regulate gene expression (10). Many pathogens, including viruses, bacteria, and protozoa, have been reported to modulate the host NFκB pathway to optimize survival in the host (12).
Mice lacking c-Rel and RelB are highly susceptible to intraperitoneally infection with T. gondii and die within 10-15 days of infection, indicating the importance of the NFκB pathway for an adequate response to T. gondii infection (13,14). C-Rel −/− mice show an early defect in the number of IL-12p40-producing cells among the peritoneal exudates cells collected at 12, 24, and 48 h post-infection, although within 2-3 days this defect is no longer apparent (14). Moreover, increased susceptibility of c-Rel −/− mice can be rescued by administration of IL-12 until 2 days post-infection, indicating that delayed production of IL-12 up to 2 days post-infection causes decreased production of IFNγ and a failure to control the parasite burden (14). Despite these findings, modulation of the NFκB pathway by T. gondii remains to be further elucidated.
In this study, we used an NFκB-luciferase assay to screen candidates for their ability to regulate NFκB activity. We found that overexpression of GRA7, 14, and 15 in a type II strain resulted in strong NFκB activity; thus, we focused on these proteins. Toxoplasma GRA15 accounts for differences in NFκB activation among different strains (15). Recombinant GRA7 protein also has potent activity against the NFκB pathway; however, it is unclear whether endogenous GRA7 is capable of affecting the NFκB pathway (16). GRA14, which is secreted into the vacuole, can be transferred to both the parasitophorous vacuole (PV) membrane (PVM) and the intravehicular network (17). However, the molecular function of this protein remains unknown. Thus, the aim of this study was to gain a comprehensive understanding of the involvement of NFκB in the pathogenesis of toxoplasmosis by comparative analysis of three molecules that modulate inflammatory cytokines and chemokines, to ultimately aid the development of strategies to control chronic Toxoplasma infections.

Ethics Statement
The use and care of animals complied with the Guide for the Care and Use of Laboratory Animals from the Ministry of Education, Culture, Sports, Science, and Technology, Japan. The experimental protocol was approved by the Committee on the Ethics of Animal Experiments at the Obihiro University of Agriculture and Veterinary Medicine (permit number: 19-50). All efforts were made to minimize animal suffering.

Experimental Design
First, we constructed 17 GRAs and 21 ROPs expressing vectors. Then they were transiently transfected into 293T cells for monitoring NFκB activity. Next, 293T cells were infected with the parental Pru, Pru gra7, Pru gra14, and Pru gra15 parasite strains and assessed their effect on NFκB activity. We evaluated the nuclear translocation of RelA in 293T cells overexpressing GRA7, GRA14, and GRA15 alone. Moreover, we also examined it in HFF cells infected with each parasite strain. Then, level of phosphorylated-IκBα in HFF cells infected with parasite strains were quantified. After that, we measured level of secreted IL-6 in Raw246.7 mouse macrophage cells infected with parasite strains, and then their RNA samples were supplied for transcriptome analysis. Lastly, we conducted survival test of both mice infected with type II T. gondii strains and mice infected with type I T. gondii strains.

Parasites and Cell Culture
Toxoplasma gondii (type II, Pru ku80 hxgprt and type I, RH hxgprt, RH hxgprt gra7, and RH hxgprt gra14) was maintained in monkey kidney adherent epithelial (Vero) cells in Eagle's minimum essential medium (MEM, Sigma, St. Louis, MO, USA) with 8% fetal bovine serum (FBS) and the appropriate antibiotics. RH hxgprt gra7 and RH hxgprt gra14 were kindly gifted by Prof. John Boothroyd (Stanford University School of Medicine) and Prof. Peter Bradley (University of California), respectively. Human embryonic kidney (293T) cells, human foreskin fibroblast (HFF) cells, and Raw264.7 mouse macrophages were cultured in Dulbecco's modified Eagle's medium (DMEM; Sigma) supplemented with 10% FBS and the appropriate antibiotics. For the purification of tachyzoites, infected cells were syringe-lysed using a 27-gauge needle to release the tachyzoite-stage parasites into the medium, which was then filtered using a 5.0-µm pore-sized filter (Millipore, Bedford, MA, USA).

Plasmid Construction
All of the plasmids and primers used in this study are listed in Tables  Generation of Pru gra7 Pru gra14, and Pru gra15 Deletion Mutants, and GRA7and GRA14-Complemented Strains The knock-out plasmid (pBS/GFP/TgGRA7KO/HX) was transfected into parental Pru strains, and selected with 25 µg/ml 3-mercaptopropionlc acid and 50 µg/ml xanthine. The electroporation of tachyzoites was performed as described previously (18). The drug-resistant parasites were cloned by limiting dilution and tested by PCR (Supplemental Figure 1). PCR-positive clones were further analyzed with western blotting and indirect fluorescent antibody test (IFAT) to confirm the protein expression. To disrupt GRA14 and GRA15 in Pru, we cotransfected the parasite with 50 µg of the CRISPR plasmid (pSAG1::CAS9-U6::sgTgGRA14 and pSAG1::CAS9-U6::sgTgGRA15), along with an amplicon containing homologous regions of GRA14 and GRA15 surrounding a pyrimethamine-resistant dihydrofolate reductase (DHFR * ) cassette (5 µg), respectively. Insert fragments were prepared by PCR amplification using the primers listed in Table 2. Selection by growth for 10 to 14 days in pyrimethamine (1 µM) was used to obtain stably resistant parasite clones that were subsequently screened by PCR to ensure the correct integration of DHFR * into the GRA14 and GRA15 gene loci (Supplemental Figure 1). PCR-positive clones were further analyzed by western blotting and IFAT to confirm the loss of GRA14 expression (Supplemental Figure 3). To complement the GRA7 and GRA14 genes, we transfected GRA7-and GRA14deficient parasites with pSAG1::CAS9-U6::sgUPRT (50 µg) to target integration to the UPRT locus, along with an amplicon containing the TgGRA7 and TgGRA14 genes containing the 5′-and 3′-untranslated regions (UTRs) (5 µg), respectively. Stably resistant clones were selected by growth on fluorouracil (10 µM) for 10 to 14 days and were subsequently screened by PCR to ensure the correct integration into the UPRT gene locus (Supplemental Figure 1). PCR-positive clones were further analyzed by western blotting and IFAT to confirm the protein expression (Supplemental Figure 2).

Cytokine ELISA
Raw246.7 mouse macrophage cells in a 12-well plate were infected with parasite lines (multiplicity of infection = 0.5) for 24 h, along with control uninfected cells. Then, supernatants were collected and IL-6 levels were determined using a cytokine enzyme-linked immunosorbent assay (ELISA) kit (Mouse OptEIA ELISA set; BD Biosciences, San Jose, CA, USA).

RNA Sequencing and KEGG Pathway Enrichment Analysis
Raw246.7 mouse macrophage cells were infected with parasite lines for 24 h, then cells were lysed and total RNA was extracted using TRI reagent (Sigma). Library preparation was performed using a TruSeq stranded mRNA sample prep kit (Illumina, San Diego, CA, USA). Sequencing was performed on an Illumina HiSeq 2500 platform in a 75-base single-end mode. Illumina Casava1.8.2 software was used for base calling and raw sequence reads were subjected to quality control, then the cleaned reads were mapped to the reference mouse genome (mm10) with CLC Genomics Workbench version 10 (GWB; CLC bio, Aarhus, Denmark) (read mapping parameters: minimum fraction length of read overlap = 0.95 and minimum sequence similarity = 0.95). Only uniquely mapped reads were retained for further analysis. We identified differentially expressed genes (DEGs) as described in detail previously (19). The expression of each gene was compared among parasite lines using the differential expression for RNA-seq function in CLC GWB. DEGs were identified as genes with a fold change in expression of >2, and a max group mean of >1. KEGG pathway analysis was also conducted as described in detail in a previous article (19). The list of DEGs was subjected to a KEGG pathway enrichment analysis using the clusterProfiler package (20) in the statistical environment R to assess their overarching function. Following CPM normalization, the expression of each gene in the enriched pathways was normalized with Z-score normalization and visualized. Normalized gene expression was visualized in a heatmap using the heatmap.2 function (21) in the gplots package in R. The genes were hierarchically clustered based on the Pearson correlation distance and the group average method.

IFAT in T. gondii-Infected Cells
HFF cells in a 12-well plate were infected with parasites (multiplicity of infection = 1) for 24 h, along with uninfected control cells. The cells were then fixed with 4% (vol/vol) paraformaldehyde in PBS for 15 min at room temperature, permeabilized with 0.1% (vol/vol) Triton X-100 and blocked in PBS with 3% (wt/vol) bovine serum albumin. Cover slips were incubated with primary antibody for 1 h at room temperature, and fluorescent secondary antibody for 1 h at room temperature.  Nuclei were counterstained with Hoechst dye. Coverslips were then mounted onto the glass slide with Mowiol 4-88 (Sigma), and photographs were taken using All-in-One microscopy (BZ-9000, Keyence, Itasca, IL, USA). Quantification of the nuclear signal was performed by randomly selecting at least 20 infected cells per T. gondii strain and measuring the mean signal intensity per nucleus using the BZ analyzer II (Keyence).

IFAT in 293T Cells With Forced Expression of GRA Proteins
The 293T cells in a collagen 1-coated 12-well plate were transiently transfected with expression vectors of GRA7, GRA14, or GRA15, or the empty p3×FLAG-cmv14 vector as a negative (empty) control, using Fugene HD. After 24 h, IFAT and quantification of the nuclear signal were performed as described above.

Western Blotting
HFF cells were infected with parasites (multiplicity of infection = 3) for 24 h, then lysed using the LysoPure TM Nuclear and Cytoplasmic Extractor Kit (Wako, Osaka, Japan) supplemented with complete mini protease inhibitors and Phos stop (Roche, Mannheim, Germany). The cell lysates were separated by SDS-polyacrylamide gel electrophoresis and transferred to a Poly Vinylidene Di-Fluoride membrane (Millipore), which was blocked in TBS/0.1% Tween-20/2% ECL Prime Blocking Reagent (GE Healthcare, Buckinghamshire, UK) and incubated with primary and secondary antibodies. The protein bands were visualized by ECL Prime Western Blotting Detection reagent (GE Healthcare), and analyzed by Versa Doc with Quantity One (Bio-Rad, Munich, Germany). Band intensity was quantified using ImageJ software developed by the US National Institutes of Health.

Survival of Mice Infected With Toxoplasma gondii
Male C57BL/6J mice, of 8 weeks of age, were obtained from Clea Japan (Tokyo, Japan). Mice were infected intraperitoneally with 500 tachyzoites of the parental strain Pru, or mutant strains Pru gra7, Pru gra14, or Pru gra15. Mice were also infected intraperitoneally with 10,000 parental Pru or Pru gra14 parasites. To determine the survival rates to the type I RH strain, 500 tachyzoites the RH gra14, RH gra7, or their parental parasites were injected into the right footpads of mice and their survival was monitored for up to 30 days.

Statistical Analyses
Statistical analyses were performed using GraphPad Prism (version 6.0) software (GraphPad Software, San Diego, CA, USA). Statistically significant differences among groups were determined using one-way ANOVA with Tukey's post-hoc test. P-values of < 0.05 represent statistically significant differences. The survival rate was compared between groups using the logrank test.

Ectopic Expression of Type II GRA14 Activates NFκB Signaling in 293T Cells
To investigate which molecules modulate the NFκB pathway in Toxoplasma, we constructed mammalian expression vectors for 17 GRAs and 21 ROPs of a Toxoplasma type II strain. Then, we assessed whether their overexpression, together with luciferase reporter plasmids carrying an element dependent on the NFκB promoter, activated the reporter. Overexpression of GRA7, GRA14, and GRA15 activated NFκB ( Figure 1A).
Overexpression of GRA14 stimulated the NFκB promoter to a similar level as that of GRA7, whereas GRA15 produced much higher levels of NFκB-dependent luciferase activity than GRA7 and GRA14 ( Figure 1B). The expression of these molecules in 293T cells was confirmed by western blotting (Supplemental Figure 3). Thus, we focused on GRA7, GRA14, and GRA15 for further analysis. Next, we generated Pru gra7, Pru gra14, and Pru gra15 parasites based on the gene-editing strategies depicted in Supplemental Figure 2. We isolated single clones of drugresistant parasites and performed diagnostic PCR to check for correct integration (Supplemental Figure 1). Moreover, we established complementation of GRA7 and GRA14. The GRA7 and GRA14 expression cassettes containing the 5′ UTR and 3'UTR were inserted into the UPRT gene locus. Drug-resistant clones were isolated and correct integration into the UPRT locus was confirmed (Supplemental Figure 1). Clones, with the exception of Pru gra15, were further analyzed by an IFAT and western blotting to confirm the protein expression (Supplemental Figure 2). The Pru gra15 mutant was excluded because of the lack of an anti-GRA15 antibody. We then assessed the physiological changes in the transgenic lines in vitro. The infection rates and egress rates of the Pru gra7 and Pru gra14 strains in Vero cells were similar to those of the parental strain (Supplemental Figures 5A-D). Whereas, the in vitro replication

Primer
Sequence (     rate of Pru gra7 parasites was significantly higher than that of the parental parasites (Supplemental Figure 5E). The replication rate of Pru gra14 was comparable to that of the parental parasites (Supplemental Figure 5F). Next, we analyzed how each GRA contributes to NFκB activation because it is known that GRA15 plays a dominant role in NFκB activation by type II T. gondii. Cells infected with the Pru gra7 and Pru gra14 mutants showed a partial decrease in luciferase activity compared with cells infected with the parental Pru ( Figure 1C). However, for Pru gra15, NFκB activity was abolished in the infected cells ( Figure 1C).

Each GRA Expression Alone Is Sufficient to Activate NFκB in 293T Cells
We assessed whether each GRA protein alone is sufficient to activate the process of NFκB signal transduction. The level of nuclear RelA in GRA7-or GRA14-expressing cells was significantly higher than the level in control cells (Figure 2).
Moreover, the level of RelA nuclear translocation in cells expressed GRA15 was even higher than that in cells expressing GRA7 and GRA14 (Figure 2). First, we performed transient expression of each GRA gene. However, it is uncertain whether the function of ectopic single parasite molecule is the same as that of its native molecule. In addition, western blotting in the Supplemental Figure 3 indicated different expression levels among the transfection with GRA genes. Thus, we conducted similar experiments using deficient parasite strains.
Cells infected with the parental Pru strain revealed a higher level of nuclear RelA than cells infected with Pru gra7 or Pru gra14, whereas the complemented strain showed a similar level of RelA signal to the parental parasite (Figures 3A,B). Moreover, GRA15 deletion almost abolished the nuclear translocation of RelA ( Figure 3C). Representative images from these experiments are shown in Figure 3D. Next, we assessed whether each GRA affects the phosphorylation of IκBα by western blotting. The phosphorylated IκBα levels were comparable between cells infected with the parental, Pru gra14 and Pru gra7 strains, whereas GRA15 deficiency obviously reduced phosphorylated IκB (Figure 4). The relative levels of phosphorylated IκBα compared with the parental Pru-infected cells were 89, 105, and 38% in cells infected with Pru gra7, Pru gra14, and Pru gra15 strains, respectively (Figure 4).

Deficiency of GRA7, GRA14, and GRA15 Predominantly Results in Downregulation of Gene Expression Mediated by NFκB in Macrophages Infected With T. gondii
We next analyzed the levels of interleukin-6 (IL-6) in mouse macrophage Raw246.7 cells infected with parasites, and found that not only GRA15 deficiency but also GRA7 and GRA14 deficiency decreased the level of secreted IL-6 in the culture supernatant ( Figure 5). This result indicated that all of these GRAs affect the induction of the host immune response. To determine the host gene expression profiles relevant to these GRAs, we conducted transcriptome analysis of Raw246.7 cells infected with each strain and the uninfected cells. In total, 49, 103, and 338 genes were downregulated and 24, 15, and 111 genes were upregulated by GRA7, GRA14, and GRA15 deficiency, respectively ( Figure 6A, the complete sets of genes are listed in Supplemental Data Sheet 1). A Venn diagram was created to illustrate the similarities and differences among the genes regulated by these three GRAs (Figure 6A). This indicated that a number of common genes were regulated by these GRAs, and that GRA15 deficiency had more diverse effects than GRA7 and GRA14 deficiency.
To gain greater insight into the pathways regulated by each GRA in host cells, we conducted Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis on the relevant genes. This analysis primarily identified immune-response-related pathways, such as the cytokine-cytokine receptor interaction pathway, the IL-17 signaling pathway, and the tumor necrosis factor (TNF) signaling pathway, that were significantly enriched in the DEGs downregulated in macrophage cultures infected with deficient parasites compared with their expression in cell cultures infected by Pru and the complemented parasites (Supplemental Data Sheet 2). A heatmap of the gene expression associated with the cytokine-cytokine receptor interactions illustrated that GRA14, similar to GRA15, regulated some cytokines and chemokines ( Figure 6B). In addition, a heatmap of the gene expression associated with the IL-17 signaling pathway defined several cytokines and chemokines as GRA7-regulated genes ( Figure 6C, the complete sets of genes are available in Supplemental Data Sheet 3. To confirm the host genes whose expression is regulated by GRAs, we quantified the expression levels of several genes: IL-1β, IL-6, Cxcl1, Cxcl5, and Ccl17 for GRA14 and GRA15; and IL-6, Ccl7, Lcn2, Csf3, and Ptgs2 for GRA7. These genes were selected because their expression appeared to be regulated by GRAs according to the heatmaps (Supplemental Data Sheet 3). In most cases, the gene expression profiles were consistent between the transcriptome and the real-time PCR data (Supplemental Figure 6). Collectively, these results indicated that GRA7, GRA14, and GRA15 deficiency robustly downregulated the immune response-related pathways induced by T. gondii infection.
To further confirm the role of GRA14 in virulence, we infected mice by intraperitoneal injection with 10,000 tachyzoites of the parental Pru and Pru gra14 parasites, and monitored mouse survival until 30 days post-infection. There was no significant difference in survival between mice infected with 10,000 parasites of the parental Pru and Pru gra14 strains ( Figure 7D).
Next, to determine the effect of GRA14 on in vivo parasite growth and the immune response at the site of infection, mice were infected with 500 tachyzoites of the parental Pru, Pru gra14, and GRA14-complemented lines. At 5 days after infection, mice were euthanized and the parasite burden and levels of cytokine secretion, including IL-12p40 and interferon-γ (IFN-γ), were examined. Mice infected with each strain showed no significant difference in parasite burden in the spleen or the peritoneal exudate cells (Supplemental Figures 7A,B). Although the differences in IL-12p40 and IFN-γ secretion by the peritoneal exudate cells were not significant, the average level of IFN-γ in Pru gra14 mutantinfected mice was higher than that in the parental Pru and complemented strains (Supplemental Figures 7C,D). Moreover, no significant difference was detected in the level of serum IFN-γ on either day 3 or day 5 among these mouse groups (Supplemental Figures 7E,F).
To investigate how the deficiency of each GRA affects host immunity at an earlier time, we conducted a time-course experiment using thioglycolate-induced peritoneal macrophages (Supplemental Figure 8). Supernatants were collected every 6 h for 24 h and measured the production of IL-12p40. IL-12p40 production was abolished in the macrophages infected with Pru gra7, Pru gra14, and Pru gra15 strains until 18 hours post-infection. However, the IL-12p40 production at 24 h post-infection decreased in the macrophages infected with the deficient parasite lines, with the highest decrease from Pru gra15, followed by Pru gra7, and then Pru gra14.
Lastly, we examined the effect of GRA7 and GRA14 on type I RH parasites. Mice were infected by intra-footpad injection with 500 parasites of the RH gra7, RH gra14, and parental RH parasite strains, respectively. Survival was monitored for 30 days and 86% (13/15) of mice infected with the parental RH strain survived, whereas all RH gra14 mutant-infected mice succumbed to the infection between 15 and 26 days after infection ( Figure 7E). By contrast, when challenged with the parental RH and RH gra7 parasites, 0% (0/8) and 25% (2/8) of mice survived after infection with the parental RH and RH gra7 parasites, respectively ( Figure 7F).

DISCUSSION
Secreted GRA15 has been identified as a major factor that contributes to the strain-specific differences in NFκB activation (15). Meanwhile, GRA7 produces a strong antibody response in the acute phase of infection (22) and has been tested as a candidate for vaccine development (23). Recent studies have revealed that GRA7 associates with ROP2 and ROP4, and functions in concert with ROP18 protein complexes that resist IFN-γ-activated host immune-related GTPase (24)(25)(26). Moreover, recombinant GRA7 interacts with inflammasomerelated molecules, such as an apoptosis-associated speck-like protein that contains a caspase recruitment domain (ASC) and phospholipase D1 (PLD1) (27). However, few studies have investigated the role of GRA7 in the pathogenesis of type II T. gondii strains. In the present study, we demonstrated that GRA14 is involved with NFκB activation by T. gondii. GRA14 seems to be implicated in the interaction with host molecules because secreted GRA14 localizes to PVs containing membranous strandlike extensions (called PVM extensions) similar to other GRA proteins such as GRA3 and GRA7 (17). Furthermore, GRA14 is anchored in the PVM with its C terminus facing the host cell cytosol (17). GRA14 has also been reported as a potential vaccine candidate against T. gondii infection. Several studies have reported the protective immunity induced by vaccination with GRA14 antigen (28)(29)(30)(31)(32). However, there have been no previous reports regarding the modification of host cell function by GRA14. Thus, we targeted GRA7, GRA14, and GRA15 in this study.
Although we focused on NFκB signaling pathway, reporter activity by GRA was also evaluated in this study using reporter  Heatmap showing that the most significantly enriched pathway associated with the GRA14 expression status was the cytokine-cytokine receptor interaction, which contained a subset of 41 genes. Rows represent samples from different processes, and columns represent genes. (C) Heatmap showing that the most significantly enriched pathway associated with the GRA7 expression status was the IL17 signaling pathway, which contained a subset of 91 genes. Rows represent samples from different processes, and columns represent genes. The cluster shown in (B,C), which was upregulated when GRA14 or GRA7 expression was restored, is enlarged on the right. The experiment for GRA14 was performed in tandem with GRA15, and the experiment for GRA7 was performed independently.
plasmids having response elements such as cAMP-responsive element, nuclear factor of activated T cells (NFAT), serum responsive element, serum responsive factor (SRF), and activated protein 1. As shown in Supplemental Figure 4, GRA14 and GRA15 activated all of them, while GRA7 activated NFAT and SRF other than NFκB. However, the main activities of GRAs were observed in NFκB activation. Interacting host factor of GRA14 is unknown, while GRA7 and GRA15 activate NFκB via TNF receptor-associated protein (TRAF). TRAF participates in the activation of the transcription factor NFκB and members of the mitogen-activated protein kinase (MAPK) family, including MAPK, c-jun N-terminal kinase, and p38. It remains possible that each GRA regulates host immunity via signaling pathway other than NFκB, but we believe that one of the primary sites of action is the NFκB pathway. Interestingly, although the levels of nuclear translocation of RelA in GRA7-and GRA14-expressing cells were significantly lower than in GRA15expressing cells, expression alone was adequate for nuclear translocation. Moreover, cells infected with Pru gra7 parasites showed no significant difference in the intensity of nuclear translocation compared with uninfected cells. In addition, GRA14 deficiency partially attenuated the intensity of nuclear RelA in cells infected with T. gondii. Collectively, these results suggest that GRA15 is the main player for NFκB activation by type II T. gondii. Additionally, GRA7 and GRA14 play a certain role in modulation of the NFκB pathway by type II T. gondii. By contrast, the levels of phosphorylated IκBα were comparable among cells infected with the parental Pru strain and mutant strains Pru gra7 and Pru gra14. It was reported that GRA15-mediated NFκB activation was dependent on TRAF6, and GRA15 deficiency caused a decrease in the levels of phosphorylated-IκBα (15), which was consistent with our results. Contrary to this, another study showed that recombinant GRA7 also interacted with TRAF6, and recombinant GRA7 protein stimulated the phosphorylation of IκBα (16). However, in the present study, GRA7 deficiency showed no clear change in the phosphorylation level of IκBα. It may be that due to the higher activity of GRA15 compared with that of GRA7 and GRA14, GRA15 compensates for the loss of GRA7 and GRA14 function.
Pru gra7, Pru gra14, and Pru gra15 strains induced significantly less cytokine secretion from infected macrophages than the parental Pru strain-infected cells. NFκB activation leads to the transcription of pro-inflammatory genes, such as those encoding IL-1β and IL-12 (15,33). In addition, our transcriptome analysis revealed that these GRAs regulated the gene expression levels of similar inflammatory cytokines and chemokines by macrophages, in turn stimulating the development of a T-helper type 1 (Th1) immune response (33). Our data suggested that either GRA7 or GRA15 deficiency is sufficient for the increase in acute virulence in infected mice. Mice infected with a type II GRA15-deficient strain had a significantly higher parasite burden than mice infected with a parental type II strain (15). GRA15 activates NFκB in host cells and induces early IL-12 secretion (15). IL-12 stimulates NK cells and T cells to secrete IFN-γ (34). On day 2 after infection, mice infected with a type II GRA15-deficient strain had significantly less IFN-γ in their intraperitoneal cavities than mice infected with a parental type II strain (15). IFN-γ is the primary cytokine of host resistance to intracellular pathogens (35). Thus, this difference in IFN-γ levels was the likely cause of the virulence differences. It has been demonstrated that GRA7 interacts with TRAF6, inducing innate immune responses via the NFκB pathway in macrophages (16). Our results suggested that the GRA7-induced reporter activity of the NFκB promotor was less than that of GRA15. However, GRA7 also interacts with a number of host cell proteins, including ASC and PLD1, revealing a new facet of the role of GRA7 in the regulation of innate immune responses (36). Thus, GRA7 deficiency might result in increased mortality comparable to that of GRA15.
GRA14 deficiency also resulted in a slight but significant increase in virulence compared with the parental strain in mice after the injection of 500 parasites. Whereas, consistent with recent research involving 2 × 10 5 parental type II gra14 parasites, such a difference was no longer detectable when 10,000 parasites were injected, which furthermore had the potential to cause lethal tissue damage (37). However, after 5 days of intraperitoneal infection, GRA14 did not affect the parasite burden or the level of cytokine secretion, including IL-12p40 and IFN-γ, from the peritoneal cavity. Moreover, no significant difference was detected in the levels of serum IFN-γ on days 3 or 5 among the groups of mice. The attenuated signal output caused by GRA14-deficiency may impair the proper immune response, resulting in an increased parasite burden in mice infected with Pru gra14 parasites at an early stage (days 1-4), explaining the slight difference in virulence of this strain. The GRA-induced protective immune response against T. gondii in mice requires activation of antigen-presenting cells such as IL-12 production in the early stages of infection. If the parasites were controlled by the protective immune response in the early stage of infection, the level of the inflammatory marker IFN-γ would be suppressed. Therefore, the increased activity of Pru gra14 at the initial stage of infection might increase IFN-γ level compared to the parental and complemented lines. Overall, our results suggest that GRA7 and GRA15 are the major contributors to in vivo virulence, whereas GRA14 has a relatively low impact on mice virulence. Furthermore, parasites deficient in GRA7 but not in GRA14 affect parasite growth in vitro. Moreover, a previous study reported that a type II gra15 mutant formed significantly larger plaques than a type II strain in HFF cells, but this was not apparent in mouse embryo fibroblast cells (15). These data indicate that growth differences in GRA7 and GRA15-deficient strains may affect their virulence in mice.
In this study, our experiments had focused on type II strains; however, we conjectured that the GRA14 proteins of type I strains are functional because there are few amino acid differences between the type I and II proteins (P43S, D323G, and S356V). The GRA15 proteins from type II and type III strains activate NFκB. Type II strains activate NFκB more strongly than type III strains, whereas the type I RH strain does not induce NFκB activation because it has a mutation in GRA15, leading to a frameshift and an early stop codon (15,38). Therefore, because GRA15 of type I strains lacks activity, it is easy to evaluate the effect of GRA14 deficiency. Thus, we hypothesized that GRA14 might be involved with the mechanisms of NFκB activation by type I T. gondii. Previous studies have shown that type I strains interfere with the host NFκB pathway to promote their survival. ROP18, a key serine/threonine kinase that phosphorylates host proteins to modulate acute virulence, is associated with phosphorylation of RelA at Ser-468 and promotes the degradation of RelA to inhibit the NFκB pathway (39). Moreover, polymorphic kinase ROP16 of type I strains is capable of suppressing the IL-12 response of infected macrophages stimulated with lipopolysaccharide, thereby inhibiting NFκB transcriptional activity (15,40). Whereas, other studies have shown that NFκB is activated by a type I strain of T. gondii, and that its activation is necessary for the inhibition of apoptosis (41)(42)(43). However, it is not known what effect GRAs have on the NFκB pathway.
Therefore, we lastly evaluated the effects of GRA7 and GRA14 deficiency in the type I RH strain on the survival of mice. Surprisingly, unlike type II parasites, all mice infected with RH gra14 parasites died within 26 days of footpad inoculation. Previous studies showed that GRA14 did not affect the growth and virulence of parasites following intraperitoneal injection of mice (17,44). Unlike intraperitoneal inoculation, which results in a rapid, acute systemic infection, intra-footpad inoculation allows us to observe the gradual spread of T. gondii in vivo (45). Generally, intraperitoneal infection by RH tachyzoites was lethal. However, intra-footpad infection led to survival or, at least, a prolonged survival time in the present study. Therefore, deleting GRA14 may result in a lethal parasitic load in mice. By contrast, mice infected with RH gra7 parasites experienced a significant delay in death compared with the parental RH strain. A previous study reported that outbred CD-1 mice infected with RH gra7 parasites exhibited a similar phenotype (25). GRA7 binds to the GTP-bound immunity-related GTPase a6 and acts synergistically with ROP18 to block immunity-related GTPases (25,26). Thus, these results suggest that GRA14 plays an important role in the control of parasite infection, creating a paradigm that protects the host animals from acute infection and death.
In conclusion, the present study demonstrated new molecular functions for GRA7 and GRA14 and confirmed their role in the induction of NFκB during a type II strain infection. NFκB activation mediated via GRA7, GRA14, and GRA15 was closely related to the Th1 response promoted by inflammatory cytokines following the activation of macrophages. This immune response limits the tissue invasion of the parasite, ensuring the survival of the host but, paradoxically, also aiding the survival of the parasite by converting it into a bradyzoite form able to persist in the muscle and brain tissues (46). GRA7 has multiple target components within the host cell that cause different virulence phenotypes dependent on the type of parasite. Whereas, the GRA14 protein has a low polymorphic phenotype and is potentially functional throughout type I, II, and III strains. Moreover, the suppressive control of virulence by early immune activation after infection, which has been regarded as a unique event to type II strains, is a conserved strategy across parasite strains. This may contribute to the high prevalence and wide distribution of this protozoan parasite. Thus, further insight into the precise role of these GRAs may help delineate the mechanism of NFκB modulation by T. gondii.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are publicly available. This data can be found in NCBI SRA: https://www.ncbi.nlm.nih.gov/sra/?term=DRA010408, accession number DRA010408.

ETHICS STATEMENT
The experimental protocol was approved by the Committee on the Ethics of Animal Experiments at the Obihiro University of Agriculture and Veterinary Medicine (permit number: 19-50). data analyses and wrote the manuscript. All authors revised the manuscript and approved the final version.
FUNDING would like to thank Kate Fox, DPhil, from Edanz Group (www. edanzediting.com/ac) and Rochelle Haidee D. Ybañez for editing a draft of this manuscript.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu. 2020.01709/full#supplementary-material Supplemental Data Sheet 1 | Detailed expression data for GRA-dependently regulated genes.
Supplemental Data Sheet 2 | Results of KEGG pathway analysis.

Supplemental Data Sheet 3 | Heat map.
Supplemental Figure 1 | GRA knockout and the complementation strategies. (A) Schematic genomic representation of the GRA7 locus and the plasmid construct used to target the GRA7 gene. The drug-resistant hypoxanthine-xanthine-guanine phosphoribosyl transferase (HXGPRT) cassette was surrounded by the 5 ′ and 3 ′ untranslated regions (UTR) of GRA7. (B) Schematic representation of the CRISPR/CAS9 strategy used to inactivate the target genes by inserting the pyrimethamine-resistance DHFR cassette (DHFR * ). Transfection of the CRISPR plasmid targeting TgGRA (G01), together with an amplicon containing the DHFR * -expressing cassette flanked by regions homologous to the target gene, was used to disrupt the corresponding target gene by insertion. (C) Transfection of the CRISPR plasmid targeting the TgUPRT gene, together with an amplicon containing 1,000 by of the 5 ′ and 3 ′ untranslated regions (5 ′ UTR and 3 ′ UTR) flanked by regions homologous to the TgUPRT gene, was used.
Supplemental Figure 2 | IFAT and western blotting to confirm the expression of TgGRA7 and TgGRA14. (A) IFAT analysis of Vero cells infected with Pm, PruAgra7 (deltaGRA7), PruAgra14 (deltaGRA14), and complemented (CompGRA7 and CompGRA14) parasites at 24 h post-infection. Cells were fixed and stained with a-TgSAG1 (green), a-TgGRA7 (red), a-TgGRA14 (red), and Hoechst dye (blue). Bars, 10 gm. (B) Western blots of the parasite strains. Into each lane, 1 × 10 6 parasites were loaded. Anti-TgGRA7 and anti-GRA14 antibodies detected 25.9 and 42 kDa proteins in the parental Pru and complemented parasites, respectively, but not in the deficient mutant parasites.
Supplemental Figure 3 | Forced expression of GRA7, GRA14, and GRA15 in 293T cells. 293T cells were transiently transfected with the expression vectors for GRA7, GRA14, and GRA15. Cells lysates were then separated by SDS-PAGE, and western blot analysis was carried out using an anti-FLAG antibody. The estimated sizes of the FLAG-tag fused to GRA7, GRA14, and GRA15 were 29.9, 48.7, and 61.8 kDa, respectively. However, the observed molecular weight of the FLAG-tag fused to GRA15 (−75 kDa), was higher than the expected predicted size (61.8 kDa). It has been shown that this is not caused by parasite-mediated modification of GRA15. Most likely, it is the particular amino acid composition of GRA15, which is enriched in Pro, Ser, and Thr, that makes it run slower than expected on an SDS-PAGE gel. Black arrows indicate the estimated band sizes of the target proteins. Supplemental Figure 6 | Expression levels of chemokines and cytokines in Raw246.7 macrophage cells. (A,B) Raw246.7 macrophage cells were infected with parasite strains for 24 h, then cells were lysed, and total RNA was extracted. RNA was used to synthesize cDNA. Real-time RT-qPCR amplification was carried out for CXCL1, CXCL5, IL-lbeta, IL-6, and Cc117. Each bar represents the mean ± the standard deviation (n = 3 for all groups), and the results are from a single experiment. Statistical analysis was performed using one-way ANOVA with post-hoc Tukey's test, * a significant difference (p < 0.05). Figure 7 | In vivo cytokine ELISA and parasite burden. Mice were infected with 500 tachyzoites of the parasites and then euthanized to examine the parasite burden and level of cytokine secretion, including IL-12p40 and IFN-y, 5 days after infection. (A,B) Parasite burden in the spleen and peritoneal cavity in mice infected with parental Pru, PruLgra14 (deltaGRA14), and the complemented parasites (Comp). (C,D) Levels of serum IFN-y and IL-12p40 in mice. (E,F) Cytokine levels in the peritoneal cavity in mice. Each plot represents data from one mouse. Statistically significant differences were analyzed by one-way ANOVA with post-hoc Tukey's test but no significant difference was found. Data were collected from one experiment.

Supplemental
Supplemental Figure 8 | Levels of interleukin-12p40 in thioglycolate-elicited macrophage cells. Macrophages were infected with the parental Pru, PruAgra7, PruAgra14, and PruAgral5 parasite strains. At every 6 h for 24 h post-infection, supernatants were collected, and IL-12p40 levels were determined by cytokine ELISA. Values are the means ± SD of four samples, * a significant difference (p < 0.05). #a significantly lower level of IL-12p40 compared with the Pm strain infected cells (p < 0.05). Differences were tested by one-way ANOVA with Turkey's post-hoc test.