Skip to main content

A comparative study of small RNAs in Toxoplasma gondii of distinct genotypes



Toxoplasma gondii is an intracellular parasite with a significant impact on human health. Inside the mammalian and avian hosts, the parasite can undergo rapid development or remain inactive in the cysts. The mechanism that regulates parasite proliferation has not been fully understood. Small noncoding RNAs (sncRNA) such as microRNAs (miRNAs) are endogenous regulatory factors that can modulate cell differentiation and development. It is anticipated that hundreds of miRNAs regulate the expression of thousands of genes in a single organism. SncRNAs have been identified in T. gondii, however the profiles of sncRNAs expression and their potential regulatory function in parasites of distinct genotypes has largely been unknown.


The transcription profiles of miRNAs in the two genetically distinct strains, RH and ME49, of T. gondii were investigated and compared by a high-through-put RNA sequencing technique and systematic bioinformatics analysis. The expression of some of the miRNAs was confirmed by Northern blot analysis.


1,083,320 unique sequences were obtained. Of which, 17 conserved miRNAs related to 2 metazoan miRNA families and 339 novel miRNAs were identified. A total of 175 miRNAs showed strain-specific expression, of which 155 miRNAs were up-regulated in RH strain and 20 miRNAs were up-regulated in ME49 strain. Strain-specific expression of miRNAs in T. gondii could be due to activation of specific genes at different genomic loci or due to arm-switching of the same pre-miRNA duplex.


Evidence for the differential expression of miRNAs in the two genetically distinct strains of T. gondii has been identified and defined. MiRNAs of T. gondii are more species-specific as compared to other organisms, which can be developed as diagnostic biomarkers for toxoplasmosis. The data also provide a framework for future studies on RNAi-dependent regulatory mechanisms in the zoonotic parasite.


Toxoplasma gondii is an obligatory intracellular parasite that infects a wide range of hosts,including humans, animals and birds. It is considered to be one of the most widely distributed protozoan parasite with a sero-prevalence in humans of up to 30% worldwide [1]. T. gondii is the etiological agent of toxoplasmosis which can be either life-threatening or long-term chronic infection.

The life cycle of Toxoplasma gondii is unusual in that the parasite is capable of indefinite proliferation in the hosts with either a sexual or an asexual cycle. The sexual cycle occurs only in the hosts of a feline species. The asexual cycle can occur in virtually any warm-blooded animals, which act as the intermediate hosts, ranging from chicken to baleen whales and humans. In the parasite’s life cycle, there are three fundamental morphotypes, named tachyzoites, bradyzoites and sporozoites. The development of T. gondii in the intermediate host involves an initial phase with a rapid proliferation of the tachyzoites, followed by the formation of tissue cysts containing slowly dividing or resting bradyzoites.

While the global population structure of T. gondii awaits further elucidation [2], three clonal lineages (Type I, II, and III) of T. gondii which comprise the majority of strains in both North America and Europe [35]. Recently, a fourth clonal lineage, designated haplogroup 12, has been identified based on isolates that are common in wild animals in the United States [6]. The virulence of T. gondii is normally defined based on the LD50 in mice. Type I strain has been regarded as a more virulent strain in mice with an LD50 <10 [4, 7]. In contrast, type II (e.g., ME49) and type III strains are less lethal in mice (LD50 > 100) and the infections are usually less severe or asymptomatic [8]. The clonal lineages also differ in a number of phenotypes such as growth rate, efficiency in migration and transmigration in tissues [7, 9].

Several studies have been carried out with aims to determine differential gene expression at various lifecycle stages of T. gondii[1013]. Likewise, efforts were taken to reveal the factors that may modulate the virulence of T. gondii. ROP5, ROP16 and ROP18, which are rhoptry-derived factors, and their expression was found to correlate with parasite virulence [1418]. Identification of strain-specific regulatory elements responsible for the distinct genotype will not only facilitate our understanding of parasite biology, but may also elucidate the association between genetically defined subpopulations (strains) and disease severity.

In recent years, the discovery of numerous small RNAs has increased our knowledge in post-transcriptional gene regulation in cell development and other biological processes. Small RNAs, such as microRNA (miRNA), small interfering RNA (siRNA), and Piwi-associated RNA (piRNA), are regulatory elements that can modulate gene expression at the post-transcriptional level. MiRNAs, a class of ≈ 22 nucleotide small-RNA sequences that participate in the post-transcriptional regulation of gene expression [19], have been known to play critical roles in diverse biological processes, including development, viral defense, metabolism, and apoptosis [2023]. The primary transcripts (pri-miRNAs) of miRNA genes are generated by either RNA polymerase II [24] or RNA polymerase III [25]. A single pri-miRNA may contain from one to six miRNA precursors. They are processed by RNase complexes (Drosha and DGCR8) into ≈ 70 nucleotide fragments with a stem-loop structure called pre-miRNA. The pre-miRNAs are then exported to the cytoplasm by Exportin-5 [26] and the hairpin is cleaved by the RNase III enzyme (Dicer) to miRNA duplexes [27]. A single miRNA can silence a number of genes while one gene can be targeted by several miRNAs [28, 29]. Previous studies have found that T. gondii has a complete machinery for small RNA generation and small RNA-mediated gene regulation [30]. There are mainly two classes of small regulatory RNAs derived from the T. gondii genome, namely miRNA and rdsRNA. And it was found that rdsRNAs were consistently more abundant in the highly virulent Toxoplasma isotype-I (RH) than in other two strains (RPU and CTG). Further, it was recently reported that the Argonaute of T. gondii (TgAgo) is methylated and its activity is Mg2+-dependent [31]. However the differential expression of miRNA in parasites of distinct genotypes has not been well investigated.

In this study, we systematically studied the expression of miRNAs in type I (RH) and type II strain (ME49) T. gondii by high-through-put RNA sequencing technology and bioinformatic analysis. A number of T. gondii specific miRNAs was identified. Meanwhile, differentially expressed miRNAs between the two strains of T. gondii were detected. The results demonstrated profound differences between the two strains of T. gondii in miRNA expression.



Tachyzoites of T. gondii RH strain and ME49 strain were propagated in mice. The tachyzoites were purified by density-gradient centrifugation on Percoll [32]. The study of using laboratory animals was reviewed and approved by the Ethical Committee of Jilin University.

RNA isolation

Total RNA of T. gondii (both RH and ME49 strain) was prepared using Trizol reagent (Invitrogen, SF, USA) according to the manufacturer’s protocol. The integrity of total RNA was examined by standard agarose gel electrophoresis, and RNA purity was reflected by the 260/280 nm absorbance ratio and the concentration was determined using a Nanodrop 1000 machine (Thermo Scientific CA, USA). The purified total RNA was stored at-80°C until use.

Construction of small RNA libraries and sequencing

For small RNA library construction and deep sequencing, RNA samples from RH and ME49 strains of T. gondii were prepared as follows: for each strain, equal quantities (10 μg) of RNA isolated from tachyzoites were pooled. Small RNA molecules in the range of 18–30 nt RNA was purified after polyacrylamide gel electrophoresis (PAGE) and ligated with proprietary adapters to the 5 and 3 termini. The samples were used as templates for cDNA synthesis. The cDNA was amplified to produce sequencing libraries which were subjected to Solexa’s sequencing-by-synthesis method. Sequencing was carried out at the Beijing Genomics Institute (BGI, China). Two separated runs with two batches of RNAs were carried out in the sequencing. Only sequences with high quality were included in the analysis.

Mapping sequence reads to the reference genome

Individual sequence reads with the base quality scores were produced by Solexa. Clean reads were obtained after removing of the low quality reads, such as adaptor null reads, insert null reads, 5 adaptor contaminants, and reads with polyA tail and ambiguous bases. Adapter sequences were trimmed from both ends of clean reads before analysis. All identical sequences were counted and eliminated from the initial data set. The resulting set of unique sequences with associated read counts was referred to as sequence tags. The unique reads were mapped onto the T. gondii genome ( using the program SOAP ( [33].

Bioinformatic analysis of T. gondii small RNAs

Sequences were, at the first step, searched against Metazoa mature miRNA of Sanger miRBase allowing two mismatches to identify homologs of known Metazoa miRNAs using the program Patscan [34]. Sequence tags of more than 5 reads that matched perfectly or near-perfectly (no more than 2 mismatches and mismatch not positioned in seed region) to metazoan mature miRNAs were assumed to be conserved miRNA candidates. Secondly, the remaining sequences were screened against the non-coding RNA database [35] and the T. gondii genome database, and some non-coding RNAs such as rRNA, tRNA and snoRNAs were identified and filtered out for further analysis. Thirdly, We searched for the inverted repeats (step loops or hairpin structure) among the remained sequences which did not match to the miRNA database by using the software Einverted of Emboss [36] with the parameters following parameters: threshold = 30, match score = 3, mismatch score = 3, gap penalty = 6, and maximum repeat length = 240 as described [37]. Each inverted repeat was extended 10 nt on each side, and the secondary structure of the inverted repeat was predicted by RNAfold [38]. Unique reads in the inverted repeats were evaluated by MirCheck [37] using modified parameters that are more suitable for organisms such as worms and protozoa. Finally, miRNA precursors that passed MirCheck were inspected manually in order to remove the false prediction (Additional file 1: Figure S1). Further, miRNAs with the same sequences derived from pre-miRNAs (the sequences may not be completely the same) located at different genomic loci are indicated with an additional dash-number suffix. Sequences of T. gondii-specific miRNAs have been deposited in the miRNA database (miRBase).

MiRNAs with statistically significant differences in relative abundance (as reflected by TPM values) between the two libraries (corresponding to the two distinct strains of the parasites) were analyzed with IDEG6 and edger [39]. We used the general Chi-square method as it is most frequently used in such analysis [40, 41]. MiRNA with a P value ≤ 0.01 were deemed to be significantly different between the samples of the two distinct strains of parasites.

Northern blot analysis of miRNA expression

Twenty micrograms of total RNA was loaded in each well of a 12.5% polyacrylamide gel containing 42% urea and run in 0.5× TBE buffer. Following the electrophoresis, RNA was transferred by capillary transferring to a Hybond-N + nylon membranes (GE Healthsystems, Uppsala, Sweden). After UV cross-linking, the membranes were baked for 1 h at 80°C. Probes complementary to small RNA sequences were end-labeled with DIG at 5 Termini (TaKaRa, Dalian, China). Pre-hybridization of the membrane was performed by overnight incubation at 53°C followed by hybridization overnight at the same temperature in Northernmax Hybridization buffer (Ambion, CA, USA). After hybridization, the membranes were washed four times for 30 min in 2×SSC, 0.05% SDS and twice for 15 min in 0.1×SSC, 0.1% SDS at RT. Chemiluminescent signal was detected using a DIG Detection Kit (Roche, Germany) following the manufacturer’s instructions. The oligonucleotide probes used for hybridization are as follows:

Tgo-novel-40 probe: TTCCTG*ATCCTA*TTTAG*CAG*GT

Tgo-novel-12-1 probe: CCA*CTTCAG*TCTTCA*AAG*TTCT


Tgo-novel-1-1-3p probe: GACAG*TGCTCG*GGTCCG*CAACA*CGCC

Tgo-novel-14-2 probe: ACA*CTT*CCCTT*CTCG*CCG

(LNA (locked nucleic acid) substitutions are indicated by a “*”).


An overview of the small RNA sequencing results

A total of 8,738,870 and 10,759,107 sequence reads were obtained from the ME49 and RH libraries, respectively. There is no major difference in total numbers of low and high quality reads between the two libraries (Additional file 2: Table S1). After removal of the low quality reads, 7,149,051 (ME49) and 8,494,754 (RH) clean reads were obtained, which contained 247,346 (ME49) and 867,853 (RH) library-specific (unique) reads, respectively (Figure 1).

Figure 1
figure 1

Work flow for profiling of small RNAs after high-through-put sequencing. Approaches and numbers of sequence tags obtained after each step were shown for both strains. The simple pie chart shows the number of strain-specific and common miRNAs identified in the two libraries. More miRNAs (155) were identified in RH strain than that in ME49. The number of miRNAs that commonly expressed in the two strains is 181.

The clean unique reads described above were mapped to the draft genome of T. gondii ( [33]. 6,068,932 and 3,088,774 non-redundant total reads from ME49 and RH were perfectly mapped onto the T. gondii genome, respectively (Additional file 3: Table S2). Among these reads, 51,832 (76.65% of total reads in ME49 library) displayed strain-specificity in ME49, while in RH the number of strain-specific reads was 117,315 (88.14% of total reads in RH library), the number of small RNA reads identified in the two strains was 15,791, which accounted for 11.86% and 23.35% of the total reads of RH and ME49, respectively (Figure 2). The length of small RNAs varied from 18 to 30 nt in the two strains. However, the length distributions of small RNAs were significantly different (Figure 3, Additional 4: Figure S2 and Additional file 5: Table S3). A majority of small RNAs in ME49 was 21 nt in length, whereas small RNAs of 26 nt were the most abundant in RH strain. The small RNAs were further categorized into, based on the sequence characteristics, rRNAs, tRNAs, small nuclear RNAs (snRNAs) and other ncRNA after BLASTN searches against the Sanger Rfam database release 9.0 (Additional file 3: Table S2). The proportions of small ncRNAs in the two libraries (ME49 and RH) were shown in Figure 2 and Figure 4A, B.

Figure 2
figure 2

Comparison of small RNAs sequences identified in the two strains of T. gondii at unique level. Of the small RNAs identified in the two strains of T. gondii, strain-specific small RNAs were much more (88.14% of RH strain and 76.65% of ME49 strain) than that shared by the two strains (11.86% of RH strain and 23.35% of ME49 strain).

Figure 3
figure 3

The length distribution of small ncRNAs in the library of ME49 (blue) and RH (red) of T. gondii . The dominant sncRNA expressed in ME49 strain is 21 nt in length, while the 26 nt sncRNAs were dominant in RH strain.

Figure 4
figure 4

Percentages of small ncRNAs in the two libraries (A. ME 49 strain, B. RH strain) of T. gondii. A majority of the sncRNAs in the two libraries is unknown. The proportion of miRNAs in ME49 and RH strain is 8.08% and 5.30% respectively.

Identification of miRNAs in T. gondii

The proportions of miRNAs in the two libraries accounted between 5-8% of the total reads of the small ncRNAs (Figure 4 A, B). In total, 17 conserved miRNAs were found based on the consensus ‘seed’ region (2–8 nt in 5 end of mature miRNA) identical to that of Homo sapiens, Mus musculus and Pongo pygmaeus (Figure 5), which were recognized by mRNAs through base pairing comparison [42]. In addition to the conserved miRNAs, we also found 339 species-specific miRNAs (novel miRNAs) of T. gondii. Moreover, 20 (5.6%) of 356 miRNAs were from the ME49 library, while 155 (43.6%) were from RH library, and 181 (50.8%) of the miRNAs were found in both libraries (Figure 1, Additional file 6: Table S4 and Additional 7: Table S5). The number of strain-specific miRNAs in RH strain was more than that in ME49 strain.

Figure 5
figure 5

Alignment of tgo-miR-574 sequence with homologues from other organisms. The seed sequences are shadowed in dark colour. Mmu-miR-574-5p, mmu-miR-1187 and mmu-Mir-466i-5p are miRNAs identified in mice. Has-miR-574-5p is a miRNA identified in humans.

Characterization of the expression of T. gondii miRNAs

Based on their genomic locations, T. gondii miRNAs were categorized into three types (Additional file 6: Table S4 and Additional file 7: Table S5) named intronic, intergenic and UTR-derived miRNAs. The numbers of miRNAs derived from the 3 genomic regions were 33, 305 and 18 respectively. No miRNA genes were found to be located in exons. Thus miRNA genes were predominantly intergenic in T. gondii. This observation was in agreement with a previous study on Schistosoma japonicum which suggested that most miRNA genes have their own control elements (or promoters) in the genome [43]. Further, it was found that a miRNA can be generated from several pre-miRNAs encoded by genes located in different genomic loci. Thus miRNAs with the same sequence but derived from different pre-miRNAs were named with an additional dash-number suffix (Additional file 6: Table S4 and Additional file 7: Table S5).

Comparison of the novel miRNAs expressed in ME49 and RH revealed strain-specific expression patterns in the two distinct strains of T. gondii. For instance, the number of reads of tgo-novel-40 in ME49 was 4000 times more than that in RH. Similarly, the expression level of tgo-novel-1-1-5p, tgo-novel-12-1, tgo-novel-41 and tgo-novel-15-1 was significantly higher in ME49 strain. On the contrary, tgo-novel-14-2, was highly expressed in RH, the reads number was about 10 times as much as that in ME49 (Additional file 7: Table 5 and Additional file 8: Table S6).

Previous studies have found that mature miRNAs can be derived from both arms of a pre-miRNA hairpin [44, 45]. In this study, we found 178 miRNAs were derived from just one hairpin arm, of which 81 mature miRNAs were located at 3 arm and 97 were located at 5 arm (Figure 6 and Additional file 7: Table S5) of the predicted hairpins. However, one miRNA, tgo-novel-1, the miRNA:miRNA* ratios showed strain-specific pattern. The reads from 3 arm were dominant in ME49 (5/3 read ratio: 443/6212), whereas in RH, the reads number from 5 arm is higher (5/3 read ratio: 9632/2838) (Additional file 7: Table S5).

Figure 6
figure 6

The sequence and the secondary stem-loop structure of tgo-novel-12-14 identified in T. gondii. Sequences and the number of reads of the mature miRNA and the complementary miR* are represented in red and blue respectively. The predicted structure of the pre-miRNA is represented on the right side.

Validation of miRNAs expression by Northern blot

Five novel miRNAs (tgo-novel-1-1-3p, tgo-novel-12-1, tgo-novel-40, tgo-novel-41, and tgo-novel-14-2) with relatively high abundance identified by sequencing were verified by Northern blot. Specific hybridization with probes of three miRNAs was observed at ≈ 23 nt (Figure 7). In addition, all probes showed hybridization signals to the pre-miRNA transcripts of about 80 nt, except for the precursor of tgo-novel-40, whose signal was detected at ≈ 100 nt. Further, the hybridization signal to the mature miRNAs of tgo-novel-40, tgo-novel-41 and their pre-miRNAs was more intense in ME49 than that in RH, which implying higher expression level in ME49. The expression of tgo-novel-1-1-3p was only detected in RH, which was similar with the results of sequencing analysis. On the contrary, the mature miRNA of tgo-novel-14-2 was only detected in ME49, though the pre-miRNAs with a similar expression level detected in the two strains. The inconsistency with the sequencing data might be due to the slow processing of pre-miRNAs in RH strain. No hybridization was seen with any probe to the mouse miRNAs.

Figure 7
figure 7

Characterization of four miRNAs by Northern-blot. Lanes from left to right are RNAs from ME49 strain (ME49), RH strain (RH) and a mouse. Total RNA isolation from ME49, RH strain and a mouse were visualized by ethidium bromide staining and served as loading controls at the top panel. LNA probes corresponding to tgo-Novel-1-1-3p, tgo-Novel-12-1, tgo-Novel-40, tgo-Novel-41, tgo-Novel-14-2 were used. Probe tgo-Novel-1-1-3p and tgo-Novel-40 only hybridized to the RNAs of RH strain and ME49, respectively. While probes of tgo-Novel-12-1, tgo-Novel-41 and tgo-Novel-14-2 hybridized to RNA of both strains. More hybridization was seen with tgo-Novel-41 probe to ME49. More hybridization was seen with pre-miRNAs except with the probe of tgo-Novel-41. No hybridization was detected with mouse RNA with any probe. Hybridization to mature miRNAs was marked with asterisks.


MiRNAs are recognized as critical regulators in gene expression at the post-transcriptional level. Previous studies have found that T. gondii possesses a complete RNA silencing pathway which suggests that small noncoding RNAs may play a critical role in the parasite development and its parasitization in the hosts [30]. In this study, the profiles of small RNA populations of the two distinct strains of T. gondii were investigated.

The distribution of small RNAs in the two small RNA libraries generated after deep sequencing was compared. A predominant number of small ncRNAs was strain-specific (Figure 1 and Figure 2) and the strain-specific small ncRNA including miRNAs were more in RH strain than in ME49 strain (Figure 1). The difference in miRNA numbers found in the two libraries could be due to the less presentation of ME49 genomic sequences in the databases which may affect small RNA identification. Further, there was a clear difference in the tendency of the length distribution in small RNAs between the two strains. In RH, the 26 nt RNAs represented the predominant species (Figure 3, Additional file 5: Table S3), while in ME49, about 18.92% of the total small RNAs were 21 nt in size, which was the most abundant class. This might be due to the genetic differences between the two parasite strains, which also suggested that the two parasite strains are biologically different. Previous study reported that the structural features of pre-miRNA hairpins might influence the efficiency of Dicer binding and specificity of precursor cleavage, which leads to the length diversity of miRNAs [46, 47]. The reasons of strain-specific length distribution of the small ncRNAs between the two strains of T. gondii remain vague, further studies are needed to dissect the mechanism in sncRNA processing that may be associated with strain-specific gene regulation.

In total, we identified 17 conserved miRNA and 339 species-specific miRNAs in the two strains of T. gondii, of which 7 miRNAs sequences were reported by Braun et al. in 2010 [30]. Interestingly, about 5% of T. gondii miRNAs were categorized as conserved and accounted for less than 1% of the read counts, whereas more than 99% of the remaining miRNAs were recognized as species-specific. This observation supports the earlier finding that T. gondii possesses a RNA-associated gene regulation machinery which is phylogenetically diverged from mammals but more similar to plants [30]. Further, let-7 and lin-4, the two most conserved miRNAs in metazoan, were not found in T. gondii indicting that the fine-tuning mechanism of miRNAs in T. gondii was distinct from other species.

During the process of the biosynthesis of miRNAs, miRNA and miRNA* (or miR*) were generated by enzymatic cleavage of the 70–80 nt precursor hairpin. The functional strand (miRNA) of the small RNA duplex is preferentially loaded into the RISC as the guide strand, while the other strand, the passenger strand (miRNA*), is degraded [44, 45]. However, recent studies suggested that mature miRNAs can be generated from both strands of the pre-miRNA duplex [4850]. In C. elegans and related nematodes, it has been reported that the diversity of miRNAs was, at least partially, due to the arm-switching and hairpin shifting [51]. We found that tgo-novel-1 changed the miRNA strand with arm switching of the same hairpin between the two strains of T. gondii. In RH strain, the dominant miRNA of tgo-novel-1 was derived from the 3 arm, while in ME49 it seemed that mature miRNA was only derived from the 5 arm. Due to the fact that the sequences of miRNAs derived from the two arms of the same hairpin were complementary, they likely regulate different target sequences. Further dissection of the function of the miRNAs derived from the two arms of the same hairpin might lead to deep understanding of the parasite biology.


In summary, 17 conserved miRNAs related to 2 metazoan miRNA families and 339 novel miRNAs were identified in the two genetically different strains of T. gondii. The majority of miRNAs were species-specific, which supports the finding that T. gondii is an evolutionarily diverged organism from other protozoana. The difference in expression abundance of certain miRNAs as well as the arm-switching in pre-miRNA processing leading to different miRNA species in the two parasite strains suggested that there was a fine-tuning mechanism of miRNA biogenesis in distinct strains of T. gondii. Understanding the genetic factors that regulate T. gondii gene expression could contribute to the development of specific tools to control the transmission of the parasite.


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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Darde ML: Toxoplasma gondii, “new” genotypes and virulence. Parasite. 2008, 15 (3): 366-371.

    Article  CAS  PubMed  Google Scholar 

  3. Lehmann T, Marcet PL, Graham DH, Dahl ER, Dubey JP: Globalization and the population structure of Toxoplasma gondii. Proc Natl Acad Sci. 2006, 103 (30): 11423-11428. 10.1073/pnas.0601438103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Sibley LD, Boothroyd JC: Virulent strains of Toxoplasma gondii comprise a single clonal lineage. Nature. 1992, 359 (6390): 82-85. 10.1038/359082a0.

    Article  CAS  PubMed  Google Scholar 

  5. Velmurugan GV, Dubey JP, Su C: Genotyping studies of Toxoplasma gondii isolates from Africa revealed that the archetypal clonal lineages predominate as in North America and Europe. Vet Parasitol. 2008, 155 (3–4): 314-318.

    Article  CAS  PubMed  Google Scholar 

  6. Khan A, Dubey JP, Su C, Ajioka JW, Rosenthal BM, Sibley LD: Genetic analyses of atypical Toxoplasma gondii strains reveal a fourth clonal lineage in North America. Int J Parasitol. 2011, 41 (6): 645-655. 10.1016/j.ijpara.2011.01.005.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Barragan A, Sibley LD: Migration of Toxoplasma gondii across biological barriers. Trends Microbiol. 2003, 11 (9): 426-430. 10.1016/S0966-842X(03)00205-1.

    Article  CAS  PubMed  Google Scholar 

  8. Boothroyd JC, Grigg ME: Population biology of Toxoplasma gondii and its relevance to human infection: do different strains cause different disease?. Curr Opin Microbiol. 2002, 5 (4): 438-442. 10.1016/S1369-5274(02)00349-1.

    Article  PubMed  Google Scholar 

  9. Kaufman HE, Remington JS, Jacobs L: Toxoplasmosis: the nature of virulence. Am J Ophthalmol. 1958, 46 (5 Part 2): 255-260. discussion 260–251

    Article  CAS  PubMed  Google Scholar 

  10. Behnke MS, Radke JB, Smith AT, Sullivan WJ, White MW: The transcription of bradyzoite genes in Toxoplasma gondii is controlled by autonomous promoter elements. Mol Microbiol. 2008, 68 (6): 1502-1518. 10.1111/j.1365-2958.2008.06249.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Cleary MD, Singh U, Blader IJ, Brewer JL, Boothroyd JC: Toxoplasma gondii asexual development: identification of developmentally regulated genes and distinct patterns of gene expression. Eukaryot Cell. 2002, 1 (3): 329-340. 10.1128/EC.1.3.329-340.2002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Dzierszinski F, Nishi M, Ouko L, Roos DS: Dynamics of Toxoplasma gondii differentiation. Eukaryot Cell. 2004, 3 (4): 992-1003. 10.1128/EC.3.4.992-1003.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Singh U, Brewer JL, Boothroyd JC: Genetic analysis of tachyzoite to bradyzoite differentiation mutants in Toxoplasma gondii reveals a hierarchy of gene induction. Mol Microbiol. 2002, 44 (3): 721-733. 10.1046/j.1365-2958.2002.02903.x.

    Article  CAS  PubMed  Google Scholar 

  14. Reese ML, Boothroyd JC: A conserved non-canonical motif in the pseudoactive site of the ROP5 pseudokinase domain mediates its effect on Toxoplasma virulence. J Biol Chem. 2011, 286 (33): 29366-29375. 10.1074/jbc.M111.253435.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Behnke MS, Khan A, Wootton JC, Dubey JP, Tang K, Sibley LD: Virulence differences in Toxoplasma mediated by amplification of a family of polymorphic pseudokinases. Proc Natl Acad Sci. 2011, 108 (23): 9631-9636. 10.1073/pnas.1015338108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Reese ML, Zeiner GM, Saeij JP, Boothroyd JC, Boyle JP: Polymorphic family of injected pseudokinases is paramount in Toxoplasma virulence. Proc Natl Acad Sci. 2011, 108 (23): 9625-9630. 10.1073/pnas.1015980108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. El Hajj H, Lebrun M, Arold ST, Vial H, Labesse G, Dubremetz JF: ROP18 is a rhoptry kinase controlling the intracellular proliferation of Toxoplasma gondii. PLoS Pathog. 2007, 3 (2): e14-10.1371/journal.ppat.0030014.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Taylor S, Barragan A, Su C, Fux B, Fentress SJ, Tang K, Beatty WL, Hajj HE, Jerome M, Behnke MS: A secreted serine-threonine kinase determines virulence in the eukaryotic pathogen Toxoplasma gondii. Science. 2006, 314 (5806): 1776-1780. 10.1126/science.1133643.

    Article  CAS  PubMed  Google Scholar 

  19. Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.

    Article  CAS  PubMed  Google Scholar 

  20. Xu P, Vernooy SY, Guo M, Hay BA: The drosophila microRNA Mir-14 suppresses cell death and is required for normal fat metabolism. Curr Biol. 2003, 13 (9): 790-795. 10.1016/S0960-9822(03)00250-1.

    Article  CAS  PubMed  Google Scholar 

  21. Ambros V: The functions of animal microRNAs. Nature. 2004, 431 (7006): 350-355. 10.1038/nature02871.

    Article  CAS  PubMed  Google Scholar 

  22. Miska EA: How microRNAs control cell division, differentiation and death. Curr Opin Genet Dev. 2005, 15 (5): 563-568. 10.1016/j.gde.2005.08.005.

    Article  CAS  PubMed  Google Scholar 

  23. Ding SW, Voinnet O: Antiviral immunity directed by small RNAs. Cell. 2007, 130 (3): 413-426. 10.1016/j.cell.2007.07.039.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Lee Y, Kim M, Han J, Yeom KH, Lee S, Baek SH, Kim VN: MicroRNA genes are transcribed by RNA polymerase II. EMBO J. 2004, 23 (20): 4051-4060. 10.1038/sj.emboj.7600385.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Borchert GM, Lanier W, Davidson BL: RNA polymerase III transcribes human microRNAs. Nat Struct Mol Biol. 2006, 13 (12): 1097-1101. 10.1038/nsmb1167.

    Article  CAS  PubMed  Google Scholar 

  26. Lund E, Guttinger S, Calado A, Dahlberg JE, Kutay U: Nuclear export of microRNA precursors. Science. 2004, 303 (5654): 95-98. 10.1126/science.1090599.

    Article  CAS  PubMed  Google Scholar 

  27. Lund E, Dahlberg JE: Substrate selectivity of exportin 5 and dicer in the biogenesis of microRNAs. Cold Spring Harb Symp Quant Biol. 2006, 71: 59-66. 10.1101/sqb.2006.71.050.

    Article  CAS  PubMed  Google Scholar 

  28. Kim VN, Nam JW: Genomics of microRNA. Trends Genet. 2006, 22 (3): 165-173. 10.1016/j.tig.2006.01.003.

    Article  CAS  PubMed  Google Scholar 

  29. Bartel DP: MicroRNAs: target recognition and regulatory functions. Cell. 2009, 136 (2): 215-233. 10.1016/j.cell.2009.01.002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Braun L, Cannella D, Ortet P, Barakat M, Sautel CF, Kieffer S, Garin J, Bastien O, Voinnet O, Hakimi MA: A complex small RNA repertoire is generated by a plant/fungal-like machinery and effected by a metazoan-like argonaute in the single-cell human parasite Toxoplasma gondii. PLoS Pathog. 2010, 6 (5): e1000920-10.1371/journal.ppat.1000920.

    Article  PubMed Central  PubMed  Google Scholar 

  31. Musiyenko A, Majumdar T, Andrews J, Adams B, Barik S: PRMT1 methylates the single argonaute of Toxoplasma gondii and is important for the recruitment of tudor nuclease for target RNA cleavage by antisense guide RNA. Cell Microbiol. 2012, 14 (6): 882-901. 10.1111/j.1462-5822.2012.01763.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Cornelissen AW, Overdulve JP, Hoenderboom JM: Separation of Isospora (Toxoplasma) gondii cysts and cystozoites from mouse brain tissue by continuous density-gradient centrifugation. Parasitology. 1981, 83 (Pt 1): 103-108.

    Article  CAS  PubMed  Google Scholar 

  33. Li R, Li Y, Kristiansen K, Wang J: SOAP: short oligonucleotide alignment program. Bioinformatics. 2008, 24 (5): 713-714. 10.1093/bioinformatics/btn025.

    Article  CAS  PubMed  Google Scholar 

  34. Dsouza M, Larsen N, Overbeek R: Searching for patterns in genomic data. Trends Genet. 1997, 13 (12): 497-498.

    Article  CAS  PubMed  Google Scholar 

  35. Griffiths-Jones S, Moxon S, Marshall M, Khanna A, Eddy SR, Bateman A: Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids Res. 2005, 33 (Database issue): D121-D124.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Rice P, Longden I, Bleasby A: EMBOSS: the European molecular biology open software suite. Trends Genet. 2000, 16 (6): 276-277. 10.1016/S0168-9525(00)02024-2.

    Article  CAS  PubMed  Google Scholar 

  37. Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14 (6): 787-799. 10.1016/j.molcel.2004.05.027.

    Article  CAS  PubMed  Google Scholar 

  38. Hofacker IL, Fontana W, Stadler PF, Bonhoeffer LS, Tacker M, Schuster P: Fast folding and comparison of RNA secondary structures. Monatshefte für Chemie/Chemical Monthly. 1994, 125 (2): 167-188. 10.1007/BF00818163.

    Article  CAS  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Vencio RZ, Brentani H, Pereira CA: Using credibility intervals instead of hypothesis tests in SAGE analysis. Bioinformatics. 2003, 19 (18): 2461-2464. 10.1093/bioinformatics/btg357.

    Article  CAS  PubMed  Google Scholar 

  41. Romualdi C, Bortoluzzi S, D’Alessi F, Danieli GA: IDEG6: a web tool for detection of differentially expressed genes in multiple tag sampling experiments. Physiol Genomics. 2003, 12 (2): 159-162.

    Article  CAS  PubMed  Google Scholar 

  42. Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120 (1): 15-20. 10.1016/j.cell.2004.12.035.

    Article  CAS  PubMed  Google Scholar 

  43. He L, Hannon GJ: MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004, 5 (7): 522-531. 10.1038/nrg1379.

    Article  CAS  PubMed  Google Scholar 

  44. Schwarz DS, Hutvagner G, Du T, Xu Z, Aronin N, Zamore PD: Asymmetry in the assembly of the RNAi enzyme complex. Cell. 2003, 115 (2): 199-208. 10.1016/S0092-8674(03)00759-1.

    Article  CAS  PubMed  Google Scholar 

  45. Khvorova A, Reynolds A, Jayasena SD: Functional siRNAs and miRNAs exhibit strand bias. Cell. 2003, 115 (2): 209-216. 10.1016/S0092-8674(03)00801-8.

    Article  CAS  PubMed  Google Scholar 

  46. Kozlowski P, Starega-Roslan J, Legacz M, Magnus M, Krzyzosiak WJ: Structures of microRNA precursors. Edited by: Shao-Yao Y. 2008, Springer, 1-16.

    Google Scholar 

  47. MacRae IJ, Zhou K, Doudna JA: Structural determinants of RNA recognition and cleavage by dicer. Nat Struct Mol Biol. 2007, 14 (10): 934-940. 10.1038/nsmb1293.

    Article  CAS  PubMed  Google Scholar 

  48. Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T: Identification of novel genes coding for small expressed RNAs. Science. 2001, 294 (5543): 853-858. 10.1126/science.1064921.

    Article  CAS  PubMed  Google Scholar 

  49. Okamura K, Phillips MD, Tyler DM, Duan H, Chou YT, Lai EC: The regulatory activity of microRNA* species has substantial influence on microRNA and 3′ UTR evolution. Nat Struct Mol Biol. 2008, 15 (4): 354-363. 10.1038/nsmb.1409.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Ruby JG, Stark A, Johnston WK, Kellis M, Bartel DP, Lai EC: Evolution, biogenesis, expression, and target predictions of a substantially expanded set of drosophila microRNAs. Genome Res. 2007, 17 (12): 1850-1864. 10.1101/gr.6597907.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. De Wit E, Linsen SE, Cuppen E, Berezikov E: Repertoire and evolution of miRNA genes in four divergent nematode species. Genome Res. 2009, 19 (11): 2064-2074. 10.1101/gr.093781.109.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


This study was supported by the National Natural Science Foundation of China (NSFC 30625029), and the national science & technology specific projects (2008ZX-10004-011).

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Ning Jiang or Qijun Chen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JW, XL, BJ, HL, SP, XP, NH, PC, JY performed the experiments. NJ and QC designed and supervised the experiments. JW, XL, NJ and QC wrote the manuscript. All authors read and approved the final version of the manuscript.

Electronic supplementary material


Additional file 1: Figure S1. The flow chart for detailed analysis of small RNAs isolated from the two strains of T. gondii. (JPEG 5 MB)


Additional file 2: Table S1. General information of the two libraries. Description: This file contains summary data from high-throughput sequencing of the two small RNA libraries. (DOC 34 KB)


Additional file 3: Table S2. Small RNA classification. Description: This file contains the reads of all small RNA transcripts identified and their relative portions in the library. (DOC 47 KB)


Additional file 4: Figure S2. Statistic analysis of sncRNAs identified in the two libraries. The reads at unique and total levels of the small RNAs in different lengths ranged from 18 to 30 nt were plotted. The difference between the lengths at unique level was significant (p < 0.0001). (TIFF 884 KB)


Additional file 5: Table S3. Length distribution of small RNAs identified in the two strains of T. gondii. Description: This file contains the reads of small RNAs with different lengths and their relative portions in the library. (DOC 57 KB)


Additional file 6: Table S4. Conserved (common) miRNAs and the genomic loci of the encoding genes identified in the two strains of T. gondii. Description: miRNAs with the same sequences could be derived from pre-miRNAs (the sequences may not be completely the same) located at different genomic loci. Their names are indicated with an additional dash-number suffix. (DOC 54 KB)


Additional file 7: Table S5. Novel (unique) miRNAs and the genomic loci of the encoding genes identified in distinct strains of T. gondii. (DOC 462 KB)


Additonal file 8: Table S6. Comparative analysis of the expression of novel miRNAs in the two libraries analyzed by software IDEG6 and Edger respectively. (XLS 12 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Wang, J., Liu, X., Jia, B. et al. A comparative study of small RNAs in Toxoplasma gondii of distinct genotypes. Parasites Vectors 5, 186 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: