Figures
Abstract
African trypanosomes, the causative agents of sleeping sickness in humans and nagana in animals, have a complex digenetic life cycle between a mammalian host and an insect vector, the blood-feeding tsetse fly. Although the importance of the insect vector to transmit the disease was first realized over a century ago, many aspects of trypanosome development in tsetse have not progressed beyond a morphological analysis, mainly due to considerable challenges to obtain sufficient material for molecular studies. Here, we used high-throughput RNA-Sequencing (RNA-Seq) to profile Trypanosoma brucei transcript levels in three distinct tissues of the tsetse fly, namely the midgut, proventriculus and salivary glands. Consistent with current knowledge and providing a proof of principle, transcripts coding for procyclin isoforms and several components of the cytochrome oxidase complex were highly up-regulated in the midgut transcriptome, whereas transcripts encoding metacyclic VSGs (mVSGs) and the surface coat protein brucei alanine rich protein or BARP were extremely up-regulated in the salivary gland transcriptome. Gene ontology analysis also supported the up-regulation of biological processes such as DNA metabolism and DNA replication in the proventriculus transcriptome and major changes in signal transduction and cyclic nucleotide metabolism in the salivary gland transcriptome. Our data highlight a small repertoire of expressed mVSGs and potential signaling pathways involving receptor-type adenylate cyclases and members of a surface carboxylate transporter family, called PADs (Proteins Associated with Differentiation), to cope with the changing environment, as well as RNA-binding proteins as a possible global regulators of gene expression.
Citation: Savage AF, Kolev NG, Franklin JB, Vigneron A, Aksoy S, Tschudi C (2016) Transcriptome Profiling of Trypanosoma brucei Development in the Tsetse Fly Vector Glossina morsitans. PLoS ONE 11(12): e0168877. https://doi.org/10.1371/journal.pone.0168877
Editor: Ben L. Kelly, Louisiana State University, UNITED STATES
Received: July 24, 2016; Accepted: December 7, 2016; Published: December 21, 2016
Copyright: © 2016 Savage et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: RNA-Seq data from this study are available from the NCBI Sequence Read Archive - SRA at http://www.ncbi.nlm.nih.gov/Traces/sra/sra.cgi - under accession numbers SRP002243 and SRR965341.
Funding: This work was supported by National Institutes of Health (http://www.nih.gov) grants AI051584 and AI076879 to S.A. and AI028798 and AI110325 to C.T., and by the Ambrose-Monell Foundation (http://www.monellfoundation.org/) to S.A. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
One of the challenges of being a successful parasitic organism cycling between different hosts is to adapt and survive in drastically different environments. For instance, the protozoan parasite Trypanosoma brucei, causing sleeping sickness in humans and nagana in animals [1, 2], is transmitted between different mammals by its insect vector, the blood-feeding tsetse fly (Diptera: Glossinidae). Successful completion of the life cycle, i.e. survival in and adaptation to new surroundings, requires activation of specific and complex developmental programs which culminate in the manifestation of at least 10 distinct morphological forms [3–5]. In the mammalian bloodstream, a pleomorphic strain exists as both a long slender form that can replicate by asexual division and a cell cycle arrested short stumpy form pre-programmed to encounter the insect host [6]. The transmission cycle begins with the tsetse fly taking a bloodmeal on an infected mammal. While in the midgut of the insect vector, slender forms are killed by the action of proteases and short stumpy forms differentiate into proliferating procyclics. Following establishment of midgut procyclics, trypanosomes must find their way to and colonize salivary glands [7]. This part of the life cycle in the insect vector consists of several morphological forms, including long trypomastigote forms which go through an asymmetric division producing short epimastigotes believed to colonize the salivary gland [8–10]. In this final stage of the T. brucei life cycle in the insect vector, epimastigotes are attached to the epithelium while they differentiate to pre-metacyclics and eventually arrest in G1, before being released as nascent metacyclics, again highlighting a complex differentiation process. The cycle closes when the fly takes another bloodmeal and transfers metacyclics into the bloodstream of the next mammalian host.
During the life cycle, as described briefly above, T. brucei needs to make crucial adaptations to new environments, including different temperatures and nutritional resources, and the parasites need to cope with the immune system in each host. In the mammalian bloodstream, T. brucei replicates extracellularly, and its cell surface is shielded from the environment by a densely packed coat of a single variant surface glycoprotein (VSG). Periodic switching to a different VSG enables the parasite population to evade the host immune response, thus establishing an infection which will be fatal, if not treated [11]. Upon entering the insect host, the VSG coat is replaced by procyclins, a family of glycoproteins characterized by internal Glu-Pro (EP) or Gly-Pro-Glu-Glu-Thr (GPEET) repeats, and at the same time trypanosomes lose their mammalian infectivity [12]. Next, following differentiation to epimastigotes, the surface coat changes to the brucei alanine-rich protein (BARP), so far the only specific molecular marker for proliferating parasites in the salivary glands [13]. The surface remodeling is complete by the acquisition of a VSG coat by the metacyclic forms, which re-establishes infectivity [14]. A second major change during the T. brucei life cycle occurs in the mitochondrial metabolism [15]. Parasites must successfully move from the glucose-rich bloodstream to the tsetse midgut, where proline appears to be an important energy source, although the evidence for this is vague. Thus, procyclics derive their energy mainly by metabolizing amino acids through pathways located in the mitochondrion as well as outside, whereas bloodstreams rely exclusively on glycolysis for energy production and have a mitochondrion with reduced morphological complexity.
This brief introduction only highlights a few of the adaptive changes that need to occur during the T. brucei life cycle. So far, most molecular investigations concentrated on stages that can be cultured in the laboratory, namely the bloodstream and procyclic forms. Unfortunately, this excludes several stages in the insect vector, in particular developmental forms that reacquire infectivity, and experimental work on parasites in the fly is extremely challenging, which resulted in the designation of this stage of the life cycle in a review of this topic as “the heart of darkness” [5]. However, recent technological advances provide a new avenue to interrogate trypanosomes in the insect vector. We hypothesized that differentiation and infection establishment processes in the tsetse are governed by differentially expressed parasite genes. Thus, to characterize the expression profiles of gene products during the T. brucei life cycle in the insect vector, we surveyed the transcriptome using RNA-Seq and Gene Ontology (GO) analysis.
Materials and Methods
Tsetse experiments
All in vivo experiments were performed with pleomorphic T. brucei brucei RUMP 503 and we used the Glossina morsitans morsitans colony maintained in the insectary at Yale University for in vivo experiments. All manipulations, including tsetse infections, tissue dissections and RNA isolations were done as previously described [16]. Flies were dissected after a minimum of 40 days post infection and 72 h after their last blood meal, and infected salivary gland, proventriculus and midgut tissues were collected from the same flies.
cDNA library preparation for RNA-Seq, read processing and differential transcript expression
RNA-Seq library preparation was done following our previously published protocols [17]. RNA samples were treated with Terminator 5′-phosphate-dependent exonuclease and then with RNase-free DNase I. First-strand cDNA synthesis was done with random primers and second-strand cDNA synthesis was initiated with an oligonucleotide complementary to the spliced leader (SL) sequence present at the 5’ end of each trypanosome mRNA, but not of the fly mRNA. The use of the SL oligonucleotide for second-strand synthesis could potentially reduce the yield of double-stranded cDNA, since only a subset of cDNAs will contain the SL sequence at the 5’ end, and thus could skew the representation of reads. However, since the same library protocol was applied to all developmental samples, the results were comparable. Two, four and three independent samples, i.e. biological replicas, were prepared from midgut, proventriculus and salivary gland tissues. Libraries were sequenced on an Illumina GAII platform at the Yale Center for Genome Analysis and the reads of 75 nt in length were pre-processed using the FASTX-toolkit on the public Galaxy webserver ([18–20]; http://galaxyproject.org/) and the SL sequence was trimmed using the Lasergene 12.1 software package from DNASTAR (http://www.dnastar.com/). All processed reads were mapped to the T. brucei 11 megabase chromosomes (GeneDB version 5) using the Lasergene 12.1 software package from DNASTAR. The SeqMan NGen layout algorithm by DNASTAR relies on unique subsequences of bases, or mers, which occur in overlapping regions of fragment reads. Mers that are common to two or more fragment reads are aligned to determine the overall layout of reads. Overlapping reads have many mers in common, but only a few mers per overlapping region are needed to identify the overlap. We chose 21 nt as the minimum length of a mer (overlapping region of a fragment read), in bases, required to be considered a match when arranging reads into contigs. By default, SeqMan NGen uses a local match percentage which requires that the match percentage threshold be met in each overlapping window of 50 bases. For the alignment, a minimum aligned length 35 nt was used and a maximum of two mismatches were allowed. For alignment of reads with multiple matches in the genome, we used two different parameters, namely once, i.e. random distribution, or never, i.e. restriction to non-repeated sequences. As a normalization method, we chose RPKM (reads assigned per kilobase of target per million mapped reads), where the signal values for each experiment are divided by the total bases of target sequence divided by one thousand; and the resulting number divided by the total number of mapped reads divided by one million. RPKM calculations were restricted to open reading frames. We used the Moderated t-Test, which is similar to the Student’s t-Test, and is used to compare the means of gene expression values for two individual replicates or two groups of replicates for a given gene. Whereas the Student’s t-Test calculates variance from the data that is available for each gene, the Moderated t-Test uses information from all of the selected genes to calculate variance. In addition, we used the FDR (Benjamini Hochberg) method as the P-value adjustment method. In this method, the P-values were first sorted and ranked. The smallest value got rank 1, the second rank 2, and the largest got rank N. Then, each P-value was multiplied by N and divided by its assigned rank to give the adjusted P-values. In order to restrict the false discovery rate to 0.05, all the genes with adjusted P-values less than 0.05 were selected. The RPKM expression data were tested by Pearson correlation analysis to evaluate sampling between biological replicates, and all correlation coefficiencies between biological pairs were over 0.81, indicating good reproducibility (S1 Fig). RNA-Seq data from this study have been submitted to the NCBI Sequence Read Archive—SRA at http://www.ncbi.nlm.nih.gov/Traces/sra/sra.cgi—under accession numbers SRP002243 and SRR965341.
Gene Ontology (GO) analysis
The functional annotation analysis of differentially transcribed genes was performed using the Gene Ontology (GO) enrichment tool on the TriTrypDB webserver (http://tritrypdb.org/). GO terms were submitted to REVIGO, a web server that summarizes and condenses long lists of GO terms by removing redundant entries [21]. The analysis was conducted on the separate lists containing the up-regulated and down-regulated transcripts.
Results and Discussion
The repertoire of expressed metacyclic VSG transcripts
One of the hallmarks of trypanosome developmental progression in the insect vector is the re-acquisition of infectivity, which requires the expression of metacyclic VSG (mVSG) genes. Metacyclics are heterogeneous, displaying a number of variants (metacyclic variable antigen types; MVATs), with each individual trypanosome expressing just one mVSG [22]. Estimates for the MVAT diversity, based on monoclonal antibodies, range between more than 14 and a maximum of 27 [23, 24]. However, there is also evidence that at least five MVATs are expressed at very low levels and, four of these, only sporadically [25]. What has been established is that mVSG genes occupy telomeric sites and that they are expressed from a monocistronic RNA polymerase I (Pol I) transcription unit [26, 27].
To begin to address the above issue and to monitor trypanosome gene expression during differentiation in the tsetse vector on a genome-wide scale, we used a modified RNA-Seq approach (Materials and Methods and ref. [17]). Tsetse flies were infected with pleomorphic T. brucei brucei RUMP 503 parasites and following verification of the infection status by microscopic examination, trypanosomes were isolated from midgut, proventriculus and salivary gland tissues [16]. Total RNA was processed for sequencing on the Illumina platform and the resulting reads were aligned to the T. brucei 11 megabase chromosomes (release 5 of the T. brucei genome strain TREU927/4 GUTat10.1). The mVSGs expressed in T. brucei brucei RUMP 503 are not known, and the mVSG repertoire is not present in the genomic sequence of the 11 T. brucei megabase chromosomes, since they are located at telomeres. Thus, to get a handle on the mVSGs expressed in RUMP 503, reads from the salivary gland transcriptome not aligning to the 11 megabase chromosomes were assembled de novo and the resulting contigs were subjected to BLAST analyses of annotated T. brucei proteins at NCBI. Although the exact lineage information on RUMP 503 is not available, it appears to be a derivative of EATRO 795 fly-transmissible ILTat 1.3 (http://tryps.rockefeller.edu/trypsru2_pedigrees.html). Consistent with this history, four metacyclic ILTat VSGs were identified, namely ILTat 1.22, ILTat 1.61, ILTat 1.63 and ILTat 1.64 (Fig 1A). Re-analysis of the transcriptome data with the mVSG sequences appended revealed the four mVSG transcripts to be up-regulated between 57- to 86-fold in salivary glands compared to the proventriculus sample (S1 Table).
(A) Four previously characterized mVSG loci are shown: ILTat 1.22 (ref. [28] and accession no. AJ012198); ILTat 1.61 (ref. [29] and accession no. AJ012199); ILTat 1.63 (ref. [30] and accession no. AJ486954); ILTat 1.64 (ref. [30] and accession no. AJ486955). Solid arrows, VSG genes; black boxes and arrows, mVSG promoters and transcription start sites, respectively; SL (black), position of trans-splices sites detected both previously by RT-PCR [30] and in the present work; SL (red), position of additional trans-splices sites detected in the present study. (B) Three independent infections with RUMP 503 were performed and the relative percentages among the four mVSGs are shown. SG1, SG2 and SG3, salivary gland transcriptomes.
It is thought that each metacyclic cell randomly initiates transcription of a single mVSG gene [22]. Our RNA-Seq data revealed four prominent mVSG transcripts, although we cannot exclude the possibility of additional member(s) being expressed, since we do not have the entire repertoire available. However, what our three independent infections with RUMP 503 clearly showed was that the four mVSG transcripts emerged with comparable representations and that mVSG ILTat 1.22 was the most abundant transcript in all experiments followed by ILTat 1.61 (Fig 1B). This result is comparable to previous investigations on the composition of metacyclic populations in trypanosome stocks isolated from an East African focus of sleeping sickness, where the 1.22 mVSG (also known as GUTat 7.1) and 1.61 mVSG (GUTat 7.15) were the most commonly detected mVSGs [25].
mVSG genes are part of the only bona fide monocistronic transcription units for protein-coding genes and the promoter and transcription start sites for ILTat 1.22, ILTat 1.61, ILTat 1.63 and ILTat 1.64 have been mapped (Fig 1A and ref [30]). Primary transcripts are matured by coupled trans-splicing and polyadenylation. In trans-splicing, the 39-nt SL sequence is added to the 5′ end of all mRNAs and the sequence signals that determine the trans-splice acceptor site appear to consist only of AG dinucleotide at the site for exon junction preceded by a poly-pyrimidine tract of varying length [31]. The RNA-Seq data enabled us to map precisely the 5’ boundaries of the mVSG mRNAs by inspecting the aligned data for regions of rapid changes in the abundance of RNA-Seq tags with the concomitant identification of an AG dinucleotide of the trans-splice acceptor site. In all four mVSGs the most prominent 5’ boundary was very close to the AUG initiation codon (Fig 1A). Nevertheless, we also mapped SL addition sites immediately downstream of the promoter in three cases and these transcripts have been detected previously by RT-PCR in RNA isolated from tsetse-derived metacyclic trypanosomes [30]. Our read coverage was high enough to determine that these transcripts extended from the promoter all the way to the coding region, but it should be noted that the abundance of these transcripts is significantly lower when compared to those covering the ORF (S2 Fig). At present the significance of the stable transcripts downstream of the promoter in mVSG expression remains to be determined.
Estimation of transcript abundance in trypanosomes colonizing different tsetse tissues
Biological replicates of trypanosomes isolated from the midgut, proventriculus and salivary gland were analysed by RNA-Seq for their global expression profile (see Materials and Methods). 10,239 T. brucei annotated and predicted genes were surveyed and transcripts were detected in at least one dataset for 9,699 genes, with 8,759 of these transcripts (90%) having at least 5 RPKM. The vast majority of genes with missing values corresponded to VSG pseudogenes and unlikely hypothetical proteins. The abundance of each transcript is summarized in S2 Table. In general, we achieved very good coverage of the overall transcriptome, i.e. we detected uniform coverage along the 11 mega-chromosomes and the majority of transcription units, with the exception of the subtelomeric regions. The procyclin isoforms (EP1, EP2 and EP3) were by far the most abundant transcripts in the midgut dataset (Table 1 and S3 Table).
Nevertheless, among the 200 most highly expressed mRNAs in the midgut transcriptome, the vast majority were coding for ribosomal proteins (115). Other highly expressed mRNAs originated from genes coding for six amino acid and nucleoside transporters, α- and β-tubulin, four cytochrome oxidase subunits, two RNA binding proteins essential for ribosomal assembly (NRBD1, Tb0927.11.14000, and NRBD2, Tb927.11.14020, aka p34 and p37, ref. [32]) and 27 hypothetical proteins of unknown function. It is worth noting that two transcripts encoding isoforms 1 (Tb927.7.5930) and 2 (Tb927.7.5940) of proteins associated with differentiation or PADs (see below) were in the top 200. Transcripts coding for ribosomal proteins (114/200) and procyclins remained dominant in the proventriculus sample, but one noticeable change was the emergence of transcripts encoding histones H1 to H4 (Tb927.11.1880, Tb927.7.2870, Tb927.1.2550 and Tb927.5.4180), and a transcript encoding a putative zinc finger protein ZC3H36 (Tb927.10.12760). In the salivary gland sample procyclin transcripts were replaced by the brucei alanine-rich protein called BARP [13] and metacyclic VSG as highly expressed transcripts encoding surface proteins. It should be noted that transcripts encoding BARP, a molecular marker for proliferating epimastigotes [13], were as abundant as mVSG transcripts, indicating that the salivary gland transcriptome most likely represented a mixture of epimastigotes and metacyclics. Transporters (5), histones (4), and α- and β-tubulin remained highly abundant, whereas cytochrome oxidase subunits diminished. Finally, PAD1 and PAD2, present in the top 200 in the midgut sample were diminished, but three transcripts encoding isoforms 4 (Tb927.7.5960), 6 (Tb927.7.5980) and 8 (Tb927.7.6000) of proteins associated with differentiation appeared in the list of highly expressed transcripts in the salivary gland data set.
Summary of differential transcript expression and Gene Ontology (GO) analysis
To begin to probe potential molecular differences between trypanosomes isolated from the three developmental stages, we defined differential expression as significant when transcript abundance changed at least two-fold and with a q-value <0.05 (FDR Bonferroni corrected) between two different stages. With this criteria and following filtering of pseudogenes, retrotransposon hot spot (RHS) protein genes and unlikely hypothetical proteins, 1,454 (15%) transcripts revealed significant differential expression in at least one developmental stage comparison (Table 2). In the proventriculus sample 889 transcripts changed abundance relative to the midgut sample, with 670 and 219 transcripts up- and down-regulated, respectively (S1 Table). Although 295 (44%) of the up-regulated transcripts encoded hypothetical proteins, using functional categorization by Gene Ontology (GO) analysis we were able to identify significantly enriched GO terms. Notably, categories involved DNA metabolism, DNA replication, transport and cillary and flagellar motility were at the top of the list (Fig 2B; see S4 Table for full list of GO terms). This is likely a reflection of proliferating cells, which is also evident in the proventriculus set of most abundant transcripts with the appearance of several histone transcripts. GO analysis of the differentially down-regulated transcripts revealed statistically significant biological processes, including carbohydrate catabolism and several transport processes (S4 Table).
(A) A venn diagram of the number of unique and common differentially expressed transcripts between stages. (B) Scatterplot of enriched GO terms in up-regulated transcripts in the proventriculus. (C) Scatterplot of enriched GO terms in down-regulated transcripts in the salivary gland.
In the salivary transcriptome 565 transcripts revealed a significant change, when compared to the proventriculus sample, with 238 transcripts up-regulated, whereas 327 transcripts were down-regulated (Table 2 and S1 Table). GO analysis revealed statistically significant up-regulation of biological processes, including plasma membrane organization, carbohydrate catabolism and several different transport processes. In contrast, in the down-regulated set of transcripts there was very a significant enrichment of categories like signal transduction, cell communication, cyclic nucleotide metabolism and DNA metabolism (Fig 2C).
Differential expression of transcripts encoding surface proteins
The most striking and best understood changes occurring during the trypanosome developmental progression in the insect vector are at the level of the parasite surface. The procyclic form, residing in the midgut is covered with a procyclin coat, whereas epimastigote forms mainly have a coat consisting of BARP [13]. Finally, the infectious metacyclics living in the salivary gland express a single mVSG. Paralleling this protein knowledge, transcripts coding for EP1, EP2 and EP3 procyclin isoforms [33] were highly up-regulated in the midgut transcriptome compared to the other two stages (Fig 3). The abundance of these transcripts decreased steadily during developmental progression with an average 12-fold down-regulation in proventriculus as compared to midgut trypanosomes. On the other hand, transcripts for the GPEET isoform, which coats early procyclic forms [33], were much less abundant and although this value decreased progressively, it was less than two-fold between the three transcriptomes. Transcripts encoding BARP, a specific marker for proliferating epimastigotes [13], were 87-fold up-regulated in the salivary gland transcriptome compared to the midgut sample (Fig 3 and S1 Table). In contrast, these transcripts were only slightly up-regulated in the proventriculus transcriptome (5.6-fold), as compared to the midgut. As described above, the four identified mVSG transcripts were up-regulated between 57- to 86-fold in salivary glands compared to the proventriculus sample. Thus, generally the transcriptome data for the major known surface coat proteins closely mirrored protein expression.
Heatmap showing RPKM changes between midgut (MG), proventriculus (PV) and salivary gland (SG). Heatmap generated with GENE-E (http://www.broadinstitute.org).
Overall, there was a clear trend for differentially expressed transcripts to encode proteins of the African trypanosome cell surface phylome or CSP [34]. Out of 1,454 transcripts 7.8% (192) encode proteins of the CSP. In particular, among the 238 transcripts up-regulated in the salivary gland transcriptome when compared to the proventriculus sample, 71 (30%) transcripts encode proteins of the CSP (S1 Table). In addition to BARP and mVSGs described above, included were transcripts encoding 12 invariant surface glycoproteins, 2 atypical VSGs, 7 transporters, including the bloodstream-specific glucose transporter 1B [35], several expression-site associated genes (ESAGs) and PAD4, 6 and 8 (see below). The invariant surface glycoproteins have been described more than 20 years ago [36] and although it was recently shown that a bloodstream stage-specific ISG75 family mediates suramin uptake [37], these invariant proteins remain poorly characterized and it is not known what role, if any, they play during the trypanosome development cycle. Atypical VSGs are another family of poorly characterized surface proteins and they were named atypical, because they lack cysteine residues in the C terminus [38]. Finally, the salivary gland CSP data set revealed four up-regulated transcripts encoding hypothetical proteins with predicted GPI anchors (Tb927.7.360, Tb927.7.380, Tb927.7.400, Tb927.7.420), that we previously showed by RT-PCR to be up-regulated in trypanosomes residing in salivary glands [16, 39].
One of the best studied differentiations in T. brucei is the transition from stumpy forms to procyclic trypomastigotes, which involves uptake of citrate/cis-aconitate (CCA) by surface-associated carboxylate-transporter family members, namely Proteins Associated with Differentiation or PADs [40]. PAD1 is only expressed at significant levels in stumpy forms and PAD2 is thermoregulated. CCA also appears to affect the activity of TbPTP1, a phosphotyrosine phosphatase, which prevents the development from stumpy to procyclic forms [41]. A downstream substrate for TbPTP1 is a glycosomal directed serine/threonine-specific phosphatase, TbPIP39, which promotes differentiation when it becomes phosphorylated as a consequence of the inhibition of TbPTP1 [42]. Intriguingly, whereas PAD1, PAD2, PAD3, PAD5 and PAD7 were down-regulated between 3.9- and 6.2-fold in the proventriculus transcriptome when compared to that of the midgut, PAD4, PAD6 and PAD8 were 5.4-, 6.0- and 4.5-fold up-regulated, respectively, in trypanosomes residing in salivary glands when compared to the proventriculus transcriptome (S1 Table).
Differential expression of transcripts encoding receptor-type adenylate cyclases (ACs)
The extracellular life cycle of T. brucei requires a mechanism to sense the environment and respond to it appropriately. It is likely that this involves a repertoire of receptors and signaling pathways and there is some evidence that cAMP may play an important role in this cascade. What has been documented is an unusual expansion of transmembrane receptor-like adenylate cyclases (ACs) with the T. brucei genome encoding more than 80 ACs [43]. This expansion is only seen in extracellular African trypanosomes, but not in related kinetoplastids, like T. cruzi and Leishmania, which live intracellularly. One well-studied subfamily of ACs is encoded by one of the expression site associated genes, namely ESAG4, is bloodstream specific, localizes to the flagellar membrane and is involved in host-parasite interaction [44–47]. Most of these ACs are referred to as genes related to ESAG4 (GRESAG4) and previous studies have revealed a spectrum of expression patterns ranging from constitutive transcription throughout the life cycle of the parasite to specific up-regulation in either bloodstream or procyclic trypanosomes [48, 49]. Inspection of the transcriptomes revealed three distinct sets of transcripts encoding ACs. The first group of 14 ACs was highly up-regulated in trypanosomes isolated from the midgut (Fig 4) and their abundance decreased steadily in the proventriculus and salivary gland transcriptome. Among the 14 ACs we identified as being up-regulated in procyclics are three (Tb927.7.7470, Tb927.10.13040, and Tb927.11.13740), which have been shown previously to be specifically expressed in this stage of the life cycle and they were also shown to be localized to the flagellar membrane [50]. Two additional midgut-specific ACs (Tb927.5.285b and Tb927.5.320) have been identified as markers for late procyclic forms [51]. A second and different set of 13 ACs was most prominently expressed in the proventriculus transcriptome with several members previously appearing in flagellum-specific proteomes (Fig 4). Lastly, two ACs (Tb927.7.7520 and Tb927.7.7530) were specifically overexpressed in the salivary gland transcriptome and previous evidence has shown that they are localized to the flagellum [47]. Taken together, our results provide evidence for the importance of distinct ACs when trypanosomes encounter different environments in the tsetse fly.
Heatmap showing RPKM changes between midgut (MG), proventriculus (PV) and salivary gland (SG). Heatmap generated with GENE-E (http://www.broadinstitute.org). Comments based on references [47, 50, 52].
Additional transcript abundance variations highlight developmental changes
Drastically different nutrient availability in the mammalian and insect host necessitates trypanosomes to adjust its energy metabolism [53]. In the midgut of the insect vector trypanosomes rely on oxidative phosphorylation, this occurring within a highly-branched mitochondrion, whereas a switch to glycolysis for ATP production takes place in bloodstream-form trypanosomes, responding to an environment rich in glucose. This adaptation is accompanied by a reduction in size and complexity of the mitochondrion and the expression of high-capacity glucose transporters on the surface of the parasite. The cytochrome oxidase (COX) complex, comprising > 10 subunits, is a key control point in respiratory activity [53, 54]. In contrast, bloodstream-form cells appear not to use a cytochrome-mediated respiratory chain and instead use an alternative pathway consisting of glycerol-3-phosphate dehydrogenase and trypanosome alternative oxidase (TAO). Among the 327 transcripts significantly down-regulated in the salivary gland transcriptome were 22 transcripts encoding known or predicted mitochondrial proteins (S1 Table), including cytochrome oxidase subunit IV (COXIV, Tb927.1.4100), cytochrome oxidase subunit VII (COXVII, Tb927.3.1410), cytochrome oxidase subunit V (COXV, Tb927.9.3170) and cytochrome oxidase subunit IX (COXIX, Tb927.10.8320). On the other hand, the transcript encoding glycerol-3-phosphate dehydrogenase (Tb927.11.7380) increased 2.9-fold from midgut to the proventriculus and then slightly decreased in the salivary gland transcriptome and TAO (Tb927.10.7090) was up-regulated 22-fold in salivary glands as compared to the proventriculus. Taken together, and considering that the salivary transcriptome represents a complex mixture of epimastigotes and metacyclics, our data nevertheless go along with observed morphological changes in mitochondrial structure and upcoming adaptations to the bloodstream environment.
In addition to changes in energy metabolism, different nutrient availability in the insect host could lead trypanosomes to modulate expression levels of various transporters. Indeed, transcripts encoding known and predicted transporters varied significantly in the three transcriptomes (Fig 5). Notably, most glucose transporter transcripts, including the bloodstream form—specific transporter 1B (Tb927.10.8440), were initially high, decreased in the proventriculus transcriptome and then resurfaced in the salivary gland transcriptome.
Heatmap showing RPKM changes between midgut (MG), proventriculus (PV) and salivary gland (SG). Heatmap generated with GENE-E (http://www.broadinstitute.org).
Comparison of the transcriptome from midgut-derived and culture adapted procyclic-form trypanosomes
With the caveat that we do not have a procyclic strain derived from RUMP 503, the strain we used for our fly infection experiments, we compared our previous RNA-Seq data from procyclic form T. brucei rhodesiense YTat 1.1 [17] to the transcriptome of trypanosomes isolated from the midgut (S5 Table). 82 transcripts were up-regulated in the midgut sample, including the procyclin isoforms (EP1, EP2 and EP3), PAD2, 6 and 8, and several receptor-type adenylate cyclases. 73 transcripts were down-regulated in the midgut sample, among them were transcripts encoding the GPEET isoform, histones H1, H2A and H4, and several components of the flagellum. It is worth noting that there was a scarcity of RBPs, with only two putative zinc finger proteins present (Tb927.10.12760 and Tb927.7.690).
How is transcript abundance regulated between life-cycle stages?
The presented catalogue of differentially expressed transcripts did not expose genes involved in the expression of protein-coding genes, i.e. components of the transcriptional apparatus. This goes along with the generally accepted view that T. brucei and related kinetoplastid organisms do not appear to employ transcriptional control to regulate the abundance level of mRNAs and adaptation to different environments. Although various post-transcriptional mechanisms come to mind, in T. brucei there is emerging evidence that RNA-binding proteins (RBPs) play a central role in the regulation of gene expression [55]. A survey of the RNA-Seq data revealed 28 RBPs that were significantly up- or down-regulated in one of the three data sets (Fig 6). Most of these RBPs are still awaiting functional assignment, but notably the list includes ZFP1 (Tb927.6.3490) and ZC3H20 (Tb927.7.2660), two CCCH-family proteins, and ALBA2 (Tb927.11.4450), which were expressed most abundantly in the midgut transcriptome. ZFP1 was one of the first trypanosome RBPs to be implicated in development [56], was shown to be transiently upregulated during the differentiation from bloodstream to procyclic forms and consistent with our data was expressed at higher levels in established procyclic cells. ALBA2 was previously identified to interact with regulatory elements in the 3’ UTR of GPEET procyclin mRNA [57] and ZC3H20 is required for growth of procyclic forms [58], which fits our transcriptome data. Another noteworthy RBP is PUF9 (Tb927.1.2600), which was most prominent in the proventriculus transcriptome. This stage of development was significantly enriched for GO terms involved in DNA metabolism and DNA replication (Fig 2), which fits with PUF9’s known function to bind and stabilize a small number of mRNAs that increase in the late G1 phase of the cell cycle [59].
Heatmap showing RPKM changes between midgut (MG), proventriculus (PV) and salivary gland (SG). Heatmap generated with GENE-E (http://www.broadinstitute.org).
To begin to address the potential function of the identified RBPs in trypanosome development in the insect vector, we chose RBP6 (Tb927.3.2930), whose transcript was up-regulated in the proventriculus compared to the midgut and the transcript level remained elevated in the salivary gland (Fig 6 and S1 Table). RBP6 is a 239 amino acid-long RNA Recognition Motif or RRM protein. Using the pLEW100 vector, we previously overexpressed RBP6 in cultured non-infectious procyclics (Lister 427 29.13.6 cells expressing the TET repressor and T7 RNA polymerase). We showed that doxycycline induction of RBP6 over the course of 9–10 days resulted in the appearance of developmental stages that have been previously described in tsetse [60]. This serendipitous discovery not only strengthens the mechanistic role RBPs play in the regulation of trypanosome gene expression, but also provides a unique system to analyze at the molecular level the events triggered by RBP6 leading to epimastigotes and then to infectious metacyclics.
Conclusions
The developmental program occurring in the insect vector ensures that trypanosomes advance in their life cycle and at the same time become a successful pathogen, i.e. they are adapted to the specific nutritional environment, as well as to develop into the infectious metacyclic form. We have provided a global transcriptome study of T. brucei during its travel through the insect vector by monitoring transcript abundance in trypanosomes isolated from three distinct tissues, i.e. the midgut, proventriculus and the salivary gland. In addition, this approach has very recently been extended to highlight the host-parasite interactions in the tsetse salivary glands [39]. Although we recognize that the salivary gland trypanosome transcriptome is rather complex and impossible to tease apart, since it contains data from proliferating epimastigotes as well as cell-cycle arrested and transmission competent metacyclics, the presented RNA-Seq datasets constitute the first resource of its kind that will facilitate further interrogation of “heart of darkness” T. brucei life-cycle stages in the tsetse vector.
Supporting Information
S1 Table. Differently Expressed Transcripts 2-fold or Higher and p<0.05.
https://doi.org/10.1371/journal.pone.0168877.s003
(XLSX)
S2 Table. All Mapped Transcripts with Biological Replicates.
https://doi.org/10.1371/journal.pone.0168877.s004
(XLSX)
S3 Table. Top 200 most Abundant Transcripts in Midgut, Proventriculus and Salivary glands.
https://doi.org/10.1371/journal.pone.0168877.s005
(XLSX)
S5 Table. Comparison of the transcriptome from midgut-derived and culture adapted procyclic-form trypanosomes.
https://doi.org/10.1371/journal.pone.0168877.s007
(XLSX)
S6 Table. Analysis of the data sets published previously by Telleria et al. 2014 [39].
https://doi.org/10.1371/journal.pone.0168877.s008
(XLSX)
Acknowledgments
We thank Emre Aksoy for generating some of the RNA-Seq libraries and Ying Yang for insectary assistance.
Author Contributions
- Conceptualization: AFS NGK SA CT.
- Formal analysis: JBF CT.
- Funding acquisition: SA CT.
- Investigation: AFS NGK JBF AV.
- Methodology: AFS NGK JBF CT.
- Resources: SA CT.
- Software: JBF CT.
- Writing – original draft: NGK CT.
- Writing – review & editing: AFS NGK JBF SA CT.
References
- 1. Simarro PP, Diarra A, Ruiz Postigo JA, Franco JR, Jannin JG. The human African trypanosomiasis control and surveillance programme of the World Health Organization 2000–2009: the way forward. PLoS Negl Trop Dis. 2011;5(2):e1007. Epub 2011/03/03. pmid:21364972
- 2. Simarro PP, Diarra A, Ruiz Postigo JA, Franco JR, Jannin JG. The human African trypanosomiasis control and surveillance programme of the World Health Organization 2000–2009: the way forward. PLoS Negl Trop Dis. 2011;5(2):e1007. pmid:21364972
- 3. Ooi CP, Bastin P. More than meets the eye: understanding Trypanosoma brucei morphology in the tsetse. Front Cell Infect Microbiol. 2013;3:71. pmid:24312899
- 4. Rico E, Rojas F, Mony BM, Szoor B, Macgregor P, Matthews KR. Bloodstream form pre-adaptation to the tsetse fly in Trypanosoma brucei. Front Cell Infect Microbiol. 2013;3:78. pmid:24294594
- 5. Sharma R, Gluenz E, Peacock L, Gibson W, Gull K, Carrington M. The heart of darkness: growth and form of Trypanosoma brucei in the tsetse fly. Trends Parasitol. 2009;25(11):517–24. pmid:19747880
- 6. MacGregor P, Szoor B, Savill NJ, Matthews KR. Trypanosomal immune evasion, chronicity and transmission: an elegant balancing act. Nat Rev Microbiol. 2012;10(6):431–8. pmid:22543519
- 7. Dyer NA, Rose C, Ejeh NO, Acosta-Serrano A. Flying tryps: survival and maturation of trypanosomes in tsetse flies. Trends Parasitol. 2013;29(4):188–96. pmid:23507033
- 8. Lewis EA, Langridge WP. Developmental forms of Trypanosoma brucei in the saliva of Glossina pallidipes and Glossina austeni. Ann Trop Med Parasitol. 1947;41(1):6–13. pmid:20249284
- 9. Sharma R, Peacock L, Gluenz E, Gull K, Gibson W, Carrington M. Asymmetric cell division as a route to reduction in cell length and change in cell morphology in trypanosomes. Protist. 2008;159(1):137–51. pmid:17931969
- 10. Van Den Abbeele J, Claes Y, van Bockstaele D, Le Ray D, Coosemans M. Trypanosoma brucei spp. development in the tsetse fly: characterization of the post-mesocyclic stages in the foregut and proboscis. Parasitology. 1999;118 (Pt 5):469–78.
- 11. Horn D. Antigenic variation in African trypanosomes. Mol Biochem Parasitol. 2014;195(2):123–9. pmid:24859277
- 12. Vassella E, Den Abbeele JV, Butikofer P, Renggli CK, Furger A, Brun R, et al. A major surface glycoprotein of Trypanosoma brucei is expressed transiently during development and can be regulated post-transcriptionally by glycerol or hypoxia. Genes Dev. 2000;14(5):615–26. pmid:10716949
- 13. Urwyler S, Studer E, Renggli CK, Roditi I. A family of stage-specific alanine-rich proteins on the surface of epimastigote forms of Trypanosoma brucei. Mol Microbiol. 2007;63(1):218–28. pmid:17229212
- 14. Vickerman K. Developmental cycles and biology of pathogenic trypanosomes. Br Med Bull. 1985;41(2):105–14. pmid:3928017
- 15. Priest JW, Hajduk SL. Developmental regulation of mitochondrial biogenesis in Trypanosoma brucei. J Bioenerg Biomembr. 1994;26(2):179–91. pmid:8056785
- 16. Savage AF, Cerqueira GC, Regmi S, Wu Y, El Sayed NM, Aksoy S. Transcript expression analysis of putative Trypanosoma brucei GPI-anchored surface proteins during development in the tsetse and mammalian hosts. PLoS Negl Trop Dis. 2012;6(6):e1708. pmid:22724039
- 17. Kolev NG, Franklin JB, Carmi S, Shi H, Michaeli S, Tschudi C. The transcriptome of the human pathogen Trypanosoma brucei at single-nucleotide resolution. PLoS Pathog. 2010;6(9).
- 18. Goecks J, Nekrutenko A, Taylor J, Galaxy T. Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010;11(8):R86. pmid:20738864
- 19. Giardine B, Riemer C, Hardison RC, Burhans R, Elnitski L, Shah P, et al. Galaxy: a platform for interactive large-scale genome analysis. Genome Res. 2005;15(10):1451–5. pmid:16169926
- 20. Blankenberg D, Gordon A, Von Kuster G, Coraor N, Taylor J, Nekrutenko A, et al. Manipulation of FASTQ data with Galaxy. Bioinformatics. 2010;26(14):1783–5. pmid:20562416
- 21. Supek F, Bosnjak M, Skunca N, Smuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6(7):e21800. pmid:21789182
- 22. Tetley L, Turner CM, Barry JD, Crowe JS, Vickerman K. Onset of expression of the variant surface glycoproteins of Trypanosoma brucei in the tsetse fly studied using immunoelectron microscopy. J Cell Sci. 1987;87(Pt 2):363–72.
- 23. Turner CM, Barry JD, Maudlin I, Vickerman K. An estimate of the size of the metacyclic variable antigen repertoire of Trypanosoma brucei rhodesiense. Parasitology. 1988;97(Pt 2):269–76.
- 24. Lenardo MJ, Rice-Ficht AC, Kelly G, Esser KM, Donelson JE. Characterization of the genes specifying two metacyclic variable antigen types in Trypanosoma brucei rhodesiense. Proc Natl Acad Sci U S A. 1984;81(21):6642–6. pmid:6593722
- 25. Barry JD, Crowe JS, Vickerman K. Instability of the Trypanosoma brucei rhodesiense metacyclic variable antigen repertoire. Nature. 1983;306(5944):699–701. pmid:6656871
- 26. Alarcon CM, Son HJ, Hall T, Donelson JE. A monocistronic transcript for a trypanosome variant surface glycoprotein. Mol Cell Biol. 1994;14(8):5579–91. pmid:8035832
- 27. Graham SV, Barry JD. Transcriptional regulation of metacyclic variant surface glycoprotein gene expression during the life cycle of Trypanosoma brucei. Mol Cell Biol. 1995;15(11):5945–56. pmid:7565747
- 28. Barry JD, Graham SV, Fotheringham M, Graham VS, Kobryn K, Wymer B. VSG gene control and infectivity strategy of metacyclic stage Trypanosoma brucei. Mol Biochem Parasitol. 1998;91(1):93–105. pmid:9574928
- 29. Graham SV, Terry S, Barry JD. A structural and transcription pattern for variant surface glycoprotein gene expression sites used in metacyclic stage Trypanosoma brucei. Mol Biochem Parasitol. 1999;103(2):141–54. pmid:10551359
- 30. Ginger ML, Blundell PA, Lewis AM, Browitt A, Gunzl A, Barry JD. Ex vivo and in vitro identification of a consensus promoter for VSG genes expressed by metacyclic-stage trypanosomes in the tsetse fly. Eukaryot Cell. 2002;1(6):1000–9. pmid:12477800
- 31. Gunzl A. The pre-mRNA splicing machinery of trypanosomes: complex or simplified? Eukaryot Cell. 2010;9(8):1159–70. pmid:20581293
- 32. Prohaska K, Williams N. Assembly of the Trypanosoma brucei 60S ribosomal subunit nuclear export complex requires trypanosome-specific proteins P34 and P37. Eukaryot Cell. 2009;8(1):77–87. pmid:18723605
- 33. Acosta-Serrano A, Vassella E, Liniger M, Kunz Renggli C, Brun R, Roditi I, et al. The surface coat of procyclic Trypanosoma brucei: programmed expression and proteolytic cleavage of procyclin in the tsetse fly. Proc Natl Acad Sci U S A. 2001;98(4):1513–8. pmid:11171982
- 34. Jackson AP, Allison HC, Barry JD, Field MC, Hertz-Fowler C, Berriman M. A cell-surface phylome for African trypanosomes. PLoS Negl Trop Dis. 2013;7(3):e2121. pmid:23556014
- 35. Bringaud F, Baltz T. A potential hexose transporter gene expressed predominantly in the bloodstream form of Trypanosoma brucei. Mol Biochem Parasitol. 1992;52(1):111–21. pmid:1625698
- 36. Ziegelbauer K, Overath P. Identification of invariant surface glycoproteins in the bloodstream stage of Trypanosoma brucei. J Biol Chem. 1992;267(15):10791–6. pmid:1587855
- 37. Alsford S, Eckert S, Baker N, Glover L, Sanchez-Flores A, Leung KF, et al. High-throughput decoding of antitrypanosomal drug efficacy and resistance. Nature. 2012;482(7384):232–6. pmid:22278056
- 38. Weirather JL, Wilson ME, Donelson JE. Mapping of VSG similarities in Trypanosoma brucei. Mol Biochem Parasitol. 2012;181(2):141–52. pmid:22079099
- 39. Telleria EL, Benoit JB, Zhao X, Savage AF, Regmi S, Alves e Silva TL, et al. Insights into the trypanosome-host interactions revealed through transcriptomic analysis of parasitized tsetse fly salivary glands. PLoS Negl Trop Dis. 2014;8(4):e2649. pmid:24763140
- 40. Dean S, Marchetti R, Kirk K, Matthews KR. A surface transporter family conveys the trypanosome differentiation signal. Nature. 2009;459(7244):213–7. pmid:19444208
- 41. Szoor B, Wilson J, McElhinney H, Tabernero L, Matthews KR. Protein tyrosine phosphatase TbPTP1: A molecular switch controlling life cycle differentiation in trypanosomes. J Cell Biol. 2006;175(2):293–303. pmid:17043136
- 42. Szoor B, Ruberto I, Burchmore R, Matthews KR. A novel phosphatase cascade regulates differentiation in Trypanosoma brucei via a glycosomal signaling pathway. Genes Dev. 2010;24(12):1306–16. pmid:20551176
- 43. Salmon D, Bachmaier S, Krumbholz C, Kador M, Gossmann JA, Uzureau P, et al. Cytokinesis of Trypanosoma brucei bloodstream forms depends on expression of adenylyl cyclases of the ESAG4 or ESAG4-like subfamily. Mol Microbiol. 2012;84(2):225–42. pmid:22340731
- 44. Pays E, Tebabi P, Pays A, Coquelet H, Revelard P, Salmon D, et al. The genes and transcripts of an antigen gene expression site from T. brucei. Cell. 1989;57(5):835–45. pmid:2720787
- 45. Paindavoine P, Rolin S, Van Assel S, Geuskens M, Jauniaux JC, Dinsart C, et al. A gene from the variant surface glycoprotein expression site encodes one of several transmembrane adenylate cyclases located on the flagellum of Trypanosoma brucei. Mol Cell Biol. 1992;12(3):1218–25. pmid:1545803
- 46. Salmon D, Vanwalleghem G, Morias Y, Denoeud J, Krumbholz C, Lhomme F, et al. Adenylate cyclases of Trypanosoma brucei inhibit the innate immune response of the host. Science. 2012;337(6093):463–6. pmid:22700656
- 47. Oberholzer M, Langousis G, Nguyen HT, Saada EA, Shimogawa MM, Jonsson ZO, et al. Independent analysis of the flagellum surface and matrix proteomes provides insight into flagellum signaling in mammalian-infectious Trypanosoma brucei. Mol Cell Proteomics. 2011;10(10):M111 010538.
- 48. Alexandre S, Paindavoine P, Tebabi P, Pays A, Halleux S, Steinert M, et al. Differential expression of a family of putative adenylate/guanylate cyclase genes in Trypanosoma brucei. Mol Biochem Parasitol. 1990;43(2):279–88. pmid:1982555
- 49. Alexandre S, Paindavoine P, Hanocq-Quertier J, Paturiaux-Hanocq F, Tebabi P, Pays E. Families of adenylate cyclase genes in Trypanosoma brucei. Mol Biochem Parasitol. 1996;77(2):173–82. pmid:8813663
- 50. Saada EA, Kabututu ZP, Lopez M, Shimogawa MM, Langousis G, Oberholzer M, et al. Insect stage-specific receptor adenylate cyclases are localized to distinct subdomains of the Trypanosoma brucei Flagellar membrane. Eukaryot Cell. 2014;13(8):1064–76. pmid:24879126
- 51. Imhof S, Knusel S, Gunasekera K, Vu XL, Roditi I. Social motility of African trypanosomes is a property of a distinct life-cycle stage that occurs early in tsetse fly transmission. PLoS Pathog. 2014;10(10):e1004493. pmid:25357194
- 52. Subota I, Julkowska D, Vincensini L, Reeg N, Buisson J, Blisnick T, et al. Proteomic analysis of intact flagella of procyclic Trypanosoma brucei cells identifies novel flagellar proteins with unique sub-localization and dynamics. Mol Cell Proteomics. 2014;13(7):1769–86. pmid:24741115
- 53. Hellemond JJ, Bakker BM, Tielens AG. Energy metabolism and its compartmentation in Trypanosoma brucei. Adv Microb Physiol. 2005;50:199–226. pmid:16221581
- 54. Capaldi RA. Structure and function of cytochrome c oxidase. Annu Rev Biochem. 1990;59:569–96. pmid:2165384
- 55. Kolev NG, Ullu E, Tschudi C. The emerging role of RNA-binding proteins in the life cycle of Trypanosoma brucei. Cell Microbiol. 2014;16(4):482–9. pmid:24438230
- 56. Hendriks EF, Robinson DR, Hinkins M, Matthews KR. A novel CCCH protein which modulates differentiation of Trypanosoma brucei to its procyclic form. EMBO J. 2001;20(23):6700–11. pmid:11726506
- 57. Mani J, Guttinger A, Schimanski B, Heller M, Acosta-Serrano A, Pescher P, et al. Alba-domain proteins of Trypanosoma brucei are cytoplasmic RNA-binding proteins that interact with the translation machinery. PLoS One. 2011;6(7):e22463. pmid:21811616
- 58. Ling AS, Trotter JR, Hendriks EF. A zinc finger protein, TbZC3H20, stabilizes two developmentally regulated mRNAs in trypanosomes. J Biol Chem. 2011;286(23):20152–62. pmid:21467035
- 59. Archer SK, Luu VD, de Queiroz RA, Brems S, Clayton C. Trypanosoma brucei PUF9 regulates mRNAs for proteins involved in replicative processes over the cell cycle. PLoS Pathog. 2009;5(8):e1000565. pmid:19714224
- 60. Kolev NG, Ramey-Butler K, Cross GA, Ullu E, Tschudi C. Developmental progression to infectivity in Trypanosoma brucei triggered by an RNA-binding protein. Science. 2012;338(6112):1352–3. pmid:23224556