Differential Brain MicroRNA Expression Profiles After Acute and Chronic Infection of Mice With Toxoplasma gondii Oocysts

Brain microRNAs (miRNAs) change in abundance in response to Toxoplasma gondii infection. However, their precise role in the pathogenesis of cerebral infection with T. gondii oocyst remains unclear. We studied the abundance of miRNAs in the brain of mice on days 11 and 33 post-infection (dpi) in order to identify miRNA pattern specific to early (11 dpi) and late (33 dpi) T. gondii infection. Mice were challenged with T. gondii oocysts (Type II strain) and on 11 and 33 dpi, the expression of miRNAs in mouse brain was investigated using small RNA (sRNA) sequencing. miRNA expression was confirmed by quantitative reverse transcription polymerase chain reaction (qRT-PCR). Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to identify the biological processes, molecular functions, and cellular components, as well as pathways involved in infection. More than 1,500 miRNAs (1,352 known and 150 novel miRNAs) were detected in the infected and control mice. The expression of miRNAs varied across time after infection; 3, 38, and 108 differentially expressed miRNAs (P < 0.05) were detected during acute infection, chronic infection and chronic vs. acute infection, respectively. GO analysis showed that chronically infected mice had more predicted targets of dysregulated miRNAs than acutely infected mice. KEGG analysis indicated that most predicted targets were involved in immune- or disease-related pathways. Our data indicate that T. gondii infection alters the abundance of miRNAs in mouse brain particularly at the chronic stage, probably to fine-tune conditions required for the establishment of a latent brain infection.

Introduction

The intracellular protozoan parasite Toxoplasma gondii is an opportunistic pathogen, which can virtually infect and replicate within any nucleated cells of warm-blooded animals and humans (Jones and Dubey, 2012). T. gondii has a complex life cycle that includes asexual propagation and the formation of tachyzoites and bradyzoites-containing cysts in the intermediate host (Black and Boothroyd, 2000) and sexual reproduction and formation of oocysts in the intestinal epithelium of felids (Tenter et al., 2000). Humans can be infected through (i) ingesting undercooked meat containing T. gondii tissue cysts, (ii) drinking water contaminated with sporulated oocysts, and (iii) transplacental (vertical) transmission (Elsheikha, 2008). This zoonotic pathogen infects approximately one-third of the world population and can cause a variety of clinical symptoms and even death in immuno-compromised patients (e.g., AIDS patients and organ transplant recipients) and in fetuses of naïve women infected during pregnancy (Dubey, 2008; Elsheikha, 2008; Zhou et al., 2011).

Previous studies have reported a correlation between a deregulated immunoinflammatory response and brain dysfunction in infected humans and animals (Cannella et al., 2014; Blanchard et al., 2015; Elsheikha and Zhu, 2016; Marra, 2018). Effective therapeutic interventions to control brain infection are therefore desirable and will be facilitated by a better understanding of the host immune response against the parasite. However, knowledge about the molecular mechanisms underlying the deregulation of immune responses observed in oocyst-induced cerebral toxoplasmosis remains limited. Oocysts have remarkable ability to endure in the environment (Lindsay and Dubey, 2009; Torrey and Yolken, 2013). Infection acquired through the consumption of water contaminated with T. gondii oocysts has been frequently reported (Dubey, 2004; Hill et al., 2011), and can be linked to even waterborne outbreaks (Boyer et al., 2011).

MicroRNA (miRNA) profiling has emerged as a useful approach to study the pathogenesis of many protozoan species, such as Cryptosporidium, Plasmodium, and Toxoplasma (Judice et al., 2016). Also, miRNAs have been promising biomarkers for the diagnosis and monitoring of progression of parasitic diseases (Hoy et al., 2014; Cai et al., 2015). miRNAs are short (20–24 nucleotides) endogenous non-coding, single stranded, RNA sequences that can control gene expression at the posttranscriptional level and mediate regulatory signals between cells in health and disease (Bartel, 2004; Triboulet et al., 2007; Hobert, 2008; Winter et al., 2009; Zeiner et al., 2010). T. gondii infection requires specific miRNA for efficient replication (Zeiner et al., 2010; Cong et al., 2017) and can manipulate host signaling pathways (Hakimi and Ménard, 2010). Infection with the parasite tissue cysts and tachyzoites can alter miRNA expression in the brain (Xu et al., 2013) and spleen (He et al., 2016) of the host, respectively. However, miRNA expression patterns in the mouse brain during acute and chronic stages of infection with T. gondii oocysts is unknown.

We previously reported, using next-generation sequencing technology, the differential expression of miRNAs in response to T. gondii infection in mouse liver (Cong et al., 2017) and brain (Xu et al., 2013). To our knowledge, no other study has yet explored the possible role of miRNAs in murine cerebral toxoplasmosis caused by oocyst infection. Expanding our earlier observations of differential expression of specific miRNAs between healthy and infected mice, we performed small RNA transcriptome sequencing analysis of the mouse brain in response to infection with T. gondii oocysts. Our data provide a platform for the design of functional studies to map the function of the differentially expressed miRNA identified in the brain of mice infected with T. gondii oocysts. Our findings indicate that specific miRNAs regulate the expression of inflammatory cytokines in the mouse brain in response to infection.

Materials and Methods

Ethics Approval

All animal experiments were approved by the Animal Administration and Ethics Committee of Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences. The animals were handled in compliance with the animal ethics requirements of the People’s Republic of China. Every effort was made to minimize animal suffering during the experiment.

Production and Purification of Oocysts

One, 10-week-old, specific-pathogen-free, kitten was infected orally with 100 freshly prepared parasite cysts obtained from brain homogenate of Kunming mice infected with T. gondii PRU strain. This strain was used because the majority of human toxoplasmosis cases have been associated with type II strains (Howe and Sibley, 1995). The cat feces were examined daily for the presence of oocysts. Once detected, T. gondii oocysts were isolated from the feces using sucrose flotation and CsCl gradient, as described previously (Staggs et al., 2009). To induce sporulation, oocysts were centrifuged at 360 × g and the oocyst’s pellet was suspended in 2% sulfuric acid and aerated on a shaker for 7 days at ambient temperature. Sporulated oocysts were washed twice with 0.85% saline and suspended in 2% sulfuric acid. Finally, the number of oocysts was determined using hemocytometer and adjusted to 100 oocysts/ml in PBS and stored at 4°C.

Infection of Mice

Female 7-week-old BALB/c mice were purchased from Lanzhou University Laboratory Animal Centre (Lanzhou, Gansu Province, China). All mice were handled according to protocols approved by the Animal Ethics Committee of Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences. Mice were housed in an Animal Biosafety Level 2 (ABSL-2) containment laboratory, with temperature-controlled room (22 ± 0.5°C) under 12 h light/dark cycles. Mice were fed a commercial rodent pellet diet and had access to water ad libitum. Mice were rested for 1 week before being infected. Twelve mice were randomly divided into four groups (three mice/group): mice infected for 11 days, mice infected for 33 days, control mice for 11 days, and control mice for 33 days. T. gondii infection was induced in each mouse via oral inoculation of 100 T. gondii oocysts in 1 ml of PBS. The uninfected (control) mice were sham-inoculated with 1 ml of PBS only. Body weight of the mice was measured daily following infection, and all mice were monitored daily for the development of clinical signs characteristics of T. gondii infection, such as ruffled hair, neurological manifestations and physical activity. At 11 and 33 dpi, mice were anesthetized by intraperitoneal injection with 100 μl xylazine (20 mg/ml) and ketamine (1 mg/ml) in PBS, brain tissues were harvested and quickly washed in PBS. All brain samples were placed separately in sterile tubes, flash frozen in liquid nitrogen, and stored frozen at -80°C, until used.

Detection of T. gondii in the Brain

TIANamp Genomic DNA kit was used to extract DNA from the brain tissues (TianGen^TM, Beijing, China) and DNA samples were stored frozen at -20°C. The presence of T. gondii in mouse brain was investigated using a semi-nested PCR assay targeting T. gondii B1 gene (Cong et al., 2016). Positive PCR amplicons were genotyped using PCR-restriction fragment length polymorphism analysis (PCR-RFLP) as previously described (Cong et al., 2015). Samples from mouse brain were collected and fixed in 10% buffered formalin (pH 7.2) for a few days before dehydration through a graded series of alcohol to xylol and embedded in paraffin wax. Sections of 5 μm thick from paraffin wax blocks were cut and stained with hematoxylin and eosin (H & E) for histopathological analysis.

RNA Extraction

Total RNA was extracted from the brain tissue of infected and non-infected mice at 11 and 33 dpi. Brain tissues were homogenized in 1 ml Trizol reagent (Invitrogen, Carlsbad, CA, United States). The concentration of RNA was determined using Nanodrop 2000 spectrophotometer (Thermo Scientific, United States). The extracted RNA samples were subjected to quality control checks in order to ensure the high quality of sRNA library construction. The purity of the RNA preparation was assessed by calculating the ratio at 260 and 280 nm. All RNA preparations had a ratio of absorbance (260/280 nm) > 1.8. The assay also confirmed that the RNA samples were free of genomic DNA contamination. The RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, United States). Only RNA samples with the RNA integrity numbers (RINs) > 7 were used for miRNA profiling analysis. The extracted RNA samples were stored frozen at -80°C, until analysis.

Small RNA Library Preparation and Sequencing

About 3 μg of total RNA per sample was used as input material to construct small RNA (sRNA) sequencing library according to the instructions of NEBNext^® Multiplex sRNA Library Prep Set for Illumina^® (NEB, United States). Index codes were added to link sequences to the respective sample. Agilent Bioanalyzer 2100 system (Agilent, Santa Clara, CA, USA) was used to assess the quality of the libraries using DNA High Sensitivity Chips. Following the library generation, the coded samples was clustered on a cBot Cluster Generation System using TruSeq SR Cluster Kit v3-cBot-HS (Illumina). After cluster generation, the libraries were sequenced on an Illumina HiSeq X Ten platform and 50 bp single-end reads were generated.

Bioinformatics Analysis

Custom Perl and Python scripts were used to process the raw reads, where ploy-N, with 5′ adapter contaminants, without 3′ adapter, or the insert tag, containing ploy A or T or G or C, and low-quality reads were removed in order to obtain clean reads. The Q20, Q30, GC content, and the error rate of the clean reads were determined. A length range of 18∼35 nt from clean reads (about 92% of total reads) was chosen to perform all subsequent analyses. The small RNA tags were mapped to the published reference Mus musculus genome sequence by Bowtie (Langmead et al., 2009), without mismatch to analyze the expression and distribution of sRNAs. Mapped sRNA tags were used to search for known miRNA. miRbase 20.0 was used as reference, and software mirdeep 2 (Friedländer et al., 2012) and srna-tools-cli were used to identify the potential miRNA.

RepeatMasker and Rfam database were used for sRNA mapping in order to remove non-coding RNAs (tRNA, rRNA, snRNA, and snoRNA), protein-coding genes and repeat gene sequences. The characteristics of hairpin structure of miRNA precursors, miREvo (Wen et al., 2012) and mirdeep2 (Friedländer et al., 2012) were integrated to predict novel miRNA. miRNA counts and base bias on the first position of the identified miRNA with certain length and on each position of all identified miRNAs were identified by custom scripts. After sRNA annotation, miFam.dat¹ was utilized to search for known miRNA’s families, and novel miRNA precursor was submitted to Rfam² to look for Rfam families, and then explore the occurrence of miRNA families identified from the samples in other families. The prediction of the target genes of miRNAs was performed using miRanda (Betel et al., 2008) and PITA³. The quantification of miRNA was evaluated by TPM (Zhou et al., 2010). The DESeq R package (1.8.3) was used for the differential expression analysis (Anders and Huber, 2010). FDR adjusted P-value (Benjamini-Hochberg method for multiple corrections) of 0.05 was considered significant.

GO and KEGG Enrichment Analyses

Gene ontology enrichment analysis of the predicted target genes of the differentially expressed miRNAs (P < 0.05) was performed using the GOseq R package, based on a Wallenius non-central hyper-geometric distribution (Young et al., 2010). KEGG⁴ pathway analysis and functional annotation of the predicted target genes were conducted using KOBAS 3.0 software (Mao et al., 2005; Kanehisa et al., 2008).

Verification of miRNA Expression by qRT-PCR

The data were validated by quantitative reverse transcription PCR (qRT-PCR) analysis of seven of the differentially expressed miRNAs (P < 0.05) in the brain sample to confirm gene expression ratios obtained by sequencing. Total RNA from infected and uninfected mouse groups was extracted using Trizol method (Invitrogen, United States) and the quality of RNA template was assessed using a NanoDrop 2000 spectrophotometer (Thermo Scientific, United States). RT-PCR reactions were carried out in biological triplicate for each RNA sample. The extracted RNA was treated with DNase I to remove any residual genomic DNA and then reverse-transcripted into single strand cDNA using Mir-X^TM miRNA First-Strand Synthesis Kit (Clontech, Mountain View, CA, United States). SYBR^® Premix Ex Taq^TM II (Takara, Shiga-ken, Japan) was used to perform qRT-PCR reaction on QIAGEN’s real-time PCR cycler (QIAGEN, Hilden, Germany). The 25 μl qRT-PCR reaction contained 9 μl ddH₂O, 12.5 μl SYBR Advantage Premix (2X), 0.5 μl ROX Dye (50 X), 0.5 μl of each miRNA-specific forward primer (Table 1), 0.5 μl mRQ 3′ Primer, and 2 μl cDNA. The qRT-PCRs were performed under the following cycling conditions: initial denaturation at 95°C for 30 s, followed by 40 cycles of 95°C for 5 s, 60°C for 30 s; melt curve analysis was performed from 60 to 95°C to confirm primer specificity by the presence of single melting curve peak, indicating a single amplicon in each qPCR reaction. No-template control and no-reverse transcriptase control were included in each plate to verify the absence of contamination. Expression of miRNAs was normalized to the level of U6 small nuclear RNA (snRNA). qRT-PCR analysis was performed using the delta-delta C_t method (ΔΔC_t) to measure the relative abundances of miRNAs (Peltier and Latham, 2008). The data were expressed as mean ± standard deviation (SD). Statistical significance was determined by Student’s t-test, with P-values of < 0.05 deemed to be statistically significant.

TABLE 1

TABLE 1. Oligonucleotides used as primers for miRNA-specific qRT-PCR analysis.

Results

Oocysts Infection in the Mouse Brain

Toxoplasma gondii infection was confirmed in the brain of infected mice at 11 and 33 dpi by positive PCR results. PCR-RFLP analysis of the positive amplicons of T. gondii B1 gene revealed a restriction fragment pattern consistent with that of T. gondii genotype II. The mouse brain of control group and negative PCR control samples revealed negative results. Histopathological analysis revealed the presence of T. gondii cyst in the brain of infected mice at 33 dpi (Supplementary Figure S1), whereas the parasite cysts were not detected in the brain of uninfected mice or the brain of mice 11 days post infection.

miRNA Expression Patterns Associated With Infection

miRNA libraries were successfully prepared from the brain of the infected and uninfected (control) mice. Key characteristics of the obtained sequencing data are summarized (Table 2). After selecting the appropriate length of sRNA, the base number distribution interval of clean reads in each sample was between 18∼35 nt; one of the highest proportions was 22 nt length of sRNA. 91.9–93.7% of the reads aligned to the reference M. musculus genome. During early and late infection, an unequal number of known and novel miRNAs in infected and uninfected mice was detected (Supplementary Tables S1, S2).

TABLE 2

TABLE 2. Characteristics of the sRNA sequences obtained in the present study.

Differentially Expressed miRNA During Acute and Chronic Infection

The pattern of global expression of miRNAs (P < 0.05) in the brain of healthy mice was compared with that of acutely and chronically infected mice. Total differentially expressed miRNAs (P < 0.05) during acute and chronic infection are summarized (Supplementary Table S3). More miRNAs were differentially expressed during chronic compared to acute stage of T. gondii infection (Figure 1). By comparing acutely infected mice with uninfected mice, 2 miRNAs were upregulated and one miRNA was downregulated (Table 3). When comparing chronically infected mice with uninfected mouse, more differentially expressed miRNAs (P < 0.05) were detected, including 25 upregulated miRNAs and 13 downregulated miRNAs (Supplementary Table S3).

FIGURE 1

FIGURE 1. Volcano plots of miRNA expression changes in the mouse brain. Differentially expressed miRNAs (P < 0.05) in (A) “acute vs uninfected,” (B) “chronic vs. uninfected” and (C) “chronic vs. acute.” Red, green and blue dots represent upregulated, downregulated and non-regulated miRNAs, respectively. The x-axis represents the log2 ratio of gene expression levels between mouse groups. The y-axis is adjusted p-value based on -log10.

TABLE 3

TABLE 3. The top differentially expressed miRNAs (P < 0.05) in mouse brains during acute and chronic infection with Toxoplasma gondii oocysts.

Next, we compared the abundance of miRNA during acute and chronic infection in order to determine the temporal changes in the expression of miRNAs during the course of infection. Out of the 108 differentially expressed miRNAs (P < 0.05), 59 were up-regulated and 49 were down-regulated. The differentially expressed miRNAs (P < 0.05) in early and late T. gondii infection are shown (Figure 2). The number of differentially expressed miRNAs (P < 0.05) increased from 3 at 11 dpi to 38 at 33 dpi as shown in Supplementary Figure S2. These data indicate that more altered expression of miRNAs characterize the brain response to chronic T. gondii infection.

FIGURE 2

FIGURE 2. Venn diagram showing the number of commonly expressed and specifically expressed miRNAs between mouse groups. The miRNAs that are significant and specific for in mice with acute infection vs. healthy control are shown in the yellow circle. The light blue circle represents the miRNA markers that discriminate chronically infected mice and control mice, while the purple circle stands for the miRNA markers that permit a distinction of chronically vs. acutely infected mice.

To validate the deep sequencing results, qRT-PCR was carried out on seven selected miRNAs showing different levels of expression (i.e., 5 upregulated and 2 downregulated) in the brain of mice. These included: mmu-miR-155-5p, mmu-miR-204-5p, mmu-miR-146a-5p, mmu-miR-142a-5p, mmu-miR-7043-3P, mmu-miR-144-3p, and mmu-miR-5114. As shown in Figure 3, qRT-PCR results of the seven examined miRNAs were consistent with those obtained by sequencing, confirming the correctness of the obtained miRNA data.

FIGURE 3

FIGURE 3. Validation of the differentially expressed miRNAs using qRT-PCR. A bar graph showing the agreement in the expression levels of a panel of seven miRNAs between qRT-PCR and the sequencing data. Columns and error bars indicate means and standard deviations of relative expression levels (n = 3), respectively. X-axis, represents the differentially expressed miRNAs (P < 0.05). The y-axis shows the relative expression levels (log₂ fold change).

Functional Annotation and KEGG Pathway Enrichment Analysis of miRNA Target Genes

To better understand the roles of the identified miRNAs in mouse brain response to infection, the target genes of miRNAs were identified using the computational prediction tools miRanda and PITA. We performed GO and KEGG enrichment analyses to identify the biological functions of the differentially expressed miRNAs (P < 0.05) during acute and chronic infection. The results revealed a total of 33 target genes of the three miRNAs of “acute vs. uninfected,” 719 target genes of the 38 miRNAs of “chronic vs. uninfected,” and 1,911 target genes of the 108 miRNAs of “chronic vs. acute.” We found that, in “acute vs. uninfected” mice, two upregulated miRNAs were successfully assigned to 14 significantly enriched GO terms, and one downregulated miRNA was significantly enriched in 21 GO terms (corrected P ≤ 0.05). Among these significantly enriched GO terms, we found that the upregulated miRNAs were significantly enriched in negative regulation of calcium ion transmembrane transporter activity, negative regulation of calcium ion transmembrane transport, negative regulation of ion transmembrane transporter activity, negative regulation of calcium ion transport, and one downregulated gene was significantly enriched in the cell part, cell, biological regulation, and intracellular membrane-bounded organelle.

Regarding “chronic infection vs. uninfected” mice, we detected 25 upregulated miRNAs that were successfully assigned to 652 significantly enriched GO terms, and 13 downregulated miRNAs that were significantly enriched in 532 GO terms (corrected P ≤ 0.05). Among these significantly enriched GO terms, many upregulated miRNAs were enriched in cell part, cell, intracellular part, intracellular, and many downregulated genes were enriched in the cell, cell part, intracellular part, and intracellular.

In “chronic vs. acute” infection, 59 upregulated miRNAs were assigned to 879 significantly enriched GO terms, and 49 downregulated miRNAs were enriched in 1,172 GO terms (corrected P ≤ 0.05). Among these significantly enriched GO terms, many upregulated miRNAs were enriched in cell part, cell, intracellular and intracellular part, and many downregulated genes were enriched in the cell part, cell, single-organism process, and intracellular.

Target genes of the differentially expressed miRNA were mapped to terms in the KEGG database to identify signaling pathways operating during acute and chronic infection. A total of 33 target genes were enriched in 52 KEGG pathways in “acute vs. uninfected” mice, 719 were enriched in 230 KEGG pathways in “chronic vs. uninfected,” and 911 were enriched in 255 KEGG pathways in “chronic vs. acute”. The top 20 highly enriched pathways in acute infection are shown in Figure 4A. such as adrenergic signaling in cardiomyocytes, calcium signaling pathway, and cardiac muscle contraction. The top 20 highly enriched pathways during chronic infection (Figure 4B), included immune-related pathways, such as chemokine signaling pathway, Rap1 signaling pathway, cAMP signaling pathway, MAPK signaling pathway, Fc gamma R-mediated phagocytosis, and Hippo signaling pathway. Other pathways associated with disease, such as pathways in cancer, chronic myeloid leukemia, HTLV-I infection were also detected. The top highly enriched pathways in “chronic vs. acute” were also involved in pathways, such as cancer, endocrine resistance and chronic myeloid leukemia (Figure 4C). These results suggest that infection with T. gondii oocysts stimulates the expression of gene targets at miRNA levels in the mouse brain especially during chronic infection and that miRNAs play key roles in mouse brain immune response to T. gondii infection.

FIGURE 4

FIGURE 4. KEGG analysis of the differentially expressed miRNAs (P < 0.05) revealed significant enrichment in immune-related pathways. The top 20 enriched pathways of the differentially expressed miRNAs are presented for (A) “Acute vs. Uninfected,” (B) “Chronic vs. Uninfected” and (C) “Chronic vs. Acute.” The degree of KEGG enrichment was determined by the enrichment factor, q-value, and miRNA gene target number. The sizes and colors of the spots represent the number of predicted gene targets of the differentially expressed miRNAs (P < 0.05) and the q-value, respectively.

Discussion

Successful outcome of T. gondii infection depends on the parasite ability to subvert host immune response, and to survive and replicate in a hostile host environment. However, the mechanisms by which this parasite interferes with host immunity in the brain remain poorly understood. MicroRNAs (miRNAs) have been established as key regulators of various biological processes with possible involvement in the pathobiology of T. gondii infection (Hakimi and Ménard, 2010). T. gondii infection is capable of modulating host microRNAs (miRNAs) in order to establish conditions that work in favor of parasite growth (Xu et al., 2013; Cong et al., 2017). In agreement with previous work, deep-sequencing analysis of miRNAs of mouse brain at 11 and 33 dpi revealed significant changes in the abundance of miRNAs attributed to infection.

Our analysis revealed > 1,500 miRNAs in both infected and uninfected mouse groups, and some of these were differentially expressed in the infected brains compared with control brain samples. For example, three differentially expressed miRNAs (P < 0.05) were detected during acute infection, whereas 38 differentially expressed miRNAs (P < 0.05) were detected in chronic infection, indicating that miRNA expression in mouse brain was influenced by the length of time elapsed after infection (Pittman and Knoll, 2015). There were 108 differentially expressed miRNAs (P < 0.05) between acutely and chronically infected mice. Although miRNAs in mouse brain tissue challenged with T. gondii type II cyst (PRU strain) at duration of infection of 14 and 21 dpi have been previously reported (Xu et al., 2013), the present study demonstrated that the infecting stage and the duration of infection can influence the expression of miRNAs in mouse brain (Supplementary Table S4). Some of the upregulated miRNAs (mmu-miR-155-5p, mmu-miR-146a-5p, mmu-miR-142a-3p, mmu-miR-142b, mmu-miR-21a-5p) and the down-regulated miRNAs (mmu-miR-409-5p, mmu-miR-127-3p, mmu-miR-493-5p) (Table 3) have been reported in previous studies (He et al., 2016; Cong et al., 2017), indicating that T. gondii effect on host miRNAs is not strictly tissue-specific.

The top 10 most differentially abundant miRNAs in chronic infection, as shown in Table 3, were the most overexpressed (3–4-fold increase) in chronic infection. The mmu-miR-155-5p was shown to be abundantly expressed in immune cells involved in mouse inflammatory responses, through inhibiting SOCS1 to enhance IFN-α (Rao et al., 2014) and decrease TNF-α (Huck et al., 2017). Also, the upregulation of mmu-miR-146a-5p is relevant to myeloid cell proliferation, inflammation and cancer (Boldin et al., 2011). The expression of the immunomodulatory mmu-miR-155-5p and mmu-miR-146a-5p was shown to be strain-specific and miR-146a ablation can affect early parasite burden, leading to significant differences in IFN-γ production and better survival in C57BL/6 mice (Cannella et al., 2014). Therefore, downregulation of miR-146a in acute infection has probably reduced parasite burden during early infection in order to limit CNS colonization and parasite cyst burden. mmu-miR-142a-3p, mmu-miR-142a-5p, and mmu-miR-142b, were also over-expressed. Of note, mmu-miR-142, plays significant roles in immune regulation (Sun et al., 2011, 2013, 2015), and in the development of nasopharynx carcinoma, cell cycling or IL-6 modulation in hematopoietic cell tissues, mature T cell proliferation and endotoxin-induced mortality in inflammatory processes; However, the roles of the miRNAs described above except mmu-miR-155-5p and mmu-miR-146a-5p in the process of T. gondii colonization remain to be elucidated.

KEGG pathway enrichment analysis based on the predicted target genes revealed that dysregulated miRNAs in the brain of infected mice was dominated by an immunological signature, in agreement with the results obtained in T. gondii brain infection with tissue cysts (Xu et al., 2013). We found that 255 biological pathways, especially immune-related and disease-related signaling pathways, were abundant among the significantly enriched mRNAs. The dominance of the immune-related signaling pathways indicates that T. gondii oocysts establish brain infection by altering host immune response in order to persist in the brain by forming a latent infection (Hunter and Sibley, 2012).

Previous studies suggested that miRNAs do not only regulate host cell development (Taganov et al., 2006; Kim et al., 2008), but were also involved in cancer or infection (O’Connell et al., 2007; Triboulet et al., 2007; Xiao and Rajewsky, 2009). Therefore, we studied the altered pathways of the predicted gene targets of the differentially expressed miRNAs (P < 0.05) in acute and chronic infection. Although the dysregulated miRNAs induced by early T. gondii infection were enriched in both disease- and signaling-related pathways (Figure 4A), the predicted targets of the upregulated or downregulated miRNAs were not significantly enriched (corrected P > 0.05). During chronic infection and “chronic vs. acute,” more disease-related pathways were involved in various physiological systems and molecular processes (Figures 4B,C). Interestingly, the targets of the upregulated miRNAs were mainly involved in cancer-related pathways (corrected P ≤ 0.05), such as pathways in cancer, chronic myeloid leukemia, and thyroid cancer; but the targets of downregulated miRNA were mainly enriched in signaling pathways (corrected P ≤ 0.05), such as Fc epsilon RI signaling pathway, cAMP signaling pathway and sphingolipid signaling pathway. These results further substantiate the important roles of miRNA in modulating brain response of mice challenged with T. gondii oocysts especially during chronic stage of infection.

Conclusion

The present study provides a comparative analysis of RNA sequencing-based miRNA differential expression patterns in mouse brain infected with T. gondii oocysts during acute and chronic infection. We identified 1,352 known miRNAs and 150 novel miRNAs. Three miRNAs were identified in mouse brain during acute infection and 38 miRNAs during chronic infection. These miRNAs are involved in the regulation of disease-related pathways and immune signaling pathways. Validation of seven differentially regulated miRNAs by qRT-PCR confirmed the validity of the alteration of their expression using RNA-seq analysis. Future studies are needed to further examine the miRNA changes in host brain and to investigate the functional roles of these differently expressed miRNAs, which may contribute to the study of regulatory signaling networks involved in the development of cerebral toxoplasmosis and may provide targets for therapeutic interventions.

Availability of Data and Materials

The RNA-seq data reported in the present study have been submitted to NCBI SRA database (accession number PRJNA418218). All other data supporting the findings, are available within the paper and its Supplementary Information Files.

Author Contributions

X-QZ and WC conceived and designed the experiments. R-SH, J-JH, F-KZ, YZ, and WC performed the experiments. R-SH, J-JH, HE, and WC contributed reagents, materials and analysis tools. R-SH and WC analyzed the data and wrote the paper. HE, WC, G-HZ, and X-QZ critically revised the manuscript. All the authors read and approved the final version of the manuscript.

Funding

Project support was kindly provided by the International Science and Technology Cooperation Project of Gansu Provincial Key Research and Development Program (Grant No. 17JR7WA031), the National Natural Science Foundation of China (Grant No. 31230073), the Elite Program of Chinese Academy of Agricultural Sciences, and the Agricultural Science and Technology Innovation Program (ASTIP) (Grant No. CAAS-ASTIP-2016-LVRI-03).

Conflict of Interest Statement

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.

Acknowledgments

We thank Beijing Novogene Biotechnology Co. Ltd. for technical assistance.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.02316/full#supplementary-material

FIGURE S1 | A representative Toxoplasma gondii cyst detected in the homogenized mouse brain at 33 dpi.

FIGURE S2 | List of the commonly expressed and specifically expressed miRNAs (P < 0.05) between the various mouse groups.

TABLE S1 | Sequencing data of the known miRNAs during acute and chronic infection with T. gondii oocysts.

TABLE S2 | Sequencing data of the novel miRNAs during acute and chronic infection with T. gondii oocysts.

TABLE S3 | The total differentially expressed miRNAs (P < 0.05) during acute and chronic infection with T. gondii oocysts.

TABLE S4 | List of the differentially expressed miRNAs in mouse brain following infection with oocysts and tissue cysts of T. gondii (PRU strain, type II).

Abbreviations

AIDS, acquired immune deficiency syndrome; cAMP, cyclic adenosine monophosphate; cDNA, complementary deoxyribonucleic acid; CNS, central nervous system; CsCl, cesium chloride; DNA, deoxyribonucleic acid; dpi, days post-infection; FDR, false discovery rate; GO, gene ontology; IFN-α, interferon alpha; IFN-γ, interferon gamma; KEGG, Kyoto encyclopedia of genes and genomes; MAPK, mitogen-activated protein kinase; miRNA, micro-ribonucleic acid; PBS, phosphate buffered saline; PCR-RFLP, polymerase chain reaction-restriction fragment length polymorphism; qRT-PCR, quantitative reverse transcription polymerase chain reaction; RINs, ribonucleic acid integrity numbers; RNA, ribonucleic acid; RNA-seq, deep-sequencing analysis of small ribonucleic acids; rRNA, ribosomal ribonucleic acid; SD, standard deviation; snoRNA, small nucleolar ribonucleic acid; snRNA, small nuclear ribonucleic acid; SOCS1, suppressor of cytokine signaling 1; sRNA, small ribonucleic acid; TNF-α, tumor necrosis factor alpha; TPM, transcript per million; tRNA, transfer ribonucleic acid.

Footnotes

1. ^http://www.mirbase.org/ftp.shtml
2. ^http://rfam.sanger.ac.uk/search/
3. ^http://genie.weizmann.ac.il
4. ^http://www.genome.jp/kegg/

References

Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11:R106. doi: 10.1186/gb-2010-11-10-r106

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartel, D. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116, 281–297. doi: 10.1016/S0092-8674(04)00045-5

CrossRef Full Text | Google Scholar

Betel, D., Wilson, M., Gabow, A., Marks, D. S., and Sander, C. (2008). The microRNA.org resource: targets and expression. Nucleic Acids Res. 36, D149–D153. doi: 10.1093/nar/gkm995

PubMed Abstract | CrossRef Full Text | Google Scholar

Black, M. W., and Boothroyd, J. C. (2000). Lytic cycle of Toxoplasma gondii. Microbiol. Mol. Biol. Rev. 64, 607–623. doi: 10.1128/MMBR.64.3.607-623.2000

CrossRef Full Text | Google Scholar

Blanchard, N., Dunay, I. R., and Schlüter, D. (2015). Persistence of Toxoplasma gondii in the central nervous system: a fine-tuned balance between the parasite, the brain and the immune system. Parasite Immunol. 37, 150–158. doi: 10.1111/pim.12173

PubMed Abstract | CrossRef Full Text | Google Scholar

Boldin, M. P., Taganov, K. D., Rao, D. S., Yang, L., Zhao, J. L., Kalwani, M., et al. (2011). miR-146a is a significant brake on autoimmunity, myeloproliferation, and cancer in mice. J. Exp. Med. 208, 1189–1201. doi: 10.1084/jem.20101823

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyer, K., Hill, D., Mui, E., Wroblewski, K., Karrison, T., Dubey, J. P., et al. (2011). Unrecognized ingestion of Toxoplasma gondii oocysts leads to congenital toxoplasmosis and causes epidemics in North America. Clin. Infect. Dis. 53, 1081–1089. doi: 10.1093/cid/cir667

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, P., Gobert, G. N., You, H., Duke, M., and McManus, D. P. (2015). Circulating miRNAs: potential novel biomarkers for hepatopathology progression and diagnosis of Schistosomiasis japonica in two murine models. PLoS Negl. Trop. Dis. 9:e0003965. doi: 10.1371/journal.pntd.0003965

PubMed Abstract | CrossRef Full Text | Google Scholar

Cannella, D., Brenier-Pinchart, M. P., Braun, L., van Rooyen, J. M., Bougdour, A., Bastien, O., et al. (2014). miR-146a and miR-155 delineate a MicroRNA fingerprint associated with Toxoplasma persistence in the host brain. Cell Rep. 6, 928–937. doi: 10.1016/j.celrep.2014.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Cong, W., Dong, X. Y., Meng, Q. F., Zhou, N., Wang, X. Y., Huang, S. Y., et al. (2015). Toxoplasma gondii infection in pregnant women: a seroprevalence and case-control study in eastern China. Biomed. Res. Int. 2015:170278. doi: 10.1155/2015/170278

PubMed Abstract | CrossRef Full Text | Google Scholar

Cong, W., Qin, S. Y., Meng, Q. F., Zou, F. C., Qian, A. D., and Zhu, X. Q. (2016). Molecular detection and genetic characterization of Toxoplasma gondii infection in sika deer (Cervus nippon) in China. Infect. Genet. Evol. 39, 9–11. doi: 10.1016/j.meegid.2016.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Cong, W., Zhang, X. X., He, J. J., Li, F. C., Elsheikha, H. M., and Zhu, X. Q. (2017). Global miRNA expression profiling of domestic cat livers following acute Toxoplasma gondii infection. Oncotarget 8, 25599–25611. doi: 10.18632/oncotarget.16108

PubMed Abstract | CrossRef Full Text | Google Scholar

Dubey, J. P. (2004). Toxoplasmosis-a waterborne zoonosis. Vet. Parasitol. 126, 57–72. doi: 10.1016/j.vetpar.2004.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Dubey, J. P. (2008). The history of Toxoplasma gondii-the first 100 years. J. Eukaryot. Microbiol. 55, 467–475. doi: 10.1111/j.1550-7408.2008.00345.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Elsheikha, H. M. (2008). Congenital toxoplasmosis: priorities for further health promotion action. Public Health 122, 335–353. doi: 10.1016/j.puhe.2007.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Elsheikha, H. M., and Zhu, X. Q. (2016). Toxoplasma gondii infection and schizophrenia: an inter-kingdom communication perspective. Curr. Opin. Infect. Dis. 29, 311–318. doi: 10.1097/QCO.0000000000000265

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedländer, M. R., Mackowiak, S. D., Li, N., Chen, W., and Rajewsky, N. (2012). miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 40, 37–52. doi: 10.1093/nar/gkr688

PubMed Abstract | CrossRef Full Text | Google Scholar

Hakimi, M. A., and Ménard, R. (2010). Do apicomplexan parasites hijack the host cell microRNA pathway for their intracellular development? F1000 Biol. Rep. 2:42. doi: 10.3410/B2-42

PubMed Abstract | CrossRef Full Text | Google Scholar

He, J. J., Ma, J., Wang, J. L., Xu, M. J., and Zhu, X. Q. (2016). Analysis of miRNA expression profiling in mouse spleen affected by acute Toxoplasma gondii infection. Infect. Genet. Evol. 37, 137–142. doi: 10.1016/j.meegid.2015.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Hill, D., Coss, C., Dubey, J. P., Wroblewski, K., Sautter, M., Hosten, T., et al. (2011). Identification of a sporozoite-specific antigen from Toxoplasma gondii. J. Parasitol. 97, 328–337. doi: 10.1645/GE-2782.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Hobert, O. (2008). Gene regulation by transcription factors and microRNAs. Science 319, 1785–1786. doi: 10.1126/science.1151651

PubMed Abstract | CrossRef Full Text | Google Scholar

Howe, D. K., and Sibley, L. D. (1995). Toxoplasma gondii comprises three clonal lineages: correlation of parasite genotype with human disease. J. Infect. Dis. 172, 1561–1566. doi: 10.1093/infdis/172.6.1561

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoy, A. M., Lundie, R. J., Ivens, A., Quintana, J. F., Nausch, N., Forster, T., et al. (2014). Parasite-derived microRNAs in host serum as novel biomarkers of helminth infection. PLoS Negl. Trop. Dis. 8:e2701. doi: 10.1371/journal.pntd.0002701

PubMed Abstract | CrossRef Full Text | Google Scholar

Huck, O., Al-Hashemi, J., Poidevin, L., Poch, O., Davideau, J. L., Tenenbaum, H., et al. (2017). Identification and characterization of microRNA differentially expressed in macrophages exposed to Porphyromonas gingivalis infection. Infect. Immun. 85:e00771-16. doi: 10.1128/IAI.00771-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Hunter, C. A., and Sibley, L. D. (2012). Modulation of innate immunity by Toxoplasma gondii virulence effectors. Nat. Rev. Microbiol. 10, 766–778. doi: 10.1038/nrmicro2858

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, J. L., and Dubey, J. P. (2012). Foodborne toxoplasmosis. Clin. Infect. Dis. 55, 845–851. doi: 10.1093/cid/cis508

PubMed Abstract | CrossRef Full Text | Google Scholar

Judice, C. C., Bourgard, C., Kayano, A. C., Albrecht, L., and Costa, F. T. (2016). MicroRNAs in the host-apicomplexan parasites interactions: a review of immunopathological aspects. Front. Cell Infect. Microbiol. 6:5. doi: 10.3389/fcimb.2016.00005

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Araki, M., Goto, S., Hattori, M., Hirakawa, M., Itoh, M., et al. (2008). KEGG for linking genomes to life and the environment. Nucleic Acids Res. 36, D480–D484. doi: 10.1093/nar/gkm882

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D. H., Saetrom, P., Snøve, O. Jr., and Rossi, J. J. (2008). MicroRNA-directed transcriptional gene silencing in mammalian cells. Proc. Natl. Acad. Sci. U.S.A. 105, 16230–16235. doi: 10.1073/pnas.0808830105

PubMed Abstract | CrossRef Full Text | Google Scholar

Langmead, B., Trapnell, C., Pop, M., and Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10:R25. doi: 10.1186/gb-2009-10-3-r25

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindsay, D. S., and Dubey, J. P. (2009). Long-term survival of Toxoplasma gondii sporulated oocysts in seawater. J. Parasitol. 95, 1019–1020. doi: 10.1645/GE-1919.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Mao, X., Cai, T., Olyarchuk, J. G., and Wei, L. (2005). Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 21, 3787–3793. doi: 10.1093/bioinformatics/bti430

PubMed Abstract | CrossRef Full Text | Google Scholar

Marra, C. M. (2018). Central nervous system infection with Toxoplasma gondii. Handb. Clin. Neurol. 152, 117–122. doi: 10.1016/B978-0-444-63849-6.00009-8

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Connell, R. M., Taganov, K. D., Boldin, M. P., Cheng, G., and Baltimore, D. (2007). MicroRNA-155 is induced during the macrophage inflammatory response. Proc. Natl. Acad. Sci. U.S.A. 104, 1604–1609. doi: 10.1073/pnas.0610731104

PubMed Abstract | CrossRef Full Text | Google Scholar

Peltier, H. J., and Latham, G. J. (2008). Normalization of microRNA expression levels in quantitative RT-PCR assays: identification of suitable reference RNA targets in normal and cancerous human solid tissues. RNA 14, 844–852. doi: 10.1261/rna.939908

PubMed Abstract | CrossRef Full Text | Google Scholar

Pittman, K. J., and Knoll, L. J. (2015). Long-term relationships: the complicated interplay between the host and the developmental stages of Toxoplasma gondii during acute and chronic infections. Microbiol. Mol. Biol. Rev. 79, 387–401. doi: 10.1128/MMBR.00027-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Rao, R., Rieder, S. A., Nagarkatti, P., and Nagarkatti, M. (2014). Staphylococcal enterotoxin B-induced microRNA-155 targets SOCS1 to promote acute inflammatory lung injury. Infect. Immun. 82, 2971–2979. doi: 10.1128/IAI.01666-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Staggs, S. E., See, M. J., Dubey, J. P., and Villegas, E. N. (2009). Obtaining highly purified Toxoplasma gondii oocysts by a discontinuous cesium chloride gradient. J. Vis. Exp. 33:e1420. doi: 10.3791/1420

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y., Oravecz-Wilson, K., Mathewson, N., Wang, Y., McEachin, R., Liu, C., et al. (2015). Mature T cell responses are controlled by microRNA-142. J. Clin. Invest. 125, 2825–2840. doi: 10.1172/JCI78753

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y., Sun, J., Tomomi, T., Nieves, E., Mathewson, N., Tamaki, H., et al. (2013). PU.1-dependent transcriptional regulation of miR-142 contributes to its hematopoietic cell-specific expression and modulation of IL-6. J. Immunol. 190, 4005–4013. doi: 10.4049/jimmunol.1202911

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y., Varambally, S., Maher, C. A., Cao, Q., Chockley, P., Toubai, T., et al. (2011). Targeting of microRNA-142-3p in dendritic cells regulates endotoxin-induced mortality. Blood 117, 6172–6183. doi: 10.1182/blood-2010-12-325647

PubMed Abstract | CrossRef Full Text | Google Scholar

Taganov, K. D., Boldin, M. P., Chang, K. J., and Baltimore, D. (2006). NF-kappaB-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses. Proc. Natl. Acad. Sci. U.S.A. 103, 12481–12486. doi: 10.1073/pnas.0605298103

PubMed Abstract | CrossRef Full Text | Google Scholar

Tenter, A. M., Heckeroth, A. R., and Weiss, L. M. (2000). Toxoplasma gondii: from animals to humans. Int. J. Parasitol. 30, 1217–1258. doi: 10.1016/S0020-7519(00)00124-7

CrossRef Full Text | Google Scholar

Torrey, E. F., and Yolken, R. H. (2013). Toxoplasma oocysts as a public health problem. Trends Parasitol. 29, 380–384. doi: 10.1016/j.pt.2013.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Triboulet, R., Mari, B., Lin, Y. L., Chable-Bessia, C., Bennasser, Y., Lebrigand, K., et al. (2007). Suppression of microRNA-silencing pathway by HIV-1 during virus replication. Science 315, 1579–1582. doi: 10.1126/science.1136319

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, M., Shen, Y., Shi, S., and Tang, T. (2012). miREvo: an integrative microRNA evolutionary analysis platform for next-generation sequencing experiments. BMC Bioinformatics 13:140. doi: 10.1186/1471-2105-13-140

PubMed Abstract | CrossRef Full Text | Google Scholar

Winter, J., Jung, S., Keller, S., Gregory, R. I., and Diederichs, S. (2009). Many roads to maturity: microRNA biogenesis pathways and their regulation. Nat. Cell Biol. 11, 228–234. doi: 10.1038/ncb0309-228

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, C., and Rajewsky, K. (2009). MicroRNA control in the immune system: basic principles. Cell 136, 26–36. doi: 10.1016/j.cell.2008.12.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, M. J., Zhou, D. H., Nisbet, A. J., Huang, S. Y., Fan, Y. F., and Zhu, X. Q. (2013). Characterization of mouse brain microRNAs after infection with cyst-forming Toxoplasma gondii. Parasit. Vect. 6:154. doi: 10.1186/1756-3305-6-154

PubMed Abstract | CrossRef Full Text | Google Scholar

Young, M. D., Wakefield, M. J., Smyth, G. K., and Oshlack, A. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11:R14. doi: 10.1186/gb-2010-11-2-r14

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeiner, G. M., Norman, K. L., Thomson, J. M., Hammond, S. M., and Boothroyd, J. C. (2010). Toxoplasma gondii infection specifically increases the levels of key host microRNAs. PLoS One 5:e8742. doi: 10.1371/journal.pone.0008742

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, L., Chen, J., Li, Z., Li, X., Hu, X., Huang, Y., et al. (2010). Integrated profiling of microRNAs and mRNAs: microRNAs located on Xq27.3 associate with clear cell renal cell carcinoma. PLoS One 5:e15224. doi: 10.1371/journal.pone.0015224

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, P., Chen, Z., Li, H. L., Zheng, H., He, S., Lin, R. Q., et al. (2011). Toxoplasma gondii infection in humans in China. Parasit. Vect. 4:165. doi: 10.1186/1756-3305-4-165

PubMed Abstract | CrossRef Full Text | Google Scholar