Disruption of the Expression of a Non-Coding RNA Significantly Impairs Cellular Differentiation in Toxoplasma gondii

The protozoan parasite Toxoplasma gondii is an important human and veterinary pathogen. Asexual replication of T. gondii in humans and intermediate hosts is characterized by two forms: rapidly growing “tachyzoites” and latent “bradyzoite” tissue cysts. Tachyzoites are responsible for acute illness and congenital neurological birth defects, while the more slowly dividing bradyzoite form can remain latent within the tissues for many years, representing a threat to immunocompromised patients. We have developed a genetic screen to identify regulatory genes that control parasite differentiation and have isolated mutants that fail to convert to bradyzoites. One of these mutants has an insertion disrupting a locus that encodes a developmentally regulated non-coding RNA transcript, named Tg-ncRNA-1. Microarray hybridizations suggest that Tg-ncRNA-1 is involved in the early steps of bradyzoite differentiation. Since Tg-ncRNA-1 does not contain an open reading frame, we used the algorithm Coding Potential Calculator (CPC) that evaluates the protein-coding potential of a transcript, to classify Tg-ncRNA-1. The CPC results strongly indicate that Tg-ncRNA-1 is a non-coding RNA (ncRNA). Interestingly, a previously generated mutant also contains an insertion in Tg-ncRNA-1. We show that both mutants have a decreased ability to form bradyzoites, and complementation of both mutants with wild-type Tg-ncRNA-1 restores the ability of the parasites to differentiate. It has been shown that an important part of bradyzoite differentiation is transcriptionally controlled, but this is the first time that a non-coding RNA is implicated in this process.

transcriptionally controlled, but this is the first time that a non-coding RNA is implicated in this process.

Introduction
Toxoplasma gondii (T. gondii) is an obligate intracellular protozoan parasite that belongs to the phylum Apicomplexa. It infects a wide variety of intermediate hosts, including humans. T. gondii causes congenital hydrocephalus and mental retardation in children, blindness in adults, and encephalitis in immune suppressed patients. This parasite is ubiquitous and the prevalence rate in humans ranges from 15% to 75%, depending on the geography [1,2].
T. gondii asexual reproduction consists of two inter-converting developmental stages: rapidly replicating tachyzoites and a slowly growing bradyzoites that form cysts. Whereas tachyzoites are responsible for disease manifestation and are susceptible to drug therapies, bradyzoites evade host immune responses and are resistant to currently available drugs [3,4]. Tachyzoites disseminate quickly throughout the body and differentiate to bradyzoites that can persist within host tissues as dormant cysts, for months or years. In immunocompromised individuals, bradyzoites de-differentiate to tachyzoites, which can cause a life-threatening disease. Hence, the interconversion process between these two developmental stages is crucial for both persistence and disease.
Tachyzoite-to-bradyzoite differentiation can be studied in vitro [5][6][7]. In vivo, the critical host molecule that drives this switch is interferon gamma [8,9]. In vitro models of differentiation mimic the "stress" of the host immune responses such as treatment with mitochondrial inhibitors, interferon gamma, high pH (8.1), high temperature, nitric oxide, CO 2 starvation. A significant part of bradyzoite differentiation has been shown to be transcriptionally controlled [9][10][11]. Despite these advances, little is known about the genes that play a central role in the regulation of this process.
We and others have previously generated differentiation mutants that fail to convert to bradyzoites [11][12][13][14][15]. However, in a few instances, the disrupted loci have been identified and linked to the mutant phenotype. In the present study, we present the characterization of two independently generated bradyzoite differentiation mutants, B7 and 29C3. These two mutants were isolated in two different genetic screens and two different labs, and importantly, they contain insertions disrupting the same transcript. This transcript is predicted to be a non-coding RNA (ncRNA). Complementing these mutants with the wild-type locus restores the ability of the parasites to differentiate. This surprising result suggests that this ncRNA plays an important role in the cellular differentiation of T. gondii.

Mutant B7 Is Severely Impaired to Form Bradyzoites
In order to understand the genetic basis and mechanisms underlying bradyzoite differentiation, we have developed a genetic screen to identify regulatory genes that control parasite differentiation, and have isolated mutants that fail to convert to bradyzoites [11,12]. In this screen, we took advantage of parasites that lack an endogenous copy of the HXGPRT gene (useful as both positive or negative selectable marker) and reintroduced this marker under the control of a bradyzoite specific promoter. In this background, following insertional mutagenesis, parasites in which a positive regulator of the differentiation signaling cascade is inactivated are resistant to negative selection, under bradyzoite growth conditions. This selection strategy has allowed us to select the population of parasites that form bradyzoites and recover a population of mutant parasites impaired in their ability to differentiate. One of these mutants (mutant B7) was selected for further analysis. Figure 1. Cellular differentiation is impaired in mutant B7. Confluent HFF cells were infected with the wild-type and mutant parasite lines and subjected to bradyzoite differentiation conditions for 72 h as described in materials and methods. (A) Immunofluorescence assay carried out with wild-type (WT) and mutant (B7) parasites after 72 h of bradyzoite growth conditions. Green: bradyzoite antigen 1 (BAG1) [16]. Red: Dolichos lectin (DL) [3]; (B) The presence/absence of marker expression (BAG1 and DL) was counted for 100 parasite vacuoles, in three independent experiments; (C) Proliferation was determined by counting the number of parasites per vacuole (in 100 vacuoles), in three independent experiments; (D) Microarray analysis reveals a global defect in the expression of bradyzoite specific genes. The expression of all the genes up-regulated more than two fold in wild-type parasites (total of 418 genes) was analyzed in the mutant B7. Out of these 418 genes, 297 genes are not up-regulated in the mutant parasites; 101 genes are up-regulated but to a lesser extent compared to the expression in wild-type parasites; and 20 genes are up-regulated to the same levels than wild-type parasites. The significance of the data was determined by Student's t-test (* p < 0.05, ** p < 0.01 and *** p < 0.001). The data are represented as mean ± SEM.
Bradyzoites can be distinguished from tachyzoites based on the expression of bradyzoite specific markers and reduced replication rate [4]. We subjected wild-type and mutant parasites to bradyzoite differentiation conditions and determined the expression of bradyzoite specific markers and growth rate ( Figure 1). The expression of the major bradyzoite antigen (BAG1) and the cyst wall marker Dolichos lectin (DL) is significantly reduced in the mutant B7. After 72 h of in vitro bradyzoite induction, ~80% of the wild-type parasites are positive for both BAG1 and DL, while only ~25% of the mutant parasites express these bradyzoite markers ( Figure 1A,B).
After invasion of a host cell, T. gondii tachyzoites establish an intracellular vacuole where they replicate synchronously with a cell cycle of ~7 h [17]. However, under bradyzoite differentiation conditions, the parasites slow down their replication rate continuously, until they stop replicating in mature cysts. Whereas the tachyzoites completely lyse the host cell monolayer after 72 h of growth, the bradyzoites never lyse the monolayer. Under tachyzoite growth conditions, there is no growth difference between the wild-type and mutant parasites (data not shown). Under bradyzoite differentiation conditions, mutant parasites grow significantly faster, and fail to stop replicating resulting in the lysis of the monolayer. After 72 h, the average number of parasites per vacuole was ~2 and ~16 for the wild-type and B7 mutant respectively ( Figure 1C). Taken together, these results show that mutant B7 has a severe defect to form bradyzoites.
In order to gain a better insight of the impact of the insertion on the global gene expression in mutant B7, we carried out microarray hybridizations with the ToxoGeneChip (details of microarray hybridization experiments can be seen in [11]). The data is archived at NCBI GEO under Series GSE23174. Looking at the whole genome, 418 genes are up-regulated more than two fold in the transition from tachyzotes to bradyzoites ( Figure 1D). Among these 418 genes up-regulated in wild-type bradyzoites, only 20 genes (4.7%) are up-regulated to the wild-type levels in the mutant B7 ( Figure 1D). Another small fraction (24.3%) of bradyzoite specific genes are up-regulated in the mutant, but the induction levels substantially lower than in wild-type parasites. This result suggests that the gene disrupted in the mutant B7, is very likely involved in the early steps of bradyzoite differentiation and, thereby, its loss prevents the expression of a large number of downstream genes.

The Locus Disrupted in Mutant B7 Encodes a Non-Coding RNA
Southern blot analysis shows that the B7 mutant has been disrupted in a single locus, and we identified the genomic DNA flanking the insertion site using a combination of plasmid rescue and inverse PCR [11]. Mutant B7 has an insertion on chromosome VI, which is 345 bp upstream of TGME49_038110, a gene predicted to be a replication factor (Figure 2A). We examined the transcripts expressed from this locus using rapid amplification of cDNA ends (RACE). Using this approach, we discovered that the insertion is directly disrupting a transcript that has no gene prediction in ToxoDB, which is not surprising since it does not contain an open reading frame. RACE also shows that this gene is alternatively spliced, producing two transcripts, one is 2601 bp long and the other is 940 bp long. This gene was named Tg-ncRNA-1 ( Figure 2A).
Real-time PCR shows that Tg-ncRNA-1 is developmentally regulated; its expression is up-regulated 24 fold in the transition from tachyzoites to bradyzoites in wild-type parasites ( Figure 2B). In the mutant bradyzoites, Tg-ncRNA-1 expression is significantly reduced, ~3 fold reduced ( Figure 2B). We examined the expression of the other genes that are close to the insertion site, and none of these genes are either affected in the mutant parasites or developmentally regulated (data not shown). These results strongly suggest that the gene responsible for the mutant phenotype is Tg-ncRNA-1, which was confirmed by complementation assays (see below, Figure 3).  Since Tg-ncRNA-1 does not contain an open reading frame, we suspected that this gene encodes a non-coding RNA (ncRNA). The basis of the classification of a transcript as ncRNA, are extensively discussed in [18]. To evaluate if Tg-ncRNA-1 is indeed a ncRNA, we used the algorithm CPC (Coding Potential Calculator), which evaluates the protein-coding potential of a transcript based on six biologically meaningful sequence features [19]. The results obtained after running CPC are shown in Table 1. In addition of running the sequence of Tg-ncRNA-1, for comparison, we ran two T. gondii protein-coding sequences (TGGT1_087710 and TGME49_000320); and two mouse sequences, one classified as ncRNA (Tsix_mus) and one mouse protein-coding gene (AF282387).
The supporting evidence for the classification of each transcript is shown in the last six columns of Table 1. These six columns show the six features used to distinguish coding from non-coding. A detailed description of each feature can be found at the CPC server [19]. Briefly, a true protein-coding transcript is likely to have more hits with known proteins than an ncRNA, and the hits are likely to have higher quality. The hit number and hit score are measures of these features. Also, for a true protein-coding transcript most of the hits are likely to reside within one frame, and the frame score measures this. The higher the hit score and frame score, the better overall quality of the hits and the more concentrated the hits are in one frame, and therefore, it is more likely that the transcript is protein-coding. The last three features in Table 1 assess the extent and quality of the ORF. A large coverage of the predicted ORF is an indicator of good ORF quality. The Log-odds score is also an indicator of the quality of a predicted ORF, and the higher the score, the higher the quality. Lastly, "type" measures the integrity of the predicted ORF, which indicates if an ORF begins with a start codon and ends with an in-frame stop codon. These six features taken together were used to calculate the coding potential score ( Table 1). The further away the coding potential score is from zero, the more reliable the prediction is. In general, scores that are smaller than −1 or bigger than 1 are significant, and transcripts with a score between −1 and 1 are classified by CPC as "weak non-coding" or "weak coding" [17]. These results strongly indicate that Tg-ncRNA-1 is a ncRNA. Table 1. The coding potential calculator (CPC) shows that Tg-ncRNA-1 is a non-coding RNA. B41(L) and B41(S) correspond to Tg-ncRNA-1 long isoform and short isoform respectively (isoforms shown in Figure 2A). The silencing machinery of T. gondii, which comprises a big diversity of plant-like and metazoan-like small RNAs (sRNAs), has been characterized [20]. Among these diverse sRNAs, there is a group named Toxoplasma REP-derived sRNAS (rdsRNAs) [20]. The rdsRNAs are generated from the repetitive elements REP1 (961 bp), REP2 (1099 bp) and REP3 (1884 bp), which are mitochondrial-like sequences dispersed throughout the nuclear genome of T. gondii [21]. Braun and co-workers suggested that the REP elements are the precursors required for rdsRNA synthesis [18]. Sequence analysis showed that Tg-ncRNA-1 (2601 bp) contains a region (positions 1868-1945 within Tg-ncRNA-1) with high similarity to the T. gondii mitochondrial-like DNA sequences REP1 and REP3. Immediately downstream of these repeats (positions 1967-1991) Tg-ncRNA-1 encodes a 24-mer that belongs to a family of small RNAs cloned and sequenced by Braun et al., family 221/ncRNA S33 (Figure 4). This 24-mer is also included in a region (positions 1945-2049) with high similarity to Plasmodium parasites mitochondrial DNA. Taken together, these findings strongly suggest that Tg-ncRNA-1 encodes a small RNA that belongs to the rdsRNA familiy, here designated Tg-rdsRNA-221 (Figure 4). More experiments need to be done to determine if Tg-rdsRNA-221 behaves functionally as its precursor Tg-ncRNA-1. These experiments will be experimentally challenging, due to it is very difficult to obtain large quantities of RNA from bradyzoites. Of note, the T. gondii silencing machinery was characterized under tachyzoite conditions only [20], but there are no reports of the bradyzoite silencing machinery. Braun et al. showed that the REP elements can encode more than one rdsRNA. Therefore, it is also possible that, in addition to Tg-rdsRNA-221, Tg-ncRNA-1 generates other bradyzoite specific rdsRNAs.  Figure 2A). Tg-ncRNA-1 contains a region with high similarity to the T. gondii mitochondrial-like DNA sequences REP1 and REP3 (REP1/3-yellow), and a region with high similarity to Plasmodium parasites mitochondrial DNA (mitochondrial DNA-red).

Sequence ID C/NC
The numbers indicate the nucleotide positions within Tg-ncRNA-1. The location of the 24-mer rdsRNA is indicated (Tg-rdsRNA-221).

Another Mutant (29C3) That Exhibits a Defect in Bradyzoite Differentiation also Has an Insertion within the Gene Tg-ncRNA-1
We obtained another insertional mutant, named 29C3 (generated by Signature Tagged Mutagenesis), and its corresponding wild-type parental line (named C3), from Laura Knoll's lab [22]. The 29C3 mutant showed a reduced number of cysts in an in vivo mouse model, and it was selected in a screen designed to isolate avirulent T. gondii mutants [20]. We were interested in analyzing 29C3 because it contains an insertion within the Tg-ncRNA-1 transcript. In 29C3, the insertion is located at a different position, compared to mutant B7 ( Figure 2A). Therefore, we examined if 29C3 also shows bradyzoite differentiation defects in vitro.
We subjected 29C3 mutant to bradyzoite differentiation condition in the same way we did for mutant B7. We determined the expression of bradyzoite specific markers and growth profile. We investigated the expression of the major bradyzoite antigen, BAG1 and the presence of Dolichos biflorus lectin staining, a marker that binds to the cyst wall. After 72 h of in vitro bradyzoite induction, the mutant 29C3 shows a reduced expression of the bradyzoite markers ( Figure 5A). The growth rate is also significantly affected in the mutant ( Figure 5B). These results suggest that the mutant 29C3 is defective in bradyzoite differentiation.
The fact that two independent mutants, B7 and 29C3, with an insertion within the gene Tg-ncRNA-1, exhibit an important defect in bradyzoite differentiation, strongly suggests that this gene plays a significant role in cyst formation.

Figure 5. Another independent insertional mutant (29C3) containing an insertion in
Tg-ncRNA-1 is defective in bradyzoite differentiation. Confluent HFF cells were infected with the wild-type (WT) and mutant parasites (29C3) and subjected to bradyzoite differentiation conditions for 72 h. (A) The presence/absence of bradyzoite marker expression (BAG1 and DL) was counted for 100 parasite vacuoles, in three independent experiments; (B) Proliferation was determined by counting the number of parasites per vacuole (in 100 vacuoles), in three independent experiments. The location of the insertion in 29C3 is shown in Figure 2A. The significance of the data was determined by Student's t-test (* p < 0.05, ** p < 0.01 and *** p < 0.001). The data are represented as mean ± SEM.

Complementation of Mutants B7 and 29C3 Restores the WT Phenotype
We used a cosmid (PSBLS72) that spans the insertion sites of both B7 and 29C3 mutants to complement both of them. A cosmid was used in order to capture all possible genomic regulatory regions flanking Tg-ncRNA-1. Stable parasite lines expressing wild-type Tg-ncRNA-1, were obtained by transfecting the cosmid PSBLS72 and selecting with phleomycin [23,24].
These complemented parasite lines were subjected to bradyzoite differentiation conditions to characterize them functionally ( Figure 3). The expression of bradyzoite markers increases significantly in the complemented B7 and 29C3 mutants ( Figure 3A,C), and the growth rate returns to wild-type levels ( Figure 3B,D). The cosmid used for complementation contains a few other genes in addition to Tg-ncRNA-1, but this gene is the only one within the cosmid that has altered expression levels in the mutant parasites; wild-type parasites up-regulate this gene ~24 fold after bradyzoite induction while it's expression is significantly reduced in the mutant bradyzoites, ~3 fold reduced (Figure 2). The other genes present in this cosmid show unaffected expression. These results indicate that Tg-ncRNA-1 is responsible for the phenotype of mutants B7 and 29C3, and complementation with the wild-type locus restores the parasite's ability to form bradyzoites.

Parasite Growth and Differentiation
All the parasites were maintained by serial passage in the confluent Human foreskin fibroblast (HFF) cells as described previously [25]. The parasites were subjected to low CO 2 condition leading to pyrimidine starvation [12] in vitro to induce bradyzoite differentiation. The conditions were achieved by using the Minimum essential medium without NaHCO 3 but supplemented with 25 mM HEPES [12].

Immunofluorescence Assays
HFF cells were grown to confluence on glass coverslips in six well plates (Becton Dickinson, Franklin Lakes, NJ, USA). The confluent HFF cells were infected with equal number of mutant and the parental wild-type (WT) parasites. After 72 h of bradyzoite induction by CO 2 starvation, the coverslips were fixed with 3.7% formaldehyde in phosphate-buffered saline (PBS) for 10 min at room temperature, permeabilized with 0.25% Triton-X-100 in PBS for 15 min at room temperature and blocked with 1× PBS 1% bovine serum albumin (BSA) for 30 min at room temperature. The coverslips were incubated for 1 h at room temperature, with antibody against BAG1 diluted 1:100 in 1× PBS 1% BSA. The coverslips were washed three times and were incubated with a secondary antibody 488 (goat anti-rabbit immunoglobulin G (IgG) conjugated to Alexa fluor 594) (Molecular Probes, Eugene, OR, USA) diluted 1:1000 and tetramethyl rhodamine isothiocyanate (TRITC)-labelled Dolichos biflorus lectin (Sigma-Aldrich, St. Louis, MO, USA) diluted 1:50, for 1 h at room temperature. The coverslips were washed three times, dipped in distilled water, blotted and mounted with Fluoromount-G (SouthernBiotech, Birmingham, AL, USA). The slides were analysed using a Leica DMIRE2 fluorescence microscope (Leica Microsystems Inc., Chatsworth, CA, USA) and images were captured using Improvision Openlab 4.0.2 software (Improvision, Coventry, UK). The samples were examined at 100× magnification and 100 vacuoles were counted. The number of vacuoles positive for BAG1 and DL staining were expressed as % positive vacuoles from duplicate experiments. The number of parasites/vacuole was recorded simultaneously and the average number of parasites/vacuole was determined.

Transfection and Bleomycin (BLE) Selection
Approximately 25-30 μg cosmid DNA with BLE cassette/transfection was resuspended in the required amount of sterile milliQ water to a final volume of 100 μL and 10 μL Sodium Acetate (3M) was added. Absolute ethanol (250 μL) was added to the sample and the sample was frozen at −80 °C for at least 10 min. DNA precipitation was carried out by centrifuging the sample at maximum speed for 15 min at RT. The DNA sample was washed with 1 mL 70% ethanol with two 10 min spins. The precipitated cosmid DNA was resuspended in 100 μL of cytomix (0.015 g ATP and 0.019 g glutathione were added fresh to 12.5 mL of Cytomix for each transfection). Freshly lysed out parasites were harvested and syringe filtered to remove the cell debris. The parasites were counted and spun down. The parasites were resuspended in appropriate volume of the cytomix, to obtain approximately 10 7 parasites/mL.
Cosmid DNA (100 μL) was mixed with parasites (300 μL) into 2 mm gap cuvettes. Parasites were electroporated at 1500 V, 25 Ω, 25 μF, and were added to fresh T25 flasks with fresh ED1 without Phleomycin and incubated at 37 °C. The parasites were harvested after complete lysis, spun down and resuspended in ED1 with 50 μg/mL phleomycin. Parasites were then subjected to extracellular BLE selection by incubating at 37 °C in a waterbath for four hours in 4 mL medium. The parasites were spun down and resuspended in fresh ED1 with 5 μg/mL phleomycin. This sample of parasites was split into two equal halves and were plated into two T25 flasks for intracellular BLE selection. After complete lysis, the parasites were subjected to two more rounds of extracellular and intracellular BLE selection. The resistant parasites were cloned out by limiting dilution and were subjected to appropriate functional analysis.

Quantitative Real Time RT PCR (TaqMan)
Indicated mutant and parental WT parasite lines were cultured under tachyzoite and bradyzoite differentiation conditions for 72 h. At the end of incubation, total RNA was extracted using RNeasy kit (Qiagen Sciences, Germantown, MD, USA). RNA (2 μg) was subjected to DNase (Invitrogen, Carlsbad, CA, USA) treatment for 45 min at room temperature. At the end of treatment, 1 μL of 25 mM EDTA was added to the reaction mixture and subjected to denaturation at 65 °C for 15 min. The primer and probe sets (Sigma Genosys, Haverhill, UK) for the candidate genes were designed using Primer Express software, version 3.0 (Applied Biosystems, Foster City, CA, USA). Real-time quantitative PCR (qPCR) was performed using a 7500 Fast Real-Time PCR system (Applied Biosystems, Foster City, CA, USA) with a reaction mixture volume of 20 μL containing TaqMan Fast Universal PCR 2× Master Mix (Applied Biosystems, Foster City, CA, USA), 900 nM of each primer, 200 nM of each probe and 1 μL of cDNA. The reaction conditions were 95 °C for 20 s, followed by 40 cycles each of 95 °C for 3 s and 60 °C for 30 s. No template controls were run with every experiment and all determinations were performed in duplicates. The relative gene expression levels were calculated as the fold change using the formula 2 −ΔΔCT , where ΔΔCT = ΔΔC T of the target genes − ΔC T of the calibrator and ΔC T = average target C T value − average endogenous control C T value. The constitutively expressed gene DHFR was used as the endogenous control.

Conclusions
Tachyzoite-bradyzoite interconversion plays a critical role in the pathogenesis of T. gondii and the genetic mechanisms that drive this process are not clearly understood [1,4]. We and others have generated bradyzoite differentiation mutants to understand the genetic basis of cellular differentiation in T. gondii [11][12][13][14][15]22]. In this report, we describe the characterization of two of these mutants with insertions disrupting the same transcript: Tg-ncRNA-1. The Coding Potential Calculator (CPC) strongly suggests that Tg-ncRNA-1 encodes a non-coding RNA. This is the first time that a ncRNA has been linked to bradyzoite differentiation in T. gondii. A growing number of studies in different species show that ncRNAs have important regulatory functions, and in many cases are associated with cell differentiation and development [18,26,27]. NcRNAs are classified as long ncRNAs (>200 nucleotides), or short ncRNAs (<200 nt and typically ~20-30 nt long). Many long ncRNAs mediate epigenetic changes recruiting chromatin remodeling complexes [18,28]. Histone modifying complexes have been linked to differentiation in T. gondii [10]. Behnke and co-workers demonstrated that conventional promoter mechanisms work in concert with the chromatin-remodeling machinery to regulate bradyzoite gene expression [10]. However, the mechanism by which the chromatin-remodeling machinery is recruited to the specific genomic loci is not clearly understood. Emerging research suggests that ncRNAs provide the specificity to the genomic loci and it is quite possible that ncRNAs function in concert with chromatin-remodeling complexes to reprogram the gene expression during tachyzoite to bradyzoite differentiation. Hence, it is tempting to implicate Tg-ncRNA-1 in this process, but more work remains to be done to test this hypothesis.
Long ncRNAs are sometimes precursors of small RNAs. The T. gondii genome has good RNAi gene candidates and yet the efficient use of RNAi for specific gene silencing is elusive [29]. Braun and coworkers reported the characterization of the tachyzoite RNA silencing machinery and they did an exhaustive analysis of the small T. gondii RNAome [20]. This RNAome comprises metazoan-like microRNAs (miRNA), and plant-like repeat associated small RNAs (rdsRNA) and satellite associated RNAs (satRNA) [20]. All these RNAs bind to the T. gondii homolog of argonaute (TgAgo). MiRNA-loaded TgAgo associates with polysomes, which likely regulates translation. TgAgo also co-purifies with chromatin-repressing complexes, which suggests a role in transcriptional silencing, most likely through association with rdsRNAs or satRNAs [20]. Our results suggest that Tg-ncRNA-1 may be the precursor of a small RNA that belongs to the rdsRNA family. Therefore, another exciting hypothesis is that rdsRNA members are involved in transcriptional gene silencing during bradyzoite differentiation.
This report provides the first evidence for a direct link between the impaired bradyzoite differentiation phenotype and a non-coding RNA gene. To date, many long and small non-coding RNAs have been sequenced and characterized in many species, but finding a phenotype associated to a ncRNA has not been straightforward. This work opens up an exciting area of research on uncovering the mechanisms by which Tg-ncRNA-1 regulates tachyzoite to bradyzoite differentiation.