Skip to main content

miRNA and circRNA expression patterns in mouse brain during toxoplasmosis development



Increasing evidence has shown that circular RNAs (circRNAs) are involved in neurodegenerative disorders, but their roles in neurological toxoplasmosis are yet to know. This study examined miRNA and circRNA expressions in mouse brain following oral infection with T. gondii Pru strain.


Total RNA extracted from acutely infected (11 days post infection (DPI)), chronically infected (35 DPI) and uninfected mouse brain samples were subjected to genome-wide small RNA sequencing. In the acutely infected mice, 9 circRNAs and 20 miRNAs were upregulated, whereas 67 circRNAs and 28 miRNAs were downregulated. In the chronically infected mice, 2 circRNAs and 42 miRNAs were upregulated, whereas 1 circRNA and 29 miRNAs were downregulated. Gene ontology analysis predicted that the host genes that produced the dysregulated circRNAs in the acutely infected brain were primarily involved in response to stimulus and ion binding activities. Furthermore, predictive interaction networks of circRNA-miRNA and miRNA-mRNA were constructed based on genome-wide transcriptome sequencing and computational analyses, which might suggest the putative functions of miRNAs and circRNAs as a large class of post-transcriptional regulators.


These findings will shed light on circRNA-miRNA interactions during the pathogenesis of toxoplasmosis, and they will lay solid foundation for studying the potential regulation roles of miRNAs and circRNAs in T. gondii induced pathogenesis.


Toxoplasma gondii, the causative agent of toxoplasmosis, is a medically important parasite infecting approximately 30% of the world’s population [1]. Immune-potent individuals usually do not show any clinical symptoms, and infection is latent and maintained inside the host as tissue cysts. When infected persons become immunocompromised, such as people suffering from AIDS or organ transplantation, parasite cyst reactivation can lead to acute toxoplasmosis, which includes eye disease, neurological problems and even death [2]. Additionally, infection of pregnant women can cause congenital infection, with dysplasia, hydrocephaly and chorioretinitis occurring in the newborns [3].

The fast proliferating T. gondii tachyzoite is fully controlled by the host potent immune response and then transforms into a slowly replicating stage (i.e. bradyzoite) enclosed in tissue cysts, and remains dormant within the central nervous system (CNS) and muscles [4]. Various kinds of brain cells, including microglia, astrocytes and neurons, can be infected. In the chronically infected mice, the parasite cysts were mainly found in the neurons [5]. In a study of congenital toxoplasmosis, cysts were also found in the neurons [6]. In toxoplasma encephalitis (TE), the rostral basal ganglion is the most commonly affected region, followed by the cerebellum and brain stem [7, 8]. In addition, it has been shown that T. gondii infection can affect levels of some neurotransmitters, such as dopamine [9]. Positive correlations between T. gondii infections and neuropsychiatric disorders like schizophrenia [10, 11], cryptogenic epilepsy [12] and Parkinson’s disease (PD) [13, 14] were revealed in the clinical and epidemiological investigations. However, the mechanisms underlying Toxoplasma-induced neuronal disorder in the brain are poorly known.

MicroRNAs (miRNAs) are endogenous small non-coding RNA molecules which regulate gene expression at the post-transcriptional level, and are now widely recognized as essential regulatory molecules involved in neuronal development and function. miRNAs are abundant in the nervous system and a link between miRNAs and the occurrence and development of neurological diseases, such as cerebral ischemia, stroke and neurodegenerative diseases, is becoming increasingly clear [15]. Notably, miRNAs are also found to be regulators of the host response to Toxoplasma infection, as miR-146a and miR-155 are highly induced in chronically infected mice brain, but worthy to mention, compared to miR-146a−/− mice, wild type mice show a more severe TE after challenge with genotype II Pru strain [16]. In addition to using host miRNAs for parasite persistence, T. gondii infection also modulates host cell miRNA profiles in a timely manner [17]. Another type of non-coding RNA, circular RNA (circRNA), is formed by exon-scrambling and was once largely neglected due to its rarity and lack of biological functions. With the development of high-throughput sequencing technology, thousands of new circRNAs have been identified in various organisms, from microorganisms to mammals [18]. circRNA is found to antagonize miRNA activity by a sponge-like mechanism and thus regulates gene expression at the post-transcriptional level [19]. Several pieces of evidence have shown that circRNAs are most highly and specifically expressed in neural tissues, and some circRNAs are expressed in a specific spatial and temporal pattern in the brain among species. More studies have shown that circRNAs may act as important regulatory factors in the pathogenesis of neurological diseases, such as epilepsy, PD and Alzheimer’s disease (AD) [20]. However, their roles in the occurrence and development processes of neurological toxoplasmosis remain largely unexplored.

To investigate the potential roles of miRNAs and circRNAs in T. gondii infection, we performed a high-throughput sequencing analysis to identify differentially expressed miRNAs and circRNAs in mice after acute and chronic T. gondii infections, respectively. Furthermore, differentially expressed miRNAs regulated by differentially expressed circRNAs were predicted, and predictive interaction analysis of miRNA-mRNA was also performed. Thus, the present study presented the first transcriptome-wide circRNA landscape and a predictive miRNA-circRNA network in T. gondii-infected mice, which provides a valuable dataset that will help elucidate the mechanisms underlying neurological disorders during T. gondii infection.


Expression of miRNAs in brains from both Toxoplasma-infected and normal mice

Infected mice with T. gondii at 11 DPI showed noticeable acute toxoplasmosis signs, including anorexia, hyperpyrexia, and messy hair. All mice survived and restored their physical status at 35 DPI. Initially, we set out to profile the miRNA expression across samples that represent typical infection states. miRNA libraries of brain samples were successfully constructed and summarized in Additional file 1: Table S1. Genome-wide sequencing identified more than 12 million clean reads in each sample. As shown in Fig. 1, the dominant small RNAs (sRNAs) were 20–24 nt in length. More than 93% sRNA sequences were mapped to the reference genome, and the majority of these mapped sRNA sequences were in the forward orientation (Additional file 2: Table S2). Distribution of sRNAs sequences showed that the number of miRNA sequences was higher in the control group (Additional file 3: Table S3). The mapped reads were divided into different categories, and the majority of small RNA reads were belonged to the annotated known miRNAs.

Fig. 1
figure 1

The average length distribution of the sRNA sequences identified in each group. Sample groups including acutely infected, chronically infected and control are labeled as AI, CI, and Con, respectively

We then compared acute infection samples against normal samples and identified 20 over-expressed and 28 under-expressed miRNAs (adjusted P value < 0.05) (Fig. 2a). Comparison between chronically infected and control samples revealed more dys-regulated miRNAs, including 42 over-expressed and 29 under-expressed miRNAs. As shown in Fig. 2b, 24 miRNAs were shared between AI vs. Con and CI vs. Con. Worthy to mention, all these 24 miRNAs showed the same variation tendency during parasite infection (Table 1). The clustering of miRNA expression profiles derived from 228 different miRNAs in control and infected samples is shown in Fig. 2c. The tree shows a very good separation between normal and infected samples. However, PCA scores plots did not clearly discriminate chronically infected mice from acutely infected group (Additional file 4: Figure S1A). A recent study showed that infection with T. gondii oocysts caused the miRNA perturbation in mice brains [21]. The high-throughput data are available in the NCBI database with the accession number PRJNA418218. As shown in Table 2, the differentially expressed miRNAs shared by both cyst and oocyst infections are listed. It is noteworthy that, mmu-miR-155-5p were up-regulated during the whole infection course in both cyst and oocyst infections, which was further confirmed by the qPCR (Additional file 5: Figure S2).

Fig. 2
figure 2

Comparison among normal and parasite infected samples reveals global changes in miRNA expression. a Statistics of differentially expressed miRNA among samples. b Venn diagram representation of differentially expressed circRNA sharing and exclusive constitutively presented among different time points vs. Control sample. c Hierarchical clustering of miRNA expression. miRNA profiles of 9 samples from different groups were clustered. Samples are in columns and miRNA in rows. Sample groups including acutely infected, chronically infected and control, are labeled as AI, CI, and Con, respectively. Brown and blue indicate higher and lower abundance, respectively

Table 1 The differentially expressed miRNAs shared by both the acute and chronic infections
Table 2 The differentially expressed miRNAs shared by both parasite cyst and oocyst infections

Mapping circRNAs from total RNA-seq data

To investigate whether T. gondii infection affects the expression of the host circRNAs, we profiled the brain circRNAs during the toxoplasmosis progression. To obtain sufficient RNA-Seq data, non- polyA-containing RNA, linear RNA and rRNA were removed prior to cDNA library preparation. Tested samples were sequenced on an Illumina HiSeq 4000 with an average data volume of 9 Gb, thereby yielding a total of more than 557 million paired-end reads with a size of 150 bp (Additional file 6: Table S4). Clean reads were mapped to the mouse reference genome (GRCm38/p6) by using Bowtie, and about 61% of these reads consisted of protein coding sequences, whereas smaller fractions aligned with ribozyme and long noncoding RNAs (Fig. 3a). circRNAs were identified and filtered by using both the find_circ and CIRI2 software as previously described [22]. A total of 16,543 circRNAs with a minimum of two reads spanning the back-splicing junction were identified. Annotations for all circRNAs identified in this study are shown in Additional file 7: Table S5. We classified these circRNAs into three groups, namely, exonic circRNAs, intergenic circRNAs and intronic circRNAs. As shown in Fig. 3b, in all of the 9 sample data sets, exonic circRNAs predominated. The size of these circRNA candidates ranged from less than 50 nt to greater than 1500 nt, most of which ranged from 100 nt to 650 nt (Fig. 3c). As shown in Fig. 3d, most genes generated one or two circRNAs, but some genes yielded multiple circRNAs.

Fig. 3
figure 3

Features of brain circRNAs. a Clean reads mapping. b The length distribution for identified circRNAs (intergenic circRNAs are excluded). c Genome origin of identified circRNAs. d Distribution of circRNAs among genes

Differentially expressed circRNAs in mice brain during T. gondii infection

We then performed expression profiling to show circRNA variations in mice brain during parasite infection. Expression abundance of circRNAs from nine samples was measured based on TPM, and no abnormal expression was found (data not shown). As a result, compared to control group, 76 circRNAs in acute infection group were detected to be differentially regulated by |log2 fold change| ≥0.5, adjusted P < 0.05, among which 9 circRNAs were up-regulated while 67 circRNAs were down-regulated. In the chronic infection group, only 3 differentially expressed circRNAs were selected with 2 up-regulated and 1 down-regulated (Fig. 4a). Among these dys-regulated circRNAs, only one was common between AI vs. Con and CI vs. Con, namely novel_circ_0057684 (Fig. 4b). Four circRNAs were selected randomly for validations. As shown in Additional file 8: Figure S3, the results are consistent with predictions. Hierarchical clustering showed that circRNA expression pattern in acute infection group was distinguishable compared to the other two groups (Fig. 4c). Meanwhile, PCA scores plots did not clearly differentiate chronically infected mice from control group (Additional file 4: Figure S1B). These data suggested that the expression of circRNAs in acutely infected brain was severely affected.

Fig. 4
figure 4

Global view of temporal circRNA expression profiles in mice brain during T. gondii infection. a Statistics of differentially expressed circRNA among samples. b Venn diagram representation of differentially expressed circRNA sharing and exclusive constitutively presented among different time points vs. Control sample. c Hierarchical clustering of circRNA expression data. Unsupervised hierarchical clustering was performed using Pearson′s correlation as distance measure and Average-linkage for linkage analysis. Brown and blue indicate higher and lower abundance, respectively. Sample groups including acutely infected, chronically infected and control are labeled as AI, CI, and Con, respectively

Functional enrichment analysis of differentially expressed circRNAs

Comparison of differentially regulated circRNAs and their corresponding host genes between the acutely infected and normal brains are shown in Additional file 9: Table S6. We performed Gene Ontology (GO) analysis on host genes. GO analysis consists of three different aspects, namely biological process, cellular component and molecular function (Fig. 5). Prediction terms with P-value less than 0.05 were selected and ranked by enrichment score (−log2 (P-value)), and the top 10 generally affected GO terms in each categories are listed. The most enriched biological process terms were related to cellular process, signal transduction and response to stimulus, such as “response to stimulus (GO:0050896)” and “single-organism cellular process (GO:0044763)”. The most enriched molecular function terms were mostly about binding activity and kinase activity, such as “ion binding (GO: 0043167)” and “transmembrane receptor protein kinase activity (GO: 0019199)”. As for cell component, the top three enriched terms were “integral to plasma membrane (GO: 0005887)”, “endosome (GO: 0005768)” and “intrinsic to plasma membrane (GO: 0031226)”.

Fig. 5
figure 5

Infection-dysregulated circRNA host gene Gene Ontology (GO) analysis. GO analysis was performed on host genes of significantly dys-regulated (P value < 0.05) in the comparison between the acute infection and normal groups. Enriched terms are grouped by GO category-biological process (a), molecular function (b) and cell component (c). The y axis on the left represents the GO term. The number of the input host genes is listed besides the bar. P values indicate the significance of enrichment of the input host genes for each GO Term

Analyses of of circRNA-miRNA and miRNA-mRNA interactions

circRNAs are a unique class of endogenous non-coding RNAs which act as miRNA decoys or sponges to regulate gene expression. Therefore, an integrated analysis of the expression profile from circRNA-miRNA predictive interactions was performed based on the high-throughput RNA sequencing results. First, using miRanda software, the binding site analysis and predictive interaction analysis of circRNAs-miRNAs identified in the comparison between acutely infected and normal brains were performed. As shown in Fig. 6, 38 dys-regulated miRNAs and 50 dys-regulated circRNAs were involved in this network. The data used to create Fig. 6 are listed in Additional file 10: Table S7. Down-regulated mmu-miR-214-3p is the most frequently targeted miRNA by 11 circRNAs, followed by up-regulated mmu-miR-21a-3p, down-regulated mmu-miR-455-3p and mmu-miR-497a-5p. It is worth noting that the novel_circ_0048152 showed the largest interaction network. All the above five RNAs were selected to validate the expression profiles obtained by RNA-Seq, and the qPCR results were consistent with the RNA-Seq findings (Additional file 11: Figure S4). Next, predictive miRNA-mRNA interacting pairs were ranked by the P value of the hypergeometric distribution. As shown in Table 3, three relationships between miRNA and targeted mRNA (|fold change|1.5, P < 0.05) in the acute infection stage and five relationships in chronic infection stage were revealed, respectively. As shown in Additional file 12: Figure S5, the expression levels of these differentially expressed mRNAs were consistent with the high-throughput sequencing data.

Fig. 6
figure 6

CircRNA-miRNA interaction network. The network consists of 50 circRNAs and 38 miRNAs identified between the acutely infected and normal mouse brains. Triangle indicates circRNA and circle indicates miRNA. The red node indicates upregulated or green node indicates downregulated RNAs (log2 fold change). The node size indicates high interaction degree (large) or low degree (small). RNAs that are associated to each other are linked by an edge. The color of the edge indicates miRanda combined interaction score for that particular circRNA-miRNA pair

Table 3 Pairwise co-expression analysis of all profiled miRNAs and mRNAs during the whole T.gondii infection


A previous study indicated that T. gondii infection could modulate the behavior of the intermediate hosts with the presence of tissue cysts in the brain, and the infected rodents showed impaired learning and memory [23]. In humans, increased rates of psychiatric disorders and suicide, and decreased psychomotor performance have been associated with persistent T. gondii infection [24, 25]. In order to reveal the mechanisms underlying neuronal disorder, we need to better understand brain gene expression and regulation induced by T. gondii infection.

Previous studies have shown that more genes were differentially expressed during chronic infection compared to acute infection after type II strain infection [26, 27]. Among these differentially expressed genes, more genes were up-regulated and were primarily involved in immune regulation and cell activation. The fast development of RNA-seq technology offers unique opportunities to discover novel regulating factors in toxoplasmosis pathogenesis. In this study, we comprehensively profiled the miRNA and circRNA expressions of mice brains in order to further our understanding of the disease pathogenesis. Using RNA-seq data, we revealed a global bias for miRNA and circRNA accumulation in T. gondii infected mouse brains.

miRNAs, widely spread in animal cells, can bind to the complementary site on the 3′ untranslated region (UTR) of the targeting mRNAs and thus facilitate mRNA degradation or translation inhibition on a post-transcriptional level. In a previous study, 637 miRNAs were identified in Kunming mouse brains at 14 and 21 days post infection with T. gondii cysts, which is a phase of rehabilitation [17]. However, the miRNAs expression has never been profiled in the acute and chronic infection stages after T. gondii infection. In this study, we identified 1314 known miRNAs and 91 novel miRNAs, of which 48 and 71 differentially expressed miRNAs were found in acute and chronic infection, respectively. Mmu-miR-155-5p is the most up-regulated miRNA during the whole infection course, which also showed significantly perturbed expression in oocyst-induced toxoplasmosis [21]. mmu-miR-155-5p showed an elevated level in chronic infection compared to acute infection, which might be resulted from the increased high chronic cyst burdens. Several studies provide evidences that miR-155 (the precursor of miR-155-5p) is involved in the innate and acquired immune responses, hematopoiesis and autoimmune disorders [28,29,30]. More evidences show that elevated levels of miR-155 occur in the formation and development of several tumors, such as leukemias and breast cancer. Both the levels of mmu-miR-142a-3p and mmu-miR-142a-5p were up-regulated in this study. miR-142a-3p has been identified as an essential player in the formation and differentiation of hematopoietic stem cells, in which miR-142a-3p directly targets transcription factor p53 which is involved in the biological process of hematopoietic stem and progenitor cells [31]. Transplantation of mesenchymal stem cells (MSCs) shows good therapeutic effects on acute lung injury, during which Beclin-1 protein is up-regulated. miR-142a-5p is found to negatively regulate Beclin-1 production, and its expression is inhibited during the transplantation of MSCs [32].

Over-expression of the miR-223 significantly reduced Plasmodium falciparum infection during the intraerythrocytic life cycle [33]. In this study, both mmu-miR-223-3p and mmu-miR-223-5p were up-regulated, which might be involved in host defense against T. gondii infection. Mmu-miR-203-3p was also constantly up-regulated during the whole infection course. Interestingly, the miR-203 functions as a tumor suppressor and is down-regulated in many kinds of cancers, such as primary prostatic tumors, cervical cancer and rhabdomyosarcoma [34,35,36]. In the present study, mmu-miR-185-3p was the most down-regulated miRNA. More evidences showed that the expression level of miR-185 decreased in carcinogenesis, in which miR-185 suppresses tumor proliferation by directly targeting DNA methyltransferase 1 (DNMT1) [37, 38].

circRNAs are produced by RNA back splicing and highly prevalent in the eukaryotic transcriptome [39]. While the functions of circRNAs remain poorly known, some circRNAs regulate gene expression by acting as miRNA sponges. For instance, a heart-related circRNA (HRCR) was found to function as an endogenous miR-223 sponge, which might be involved in the protection of the heart from pathological hypertrophy and heart failure [40]. In the present study, a large number of circRNAs were identified in both infected and normal mouse brains. Most circRNAs have low abundance, and certain circRNAs are expressed in one gene locus. Interestingly, 76 and 3 differentially expressed circRNAs were identified in the acute and chronic infection stages, respectively, which suggests that these circRNAs may play key roles in toxoplasmosis pathogenesis. In order to explore the roles of these differential circRNAs in mouse toxoplasmosis, GO analysis was used to annotate the biological functions of host linear transcripts in the acute infection stage, and the results showed that a significant amount of GO terms were related with the cellular process and binding activity. circRNA-miRNA predictive interaction networks were also established to predict the relationships between circRNAs and miRNAs. Generally, a circRNA harbors more binding sites for different miRNAs, which means that the circRNA may play multiple functional roles in the regulation of miRNA target gene expression. In this study, it was found that the novel_circ_0048152 showed the largest interaction network and might potentially interact with seven miRNAs including mmu-miR-206-3p, mmu-miR-1a-3p, mmu-miR-598-3p, mmu-miR-147-3p, mmu-miR-484, mmu-miR-34a-5p and mmu-miR-21a-3p. Further biological roles of the novel_circ_0048152 still need to be explored. Moreover, we found that mmu-miR-214-3p, mmu-miR-21a-3p, mmu-miR-455-3p and mmu-miR-497a-5p were potentially co-expressed with multiple circRNAs. However, their functional roles during toxoplasmosis have never been reported.


The present study revealed unique sets of miRNAs and circRNAs along with their expression profiles. Their potential roles were predicted by bioinformatics softwares. Predictive interaction networks were constructed for both miRNA-mRNA and circRNA-miRNA. Our data lay a solid foundation for studying the functions of miRNAs and circRNAs in toxoplasmosis pathogenesis, which will be performed in further experimental studies.


Mice and parasite

Six to eight-week-old female BALB/c mice were purchased from the Experimental Animal Center of Lanzhou Veterinary Research Institute, Chinese Academy of Agriculture Sciences, PR China. Strict animal handling guidelines were followed in this study. All animals had free access to food and water. They were kept at 22 °C in a 12/12 h light/dark cycle.T. gondii type II Pru strain was used in this study, and was maintained throughout the study via oral inoculation of cysts in mice. Cysts were obtained from brain tissues of infected Kunming mice after 60 days post infection (DPI).

Mouse infection and sample collection

Nine BALB/c mice were randomly assigned to three experimental groups and each group contains three individuals. Six mice were infected orally with ~ 10 freshly prepared T. gondii cysts by gavage. Meanwhile, three mock-infected (control) mice received 100 μL saline alone. All infected mice were checked daily for the development of clinical signs of toxoplasmosis. At 11 and 35 DPI with T. gondii cysts, infected (n = 3 at each time point) and control mice (n = 3) were anesthetized with 5% isoflurane gas at around 10 a.m., and then sacrificed by cervical dislocation. The brain tissues were collected and immediately snap-frozen in liquid nitrogen and kept at − 80 °C until use.

RNA extraction

Total RNA was extracted from brain tissues using Trizol method (Invitrogen, Carlsbad, CA, USA). The purity of RNA was determined using the NanoPhotometer® spectrophotometer (IMPLEN, CA, USA). Integrity of RNA was confirmed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, USA). The concentration of RNA was quantified using Qubit® RNA Assay Kit in Qubit® 2.0 Flurometer (Life Technologies, CA, USA). RNA samples were stored at − 80 °C until used.

RNA sequencing

The detailed method for mRNA sequencing was described previously [41]. Three micrograms of the total RNA per sample was used for the construction of the small RNA library. Sequencing libraries were constructed using NEBNext® Multiplex Small RNA Library Prep Set for Illumina® (NEB, USA.) according to the manufacturer’s recommendations. The quality of sequencing library was evaluated on the Agilent Bioanalyzer 2100 system using DNA High Sensitivity Chips. Small RNA sequencing libraries were constructed with HiSeq 2500 sequencing platform and 50 bp single-end reads were generated. Raw data of fastq format were firstly processed through in-house perlscripts. Q20, Q30, and GC-content of the clean data were calculated. In order to analyze the small RNA expression and distribution on the genome, small RNA tags were mapped to reference sequence by Bowtie [42]. MiRBase 20.0 was used to align small RNA to the miRNA precursor to obtain the known miRNA. By exploring the secondary structure, the Dicer cleavage site and the minimum free energy of the unannotated small RNA tags, softwares miREvo and mirdeep2 were integrated to predict novel miRNAs [43, 44].

For the construction of circRNA library, 5 μg RNA per sample was used as input material. Ribosomal RNA was removed by Epicentre Ribozero™ rRNA Removal Kit (Epicentre, USA) followed by the linear RNA digestion with RNase R (Epicentre, USA). The sequencing libraries were constructed by NEBNext® Ultra™ Directional RNA Library Prep Kit for Illumina® (NEB, USA) following the manufacturer’s instructions. According to the manufacturer’s instructions, circRNA sequencing libraries were constructed with HiSeq 4000 sequencing platform and 150 bp paired-end reads were generated. Reference mouse genome was downloaded from the NCBI genome website directly. Index of the reference genome was built using bowtie2 v2.2.8 and paired-end clean reads were aligned to the reference genome using Bowtie. The mapped reads were removed. The unmapped reads, including splice junctions-aligned reads were further processed by find_circ and CIRI2 [45, 46]. At last, the circRNAs that contain splice junctions-aligned reads were identified. The intersection between the two algorithms for circRNA prediction was selected.

RNA quantification and differential expression analysis

The FPKM (fragments per kilo-base of exon per million fragments) of the coding genes was calculated using Cuffdiff (v.2.1.1) to evaluate the mRNA expression level in each sample [47]. The expression levels of miRNAs and circRNAs were estimated by TPM (transcripts per million reads) through the following criteria: Normalized expression = actual miRNA readcount or circRNA junction readcount× 1,000,000/libsize (libsize is the sum of miRNA or circRNA readcount) [48].

Differential expression analysis of two groups was conducted using the DESeq 2(1.14.1). The P-values were adjusted using the Benjamini & Hochberg method. Corrected P-value of 0.05 was set as the threshold for significantly differential expression by default [49]. By using metaboanalyst 4.0 online software, these differentially expressed RNAs were selected to perform hierarchical cluster analysis and principal component analysis (PCA) [50].

Quantitative PCR

Total RNA was extracted using Trizol reagent (Life Technologies). To quantify the amount of mature miRNA, we used the miRcute Plus miRNA qPCR Detection Kit (SYBR Green) (FP411, Tiangen Biotech, Beijing, China) and small nuclear U6 as an internal standard. The reactions were performed at 95 °C for 15 min, 94 °C for 20 s, 60 °C for 30 s, 72 °C for 34 s. All reactions were run in triplicate. To quantify the amount of circRNAs, cDNA was synthesized with the HiScript II Q RT SuperMix (Vazyme, Nanjing, China). Quantitative real-time PCR (qRT-PCR) was performed with a ChamQ SYBR qPCR Master Mix kit (Vazyme, Nanjing, China) and GAPDH was used as a normalization control. The reactive cycle consisted of 95 °C for 30 s, then 40 cycles of 95 °C for 5 s and 60 °C for 34 s. A set of divergent primers (the forward primer being located downstream of the reverse primer) were designed and used for backspliced circRNA detection. All PCR primers were designed using oligo7 software and are listed in Additional file 13: Table S8. The qRT-PCR was run on a Bio-Rad CFX connect Real-Time PCR Detection System (BIORAD, USA) and the qRT-PCR data were analyzed using the comparative 2−ΔΔCT method.

miRNA target gene prediction and functional analysis

The miRanda, PITA and RNAhybrid softwares were used to predict potential target genes of all the differentially expressed miRNAs. An intersection among the three combined predictions were selected. Threshold parameters for PITA were set as follows: max target length: 50,000; energy cutoff: − 10; p-value cutoff: 0.05. Parameters for RNAhybrid were set as follows: utr: 3utr.fa; mir: mature.fa. Functional annotation of the predicted microRNA targets were conducted using the GO database. GOseq package was downloaded ( and used in this study [51]. GOseq based on Wallenius non-central hyper-geometric distribution was implemented for GO enrichment analysis.

Prediction of circRNA-miRNA interactions

circRNAs function as miRNA sponges and affect miRNA-mediated regulation of target gene expression via natural sequestration or competitive suppression of miRNA activity [52]. circRNA-miRNA interactions were explored with these obtained differentially expressed miRNAs and circRNAs. By using miRanda software, circRNA-miRNA predictive interaction network was built based on the prediction and analysis of miRNA binding sites and the correlations between circRNA and miRNA that were ranked as previously described [53]. In this study, we set the Smith-Waterman hybridization alignments match score higher than 140 and the minimum free energy of the duplex structure less than − 10 to guarantee the reliability of our results. Other setting parameters were as follows: gap open penalty:-9; gap extend penalty: − 4; scaling parameter: 4. Cytoscape 3.71 software was used to diagram the map of the circRNA-miRNA interaction network [54].

Statistical analysis

Statistical analyses were performed using SPSS 20.0 software. All data were expressed as the mean ± SEM. P value less than 0.05 was considered statistically significant. Student’s t test was used to compare the qPCR results.

Availability of data and materials

The high-throughput data are available in the NCBI database. The bioproject accession numbers for the mRNA-seq, miRNA-seq and circRNA-seq data reported in this paper are PRJNA516994(, PRJNA516517( and PRJNA516718(, respectively. The full dataset is also available from Chun-Xue Zhou upon request (



Acute infection


Acquired immune deficiency syndrome


Chronic infection


DiGeorge syndrome critical region gene 6


Diacylglycerol kinase zeta


DNA methyltransferase 1


Days post infection


Fold change


Fragments per kilo-base of exon per million fragments


Glyceraldehyde-3-phosphate dehydrogenase


Gene Ontology


Interferon-induced protein 44 like


Interferon regulatory factor 1;MSCs

Mesenchymal stem cells




Serine (or cysteine) peptidase inhibitor, clade G, member 1




Transforming growth factor, beta receptor I


Tripartite motif-containing 21


  1. Hide G. Role of vertical transmission of Toxoplasma gondii in prevalence of infection. Expert Rev Anti Infect Ther. 2016;14(3):335–44.

    Article  CAS  PubMed  Google Scholar 

  2. Tenter AM, Heckeroth AR, Weiss LM. Toxoplasma gondii: from animals to humans. Int J Parasitol. 2000;30(12–13):1217–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Pappas G, Roussos N, Falagas ME. Toxoplasmosis snapshots: global status of Toxoplasma gondii seroprevalence and implications for pregnancy and congenital toxoplasmosis. Int J Parasitol. 2009;39(12):1385–94.

    Article  PubMed  Google Scholar 

  4. Sullivan WJ Jr, Jeffers V. Mechanisms of Toxoplasma gondii persistence and latency. FEMS Microbiol Rev. 2012;36(3):717–33.

    Article  CAS  PubMed  Google Scholar 

  5. Ferguson DJ, Hutchison WM. The host-parasite relationship of Toxoplasma gondii in the brains of chronically infected mice. Virchows Arch A Pathol Anat Histopathol. 1987;411(1):39–43.

    Article  CAS  PubMed  Google Scholar 

  6. Sims TA, Hay J, Talbot IC. An electron microscope and immunohistochemical study of the intracellular location of Toxoplasma tissue cysts within the brains of mice with congenital toxoplasmosis. Br J Exp Pathol. 1989;70(3):317.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. Ho YC, Sun HY, Chen MY, Hsieh SM, Chang SC. Clinical presentation and outcome of toxoplasmic encephalitis in patients with human immunodeficiency virus type 1 infection. J Microbiol Immunol Infect. 2008;41(5):386–92.

    PubMed  Google Scholar 

  8. Strittmatter C, Lang W, Wiestler OD, Kleihues P. The changing pattern of human immunodeficiency virus-associated cerebral toxoplasmosis: a study of 46 postmortem cases. Acta Neuropathol. 1992;83(5):475–81.

    Article  CAS  PubMed  Google Scholar 

  9. Prandovszky E, Gaskell E, Martin H, Dubey JP, Webster JP, McConkey GA. The neurotropic parasite Toxoplasma gondii increases dopamine metabolism. PLoS One. 2011;6(9):e23866.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Torrey EF, Bartko JJ, Lun ZR, Yolken RH. Antibodies to Toxoplasma gondii in patients with schizophrenia: a meta-analysis. Schizophr Bull. 2007;33(3):729–36.

    Article  PubMed  Google Scholar 

  11. Mortensen PB, Nørgaard-Pedersen B, Waltoft BL, Sørensen TL, Hougaard D, Torrey EF, et al. Toxoplasma gondii as a risk factor for early-onset schizophrenia: analysis of filter paper blood samples obtained at birth. Biol Psychiatry. 2007;61(5):688–93.

    Article  PubMed  Google Scholar 

  12. Palmer BS. Meta-analysis of three case controlled studies and an ecological study into the link between cryptogenic epilepsy and chronic toxoplasmosis infection. Seizure. 2007;16(8):657–63.

    Article  PubMed  Google Scholar 

  13. Miman O, Kusbeci OY, Aktepe OC, Cetinkaya Z. The probable relation between Toxoplasma gondii and Parkinson's disease. Neurosci Lett. 2010;475(3):129–31.

    Article  CAS  PubMed  Google Scholar 

  14. Çelik T, Kamişli Ö, Babür C, Cevik MO, Oztuna D, Altinayar S. Is there a relationship between Toxoplasma gondii infection and idiopathic Parkinson's disease? Scand J Infect Dis. 2010;42(8):604–8.

    Article  PubMed  Google Scholar 

  15. Christensen M, Schratt GM. microRNA involvement in developmental and functional aspects of the nervous system and in neurological diseases. Neurosci Lett. 2009;466(2):55–62.

    Article  CAS  PubMed  Google Scholar 

  16. Cannella D, Brenier-Pinchart MP, Braun L, van Rooyen JM, Bougdour A, Bastien O, et al. miR-146a and miR-155 delineate a microRNA fingerprint associated with Toxoplasma persistence in the host brain. Cell Rep. 2014;6(5):928–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Xu MJ, Zhou DH, Nisbet AJ, Huang SY, Fan YF, Zhu XQ. Characterization of mouse brain microRNAs after infection with cyst-forming Toxoplasma gondii. Parasit Vectors. 2013;6(1):154.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Conn SJ, Pillman KA, Toubia J, Conn VM, Salmanidis M, Phillips CA, et al. The RNA binding protein quaking regulates formation of circRNAs. Cell. 2015;160(6):1125–34.

    Article  CAS  PubMed  Google Scholar 

  19. Granados-Riveron JT, Aquino-Jarquin G. The complexity of the translation ability of circRNAs. Biochim Biophys Acta. 2016;1859(10):1245–51.

    Article  CAS  PubMed  Google Scholar 

  20. Shao Y, Chen Y. Roles of circular RNAs in neurologic disease. Front Mol Neurosci. 2016;9:25.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Hu RS, He JJ, Elsheikha HM, Zhang FK, Zou Y, Zhao GH, et al. Differential brain microRNA expression profiles after acute and chronic infection of mice with Toxoplasma gondii oocysts. Front Microbiol. 2018;9:2316.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Shen M, Li T, Zhang G, Wu P, Chen F, Lou Q, et al. Dynamic expression and functional analysis of circRNA in granulosa cells during follicular development in chicken. BMC Genomics. 2019;20(1):96.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Kannan G, Moldovan K, Xiao JC, Yolken RH, Jones-Brando L, Pletnikov MV. Toxoplasma gondii strain-dependent effects on mouse behaviour. Folia Parasitol (Praha). 2010;57(2):151–5.

    Article  Google Scholar 

  24. Torrey EF, Yolken RH. Toxoplasma gondii and schizophrenia. Emerg Infect Dis. 2003;9(11):1375.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Carruthers VB, Suzuki Y. Effects of Toxoplasma gondii infection on the brain. Schizophr Bull. 2007;33(3):745–51.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Pittman KJ, Aliota MT, Knoll LJ. Dual transcriptional profiling of mice and Toxoplasma gondii during acute and chronic infection. BMC Genomics. 2014;15:806.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Tanaka S, Nishimura M, Ihara F, Yamagishi J, Suzuki Y, Nishikawa Y. Transcriptome analysis of mouse brain infected with Toxoplasma gondii. Infect Immun. 2013;81(10):3609–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Tili E, Croce CM, Michaille JJ. miR-155: on the crosstalk between inflammation and cancer. Int Rev Immunol. 2009;28(5):264–84.

    Article  CAS  PubMed  Google Scholar 

  29. O'connell RM, Rao DS, Chaudhuri AA, Baltimore D. Physiological and pathological roles for microRNAs in the immune system. Nat Rev Immunol. 2010;10(2):111–22.

    Article  CAS  PubMed  Google Scholar 

  30. O'Connell RM, Taganov KD, Boldin MP, Cheng G, Baltimore D. MicroRNA-155 is induced during the macrophage inflammatory response. Proc Natl Acad Sci U S A. 2007;104(5):1604–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Lu X, Wei Y, Liu F. Direct regulation of p53 by miR-142a-3p mediates the survival of hematopoietic stem and progenitor cells in zebrafish. Cell Discov. 2015;1:15027.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Zhou Z, You Z. Mesenchymal stem cells alleviate LPS-induced acute lung injury in mice by MiR-142a-5p-controlled pulmonary endothelial cell autophagy. Cell Physiol Biochem. 2016;38(1):258–66.

    Article  CAS  PubMed  Google Scholar 

  33. LaMonte G, Philip N, Reardon J, Lacsina JR, Majoros W, Chapman L, et al. Translocation of sickle cell erythrocyte microRNAs into Plasmodium falciparum inhibits parasite translation and contributes to malaria resistance. Cell Host Microbe. 2012;12(2):187–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Viticchiè G, Lena AM, Latina A, Formosa A, Gregersen LH, Lund AH, et al. MiR-203 controls proliferation, migration and invasive potential of prostate cancer cell lines. Cell Cycle. 2011;10(7):1121–31.

    Article  PubMed  CAS  Google Scholar 

  35. Zhu X, Er K, Mao C, Yan Q, Xu H, Zhang Y, et al. miR-203 suppresses tumor growth and angiogenesis by targeting VEGFA in cervical cancer. Cell Physiol Biochem. 2013;32(1):64–73.

    Article  CAS  PubMed  Google Scholar 

  36. Diao Y, Guo X, Jiang L, Wang G, Zhang C, Wan J, et al. miR-203, a tumor suppressor frequently down-regulated by promoter hypermethylation in rhabdomyosarcoma. J Biol Chem. 2014;289(1):529–39.

    Article  CAS  PubMed  Google Scholar 

  37. Qadir XV, Han C, Lu D, Zhang J, Wu T. miR-185 inhibits hepatocellular carcinoma growth by targeting the DNMT1/PTEN/Akt pathway. Am J Pathol. 2014;184(8):2355–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Tang H, Liu P, Yang L, Xie X, Ye F, Wu M, et al. miR-185 suppresses tumor proliferation by directly targeting E2F6 and DNMT1 and indirectly upregulating BRCA1 in triple-negative breast cancer. Mol Cancer Ther. 2014;13(12):3185–97.

    Article  CAS  PubMed  Google Scholar 

  39. Wilusz JE, Sharp PA. A circuitous route to noncoding RNA. Science. 2013;340(6131):440–1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Wang K, Long B, Liu F, Wang JX, Liu CY, Zhao B, et al. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223. Eur Heart J. 2016;37(33):2602–11.

    CAS  PubMed  Google Scholar 

  41. Zhou CX, Elsheikha HM, Zhou DH, Liu Q, Zhu XQ, Suo X. Dual identification and analysis of differentially expressed transcripts of porcine PK-15 cells and Toxoplasma gondii during in vitro infection. Front Microbiol. 2016;7:721.

    PubMed  PubMed Central  Google Scholar 

  42. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Wen M, Shen Y, Shi S, Tang T. miREvo: an integrative microRNA evolutionary analysis platform for next-generation sequencing experiments. BMC Bioinform. 2012;13(1):140.

    Article  CAS  Google Scholar 

  44. Friedländer MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2011;40(1):37–52.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–8.

    Article  CAS  PubMed  Google Scholar 

  46. Gao Y, Wang J, Zhao F. CIRI: an efficient and unbiased algorithm for de novo circular RNA identification. Genome Biol. 2015;16(1):4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Zhou L, Chen J, Li Z, Li X, Hu X, Huang Y, et al. Integrated profiling of microRNAs and mRNAs: microRNAs located on Xq27. 3 associate with clear cell renal cell carcinoma. PLoS One. 2010;5:e15224.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodology. 1995;57(1):289–300.

    Google Scholar 

  50. Xia JG, Sinelnikov IV, Han B, Wishart DS. MetaboAnalyst 3.0-making metabolomics more meaningful. Nucleic Acids Res. 2015;43(W1):W251–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Sand M, Bechara FG, Gambichler T, Sand D, Bromba M, Hahn SA, et al. Circular RNA expression in cutaneous squamous cell carcinoma. J Dermatol Sci. 2016;83(3):210–8.

    Article  CAS  PubMed  Google Scholar 

  53. Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in drosophila. Genome Biol. 2003;5(1):R1.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We are grateful to Beijing Novogene biotechnology co. ltd. for technical assistance.


Project support was provided by the Natural Science Foundation of Shandong Province of China (Grant No. ZR2018BH30), Basal Research Fund of Shandong University (Grant No. 2017GN0015), China Postdoctoral Science Foundation (Grant No. 2018 M642665), the Fund of Fujian Provincial Key Laboratory for the Prevention and Control of Animal infectious Diseases and Biotechnology (Grant No. 2018KF04) and the Agricultural Science and Technology Innovation Program (ASTIP) (Grant No. CAAS-ASTIP-2016-LVRI-03).The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



XQZ and SYH conceived and designed the experiments. CXZ and KA performed the experiments. CQH, JJG and HC contributed reagents/materials/analysis tools. CXZ analyzed the data and wrote the paper. SYH and XQZ critically revised the manuscript. All authors read and approved the final version of the manuscript.

Corresponding authors

Correspondence to Chun-Xue Zhou, Shen-Yi He or Xing-Quan Zhu.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Animal Ethics Committee of Lanzhou Veterinary Research Institute, Chinese Academy of Agriculture Sciences (Permit No. LVRIAEC2016–007). All experiments were performed strictly according to the requirements of the Animal Ethics Procedures and Guidelines of the People’s Republic of China. We adhered to the ARRIVE guidelines wherever possible.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Reads quality of miRNA libraries.

Additional file 2: Table S2.

Statistics of small RNA reads mapped to the reference sequence.

Additional file 3: Table S3.

RNA reads counts in the infected and control brains.

Additional file 4: Figure S1.

Principal Component Analysis (PCA) score scatter plots of miRNA (A) and circRNA (B) obtained from high-throughput RNA sequencing. Ellipses enclose the 95% confidence intervals estimated by the sample means and covariances of each group.

Additional file 5: Figure S2.

Validation of the expression of mmu-miR-155-5p using qPCR analysis. Y axis indicates the relative expression with log2 fold change.

Additional file 6: Table S4.

Reads quality of circRNA libraries.

Additional file 7: Table S5.

Annotation for all circRNAs identified in this study.

Additional file 8: Figure S3.

Validation of circRNA expression by qPCR. Y axis indicates the relative expression with log2 fold change.

Additional file 9: Table S6.

Differentially expressed circRNAs and their corresponding host genes.

Additional file 10: Table S7.

The data used to construct Fig. 6.

Additional file 11: Figure S4.

qPCR verification.

Additional file 12: Figure S5.

Validation of mRNA expression by qPCR. Y axis indicates the relative expression with log2 fold change.

Additional file 13: Table S8.

Primers for miRNA and circRNA q-PCR analysis.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhou, CX., Ai, K., Huang, CQ. et al. miRNA and circRNA expression patterns in mouse brain during toxoplasmosis development. BMC Genomics 21, 46 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Toxoplasma gondii
  • Brain
  • miRNA
  • circRNAs
  • Interaction