Skip to main content

Advertisement

Unveiling the partners of the DRBD2-mRNP complex, an RBP in Trypanosoma cruzi and ortholog to the yeast SR-protein Gbp2

Abstract

Background

RNA-binding proteins (RBPs) are well known as key factors in gene expression regulation in eukaryotes. These proteins associate with mRNAs and other proteins to form mRNP complexes that ultimately determine the fate of target transcripts in the cell. This association is usually mediated by an RNA-recognition motif (RRM). In the case of trypanosomatids, these proteins play a paramount role, as gene expression regulation is mostly posttranscriptional. Despite their relevance in the life cycle of Trypanosoma cruzi, the causative agent of Chagas’ disease, to date, few RBPs have been characterized in this parasite.

Results

We investigated the role of DRBD2 in T. cruzi, an RBP with two RRM domains that is associated with cytoplasmic translational complexes. We show that DRBD2 is an ortholog of the Gbp2 in yeast, an SR-rich protein involved in mRNA quality control and export. We used an immunoprecipitation assay followed by shotgun proteomics and RNA-seq to assess the interaction partners of the DRBD2-mRNP complex in epimastigotes. The analysis identified mostly proteins involved in RNA metabolism and regulation, such as ALBA1, ALBA3, ALBA4, UBP1, UBP2, DRBD3, and PABP2. The RNA-seq results showed that most of the transcripts regulated by the DRBD2 complex mapped to hypothetical proteins related to multiple processes, such as to biosynthetic process, DNA metabolic process, protein modification, and response to stress.

Conclusions

The identification of regulatory proteins in the DRBD2-mRNP complex corroborates the important role of DRBD2 in gene expression regulation in T. cruzi. We consider these results an important contribution to future studies regarding gene expression regulation in T. cruzi, especially in the field of RNA-binding proteins.

Background

Trypanosoma cruzi is the causative agent of Chagas’ disease, which affects approximately 7 million people worldwide and is endemic in Latin America [1]. In trypanosomes, the control of gene expression occurs mainly at the posttranscriptional level through messenger RNA (mRNA) processing in the nucleus, its transport to the cytoplasm, translation, and degradation [2, 3].

T. cruzi is exposed to diverse biological conditions in nature because the life cycle of the parasite occurs in two different hosts (mammalian and insect) [4]. The environmental and morphological changes throughout the life cycle of T. cruzi lead to changes in gene expression and metabolic pathways in this organism [5]. RNA-binding proteins (RBPs) are well described as key players in gene expression regulation in eukaryotes [6]. These proteins associate with mRNAs and other proteins to form ribonucleoprotein (mRNP) complexes that ultimately determine the fate of target transcripts in the cell [7]. This association is usually mediated by the RNA-recognition motif (RRM), the most frequent domain in RBPs in eukaryotes [8,9,10].

Despite their importance in mRNA regulation, few RBPs have been characterized in T. cruzi [11,12,13]. Some examples are TcUBP1 and TcUBP2, two RBPs that act in the destabilization and degradation of target transcripts by binding to AU-rich elements (ARE) present in the 3′-UTRs of target mRNAs [14, 15]. TcDHH1 is a cytoplasmic DEAD box helicase involved in mRNA metabolism [16]. TcNRBD1 is associated with the translation machinery; conversely, its ortholog in T. brucei is involved in ribosomal RNA metabolism [11]. The zinc finger protein ZC3H39 in T. cruzi regulates specific sets of transcripts depending on the physiological condition of the parasite [17]; when subjected to nutritional stress, this protein binds to highly expressed mRNA transcripts, such as ribosomal proteins and cytochrome C oxidase. TcRBP9 is an RBP involved in translation regulation in T. cruzi [17]; TcRBP9 associates with proteins involved in RNA metabolism in epimastigotes and with several translation initiation factors in nutritionally stressed epimastigotes. RBSR1 is a recently characterized nuclear RBP [18]; this protein is an ortholog to the human splicing factor SRSF7 and regulates a set of small nucleolar RNAs and small nuclear RNAs, indicating a role in RNA processing in the nucleus.

In T. brucei, a considerable number of RBPs have previously been described [19, 20]. We highlight the ALBA proteins, which are RBPs associated with the translation machinery in T. brucei [21, 22]; ALBAs 1, 2, 3, and 4 are recruited to mRNP complexes under nutritional stress and thus are likely linked to a role in controlling translation. Interestingly, ALBAs 3 and 4 associate with DHH1 to form stress granules [21, 22]. In Leishmania, ALBAs 1 and 3 associate with other RBPs and ribosomal subunits to form a cytoplasmic mRNP that ultimately represses translation [23].

Here, we investigated the role of DRBD2 in T. cruzi, a cytoplasmic RBP with two RRM domains. We show that DRBD2 is an ortholog of Gbp2 in yeast, a shuttling serine/arginine (SR)-rich protein involved in mRNA quality control and export [24]. In T. brucei, previous reports have indicated that DRBD2 destabilizes its target mRNAs, as described in yeast [25]. In our study, a DRBD2 sedimentation profile on a sucrose density gradient demonstrated that the protein is present in translational complexes. With this result as motivation, we used immunoprecipitation assays followed by shotgun proteomics and RNA-seq to assess the partners of the DRBD2-mRNP complex; our analysis identified mostly proteins involved in RNA regulation and metabolism. We highlight the RBPs ALBA1 [22, 23], ALBA3 and ALBA4 [21, 22], UBP1 and UBP2 [14, 15], DRBD3 [19, 26], and the poly(A) binding protein PABP2 [27]. Interestingly, the RNA-seq analysis identified mostly transcripts mapped to hypothetical proteins. Taken together, our results help elucidate the role of DRBD2 in gene expression regulation in T. cruzi through its association to well-known RNA binding proteins.

Results

DRBD2 is an ortholog of the SR-protein Gbp2 in yeast

The results from BLAST together with those from the multiple sequence analysis provided convincing evidence that the T. cruzi DRBD2 (TcCLB.510755.120) protein is an ortholog of Gbp2 proteins in other organisms, such as T. brucei and S. cerevisiae, with approximately 40% similarity between the sequences. However, domain conservation and the protein sizes show that these proteins acquired - or lost - different functions during evolution (Additional file 1: Figure S1).

DRBD2 is an RBP with two RRM domains Fig. 1. Its ortholog in yeast, Gbp2, is an SR-rich protein with three RRM domains (Fig. 1) [28]. The SR-protein Gbp2 forms a mRNP complex that acts in the quality control and export of recently spliced mRNAs [24]; after mRNA transport, the protein remains associated with the transcripts, suggesting a role in translation regulation [29].

Fig. 1
figure1

DRBD2 and yeast Gbp2 alignment and secondary structure maps. In green is the first RRM domain, in blue the second RRM domain, present only in Gbp2 and in yellow the third RRM domain. In bold and underlined the amino acids that are part of the secondary structure (S - β-sheets and H - α-helixes). The Serine residues in red are the phosphorylation sites of DRBD2. Symbols - “*” (asterisk) fully conserved amino acid residues; “:“ (colon) conservation between groups of strongly similar properties; “.” (period) conservation between groups of weakly similar properties

Although DRBD2 does not have SR-rich motifs as in Gbp2, DRBD2 contains several predicted phosphorylation sites in its sequence (Additional file 2: Figure S2). Phosphorylation is the main posttranslational modification responsible for nuclear transport and other cellular processes [30]. The phosphoproteome analysis of T. cruzi identified phosphorylated cysteine residues in DRBD2 from exponentially growing epimastigotes (Fig. 1a) [31]. Based on this result, we hypothesize that DRBD2 is involved in (i) mRNA transport, such as the Gbp2 modus operandi [29], and/or (ii) binding transcripts being exported from the nucleus, as DRBD2 is an RBP with two RRM domains [8].

TAP-tagged DRBD2 is detected throughout metacyclogenesis

TAP-tagged DRBD2 is detected throughout the metacyclogenesis process (Fig. 2a and b). It is a cytoplasmic protein, slightly concentrated around the nucleus in unstressed epimastigotes and metacyclic trypomastigotes (Fig. 2a). This pattern of cellular localization would corroborate the hypothesis raised above, that DRBD2 could have a role in the nuclear transport of its target transcripts. In stressed parasites, however, DRBD2 altered its localization, displaying a granular formation scattered in the cytoplasm (Fig. 2a); this shift in the protein localization suggests an alteration of DRBD2 function during stress. In addition, RNA-seq data available at tritryp database, showed that the drbd2 gene is expressed in all developmental forms of the parasite (Additional file 3: Figure S3 A and B) [32,33,34].

Fig. 2
figure2

DRBD2 cellular localization and expression pattern throughout metacyclogenesis. a Detection of TAP-tagged DRBD2 by indirect immunofluorescence. DIC: differential interference contrast; DAPI: DNA stained with 4′,6-diamidino-2-phenylindole (DAPI); TcDRBD2: endogenous TAP-tagged DRBD2 detected by anti-protein A primary antibody (1:40,000 dilution), and secondary antibody conjugated to Alexa 488 (1:400 dilution); Merge: merge between TcDRBD2 and DAPI images. b Detection of DRBD2 throughout metacyclogenesis by immunoblotting. Epi exp. - Epimastigotes in exponential growth; Epi late - Epimastigotes in late exponential growth; Str - Epimastigotes subjected to nutritional stress; Meta - Metacyclic trypomastigotes. Estimated molecular weight of TAP-tagged DRBD2: 52 kDa; TcDHH1 (control): 47 kDa. M: BenchMark™ Protein Ladder (Thermo Fisher Scientific). Bound antibodies were detected by the Odyssey® imager. Approximately 5,0 × 106 parasites per well were used

DRBD2 is present in translational complexes in unstressed epimastigotes

We investigated if DRBD2 was present in cytoplasmic complexes involved in translation by performing a polysome fractionation of epimastigotes in exponential growth (Fig. 3a) and epimastigotes subjected to a nutritional stress (Fig. 3b). Nutritional stress is one of the key steps that triggers T. cruzi metacyclogenesis [35]. DRBD2 was detected in polysome-enriched and polysome-independent fractions in both conditions analyzed - unstressed and stressed epimastigotes (Fig. 3a and b, respectively).

Fig. 3
figure3

DRBD2 is part of translational complexes in unstressed epimastigotes. Polysome fractionation by sucrose density gradient. The fractions (1–22) were collected after the sedimentation of cytoplasmic lysates from epimastigotes (a) in exponential growth and (b) subjected to nutritional stress, both treated with 100 mg/ml cycloheximide or 2 mM puromycin. The 40S and 60S ribosomal subunits, the 80S ribosome monomer, and polysomes are indicated. Western blotting of alternate fractions was performed with anti-protein A antibody (1:40,000 dilution) for TAP-tagged DRBD2 (52 kDa) detection, and mouse anti-S7 (25 kDa) antibody (1:1000 dilution) for S7 detection, a ribosomal protein commonly used as a control

The sedimentation profile of DRBD2 on the sucrose density gradients of unstressed epimastigotes showed that its migration was similar to the S7 ribosomal protein - used as a control of polysomal fractions - after treatment with both cycloheximide and puromycin (Fig. 3a). These results indicate that DRBD2 is, in fact, associated with the translation machinery in unstressed epimastigotes. In stressed parasites, however, DRBD2 was associated with polysome-independent heavy complexes after puromycin treatment, differently from the control S7 (Fig. 3b), indicating that DRBD2 is not associated with the translation machinery in epimastigotes subjected to nutritional stress.

In yeast, Gbp2 is also present in cytoplasmic mRNP complexes formed by Gbp2, Npl3, and Hrbd1 [29]; this complex represses the translation of the regulated transcripts. During glucose deprivation, Gbp2, Npl3, and Hrbd1 were found in P-bodies and stress granules in yeast [36, 37], suggesting a translation arrest of the associated transcripts in this condition.

DRBD2 associates with proteins involved in RNA metabolism in epimastigotes

The association of DRBD2 with translation complexes in unstressed epimastigotes motivated us to investigate the proteins partners of its mRNP complex. For this, we used immunoprecipitation (IP) assays and analyzed them by shotgun proteomics. Initially, we performed an SDS-PAGE of the IP assay eluates to compare DRBD2 immunoprecipitated complexes with the control, i.e., parasites expressing only TAP-tagging empty vector (Fig. 4). Figure 4a shows the detection of differentially abundant proteins in the DRBD2 immunoprecipitated complex compared to the control IP assay (lanes 2 and 1, respectively).

Fig. 4
figure4

DRBD2 immunoprecipitation assay. a SDS-PAGE stained by silver nitrate of the eluted fractions of (1) the TAP-tagged control and (2) TAP-tagged DRBD2 IP assays. b Detection by western blotting of (1) TAP-tagged control and (2) TAP-tagged DRBD2 in the eluted fractions of the IP assays. Protein detection is indicated in boxes. TAP-tagged DRBD2 (52 kDa) and the TAP-tagged control (21 kDa) were detected by primary rabbit anti-protein A antibody (1:40,000 dilution). Bound antibodies were detected by the alkaline phosphatase method. M: PageRuler prestained Protein Ladder (Thermo Fisher Scientific). The additional protein bands detected in the membrane, at approximately 55 kDa, refer to the IgG chains from the magnetic beads used in the IP assay

The complete proteomic results are presented in Additional file 4: Table S1. Tables 1 and 2 shortlist manually selected results based on abundancy (approximated by NIAF) and having a counterpart with an annotation in T. brucei. These proteins are involved in different processes of gene expression regulation and RNA metabolism in trypanosomatids, such as the RBPs ALBA1, ALBA3 and ALBA4, UBP1 and UBP2, DRBD3, and the poly(A) binding protein PABP2. In addition, we also identified a number of ribosomal proteins, reinforcing the previous observation, from the polysome fractionation, indicating that DRBD2 might be associated to translation regulation (Additional file 4: Table S1).

Table 1 Highlighted unique proteins identified in the DRBD2-mRNP complex. Unique proteins identified only in the TAP-tagged DRBD2 IP assay and not in the control. TriTrypDB ID and protein description: identification number and T. cruzi protein description from TriTryp database. Signal: sum of PatternLab’s normalized label-free quantitation derived from the extracted ion chromatograms - XIC; T. brucei protein description: the corresponding protein description of T. brucei from TriTrypDB
Table 2 Highlighted enriched proteins in the DRBD2-mRNP complex. Enriched proteins identified in TAP-tagged DRBD2 and control IP assays with differential abundance. TriTrypDB ID and protein description: TriTryp database identification number and description for a given protein; Signal control and Signal DRBD2: PatternLab’s Normalized Ion Abundance Factors (NIAF) from control and DRBD2 replicates; Proteins with fold change > 2 were considered; Product: a given name for the corresponding protein

Assessment of transcripts regulated by the DRBD2-mRNP complex

We performed RNA-seq analysis of the immunoprecipitated DRBD2-mRNP complex to assess the regulated transcripts. The analysis pinpointed 86 transcripts as more abundant in the DRBD2-mRNP in epimastigotes when compared to the TAP-tag strain control, with an absolute fold-change of at least two (Additional file 5: Table S2 and Additional file 7: Table S4); however, most transcripts mapped to hypothetical proteins. To obtain further information regarding the transcripts, we performed gene ontology analysis (Fig. 5). Most of the mRNAs code for proteins related to biosynthetic process, DNA metabolic process, protein modification, and response to stress (Fig. 5). The complete information of the RNA-seq analysis is found in Additional file 5: Table S2, Additional file 6: Table S3, Additional file 7: Table S4. Table 3 lists the most abundant transcripts identified in this work.

Fig. 5
figure5

Gene ontology classification for biological processes of the identified transcripts of the DRBD2-mRNP complex. The numbers indicate the score of each gene ontology term and the percentage of each term in the group of transcripts being analyzed

Table 3 The most abundant transcripts associated with the DRBD2-mRNP complexes in epimastigotes. Gene ID: T. cruzi Dm28c strain identification number from the TriTryp database; Product description: resulting protein description provided by TriTrypDB

We then used the results from the T. cruzi ribosome profiling by Smircich et al. [34] to verify if the transcripts identified in this work are associated with polysomes (Additional file 5: Table S2 and Additional file 7: Table S4). Interestingly, some of the most abundant transcripts associated with the DRBD2-mRNP complex showed high levels of expression in the total RNA fraction rather than in the polysomes, in both epimastigotes and metacyclic trypomastigotes (Additional file 3: Figure S3C), indicating a possible role of DRBD2 in negatively regulating the translation of its target transcripts, as described in T. brucei [25].

Discussion

RBPs play lead roles in gene expression regulation in trypanosomatids [10]. The increasing knowledge about the molecular functioning of these proteins in T. cruzi has helped us to gradually understand how the parasite copes to survive facing various biological diversities throughout its life cycle [38].

Here, we studied the RBP DRBD2 in T. cruzi, an ortholog of the SR-rich protein Gbp2 in yeast, which forms a complex with Hrbd1 and Npl3 to act in the quality control and export of nascent mRNAs [24, 29]. DRBD2 is a cytoplasmic protein that shifts its cellular localization to a granular pattern when subjected to a nutritional stress. Additionally, we showed that DRBD2 participates in translational complexes in epimastigotes, according to its sedimentation profile on a sucrose gradient. The shotgun proteomic analysis identified mostly proteins involved in RNA metabolism and regulation; some examples are ZC3H39, UBP1, UBP2, ALBA3, ALBA4, and NRBD1.

In Leishmania, ALBA1, 3 and 4 coimmunoprecipitate with DRBD2 in amastigotes [23]; Ntf2-like and PABP2 also associate with ALBA proteins, suggesting that a protein complex is formed by these partners. These results are consistent with our results; PABP2, Ntf2-like, ALBA1, 3, and 4 were also detected in the DRBD2-mRNP complex in T. cruzi. Moreover, ALBA proteins are associated with the translation machinery in T. brucei [21, 22], such as DRBD2 in T. cruzi. When T. brucei parasites were subjected to a nutritional stress, ALBA3 and ALBA4 associated with DHH1 to form cytoplasmic stress granules, resulting in a reduced translation rate [22]. Interestingly, we observed that DHH1, ALBA3, and ALBA4 associate with the DRBD2-mRNP complex.

In T. brucei, DRBD3 is an RBP that regulates transcripts of ribosomal proteins and translation factors and is thus involved in translation control [19]. Furthermore, TbDRBD3 is part of a cytoplasmic mRNP complex with PABP2 and Ntf2-like [26], which is consistent with our results (Tables 1 and 2). Interestingly, Hel67 was also associated with the TcDRBD2 complex; in Leishmania, Hel67 is a cytoplasmic ATP-dependent DEAD-box helicase that binds to specific ribosomal RNAs, preventing its degradation and ultimately inhibiting cell death [39].

The poly(A)-binding protein PABP2 is highly abundant in the DRBD2 complex (Additional file 4: Table S1). In T. brucei, PABP2 may accumulate in nuclear periphery granules, most likely regulating the nuclear export of transcripts [40]. TbDRBD2 is also present in the TbPABP2 complex, together with TbZC3h40, TbALBA1, and TbALBA3 [41]; these results are consistent with ours (Tables 1 and 2).

Finally, the elongation factor 1 alpha (EF-1α) is also present in the DRBD2-mRNP complex. This protein is conserved in eukaryotic organisms [42]; it is a typical moonlighting protein and is involved in cellular processes such as translation and RNA transport [43]. In T. cruzi, EF-1α is associated with both polysome and polysome-free fractions in epimastigotes 6; also, the EF-1α-mRNP complex seems to play a role in stress response by regulating different sets of transcripts in this condition [42]. In addition, ribosomal proteins are abundantly present in ribonucleoprotein (mRNP) complexes because they are usually involved in translation, hence they are considered an indicative of translation regulation by these complexes, as described for NRBD1 from T. cruzi [11].

The RNA-seq results showed that most of the transcripts regulated by the DRBD2 complex mapped to hypothetical proteins related to multiple processes, such as to biosynthetic process, DNA metabolic process, protein modification, and response to stress.

Conclusion

RBPs are key proteins in gene expression regulation of T. cruzi. The information regarding the actors and mechanisms in this field are continually increasing. DRBD2’s orthology to yeast Gbp2 is a promising indicative of the similarity between these proteins roles in negatively regulating their targets in mRNP complexes. Moreover, DRBD2 seems to be involved in RNA metabolism regulation in T. cruzi once several regulatory proteins were identified by the shotgun proteomic analysis of the immunoprecipitated complexes.

We consider these results an important contribution to future studies regarding gene expression regulation in T. cruzi, especially in the field of RNA-binding proteins.

Methods

In silico analysis

Similarities between the DRBD2 (TcCLB.510755.120) sequence from T. cruzi and the Gbp2 sequences from several other organisms from NCBI Refseq and UniProt SwissProt (accessed in September, 2018), including Trypanosoma brucei, Saccharomyces cerevisiae, and Homo sapiens, were assessed using the BLAST algorithm [44]. This analysis was achieved by querying the TcDRBD2 sequence against other organism genomes, followed by inverse alignment using the best candidates from the first search, using the best reciprocal hit (BRH) method [45]. Plus, we used Constraint-based Multiple Alignment Tool (COBALT), contained in NCBI platform, to perform a multiple protein sequence alignment using conserved domain and local sequence similarity information. To predict phosphorylation sites in the DRBD2 sequence, we used the NetPhos 2.0 Server [46] (http://www.cbs.dtu.dk/services/NetPhos/). Secondary structure mapping was obtained with PSIPRED [47].

Parasite cultures and drug treatments

T. cruzi Dm28c strain was originally isolated from Didelphis marsupialis [35, 48] and epimastigotes forms were maintained in axenic culture in liver infusion tryptose (LIT) medium at 28 °C [48]. To prepare 1.0 L of LIT medium, we used 9.0 g of liver infusion broth, 5.0 g of tryptose, 1.0 g of NaCl, 8.0 g of Na2HPO4, 0.4 g of KCl, 1.0 g of glucose, 100 ml of fetal bovine serum (heat-inactivated), 10 mg of hemin, and 1.0 L of distilled water. For nutritional stress induction, epimastigotes in the late logarithmic growth phase were harvested by centrifugation at 7000 x g for 5 min at 10 °C and incubated for 2 h at 28 °C in TAU medium (190 mM NaCl, 17 mM KCl, 2 mM MgCl2, 2 mM CaCl2, 8 mM phosphate buffer pH 6.0) at a density of 5.0 × 108 parasites/ml. Metacyclic trypomastigotes were prepared as previously described [48].

Recombinant protein expression

The T. cruzi Dm28c strain was used as the DNA template. The DRBD2 (TriTrypDB ID TcCLB.510755.120) open reading frame polymerase chain reaction (PCR) product was cloned into the pDONR™221 vector from Gateway Technology (Invitrogen) and was then recombined into a pTcTAPN vector that adds the Tap-tag in the N-terminus portion of the protein [49]. T. cruzi Dm28c epimastigotes were transfected with these plasmids as described by Lu et al. [50]. Stable lines were selected by adding 500 mg/ml G418 to the culture medium. It is episomally expressed, not recombining into a genomic locus (DRBD2 molecular weight: 31 kDa; TAP-tagging: 21 kDa).

Immunoblotting and immunofluorescence

Immunoblotting and immunofluorescence were performed as previously described [51]. A rabbit anti-protein A antibody (Sigma-Aldrich) was used to detect TAP-tagging (1:40,000 dilution). TcDHH1 (47 kDa) was used as a control of the metacyclogenesis western blotting assay. Anti-TcDHH1 (1:100 dilution) was kindly provided by Jimena Ferreira da Costa. Bound antibodies were detected by the alkaline phosphatase method or by the Odyssey® imager (LI-COR) using the secondary antibodies IRDye 800CW or IRDye 680RD (1:15,000 dilution) [52]. Immunofluorescence slides were analyzed by inverted microscopy (Leica DMI6000B) associated with deconvolution software (Leica AF6000, microscope facility RPT07C PDTIS/ Carlos Chagas Institute - Fiocruz Paraná).

Polysome sedimentation profiles

Polysomes were separated on a sucrose density gradient as previously described, performed once for each condition [53]. The drugs used here were puromycin, which causes premature chain termination during translation, thus disassembling polysomes [54]; and cycloheximide, that interferes with the translocation step in protein synthesis, thus blocking eukaryotic translational elongation and stabilizing polysomes [54]. The resulting fractions were analyzed by western blotting using rabbit anti-protein A (1:40,000 dilution) for DRBD2 detection and mouse anti-S7 (1:1000 dilution), kindly provided by Dr. Stênio P. Fragoso. S7 (25 kDa) is a small ribosomal protein commonly used as a control. The fractionation was performed twice for each treatment and condition.

Immunoprecipitation assays

Immunoprecipitation (IP) assays were performed with cytoplasmic extracts of exponentially growing epimastigotes. Primary rabbit anti-protein A antibodies were coupled to 30 mg of Dynabeads® M-270 Epoxy (Thermo Fisher), following the manufacturer’s instructions. The equivalent of 5.0 × 109 parasites expressing TAP-tagged DRBD2 and control parasites expressing only the TAP-tag were cryo-milled using the Mixer Mill MM 400 (Retsch), resulting in 1 g of cell powder. Approximately 5.0 × 108 parasites were resuspended in 1 ml of solubilization buffer [25 mM and 50 mM sodium citrate; 20 mM Tris-HCl pH 7.5; 5 mM MgCl2; 0.1% CHAPS; protease inhibitor (Complete Mini Protease inhibitor cocktail tablet, Roche, 1:100) and 40 U/ml RNase OUT (Invitrogen)], followed by centrifugation at 20,000 x g for 10 min at 4 °C. The supernatant was incubated with 3 μl of the magnetic beads under constant stirring at 4 °C for 2 h. The magnetic beads were washed 3 times with the solubilization buffer, and the proteins linked to the beads were eluted with 30 μl of elution buffer (2% SDS and 20 mM Tris-HCl pH 8.0) at 72 °C for 20 min. The eluted fractions were visualized using 13% SDS-PAGE with silver nitrate staining, as previously described [55], and then, TAP-tagged DRBD2 was detected by immunoblotting.

Shotgun proteomic data acquisition and analysis

Digested extracts of the immunoprecipitated lysates were fractionated by a NanoLC Easy 1000 system (Thermo Fisher Scientific). Ionized peptides were analyzed by an LTQ Orbitrap XL ETD (Thermo Fisher Scientific). Sequences of T. cruzi were downloaded from TriTrypDB on March 16, 2016. PatternLab for proteomics 4.0 was used for data analysis, following the provided instructions [56]. Briefly, the search was performed using PatternLab’s integrated Comet [57] search engine. The validity of the Peptide Sequence Matches (PSM) was assessed using PatternLab’s Search Engine Processor module (SEPro) [58]. A cutoff score was established to accept a false discovery rate (FDR) of 1% based on the number of decoys [59]. The results were postprocessed to accept PSM with less than 5 ppm from the average ppm of the pre-validated results. Relative quantitation was performed according to PatternLab’s Normalized Ion Abundance Factors (NIAF) as a relative quantitation strategy [56]. For proteins present in both conditions (DRBD2 assay and control), we used PatternLab’s TFold module [60], with an FDR of 0.05. Only proteins having an absolute fold change ≥ 2 were considered for the TFold analysis.

Ribonomic data acquisition and analysis

The eluted fractions from the immunoprecipitation assays were purified with the miRCURY™ RNA Isolation Kit (Exiqon), with the “Cell & Plant” protocol from the manufacturer’s manual, initiating with the “Cell Lysis” step in the manual. The purified RNAs were quantified and analyzed using the Agilent DNA 1000 kit on an Agilent Technologies 2100 Bioanalyzer and then subjected to sequencing on the Illumina® platform using paired-end configuration on the MiSeq Desktop Sequencer (Illumina®). The sequencing data obtained from 3 biological replicates were analyzed using CLC Genomics Workbench© 7.5.5. The T. cruzi Dm28c reference genome used for mapping was obtained from the NCBI database (AYLP01), and the alignment was performed as follows: additional upstream and downstream sequences of 75 bases, minimum number of reads: 10, maximum number of mismatches: 2, and nonspecific match limit: - 2. We selected possible targets of the DRBD2-mRNP with the β binomial statistical test (Baggerly’s test [61]) with a q-value ≤1%, a minimum RPKM value of 50 and a fold change > 5 when compared to the control (parasites expressing TAP-tag) for significance. The RNA-seq data were deposited in the NCBI Sequence Read Archive (SRA) database under the accession number SRX4560510.

Availability of data and materials

The RNA-seq data were deposited in the NCBI Sequence Read Archive (SRA) database under the accession number SRX4560510. The proteomic raw data is available at http://proteomics.fiocruz.br/supplementaryfiles/wippel2018/.

Abbreviations

ARE:

AU-rich elements

FDR:

False discovery rate

IP:

Immunoprecipitation

LIT:

Liver infusion tryptose

mRNA:

messenger RNA

mRNP:

ribonucleoprotein

NIAF:

Normalized Ion Abundance Factors

PCR:

Polymerase chain reaction

PSM:

Peptide Sequence Matches

RBP:

RNA-binding proteins

RRM:

RNA Recognition Motif

SEPro:

Search Engine Processor

SR:

Serine-arginine

SRA:

Sequence Read Archive

References

  1. 1.

    Rassi A Jr, Rassi A, Marin-Neto JA. Chagas disease. Lancet. 2010;375:1388–402.

  2. 2.

    Clayton CE. Life without transcriptional control? From fly to man and back again. EMBO J. 2002;21:1881–8. https://doi.org/10.1093/emboj/21.8.1881.

  3. 3.

    Alves LR, Guerra-Slompo EP, De Oliveira AV, Malgarin JS, Goldenberg S, Dallagiovanna B. MRNA localization mechanisms in Trypanosoma cruzi. PLoS One. 2013;8:1106–18.

  4. 4.

    Bartholomeu DC, de Paiva RMC, Mendes TAO, DaRocha WD, Teixeira SMR. Unveiling the intracellular survival gene kit of Trypanosomatid parasites. PLoS Pathog. 2014;10:e1004399.

  5. 5.

    Glisovic T, Bachorik JL, Yong J, Dreyfuss G. RNA-binding proteins and post-transcriptional gene regulation. FEBS Lett. 2008;582:1977–86. https://doi.org/10.1016/j.febslet.2008.03.004.

  6. 6.

    Alves LR, Ávila AR, Correa A, Holetz FB, Mansur FCB, Manque PA, et al. Proteomic analysis reveals the dynamic association of proteins with translated mRNAs in Trypanosoma cruzi. Gene. 2010;452:72–8.

  7. 7.

    Keene JD. RNA regulons: coordination of post-transcriptional events. Nat Rev Genet. 2007;8:533–43. https://doi.org/10.1038/nrg2111.

  8. 8.

    Cléry A, Blatter M, Allain FH. RNA recognition motifs: boring? Not quite. Curr Opin Struct Biol. 2008;18:290–8. https://doi.org/10.1016/j.sbi.2008.04.002.

  9. 9.

    Dreyfuss G, Kim VN, Kataoka N. Messenger-RNA-binding proteins and the messages they carry. Nat Rev Mol Cell Biol. 2002;3:195–205.

  10. 10.

    Kolev NG, Ullu E, Tschudi C. The emerging role of RNA-binding proteins in the life cycle of Trypanosoma brucei. Cell Microbiol. 2014;16:482–9.

  11. 11.

    Oliveira C, Carvalho PC, Alves LR, Goldenberg S. The role of the trypanosoma cruzi TcNRBD1 protein in translation. PLoS One. 2016;11:e0164650.

  12. 12.

    Oliveira C, Faoro H, Alves LR, Goldenberg S. RNA-binding proteins and their role in the regulation of gene expression in trypanosoma cruzi and saccharomyces cerevisiae. Genet Mol Biol. 2017;40:22–30.

  13. 13.

    Tyler Weisbarth R, Das A, Castellano P, Fisher MA, Wu H, Bellofatto V. The Trypanosoma cruzi RNA-binding protein RBP42 is expressed in the cytoplasm throughout the life cycle of the parasite. Parasitol Res. 2018;117:1095–104.

  14. 14.

    D’Orso I, Frasch ACC. TcUBP-1, an mRNA destabilizing factor from trypanosomes, homodimerizes and interacts with novel AU-rich element- and poly(a)-binding proteins forming a ribonucleoprotein complex. J Biol Chem. 2002;277:50520–8.

  15. 15.

    Noé G, De Gaudenzi JG, Frasch AC. Functionally related transcripts have common RNA motifs for specific RNA-binding proteins in trypanosomes. BMC Mol Biol. 2008;9:107.

  16. 16.

    Holetz FB, Alves LR, Probst CM, Dallagiovanna B, Marchini FK, Manque P, et al. Protein and mRNA content of TcDHH1-containing mRNPs in Trypanosoma cruzi. FEBS J. 2010;277:3415–26.

  17. 17.

    Wippel HH, Inoue AH, Vidal NM, Costa JFD, Marcon BH, BAA R, et al. Assessing the partners of the RBP9-mRNP complex in Trypanosoma cruzi using shotgun proteomics and RNA-seq. RNA Biol. 2018;15(8):1106–18.

  18. 18.

    Wippel HH, Malgarin JS, Martins SDT, Vidal NM, Marcon BH, Miot HT, et al. The Nuclear RNA-binding Protein RBSR1 Interactome in Trypanosoma cruzi. J Eukaryot Microbiol. 2018;66(2):244–53.

  19. 19.

    Das A, Bellofatto V, Rosenfeld J, Carrington M, Romero-Zaliz R, Del Val C, et al. High throughput sequencing analysis of Trypanosoma brucei DRBD3/PTB1-bound mRNAs. Mol Biochem Parasitol. 2015;199:1–4. https://doi.org/10.1016/j.molbiopara.2015.02.003.

  20. 20.

    Clayton C. The regulation of trypanosome gene expression by RNA-binding proteins. PLoS Pathog. 2013;9(11):e1003680.

  21. 21.

    Mani J, Güttinger 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:e22463. https://doi.org/10.1371/journal.pone.0022463.

  22. 22.

    Subota I, Rotureau B, Blisnick T, Ngwabyt S, Durand-Dubief M, Engstler M, et al. ALBA proteins are stage regulated during trypanosome development in the tsetse fly and participate in differentiation. Mol Biol Cell. 2011;22:4205–19.

  23. 23.

    Dupé A, Dumas C, Papadopoulou B. Differential subcellular localization of Leishmania Alba-domain proteins throughout the parasite development. PLoS One. 2015;10:e0137243.

  24. 24.

    Hackmann A, Wu H, Schneider U-M, Meyer K, Jung K, Krebber H. Quality control of spliced mRNAs requires the shuttling SR proteins Gbp2 and Hrb1. Nat Commun. 2014;5:3123.

  25. 25.

    Lueong S, Merce C, Fischer B, Hoheisel JD, Erben ED. Gene expression regulatory networks in T rypanosoma brucei : insights into the role of the mRNA-binding proteome: gene expression regulatory networks in Trypanosoma brucei. Mol Microbiol. 2016;100:457–71.

  26. 26.

    Fernández-Moya SM, García-Pérez A, Kramer S, Carrington M, Estévez AM. Alterations in DRBD3 ribonucleoprotein complexes in response to stress in Trypanosoma brucei. PLoS One. 2012;7:e48870. https://doi.org/10.1371/journal.pone.0048870.

  27. 27.

    Anderson P, Kedersha N. RNA granules. J Cell Biol. 2006;172:803–8. https://doi.org/10.1083/jcb.200512082.

  28. 28.

    Windgassen M, Krebber H. Identification of Gbp2 as a novel poly(a)+ RNA-binding protein involved in the cytoplasmic delivery of messenger RNAs in yeast. EMBO Rep. 2003;4:278–83.

  29. 29.

    Windgassen M, Sturm D, Cajigas IJ, Gonzalez CI, Seedorf M, Bastians H, et al. Yeast shuttling SR proteins Npl3p, Gbp2p, and Hrb1p are part of the translating mRNPs, and Npl3p can function as a translational repressor. Mol Cell Biol. 2004;24:10479–91.

  30. 30.

    Nardozzi JD, Lott K, Cingolani G. Phosphorylation meets nuclear import: a review. Cell Commun Signal. 2010;8:32.

  31. 31.

    Marchini FK, de Godoy LMF, Rampazzo RCP, Pavoni DP, Probst CM, Gnad F, et al. Profiling the Trypanosoma cruzi Phosphoproteome. PLoS One. 2011;6:e25381.

  32. 32.

    Belew AT, Junqueira C, Rodrigues-Luiz GF, Valente BM, Oliveira AER, Polidoro RB, et al. Comparative transcriptome profiling of virulent and non-virulent Trypanosoma cruzi underlines the role of surface proteins during infection. PLoS Pathog. 2017;13:e1006767.

  33. 33.

    Li Y, Shah-Simpson S, Okrah K, Belew AT, Choi J, Caradonna KL, et al. Transcriptome remodeling in Trypanosoma cruzi and human cells during intracellular infection. PLoS Pathog. 2016;12:1005511.

  34. 34.

    Smircich P, Eastman G, Bispo S, Duhagon MA, Guerra-Slompo EP, Garat B, et al. Ribosome profiling reveals translation control as a key mechanism generating differential gene expression in Trypanosoma cruzi. BMC Genomics. 2015;16:443.

  35. 35.

    Contreras VT, Araujo-Jorge TC, Bonaldo MC, Thomaz N, Barbosa HS, Meirelles MN, et al. Biological aspects of the Dm 28c clone of Trypanosoma cruzi after metacyclogenesis in chemically defined media. Mem Inst Oswaldo Cruz. 1988;83:123–33.

  36. 36.

    Mitchell SF, Jain S, She M, Parker R. Global analysis of yeast mRNPs. Nat Struct Mol Biol. 2013;20:127–33.

  37. 37.

    Buchan JR, Muhlrad D, Parker R. P bodies promote stress granule assembly in Saccharomyces cerevisiae. J Cell Biol. 2008;183:441–55.

  38. 38.

    Alves LR, Goldenberg S. RNA-binding proteins related to stress response and differentiation in protozoa. World J Biol Chem. 2016;7:78–87. https://doi.org/10.4331/wjbc.v7.i1.78.

  39. 39.

    Padmanabhan PK, Samant M, Cloutier S, Simard MJ, Papadopoulou B. Apoptosis-like programmed cell death induces antisense ribosomal RNA (rRNA) fragmentation and rRNA degradation in Leishmania. Cell Death Differ. 2012;19(12): 1972–82.

  40. 40.

    Kramer S, Bannerman-Chukualim B, Ellis L, Boulden EA, Kelly S, Field MC, et al. Differential localization of the two T. brucei poly(a) binding proteins to the nucleus and RNP granules suggests binding to distinct mRNA pools. PLoS One. 2013;8:e54004. https://doi.org/10.1371/journal.pone.0054004.

  41. 41.

    Zoltner M, Krienitz N, Field MC, Kramer S. Comparative proteomics of the two T. brucei PABPs suggests that PABP2 controls bulk mRNA. PLoS Negl Trop Dis. 2018;12:e0006679.

  42. 42.

    Alves LR, Oliveira C, Goldenberg S. Eukaryotic translation elongation factor-1 alpha is associated with a specific subset of mRNAs in Trypanosoma cruzi. BMC Microbiol. 2015;15:104.

  43. 43.

    Murray JW, Edmonds BT, Liu G, Condeelis J. Bundling of actin filaments by elongation factor 1 alpha inhibits polymerization at filament ends. J Cell Biol. 1996;135:1309–21 http://www.ncbi.nlm.nih.gov/pubmed/8947553.

  44. 44.

    Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.

  45. 45.

    Ward N, Moreno-Hagelsieb G. Quickly finding Orthologs as reciprocal best hits with BLAT, LAST, and UBLAST: how much do we miss? PLoS One. 2014;9:e101850.

  46. 46.

    Blom N, Gammeltoft S, Brunak S. Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. J Mol Biol. 1999;294:1351–62.

  47. 47.

    Buchan DWA, Minneci F, Nugent TCO, Bryson K, Jones DT. Scalable web services for the PSIPRED protein analysis workbench. Nucleic Acids Res. 2013;41(Web Server issue):W349–57.

  48. 48.

    Contreras VT, Salles JM, Thomas N, Morel CM, Goldenberg S. In vitro differentiation of Trypanosoma cruzi under chemically defined conditions. Mol Biochem Parasitol. 1985;16:315–27.

  49. 49.

    Batista M, Marchini FK, Celedon PA, Fragoso SP, Probst CM, Preti H, et al. A high-throughput cloning system for reverse genetics in Trypanosoma cruzi. BMC Microbiol. 2010;10:259. https://doi.org/10.1186/1471-2180-10-259.

  50. 50.

    Lu HY, Buck GA. Expression of an exogenous gene in Trypanosoma cruzi epimastigotes. Mol Biochem Parasitol. 1991;44:109–14.

  51. 51.

    Inoue AH, Serpeloni M, Hiraiwa PM, Yamada-Ogatta SF, Muniz JRC, Motta MCM, et al. Identification of a novel nucleocytoplasmic shuttling RNA helicase of trypanosomes. PLoS One. 2014;9:e109521.

  52. 52.

    Towbin H, Staehelin T, Gordon J. Electrophoretic transfer of proteins from polyacrylamide gels to nitrocellulose sheets: procedure and some applications. Proc Natl Acad Sci U S A. 1979;76:4350–4.

  53. 53.

    Holetz FB, Correa A, Avila AR, Nakamura CV, Krieger MA, Goldenberg S. Evidence of P-body-like structures in Trypanosoma cruzi. Biochem Biophys Res Commun. 2007;356:1062–7. https://doi.org/10.1016/j.bbrc.2007.03.104.

  54. 54.

    Chassé H, Boulben S, Costache V, Cormier P, Morales J. Analysis of translation using polysome profiling. Nucleic Acids Res. 2017;45(3):e15.

  55. 55.

    Chevallet M, Luche S, Rabilloud T. Silver staining of proteins in polyacrylamide gels. Nat Protoc. 2006;1:1852–8.

  56. 56.

    Carvalho PC, Lima DB, Leprevost FV, Santos MDM, Fischer JSG, Aquino PF, et al. Integrated analysis of shotgun proteomic data with PatternLab for proteomics 4.0. Nat Protoc. 2015;11:102–17.

  57. 57.

    Eng JK, Jahan TA, Hoopmann MR. Comet: an open-source MS/MS sequence database search tool. Proteomics. 2013;13:22–4.

  58. 58.

    Carvalho PC, Fischer JSG, Xu T, Cociorva D, Balbuena TS, Valente RH, et al. Search engine processor: filtering and organizing peptide spectrum matches. Proteomics. 2012;12:944–9.

  59. 59.

    Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57:289–300.

  60. 60.

    Carvalho PC, Yates JR, Barbosa VC. Improving the TFold test for differential shotgun proteomics. Bioinformatics. 2012;28:1652–4.

  61. 61.

    Baggerly KA, Deng L, Morris JS, Aldaz CM. Differential expression in SAGE: accounting for normal between-library variation. Bioinformatics. 2003;19:1477–83 http://www.ncbi.nlm.nih.gov/pubmed/12912827.

Download references

Acknowledgments

We would like to thank Nilson Fidêncio and Vanessa M. dos Santos for their technical assistance and the Program for Technological Development in Tools for Health-PDTIS-FIOCRUZ for the use of its facilities (RPT02H, mass spectrometry facility; and RPT07C, microscopy facility at the Carlos Chagas Institute/Fiocruz-PR, Brazil).

Funding

This work was supported by the CNPq Conselho Nacional de Desenvolvimento Científico e Tecnológico and Fiocruz-PAPES. PCC and SG are research fellows from CNPq.

Author information

JSM: performed experiments and analyzed the data; HHW: performed experiments, analyzed the data, and prepared the manuscript; AHI: performed experiments, analyzed the data, and critically discussed the results; FVL: analyzed the data and discussed the results; PCC: helped in the proteomic data analysis; SG: critically discussed the results; and LRA: analyzed the RNA-seq data, prepared the manuscript and critically discussed the results. All authors read and approved the final manuscript.

Correspondence to Lysangela Ronalte Alves.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Additional files

Additional file 1:

Figure S1. DRBD2 is an ortholog of yeast Gbp2. Protein alignment based on NCBI Constraint-based Multiple Alignment Tool (COBALT); this tool performs multiple protein sequence alignment using conserved domain and local sequence similarity information. Red bars indicate the regions of proteins sequences similarities. (TIFF 781 kb)

Additional file 2:

Figure S2. Phosphorylation sites predicted by NetPhos 2.0 server. DRBD2 sequence provided by TriTrypDB with 276 amino acids. DRBD2 amino acid sequence and DRBD2 predicted phosphorylation sites. Ser: serine; Thr: threonine; Tyr: tyrosine. (TIF 54 kb)

Additional file 3:

Figure S3. Expression levels of drbd2 transcript and the targets of DRBD2-mRNP. (A) drbd2 expression level in the developmental forms (data from Li et al., 2016 [33] available at tritrypdb). (B) drbd2 expression level in the epimastigotes and trypomastigote forms - RNA-seq (total RNA) and ribosome profiling sequencing (Ribo-seq) (data from Smircich et al., 2015 [34] aviailable at tritrypdb). (C) Expression levels of the most abundant transcripts associated to DRBD2-mRNP complex in the epimastigotes and trypomastigote forms - RNA-seq (Total RNA) and ribosome profiling sequencing (Ribo-seq) (data from Smircich et al., 2015 aviailable at tritrypdb). In the X-axis the expression levels in FPKM (Fragments Per Kilobase Million). (TIFF 385 kb)

Additional file 4:

Table S1. Complete information of the shotgun proteomic analysis of the DRBD2-mRNP complex. Unique proteins: proteins identified only in DRBD2 IP assays. Signal: sum of PatternLab’s normalized label-free quantitation derived from the extracted ion chromatograms – XIC; Protein description and ID from TriTrypDB: TriTryp database description and identification number for a given protein. Product: the protein given name from TriTrypDB. Enriched proteins: proteins identified in DRBD2 IP and control IP replicates with differential abundance. Proteins with fold change ≥2 were considered. Signal DRBD2 and Signal control: signal intensity of the proteins in DRBD2 IP and control IP replicates, respectively. (XLSX 38 kb)

Additional file 5:

Table S2. Transcripts identified by the ribonomic analysis of the DRBD2-mRNP complex. Transcripts identified by RNA-seq from the DRBD2-mRNP complex from epimastigotes in exponential growth. The table also shows the transcripts regulated by DRBD2 that are associated with polysomes in epimastigotes and metacyclic trypomastigotes, information provided by Smircich et al. (XLSX 224 kb)

Additional file 6:

Table S3. DRBD2-mRNP complex RNA-seq mapping statistics information. (XLSX 5 kb)

Additional file 7:

Table S4. Complete information of the ribonomic analysis of the DRBD2-mRNP complex. (XLSX 4213 kb)

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • RNA-binding protein
  • Gene expression regulation
  • RNA metabolism
  • Trypanosoma cruzi