Abstract
The obligate intracellular parasite, Toxoplasma gondii, is highly prevalent among livestock species. Although cattle are generally resistant to Toxoplasma strains circulating in Europe and North America, the underlying mechanisms are largely unknown. Here, we report that bovine bone marrow-derived macrophage (BMDM) pre-stimulated with interferon gamma (IFNγ) restricts intracellular Toxoplasma growth independently of nitric oxide. While Toxoplasma promoted the expression of genes associated with alternative macrophage activation and lipid metabolism, IFNγ abrogated parasite-induced transcriptional responses and promoted the expression of genes linked to the classical macrophage activation phenotype. Additionally, several chemokines, including CCL22, that are linked to parasite-induced activation of the Wnt/β-catenin signaling were highly expressed in Toxoplasma-exposed naïve BMDMs. A chemical Wnt/β-catenin signaling pathway antagonist (IWR-1-endo) significantly reduced intracellular parasite burden in naïve BMDMs, suggesting that Toxoplasma activates this pathway to evade bovine macrophage anti-parasitic responses. Congruently, intracellular burden of a mutant Toxoplasma strain (RHΔASP5) that does not secrete dense granule proteins into the host cell, which is an essential requirement for parasite-induced activation of the Wnt/β-catenin pathway, was significantly reduced in naïve BMDMs. However, both the Wnt/β-catenin antagonist and RHASPΔ5 did not abolish parasite burden differences in naïve and IFNγ-stimulated BMDMs. Finally, we observed that parasites infecting IFNγ-stimulated BMDMs largely express genes associated with the slow dividing bradyzoite stage. Overall, this study provides novel insights into bovine macrophage transcriptional response to Toxoplasma. It establishes a foundation for a mechanistic analysis IFNγ-induced bovine anti-Toxoplasma responses and the counteracting Toxoplasma survival strategies.
Introduction
Toxoplasma is a zoonotic protozoan parasite that infects virtually all warm-blooded vertebrates and is perhaps the most successful and widespread human pathogen. It is the leading cause of encephalitis and death in HIV/AIDS patients (Basavaraju, ) and is ranked 4th among foodborne parasites with the greatest global impact (FAO/WHO, ) and contributors to years lived with disability and disability-adjusted life years per million persons (Gkogka et al., ). Toxoplasma undergoes sexual reproduction exclusively in cats while asexual reproduction can occur in any intermediate host, including humans and farmed animals. Definitive-to-intermediate host parasite transmission occurs via ingestion of oocysts from infected cat feces, while transmission between intermediate hosts can occur through the consumption of contaminated food products or vertically from mother to fetus (Hill and Dubey, ). Although Toxoplasma infections in healthy individuals are mostly asymptomatic, severe disease or even death, often caused by damage to the brain or other organs, can occur in immunocompromised or congenitally infected individuals (Hill and Dubey, ). Toxoplasma typically establishes lifelong chronic infection in healthy intermediate hosts by encysting, after an initial phase of rapid intracellular proliferation and cell–cell spread, in the central nervous system and muscle tissues (Hill and Dubey, ). To do this, the parasite must carefully regulate immune activation and host cell anti-Toxoplasma effector mechanisms.
Toxoplasma invades host cells mostly through an active process powered by parasite-generated actin motor activity to form a specialized non-fusogenic parasitophorous vacuole (PV) (Fleckenstein et al., ) that helps the parasite avoid host cell immune surveillance (Morisaki et al., ). In phagocytic cells, such as macrophages, the parasite can also be taken up by phagocytosis and, sometimes the phagocytosed parasites can exit the phagosome to form a PV (Zhao et al., 2014). Besides hiding within the PV, the parasite also evades host cell immunity by sequentially discharging several effector proteins that modulate a variety of host immune and metabolic processes, including the inflammatory pathway (Hunter and Sibley, ). Generally, effective host anti-Toxoplasma responses are dependent on the production of interleukin (IL)-12 (Gazzinelli et al., ) by macrophages and dendritic cells (Gazzinelli et al., ). IL-12 in turn activates natural killer (NK) and T cells to secrete interferon gamma (IFNγ) (Gazzinelli et al., , ), a pro-inflammatory cytokine that activates several anti-Toxoplasma effector mechanisms such as the interferon-regulated GTPases (IRGs) in mice (Zhao et al., 2009), reactive nitrogen/oxygen intermediates (Scharton-Kersten et al., ), tryptophan degradation and cell death in human cells (Pfefferkorn, ; Niedelman et al., ), and inflammasome activation (Cirelli et al., ). In return, Toxoplasma has evolved several mechanisms to counteract the IFNγ-induced host defenses, including the secretion of effectors proteins from specialized apical organelles that co-opt host transcription and signaling pathways to control host cell responsiveness to inflammatory signals (Jensen et al., ; Koshy et al., ; Bougdour et al., ; He et al., ).
Besides producing the IL12 that primes IFNγ secretion by NK and T cells, macrophages are also the preferred intracellular niche for the fast-dividing parasite stage (Jensen et al., ). Thus, the innate defenses triggered by monocytes/macrophages are key to toxoplasmosis pathogenesis in humans and rodents (Channon et al., ; Dunay and Sibley, ; Gregg et al., ; Tosh et al., 2016; Song et al., ). Besides killing microbes, including Toxoplasma, macrophages can also initiate adaptive immune responses (van de Vosse et al., 2009; Thi et al., 2012). Human and mouse models show that upon the engagement of surface signaling receptors or pattern recognition receptors (PRRs) such as, toll-like (TLRs), RIG-I-like (RLRs), and the cytosolic NOD-like (NLRs) receptors by conserved pathogen-associated molecular patterns (PAMPs) such as, lipopolysaccharide (LPS) or immune factors, including cytokines, macrophages assume different activation phenotypes. The most extreme classical [M1, M(IFNγ)] and the alternative [M2, M (IL-4)] phenotypes are separated by several intermediate activation states (Murray et al., ). The M1 phenotype, which can be induced by IFNγ, is highly microbicidal and characterized by the production of reactive oxygen and nitrogen intermediates such as nitric oxide (NO), and a range of pro-inflammatory cytokines and chemokines, such as Tumor necrosis factor (TNF) alpha (De Paoli et al., ; Murray et al., ). In contrast, the M2 phenotype, which is induced by IL-4 and IL-13 and is important for regulating inflammation, is characterized by the production of anti-inflammatory cytokines and growth factors (Sindrilaru and Scharffetter-Kochanek, ). The general hypothesis is that macrophage activation phenotypes, which are underpinned by discrete transcriptional programs (Hassan et al., ), provide a high degree of plasticity that is exploited by some intracellular pathogens, including Toxoplasma and Mycobacteria, to turn this potentially hostile host cell into a favorable replication niche (Price and Vance, ). Although effective host response to Toxoplasma require the induction of inflammation, characterized by classical macrophage activation, Toxoplasma promotes its survival in macrophages by secreting effector proteins to dampen inflammatory responses and favor alternative macrophage activation (Jensen et al., ). Indeed, Toxoplasma strain differences in virulence in mice is partly due to strain differences in inducing alternative macrophage activation (Jensen et al., ), while host differences in susceptibility to Toxoplasma is due in part to differences in macrophage activation phenotypes after infection (Jensen et al., , ). Therefore, to devise strategies to improve the early defense against Toxoplasma and a variety of intracellular pathogens, it is important to understand the molecular mechanisms underpinning macrophage response to Toxoplasma and/or immune effector proteins, such as cytokines.
Toxoplasma is common in many species of livestock, including cattle. Compared to other livestock species, such as sheep and pigs, cattle are highly resistant to Toxoplasma and rarely transmit the parasite to other intermediate hosts (Dubey, ; Esteban-Redondo and Innes, ). Unlike sheep and pigs, natural Toxoplasma infection in cattle is mostly asymptomatic and does not appear to result in abortion. However, the molecular factors and mechanisms that modulate bovine-Toxoplasma interactions, which can be exploited to enhance resistance in other ruminants, are ambiguous. Although IFNγ is central to anti-Toxoplasma responses in virtually all vertebrates, the role of this cytokine, and the mechanism underpinning, bovine resistance to Toxoplasma is equivocal. In the present study we sought to determine the role of IFNγ, and host genes associated with effective bovine macrophage response to Toxoplasma. We performed RNA-sequencing on naïve or IFNγ-stimulated bovine bone marrow derived macrophages (BMDMs) that were either unexposed or exposed to a Toxoplasma strain that is highly virulent in laboratory inbred mice (RH) for 24 h. Analysis of the datasets provides novel insights into the Toxoplasma-induced transcriptional responses in naïve and pre-stimulated BMDMs. We report that IFNγ enhances bovine BMDMs anti-Toxoplasma responses and that, despite producing large amounts of nitric oxide, bovine macrophages restrict Toxoplasma independently of nitric oxide. Toxoplasma induces the Wnt/β-catenin signaling pathway and the expression of several anti-inflammatory chemokines and arginine metabolism in naïve BMDMs. On the other hand, to survive in IFNγ-primed BMDMs, the parasite expresses mostly genes that are associated with its slow dividing bradyzoite stage.
Materials and Methods
Parasites
The type I Toxoplasma strain (RH) engineered to express green fluorescent protein (GFP) and firefly luciferase has previously been described (Jensen et al., ). The RHΔASP5 (Hammoudi et al., ) was a generous gift from Dr. Mohamed-Ali Hakimi (INSERM). All parasite strains were maintained by serial passage on confluent human foreskin fibroblast (HFF) monolayer.
Primary Bone Marrow Derived Macrophages (BMDMs)
Marrow cells were flushed from the ribs of three (n = 3) ~2 year old calves using phosphate buffered saline (PBS, Invitrogen). The cells were centrifuged at 500 × g for 5 min at 4°C and re-suspended in red cell lysis buffer (Sigma) and incubated on ice for 5 min. Next, the cells were passed through a 70 μm cell strainer (BD Biosciences) and centrifuged at 500 × g for 5 min at 4°C. Cells from each calf were subsequently differentiated into macrophages in 10 cm non-tissue culture petri dishes (Corning) in RPMI 1640 (Sigma-Aldrich) supplemented with heat-inactivated 20% fetal bovine serum (FBS, Thermo Fisher Scientific), GlutaMAX (Thermo Fisher Scientific), penicillin/streptomycin (Thermo Fisher Scientific), and recombinant human CSF1 (104 U/ml; a gift from Chiron, Emeryville, CA) for 10 days as previously described (Young et al., 2018).
In vitro Measurements
BMDMs from each calf were detached using a cell scraper, washed, counted, and seeded separately in triplicates in 96-well plates at a density of 105 cells/well. The BMDMs were left unstimulated (naïve) or stimulated with: recombinant bovine IFNγ (100 ng/mL), lipopolysaccharide (LPS) from Salmonella enterica serotype Minnesota Re 595 (Sigma-Aldrich, 100 ng/mL), or a combination of IFNγ (100 ng/mL) and LPS (100 ng/mL) with or without aminoguanidine (final concentration of 500 μM) and incubated at 37°C in 5% CO2. 24 h post-stimulation, cell free supernatants were collected for Griess reagent-based nitric oxide assay, as previously described (Young et al., 2018), and cell viability, using the CellTiter 96® Aqueous One Solution Cell Proliferation kit (Promega) according to the manufacturer's recommendations. Freshly lysed parasites (by sequential passage through 25G and 27G needles) were passed through a 5 μm filter to remove HFFs, counted, diluted in RPMI, and added to the BMDMs at a multiplicity of infection (MOI) 1 for 24 h before luciferase activity was measured using a luciferase assay kit (Promega) according to the manufacturer recommendations.
RNA Sequencing and Analysis
2 × 106 BMDMs from each calf were seeded in 6-well plates and left unstimulated or stimulated with IFNγ (100 ng/mL) for ~18 h at 37°C in 5% CO2. Next, the media was replaced with fresh media containing Toxoplasma at a MOI 1 and incubated further for ~18 h. Fresh media was added to the non-infected control BMDMs. At the end of the incubation period, total RNA was isolated from each BMDM sample using QIAzol® Lysis Reagent and miRNeasy Mini Kit (Qiagen) according to the manufacturer recommendations. RNA quality and integrity were assessed on the Agilent 2200 TapeStation using an Agilent RNA ScreenTape and quantified using a Qubit RNA BR Assay Kit and Fluorometer. RNA-seq libraries were generated and sequenced by Edinburgh Genomics. All libraries were prepared using the Illumina TruSeq Stranded library protocol for total RNA libraries (Part: 15031048, Revision E). Briefly, polyA-tail enrichment (Dynabeads mRNA Purification Kit; Invitrogen) was performed on the total RNA and the mRNA fragmented into 200–400 base-pairs, and reverse transcribed into cDNA before Illumina sequencing adapters were added to each end. Twelve barcoded libraries were multiplexed and sequenced on a single S2 sequencing lane on the Illumina NovaSeq 6000 machine to yield ~60 million 50 bp high quality strand-specific paired reads per sample. Reads were pseudoaligned to the Ensembl bovine genome (ARS-UCD1.2) using Kallisto v.0.46.0 (Bray et al., ) with 100 bootstraps to generate transcript-level expression estimates as transcripts per million (TPM) as previously described (Young et al., 2018). Gene-level differential expression analysis was performed in sleuth as previously described (Pimentel et al., ). TPM values derived from pseudoalignment of RNA-seq reads to the GT1 genome (ToxoDB.org) were used as input in sleuth to identify differentially expressed Toxoplasma genes.
Quantitative Real-Time PCR
One microgram of total RNA from each BMDMs sample was reverse transcribed using Superscript III (Thermo Fisher Scientific), 10 mM dNTP mix (Thermo Fisher Scientific) and Oligo (dT)12−18 Primer (Thermo Fisher Scientific) according to the manufacturer's instructions. Gene-specific primers (Table 1) were designed using Primer-BLAST (Ye et al., 2012). All primers were commercially synthesized by Thermo Fischer Scientific. Quantitative real-time RT-PCR (qPCR) was performed in a Rotor-Gene Q real-time PCR cycler (Qiagen) using FastStart Universal SYBR Green Master (Rox) in final volume of 20 μL. The linearity and efficiency of qPCR amplification was determined for each primer pair using a standard curve generated by a serial dilution of cDNA pooled from all the samples. All reactions were performed in duplicate and “no template” controls included for each gene. Agarose gel electrophoresis was used to confirm product sizes and melt curve analysis confirmed specificity of amplification. Data were analyzed using the 2−ΔΔCt method and statistical analyzed performed on the ΔCt values. qPCR was performed on parasite DNA as previously described (Li et al., ).
Table 1
| Symbol | Primer 1 forward | Primer 1 reverse | Amplicon length |
|---|---|---|---|
| C1R | GGTGCAGGATCAAGGACTGC | GTGTGCATCTTGTAGAAGGGCT | 136 |
| ISG15 | GACCTGACGGTGAAGATGCTA | ATCTTCTGGGCGATGAACTGC | 98 |
| PSMB8 | TGTCAATATGTACCACATGAAGGAG | CACCATCACTGACTGGCCTC | 102 |
| SGK1 | GCCAAGGATGACTTTATGGAGA | AGGATCAAAGTGTCGCAGG | 138 |
| IRF1 | AGGACATCATGAAGCTCTTGGA | GCTCCTCCTTGCAACTGAACT | 129 |
| GBP5 | CCAGGAAAGGAATACAGGCTGA | TTCCATTGCTGTGAGAGCCAG | 107 |
| RSAD2 | GTGGTTCCAGAAGTACGGTGA | CTTCTTTCCTTGACCACGGC | 103 |
| CFB | CTTGCAAAGGTGATTCTGGTGG | CGCTTGCAAACATCCACGAC | 100 |
| GBP1 | CTCTCAAACTGCAGGAACAGTC | TGCTTTGGATAAGAGTGACCAG | 175 |
| CXCL10 | TCCTCGAACACGGAAAGAGGCATA | AGCTGATATGGTGACTGGCTTGGT | 164 |
| CCL22 | CGGGACTACATCCGTTACCC | CAGCACAGATCTCTCGGTCC | 121 |
| CCl24 | GCAGGAGTGATCTTCACCACC | TAGCGGAGGCTTTCTTCTGC | 115 |
| SDHA | ACCTGATGCTTTGTGCTCTGC | CCTGGATGGGCTTGGAGTAA | 126 |
| GAPDH | GGTGATGCTGGTGCTGAGTA | TCATAAGTCCCTCCACGATG | 265 |
Primer sets for qPCR validation of representative differentially expressed genes.
Results
IFNγ Induces Restriction of Toxoplasma Growth in Bovine Macrophages Independently of Nitric Oxide
To determine whether IFNγ is essential for the control of Toxoplasma growth in bovine macrophages, we exposed naïve or IFNγ-stimulated bovine bone marrow-derived macrophages (BMDMs) to a luciferase-expressing type I Toxoplasma strain (RH) for 24 h and measured luciferase activity, a proxy for parasite burden (Hassan et al., ). In murine macrophages, IFNγ is known to require a second stimulant, such as LPS or TNF, to effectively restrict Toxoplasma (Sibley et al., ; Hassan et al., ). Therefore, we included in our experiments, BMDMs that were pre-stimulated with IFNγ+LPS. As expected, there was reduced intracellular parasite burden in BMDMs pre-stimulated with IFNγ or IFNγ+LPS, relative to naïve BMDMs (Figure 1A). However, there were no significant differences in parasite burden between IFNγ- and IFNγ+LPS-stimulated BMDMs, suggesting that IFNγ alone can sufficiently induce anti-Toxoplasma activities in bovine BMDMs. A quantitative real-time polymerase chain reaction (qRT-PCR) analysis of parasite DNA (Li et al., ) confirmed the intracellular parasite burden differences between naïve and IFNγ-stimulated BMDMs (Figure 1B).
Figure 1
Previously, others and we reported that Nitric oxide (NO), a product of the inducible nitric oxide synthase (iNOS or NOS2)-catalyzed L-arginine metabolism (MacMicking et al.,
Toxoplasma Induces a Robust Bovine Macrophage Transcriptional Response
To gain mechanistic insight into bovine innate immune response to Toxoplasma, we leveraged high throughput RNA-sequencing to profile the transcriptional landscape of bovine BMDMs that were left unstimulated (naïve) or pre-stimulated with IFNγ before being exposed to Toxoplasma for 24 h. Transcriptional analysis was also performed on uninfected controls. Downstream analysis was restricted to genes that were differentially expressed by more than 2-fold change and had at least 10 reads that uniquely aligned to the genome in at least two samples when compared the uninfected naïve BMDMs. In total, 1,349 unique genes were differential expressed in at least one condition, of which 887 were differentially expressed in a condition-specific manner (Figure 2A and Supplementary Table 1). Unlike IFNγ that interacts mainly with the IFNγ receptors to induce the expression of several IFNγ-specific genes (ISGs), Toxoplasma is likely to interact with many, yet to be defined, bovine BMDM pattern recognition receptors to induce bovine gene expression. Indeed, there were more differentially expressed genes that were unique to naïve BMDMs exposed to Toxoplasma (RH-specific) than in the uninfected IFNγ-stimulated BMDMs, (657 vs. 134, respectively; Hypergeometric P-value ≤ 0.05). Interestingly, only 96 genes were differentially expressed exclusively in IFNγ-stimulated BMDMs exposed to Toxoplasma (IRH-specific), suggesting that pre-stimulation with IFNγ abrogates a majority of Toxoplasma-induced transcriptional changes in the BMDMs (Figure 2A). The variable expression of most genes (579/657, ~88%), including immunoregulatory genes such as Krüppel-like factor 4 (KLF4) and Ornithine decarboxylase (ODC1), in the RH BMDMs was due to increased transcript abundance in the Toxoplasma-exposed BMDMs, rather than a downregulation of innately expressed genes, as revealed by hierarchical clustering (Figure 2B), suggesting that effective response to the parasite is inducible. Three hundred and forty genes, including several inflammatory cytokines such C-X-C motif chemokine ligand 10 (CXCL10) that were differentially expressed in unexposed IFNγ-stimulated BMDMs, remained dysregulated by a similar fold-change magnitude when the IFNγ-stimulated BMDMs were exposed to Toxoplasma (IRH), suggesting an inability by the parasite to overcome most of the IFNγ-induced transcriptional changes. Arginase 2 (ARG2) was among 13 genes that were highly expressed in RH but downregulated in IFNγ-stimulated BMDMs. Similarly, 60 genes, including the C-C motif chemokine 22 (CCL22), were upregulated in RH but downregulated by a difference of more than 2-fold change between in RH and IRH BMDMs.
Figure 2

Distinct changes in gene expression underpin bovine macrophage response to IFNγ and/or Toxoplasma. (A) RNA-seq analysis of bovine BMDMs that were either left unstimulated and uninfected (Control), unstimulated and infected (RH), or stimulated with IFNγ (100 ng/mL) and uninfected (IFNγ) or stimulated with IFNγ and infected (IRH). Venn diagram representation of the number of dysregulated genes (≥2-fold; ≥10 uniquely aligned reads) relative to the control BMDMs. (B) Heatmap representation of genes dysregulated in RH BMDMs shows most of the stochastic changes in gene expression was due was due to Toxoplasma-induced increase in transcript abundance. (C) Enrichment analysis in the HALLMARK gene sets of the differentially expressed genes shows that Toxoplasma induced a unique gene sets, including the Myc targets V1 and V2.
To gain greater insights into the biological processes underpinning bovine BMDMs transcriptional response to Toxoplasma, and to determine whether the transcriptional changes in the different BMDM conditions was due to differences in the induction of distinct gene sets or the magnitude of induction of the same gene sets, we performed a pre-ranked gene set enrichment analysis (GSEA) (Subramanian et al., 2005) on the differentially expressed genes in IFNγ-stimulated, IRH, and RH BMDMs using the curated “HALLMARK” gene set database (Liberzon et al.,
IFNγ Reverses Toxoplasma-Induced Dysregulation of Genes Linked to Inflammatory and Metabolic Pathways to Restrict Parasite Replication in Bovine Macrophages
As demonstrated above, Toxoplasma significantly dysregulate several bovine genes, potentially to support its intracellular survival. We reasoned that since IFNγ-stimulated BMDMs are refractory to the parasite, Toxoplasma-induced genes that support the intracellular parasite lifestyle are likely to be targeted and reversed when BMDMs are pre-stimulated with IFNγ prior to infection. Several inflammatory genes were dysregulated in IRH but not in RH BMDMs. For example, CXCL9, 10, and 11, were 22-, 10-, and 12-fold, respectively, upregulated in IRH, but not dysregulated in RH, BMDMs. Markers for the Toxoplasma-refractory classically activated macrophages, such as CD180, CD74, and CD1D were highly expressed in IRH but not in RH BMDMs. Similarly, c-type lectins including CLEC6A, which are associated with classically activated macrophages (Jensen et al.,
Figure 3

Dysregulated genes in RH and IRH are enriched in unique and common functional terms. Enrichment analysis in functional annotation and gene ontology (GO) analysis of the differentially expressed genes using metascape show that enrichment for unique functional terms in (A) Genes upregulated and/or unchanged IRH BMDMs, and (B) Genes upregulated in RH but not IRH BMDMs. (C) A subset of differentially expressed genes were validated in quantitative real-time polymerase chain reaction (qPCR) in naïve infected (RH), IFNγ-stimulated and infected (IRH), and IFNγ-stimulated BMDMs. Data are average value ±s.d. of three replicates. P-values of two-tailed unpaired Student's t-test; ***p < 0.001. Data are representative of two independent experiments.
Genes implicated in chemotaxis or anti-inflammatory responses were differentially expressed by a margin of more than 2-fold change between RH and IRH BMDMs. For example, chemotaxis-related CCL22 and CCL24 chemokines were upregulated by 20- and 12-fold in RH but upregulated by 14- and 8-fold in IRH BMDMs. Similarly, the chemokine receptor (CCR7) was upregulated by 26-fold in RH but not differentially expressed in IRH BMDMs. The anti-inflammatory-related suppressor of cytokine signaling (SOCS) 3, and TNF receptor-associated factor 1 (TRAF1) were also upregulated in RH but downregulated in IRH BMDMs. Mitogen-activated protein kinases (MAPK), including MAP3K14 and MAP4K1, which induce non-canonical (Jin et al.,
IFNγ Abrogates Toxoplasma-Induced Wnt/β-Catenin Signaling in Bovine Macrophages
Recent studies have documented Toxoplasma-induced expression of several chemokines, including CCL17, CCL22, and CCL24 (He et al.,
Unlike Toxoplasma, IFNγ is known to suppress Wnt/β-catenin signaling (Nava et al.,
Figure 4

A role for Wnt/β-catenin signaling in the IFNγ-induced response to Toxoplasma in bovine macrophages. (A) Naïve (white bars) or IFNγ-stimulated (black bars) were cultured overnight with or without increasing concentrations of the Wnt/β-catenin signaling antagonist IWR-1-endo and infected with a parental wildtype (RHTy3) Toxoplasma strains. IWR-1-endo did not abolish differences in intracellular parasite burden between naïve and IFNγ-stimulated BMDMs. (B) The growth of the knockout RHΔASP5 parasite strain was significantly inhibited in naïve BMDMs infected, relative to the wildtype. RHΔASP5 parasites do not secrete most dense granule proteins, including GRA18, beyond the parasitophorous vacuole membrane. (C) Unlike wildtype parasites, IWR-1-endo did not significantly alter parasite burden in naïve or IFNγ-stimulated BMDMs infected with RHΔASP5. (D) The parental RHTy3 strain induced significantly higher expression of CCL22 in naïve BMDMs when compared to the knockout RHΔASP5 strain. Data are average value ±s.d. of three replicates. P-values of two-tailed unpaired Student's t-test; *p < 0.05, ***p < 0.001 and ns = not significant. Data are representative of three independent experiments.
Parasites in IFNγ-Stimulated Bovine Macrophages Transcribe Mostly Bradyzoite-Related Genes
To evaluate whether the infection of bovine BMDMs is underpinned by distinct parasite expression signatures, we used the RNA-sequencing reads uniquely aligning to the parasite genome (GT1 v.46) to evaluate stochastic changes in Toxoplasma transcript abundance in naïve and IFNγ-stimulated BMDMs. We focused our analysis on genes that were modulated by more than 3-fold change and had more than 10 uniquely aligned RNA-sequencing reads in at least two samples when comparing parasites infecting naïve and IFNγ-stimulated BMDMs. One hundred and eighty one genes were differentially expressed, of which 94 corresponded to genes upregulated in IFNγ-stimulated BMDMs (Figure 5 and Supplementary Table 1). Parasite genes that were highly expressed when the parasite is in IFNγ-stimulated BMDMs included DnaK-tetratricopeptide repeat (DnaK-TPR) and Cyclic AMP-Dependent Protein Kinase Subunit 3 (cAMPK3), which are associated with or involved in stress-induced Toxoplasma stage conversion from the fast dividing tachyzoite to the semi-dominant bradyzoite parasite stage (Ueno et al., 2011; Sugi et al., 2016). Thus, consistent with the well-defined immune pressure-induced parasite stage conversion, Toxoplasma potentially responds to the IFNγ-induced immune pressure by slowing down its rate of replication.
Figure 5

Toxoplasma gene expression varies in naïve and stimulate bovine macrophage. Heatmap representation of the differentially expressed genes (≥3-fold; ≥10 uniquely aligned reads in at least two samples) between parasites in naïve (RH) and IFNγ-stimulated (IRH) BMDMs.
Discussion
The outcome of Toxoplasma infections in cattle is fundamentally different from that of rodents, small ruminants, and pigs: clinical disease in naturally infected cattle are rare and many large-scale studies have failed to detect viable parasites in bovine tissues (Stelzer et al., 2019), suggesting that, like humans, cattle are largely an end-stage host for the parasite. However, the mechanisms underlying bovine resistance to Toxoplasma are largely unknown. In this study, we performed high throughput functional genomics and parasitological assays to investigate the molecular factors that modulate Toxoplasma interactions with bovine bone marrow-derived macrophages (BMDMs). Similar to observations in rodent and human macrophages, interferon gamma (IFNγ) significantly enhanced intracellular parasite growth restriction in bovine BMDMs. Interestingly, despite the ability to restrict intracellular parasite growth, IFNγ on its own was not able to induce nitric oxide (NO) secretion in bovine BMDMs, indicating that, unlike murine BMDMs, IFNγ-induced toxoplasmastatic activities in bovine BMDMs are independent of NO. Additionally, we found that secreted parasite dense granule proteins potentially induce the expression of chemokines, including CCL22, to enhance parasite growth in naïve bovine BMDMs.
Understanding the mechanisms of innate recognition and response to Toxoplasma in mice has been invaluable in advancing knowledge on response to Toxoplasma in other vertebrates, including cattle. The IFNγ-induced murine anti-Toxoplasma effectors, which are largely dependent on the activation of immunity-related GTPases (IRGs) and guanylate binding proteins (GBPs), are well-characterized. However, cattle, like humans, lack functional toll-like receptor (TLR) 11 and 12 (O'Neill et al.,
Recent studies have documented a role for Wnt/β-catenin signaling in host responses to Toxoplasma. Enhanced intracellular Toxoplasma growth was found to co-occur with increased β-catenin (CTNNB1) gene expression in human cells (Majumdar et al.,
The mechanisms regulating the impact of Wnt/β-catenin on intracellular Toxoplasma replication are largely equivocal. Recent studies indicate that β-catenin alters intracellular parasite growth dynamics by interacting with indoleamine 2,3-dioxygenase 1 (IDO1): IDO1 promoter activity is supported by β-catenin. IDO1 is an IFNγ-inducible protein that degrades tryptophan to kynurenine and is known to impede Toxoplasma growth in human fibroblasts (Pfefferkorn,
Statements
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/ PRJNA646376.
Ethics statement
The animal study was reviewed and approved by the Protocols and Ethics Committees of The Roslin Institute, The University of Edinburgh, and the Royal (Dick) School of Veterinary Medicine. In accordance with the United Kingdom Animal (Scientific Procedures) Act 1986, this study did not require a Home Office project license as no regulated procedures were carried out.
Author contributions
MH conceived and designed the experiments and wrote the manuscript. MH and AG performed the experiments and analyzed the data. All authors contributed to the article and approved the submitted version.
Funding
MH was funded by a University of Edinburgh Chancellor's Fellowship and a Bill and Melinda Gates Foundation award to the Centre for Tropical Livestock Genetics and Health (OPP1127286). The Roslin Institute receives strategic investment funding from the Biotechnology and Biological Sciences.
Acknowledgments
The authors thank the University of Edinburgh Genomics Core for providing facilities and services.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2020.00437/full#supplementary-material
Supplementary Figure 1Quantitative real-time polymerase chain reaction (qPCR) of some differentially expressed genes in naïve infected (RH; black bars), IFNγ-stimulated and infected (IRH; blue bars), and IFNγ-stimulated (red bars), unstimulated-uninfected control (open bars) BMDMs. Data are average value ± s.d. of three replicates. Data are representative of two independent experiments.
References
1
BaiM.LiW.YuN.ZhangH.LongF.ZengA. (2017). The crosstalk between β-catenin signaling and type I, type II and type III interferons in lung cancer cells. Am. J. Transl. Res.9, 2788–2797.
2
BasavarajuA. (2016). Toxoplasmosis in HIV infection: an overview. Trop. Parasitol.6, 129–135. 10.4103/2229-5070.190817
3
BlumeM.SeeberF. (2018). Metabolic interactions between Toxoplasma gondii and its host. F1000Res. 7:F1000. 10.12688/f1000research.16021.1
4
BougdourA.DurandauE.Brenier-PinchartM. P.OrtetP.BarakatM.KiefferS.et al. (2013). Host cell subversion by Toxoplasma GRA16, an exported dense granule protein that targets the host cell nucleus and alters gene expression. Cell Host Microbe13, 489–500. 10.1016/j.chom.2013.03.002
5
BrayN. L.PimentelH.MelstedP.PachterL. (2016). Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol.34, 525–527. 10.1038/nbt.3519
6
BrennerD.GolksA.KieferF.KrammerP. H.ArnoldR. (2005). Activation or suppression of NFκB by HPK1 determines sensitivity to activation-induced cell death. EMBO J.24, 4279–4290. 10.1038/sj.emboj.7600894
7
ChannonJ. Y.SeguinR. M.KasperL. H. (2000). Differential infectivity and division of Toxoplasma gondii in human peripheral blood leukocytes. Infect. Immun.68, 4822–4826. 10.1128/IAI.68.8.4822-4826.2000
8
CirelliK. M.GorfuG.HassanM. A.PrintzM.CrownD.LepplaS. H.et al. (2014). Inflammasome sensor NLRP1 controls rat macrophage susceptibility to Toxoplasma gondii. PLoS Pathog.10:e1003927. 10.1371/journal.ppat.1003927
9
CoffeyM. J.SleebsB. E.UboldiA. D.GarnhamA.FrancoM.MarinoN. D.et al. (2015). An aspartyl protease defines a novel pathway for export of Toxoplasma proteins into the host cell. elife4:e10809. 10.7554/eLife.10809.025
10
Curt-VaresanoA.BraunL.RanquetC.HakimiM. A.BougdourA. (2016). The aspartyl protease TgASP5 mediates the export of the Toxoplasma GRA16 and GRA24 effectors into host cells. Cell Microbiol.18, 151–167. 10.1111/cmi.12498
11
De PaoliF.StaelsB.Chinetti-GbaguidiG. (2014). Macrophage phenotypes and their modulation in atherosclerosis. Circ. J.78, 1775–1781. 10.1253/circj.CJ-14-0621
12
DenisM.WedlockD. N.BuddleB. M. (2005). IFN-γ enhances bovine macrophage responsiveness to Mycobacterium bovis: impact on bacterial replication, cytokine release and macrophage apoptosis. Immunol. Cell Biol.83, 643–650. 10.1111/j.1440-1711.2005.01386.x
13
DoladoI.NebredaA. R. (2008). AKT and oxidative stress team up to kill cancer cells. Cancer Cell14, 427–429. 10.1016/j.ccr.2008.11.006
14
DubeyJ. P. (1986). A review of toxoplasmosis in cattle. Vet. Parasitol.22, 177–202. 10.1016/0304-4017(86)90106-8
15
DunayI. R.SibleyL. D. (2010). Monocytes mediate mucosal immunity to Toxoplasma gondii. Curr. Opin. Immunol.22, 461–466. 10.1016/j.coi.2010.04.008
16
Esteban-RedondoI.InnesE. A. (1997). Toxoplasma gondii infection in sheep and cattle. Comp. Immunol. Microbiol. Infect. Dis.20, 191–196. 10.1016/S.0147-9571(96)00039-2
17
FAO/WHO (2014). Multicriteria-Based Ranking for Risk Management of Food-Borne Parasites. Available online at: http://www.fao.org/3/a-i3649e.pdf%5Cnfiles/540/a-i3649e.pdf (accessed August 9, 2020).
18
FleckensteinM. C.ReeseM. L.Konen-WaismanS.BoothroydJ. C.HowardJ. C.SteinfeldtT. (2012). A Toxoplasma gondii pseudokinase inhibits host IRG resistance proteins. PLoS Biol.10:e1001358. 10.1371/journal.pbio.1001358
19
GazzinelliR. T.HienyS.WynnT. A.WolfS.SherA. (1993). Interleukin 12 is required for the T-lymphocyte-independent induction of interferon γ by an intracellular parasite and induces resistance in T-cell- deficient hosts. Proc. Natl. Acad. Sci. U.S.A.90, 6115–6119. 10.1073/pnas.90.13.6115
20
GazzinelliR. T.WysockaM.HayashiS.DenkersE. Y.HienyS.CasparP.et al. (1994). Parasite-induced IL-12 stimulates early IFN-gamma synthesis and resistance during acute infection with Toxoplasma gondii. J. Immunol.153, 2533–2543.
21
GkogkaE.ReijM. W.HavelaarA. H.ZwieteringM. H.GorrisL. G. M. (2011). Risk-based estimate of effect of foodborne diseases on public health, Greece. Emerg. Infect. Dis.17, 1581–1590. 10.3201/eid1709.101766
22
GomesA. F.MagalhãesK. G.RodriguesR. M.De CarvalhoL.MolinaroR.BozzaP. T.et al. (2014). Toxoplasma gondii-skeletal muscle cells interaction increases lipid droplet biogenesis and positively modulates the production of IL-12, IFN-g and PGE 2. Parasit. Vectors7:47. 10.1186/1756-3305-7-47
23
GreggB.TaylorB. C.JohnB.Tait-WojnoE. D.GirgisN. M.MillerN.et al. (2013). Replication and distribution of Toxoplasma gondii in the small intestine after oral infection with tissue cysts. Infect. Immun.81, 1635–1643. 10.1128/IAI.01126-12
24
HammoudiP. M.JacotD.MuellerC.Di CristinaM.DoggaS. K.MarqJ. B.et al. (2015). Fundamental roles of the golgi-associated Toxoplasma aspartyl protease, ASP5, at the host-parasite interface. PLoS Pathog.11:e1005211. 10.1371/journal.ppat.1005211
25
HardbowerD. M.AsimM.LuisP. B.SinghK.BarryD. P.YangC.et al. (2017). Ornithine decarboxylase regulates M1 macrophage activation and mucosal inflammation via histone modifications. Proc. Natl. Acad. Sci. U.S.A.114, E751–E760. 10.1073/pnas.1614958114
26
HassanM. A.JensenK. D.ButtyV.HuK.BoedecE.PrinsP.et al. (2015). Transcriptional and linkage analyses identify loci that mediate the differential macrophage response to inflammatory stimuli and infection. PLOS Genet.11:e1005619. 10.1371/journal.pgen.1005619
27
HeH.Brenier-PinchartM.-P. P.BraunL.KrautA.TouquetB.CoutéY.et al. (2018). Characterization of a Toxoplasma effector uncovers an alternative GSK3/β-catenin-regulatory pathway of inflammation. elife7:e39887. 10.7554/eLife.39887.027
28
HillD.DubeyJ. P. (2002). Toxoplasma gondii: transmission, diagnosis, and prevention. Clin. Microbiol. Infect.8, 634–640. 10.1046/j.1469-0691.2002.00485.x
29
HunterC. A.SibleyL. D. (2012). Modulation of innate immunity by Toxoplasma gondii virulence effectors. Nat. Rev. Microbiol.10, 766–778. 10.1038/nrmicro2858
30
ImrieH.WilliamsD. J. L. (2019). Stimulation of bovine monocyte-derived macrophages with lipopolysaccharide, interferon-?, Interleukin-4 or Interleukin-13 does not induce detectable changes in nitric oxide or arginase activity. BMC Vet. Res.15:45. 10.1186/s12917-019-1785-0
31
JensenK.KennethH.WhitmarshR.HassanM. A.JulienL.LuD.et al. (2013). Toxoplasma gondii rhoptry 16 kinase promotes host resistance to oral infection and intestinal inflammation only in the context of the dense granule protein GRA15. Infect. Immun.81, 2156–2167. 10.1128/IAI.01185-12
32
JensenK. D. C.WangY.WojnoE. D. T.ShastriA. J.HuK.CornelL.et al. (2011). Toxoplasma polymorphic effectors determine macrophage polarization and intestinal inflammation. Cell Host Microbe9, 472–483. 10.1016/j.chom.2011.04.015
33
JinJ.HuH.LiH. S.YuJ.XiaoY.BrittainG. C.et al. (2014). Noncanonical NF-κB pathway controls the production of type I interferons in antiviral innate immunity. Immunity40, 342–354. 10.1016/j.immuni.2014.02.006
34
JungiT. W.AdlerH.AdlerB.ThönyM.KrampeM.PeterhansE. (1996). Inducible nitric oxide synthase of macrophages. Present knowledge and evidence for species-specific regulation. Vet. Immunol. Immunopathol.54, 323–330. 10.1016/S0165-2427(96)05690-5
35
KoshyA. A.DietrichH. K.ChristianD. A.MelehaniJ. H.ShastriA. J.HunterC. A.et al. (2012). Toxoplasma co-opts host cells it does not invade. PLoS Pathog.8:e1002825. 10.1371/journal.ppat.1002825
36
LiM. H.CheT. C.KuoT. T.TsengC. C. P.TsengC. C. P. (2000). Real-time PCR for quantitative detection of Toxoplasma gondii. J. Clin. Microbiol.38, 4121–4125. 10.1128/JCM.38.11.4121-4125.2000
37
LiberzonA.BirgerC.ThorvaldsdóttirH.GhandiM.MesirovJ. P.TamayoP. (2015). The molecular signatures database hallmark gene set collection. Cell Syst.1, 417–425. 10.1016/j.cels.2015.12.004
38
LiuW.JiangJ.YanD.LiD.LiW.MaY.et al. (2014). Pentraxin 3 promotes oxLDL uptake and inhibits cholesterol efflux from macrophage-derived foam cells. Exp. Mol. Pathol.96, 292–299. 10.1016/j.yexmp.2014.03.007
39
MacKenzieC.HeselerK.MullerA.DaubenerW. (2007). Role of indoleamine 2,3-dioxygenase in antimicrobial defence and immuno-regulation: tryptophan depletion versus production of toxic kynurenines. Curr. Drug Metab.8, 237–244. 10.2174/138920007780362518
40
MacMickingJ.XieQ. W.NathanC. (1997). Nitric oxide and macrophage function. Annu. Rev. Immunol.15, 323–350. 10.1146/annurev.immunol.15.1.323
41
MajumdarT.SharmaS.KumarM.HussainM. A.ChauhanN.KaliaI.et al. (2019). Tryptophan-kynurenine pathway attenuates β-catenin-dependent pro-parasitic role of STING-TICAM2-IRF3-IDO1 signalosome in Toxoplasma gondii infection. Cell Death Dis.10:161. 10.1038/s41419-019-1420-9
42
MeloM.NguyenQ.CordeiroC.HassanM.YangN.McKellR.et al. (2013). Transcriptional analysis of murine macrophages infected with different Toxoplasma strains identifies novel regulation of host signaling pathways. PLoS Pathog.9:e1003779. 10.1371/journal.ppat.1003779
43
MeloM. B.JensenK. D.SaeijJ. P. (2011). Toxoplasma gondii effectors are master regulators of the inflammatory response. Trends Parasitol.27, 487–495. 10.1016/j.pt.2011.08.001
44
MeloM. B.KasperkovitzP.CernyA.Konen-WaismanS.Kurt-JonesE. A.LienE.et al. (2010). UNC93B1 mediates host resistance to infection with Toxoplasma gondii. PLoS Pathog.6:e1001071. 10.1371/journal.ppat.1001071
45
MorisakiJ. H.HeuserJ. E.SibleyL. D. (1995). Invasion of Toxoplasma gondii occurs by active penetration of the host cell. J. Cell Sci.108, 2457–2464.
46
MurrayP. J.AllenJ. E.BiswasS. K.FisherE. A.GilroyD. W.GoerdtS.et al. (2014). Macrophage activation and polarization: nomenclature and experimental guidelines. Immunity41, 14–20. 10.1016/j.immuni.2014.06.008
47
NavaP.KamekuraR.QuirósM.Medina-ContrerasO.HamiltonR. W.KolegraffK. N.et al. (2014). IFNγ-induced suppression of β-catenin signaling: evidence for roles of Akt and 14.3.3ζ. Mol. Biol. Cell25, 2894–2904. 10.1091/mbc.e13-09-0512
48
NiedelmanW.SprokholtJ. K.CloughB.FrickelE. M.SaeijJ. P. (2013). Cell death of gamma interferon-stimulated human fibroblasts upon Toxoplasma gondii infection induces early parasite egress and limits parasite replication. Infect. Immun.81, 4341–4349. 10.1128/IAI.00416-13
49
O'NeillL. A. J.GolenbockD.BowieA. G. (2013). The history of toll-like receptors-redefining innate immunity. Nat. Rev. Immunol.13, 453–460. 10.1038/nri3446
50
PfefferkornE. R. (1984). Interferon gamma blocks the growth of Toxoplasma gondii in human fibroblasts by inducing the host cells to degrade tryptophan. Proc. Natl. Acad. Sci. U.S.A.81, 908–912. 10.1073/pnas.81.3.908
51
PimentelH.BrayN. L.PuenteS.MelstedP.PachterL. (2017). Differential analysis of RNA-seq incorporating quantification uncertainty. Nat. Methods14, 687–690. 10.1038/nmeth.4324
52
PriceJ. V.VanceR. E. (2014). The macrophage paradox. Immunity41, 685–693. 10.1016/j.immuni.2014.10.015
53
ReyaT.DuncanA. W.AillesL.DomenJ.SchererD. C.WillertK.et al. (2003). A role for wnt signalling in self-renewal of haematopoietic stem cells. Nature423, 409–415. 10.1038/nature01593
54
Scharton-KerstenT. M.YapG.MagramJ.SherA. (1997). Inducible nitric oxide is essential for host control of persistent but not acute infection with the intracellular pathogen Toxoplasma gondii. J. Exp. Med.185, 1261–1273. 10.1084/jem.185.7.1261
55
SchneemannM.SchoedonG. (2002). Species differences in macrophage NO production are important. Nat. Immunol.3:102. 10.1038/ni0202-102a
56
SibleyL. D.AdamsL. B.FukutomiY.KrahenbuhlJ. L. (1991). Tumor necrosis factor-alpha triggers antitoxoplasmal activity of IFN-gamma primed macrophages. J. Immunol.147, 2340–2345.
57
SindrilaruA.Scharffetter-KochanekK. (2013). Disclosure of the culprits: macrophages—versatile regulators of wound healing. Adv. Wound Care2, 357–368. 10.1089/wound.2012.0407
58
SongH. B.JunH. O.KimJ. H. J. H.LeeY. H.ChoiM. H.KimJ. H. J. H. (2017). Disruption of outer blood-retinal barrier by Toxoplasma gondii-infected monocytes is mediated by paracrinely activated FAK signaling. PLoS ONE12:e0175159. 10.1371/journal.pone.0175159
59
StaalF. J. T.LuisT. C.TiemessenM. M. (2008). WNT signalling in the immune system: WNT is spreading its wings. Nat. Rev. Immunol.8, 581–593. 10.1038/nri2360
60
StelzerS.BassoW.Benavides SilvánJ.Ortega-MoraL. M.MaksimovP.GethmannJ.et al. (2019). Toxoplasma gondii infection and toxoplasmosis in farm animals: risk factors and economic impact. Food Waterborne Parasitol.15:e00037. 10.1016/j.fawpar.2019.e00037
61
SubramanianA.TamayoP.MoothaV. K.MukherjeeS.EbertB. L.GilletteM. A.et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U.S.A.102, 15545–15550. 10.1073/pnas.0506580102
62
SugiT.MaY. F.TomitaT.MurakoshiF.EatonM. S.YakubuR.et al. (2016). Toxoplasma gondii cyclic AMP-dependent protein kinase subunit 3 is involved in the switch from tachyzoite to bradyzoite development. mBio7, e00755–16. 10.1128/mBio.00755-16
63
ThiE. P.LambertzU.ReinerN. E. (2012). Sleeping with the enemy: how intracellular pathogens cope with a macrophage lifestyle. PLoS Pathog.8:e1002551. 10.1371/journal.ppat.1002551
64
TombáczK.BurgessG.HolderA.WernersA.WerlingD. (2018). Toxoplasma gondii profilin does not stimulate an innate immune response through bovine or human TLR5. Innate Immun.24, 422–429. 10.1177/1753425918798450
65
ToshK. W.MitterederL.Bonne-AnneeS.HienyS.NutmanT. B.SingerS. M.et al. (2016). The IL-12 response of primary human dendritic cells and monocytes to Toxoplasma gondii is stimulated by phagocytosis of live parasites rather than host cell invasion. J. Immunol.196, 345–356. 10.4049/jimmunol.1501558
66
UenoA.DautuG.HagaK.MunyakaB.CarmenG.KobayashiY.et al. (2011). Toxoplasma gondii: a bradyzoite-specific DnaK-tetratricopeptide repeat (DnaK-TPR) protein interacts with p23 co-chaperone protein. Exp. Parasitol.127, 795–803. 10.1016/j.exppara.2011.01.015
67
van de VosseE.van DisselJ. T.OttenhoffT. H. (2009). Genetic deficiencies of innate immune signalling in human infectious disease. Lancet Infect. Dis.9, 688–698. 10.1016/S1473-3099(09)70255-5
68
WileyM.SweeneyK. R.ChanD. A.BrownK. M.McMurtreyC.HowardE. W.et al. (2010). Toxoplasma gondii activates hypoxia-inducible factor (HIF) by stabilizing the HIF-1α subunit via type I activin-like receptor kinase receptor signaling. J. Biol. Chem.285, 26852–26860. 10.1074/jbc.M110.147041
69
YamamotoM.MaJ. S.MuellerC.KamiyamaN.SaigaH.KuboE.et al. (2011). ATF6beta is a host cellular target of the Toxoplasma gondii virulence factor ROP18. J. Exp. Med.208, 1533–1546. 10.1084/jem.20101660
70
YeJ.CoulourisG.ZaretskayaI.CutcutacheI.RozenS.MaddenT. L. (2012). Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinf.13:134. 10.1186/1471-2105-13-134
71
YoungR.BushS. J.LefevreL.McCullochM. E. B.LisowskiZ. M.MuriukiC.et al. (2018). Species-specific transcriptional regulation of genes involved in nitric oxide production and arginine metabolism in macrophages. Immunohorizons2, 27–37. 10.4049/immunohorizons.1700073
72
ZhaoY.MarpleA. H.FergusonD. J. P.BzikD. J.YapG. S. (2014). Avirulent strains of Toxoplasma gondii infect macrophages by active invasion from the phagosome. Proc. Natl. Acad. Sci. U.S.A.111, 6437–6442. 10.1073/pnas.1316841111
73
ZhaoY. O.KhaminetsA.HunnJ. P.HowardJ. C. (2009). Disruption of the Toxoplasma gondii parasitophorous vacuole by IFNgamma-inducible immunity-related GTPases (IRG proteins) triggers necrotic cell death. PLoS Pathog.5:e1000288. 10.1371/journal.ppat.1000288
Summary
Keywords
bovine toxoplasmosis, macrophages, RNA-sequencing, Toxoplasma, chemokine, Wnt/β-catenin-signaling
Citation
Gossner A and Hassan MA (2020) Transcriptional Analyses Identify Genes That Modulate Bovine Macrophage Response to Toxoplasma Infection and Immune Stimulation. Front. Cell. Infect. Microbiol. 10:437. doi: 10.3389/fcimb.2020.00437
Received
27 April 2020
Accepted
16 July 2020
Published
20 August 2020
Volume
10 - 2020
Edited by
Prajwal Gurung, The University of Iowa, United States
Reviewed by
David Arranz-Solís, University of California, Davis, United States; Xiaoyu Hu, University of Texas Southwestern Medical Center, United States
Updates

Check for updates
Copyright
© 2020 Gossner and Hassan.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Musa A. Hassan musa.hassan@roslin.ed.ac.uk
This article was submitted to Microbes and Innate Immunity, a section of the journal Frontiers in Cellular and Infection Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.