TbRAP1 interacts with the active VSG RNA in vivo
RAP1 homologs have been identified from kinetoplastids to mammals21. None of the known RAP1 homologs has been reported to have any RNA binding activity. We previously found that TbRAP1 does not bind the telomeric repeat-containing RNA (TERRA) in RNA IP experiments22,23,24. However, to our great surprise, we found that TbRAP1 interacts with the active VSG RNA (Fig. 1a, b). RNA crosslinking immunoprecipitation (RNA CLIP) assays were performed in TbRAP1F2H+/- cells that express VSG2 as the major surface antigen (Table 1), in which one TbRAP1 allele is deleted and the other has an N-terminal FLAG-HA-HA (F2H) tag17. Quantitative RT-PCR (qRT-PCR) analysis detected significantly more active VSG2 RNA in the TbRAP1 CLIP product than in the negative control IgG CLIP product (Fig. 1a). RNAs of the telomerase reverse transcriptase (TbTERT25), small nuclear RNA gene activation protein 50 (SNAP50), and Protein Kinase A catalytic subunit (PKAC1) were also examined in CLIP products. Approximately the same amount of TbTERT, SNAP50, and PKAC1 RNAs were detected in both TbRAP1 and IgG CLIP products (Fig. 1a). Therefore, TbRAP1 interacts with the active VSG RNA but not TbTERT, SNAP50, or PKAC1 RNAs. We also performed RNA CLIP in PVS3-2/OD1-1 cells (Table 1) that express VSG9 as the major surface antigen15. Again, qRT-PCR detected significantly more VSG9 RNA in the TbRAP1 CLIP product than in the IgG control (Fig. 1b), indicating that TbRAP1 can interact with the active VSG RNA regardless of which VSG is expressed. Therefore, we report for the first time that TbRAP1 is associated with the active VSG RNA in vivo, an unprecedented finding for RAP1 homologs.
a RNA CLIP experiments were performed in TbRAP1F2H+/- cells that express VSG2. qRT-PCR was performed to estimate the amount of the VSG2 RNA and the TbTERT, SNAP50, and PKAC1 RNAs in the RNA CLIP product. Enrichment of the VSG2, TbTERT, SNAP50, and PKAC1 RNAs (CLIP/Input) was calculated for the CLIP experiment using the HA antibody 12CA5 and that using IgG. Relative enrichment was calculated using the enrichment of IgG CLIP as a reference. Average and standard deviation were calculated from three (SNAP50 & PKAC1), five (TbTERT), and seventeen (VSG2) independent experiments. P values of two-sided unpaired t-tests (compared to VSG2 RNA enrichment) are shown. b RNA CLIP was performed in VSG9-expressing PVS3-2/OD1-1 cells using a TbRAP1 rabbit antibody15 and IgG and the enrichment of the VSG9 RNA in the CLIP product was calculated. Average and standard deviation were calculated from three independent experiments. Error bars represent standard deviation. Source data are provided as a Source Data file. c Domain structure of TbRAP1. Inset, an enlarged diagram of the TbRAP1 MybLike domain (aa 639–761)15, which contains an RRM (aa 653–727) and the DNA Binding (DB) domain (aa 734–761)18. Arrowheads mark the conserved F655 and F694 residues. d Superposition of TbRAP1653-727 (green) with the RRM1 domain of hnRNP A1 (orange) bound with a short RNA oligo (golden) [doi.org/10.2210/pdb5MPG/pdb]28. Inset highlights that F655 and F694 in TbRAP1 superimpose well with F17 and F59 in hnRNP A1 that form stacking interactions with the RNA substrate.
The TbRAP1 MybLike domain contains an RRM module
To investigate whether the TbRAP1 MybLike domain (aa 639–761) is responsible for binding to the active VSG RNA, we first determined the solution structure of TbRAP1639-761 by NMR spectroscopy (Fig. 1c, Supplementary Fig. 1a, Supplementary Table 1). The N-terminal region of TbRAP1639-761 does not adopt a typical Myb fold but forms a canonical RRM (aa 653–727)26 with the characteristic topology of a four-stranded anti-parallel β-sheet and two α-helices packed behind the β-sheet (Fig. 1d; Supplementary Fig. 1b). The DB domain (aa 734–761) at the C-terminus of TbRAP1639-761 forms a long and flexible loop (Supplementary Fig. 1b, left). In contrast, none of the known RAP1 homologs has been reported to have an RRM domain.
RRM is a conserved structural platform that binds to diverse RNAs and ssDNAs26,27. Sequence analysis shows that TbRAP1 RRM contains the signature RNP1 and RNP2 sequence motifs, with F655 in RNP2 and F694 in RNP1 representing the two conserved aromatic residues critical for substrate binding (Supplementary Fig. 1c)26. TbRAP1 RRM superimposes well with RRM1 of heterogeneous nuclear ribonucleoprotein (hnRNP) A1 bound with an RNA oligo [doi.org/10.2210/pdb5MPG/pdb]28 (Fig. 1d), with a Root Mean Square Deviation (RMSD) of ~3.3-3.5 Å for the main chain atoms. Notably, F655 and F694 of TbRAP1 match exactly to F17 and F59 of hnRNP A1 that form stacking interactions with RNA (Fig. 1d). In addition, sequence alignment and structural prediction by AlphaFold229 confirm that RAP1 homologs in representative Trypanosomatida organisms all have a highly conserved RRM (Supplementary Fig. 1d), while vertebrate and fungal RAP1s do not seem to have any RRM domain (Supplementary Fig. 1e). Thus, the RRM domain is uniquely conserved in RAP1 homologs of these microbial parasites but absent in RAP1s from higher eukaryotes.
TbRAP1 RRM binds to the consensus VSG 3’UTR region in vitro
Since TbRAP1639-761 contains an RRM domain plus a flexible DB domain, we then used NMR titration to test whether it binds to the active VSG RNA. We used 34-VSG-UTR, a 34 nt RNA from the 3’UTR of VSG2 that contains the consensus 16-mer found in all VSG 3’UTRs (Supplementary Table 2)19,20. We titrated 34-VSG-UTR into 15N-labeled TbRAP1639-761 (Supplementary Table 3) and observed significant concentration-dependent chemical shifts for RNP1 and RNP2 residues, particularly F655 and F694, in heteronuclear single quantum correlation (HSQC) spectra (Fig. 2a–c). A few residues in the DB domain also showed noticeable chemical shifts, although at much lower magnitudes compared to RNP1 and RNP2 residues (Fig. 2a, b). These results suggest that both the RRM module and the DB domain interact with the 34-VSG-UTR, with RRM playing a major role. Compared to other RRM domains, the chemical shifts induced by 34-VSG-UTR in TbRAP1639-761 are mostly in the slow-to-intermediate exchange region, indicative of a moderate micromolar binding affinity28.
a, d, g 1H-15N HSQC NMR spectra of 15N-labeled TbRAP1639–761 (a), TbRAP1639–733 (d) and TbRAP1639–7332FL (g) in the absence (black) and presence of 34-VSG-UTR in 3× molar excess (red). In (a), residues located in the RRM domain (labeled in green) showed noticeable chemical shifts (arrow) while residues in the DB domain (labeled in blue) did not (underline). In (g), no chemical shifts were observed. b, e Chemical shift differences of individual TbRAP1 residues in NMR titration when TbRAP1639-761 (b) or TbRAP1639-733 (e) was used. Source data are provided as a Source Data file. c Inset of overlaid 1H-15N-HSQC spectra in (a) highlighting chemical shift perturbations for key residues in RRM in the absence (black) and presence of 34-VSG-UTR in 1x (blue), 2x (green) and 3× (red) molar excess. Residues located on RNP1 and RNP2 of RRM, including the conserved F655 and F694 are highlighted in (c). f Insets of overlaid 1H-15N-HSQC spectra of 15N-labeled TbRAP1639–733 in the absence (black) or presence of 34-VSG-UTR (top), (UUAGGG)2 (middle), and 35-random (bottom) in 1x (blue), 2x (green) and 3× (red) molar excess. Highlighted residues are the same as in (c). Only 34-VSG-UTR induced noticeable chemical shifts in the RRM domain. PPM, parts per million.
To further characterize how the RRM and DB domains bind RNA, we did similar NMR titration studies using the RRM-containing TbRAP1639-733, TbRAP1639-7332FL with the two key aromatic residues F655 and F694 of the RRM domain mutated to leucine residues, and TbRAP1639-7615A with the R/K patch in the DB domain mutated to five alanines18 (Supplementary Table 3). For both TbRAP1639-733 and TbRAP1639-761, 34-VSG-UTR induced similar patterns of chemical shift in RNP1 and RNP2 (Fig. 2, a, b, d, e), but the magnitude was smaller for TbRAP1639-733 than for TbRAP1639-761 (Fig. 2, b, c, e, f). NMR titration using TbRAP1639-7615A also showed similar results as TbRAP1639-733 (Supplementary Fig. 1f). However, no chemical shifts were observed for TbRAP1639-7332FL even when 34-VSG-UTR was in 3-fold molar excess (Fig. 2g). These results indicate that RRM alone can bind 34-VSG-UTR, which requires the two conserved aromatic residues F655 and F694, while the DB domain helps to strengthen this binding.
To explore the sequence specificity of TbRAP1 RRM, we tested TbRAP1639-733’s binding to (UUAGGG)2, an oligo that contains the TERRA sequence22,23,24. (UUAGGG)2 did not induce any noticeable chemical shifts when titrated to TbRAP1639-733 (Fig. 2f; Supplementary Fig. 2a), which is consistent with our previous observation that TbRAP1 does not bind TERRA24. We also tested TbRAP1639-733’s binding to 35-random, a 35 nt RNA with a random sequence (Supplementary Table 2) by NMR titration. No noticeable chemical shifts were detected, either (Fig. 2f; Supplementary Fig. 2b). These data suggest that TbRAP1639-733 does bind RNA with certain sequence specificity.
RRM domains are known to recognize short RNA motifs of 2-8 nucleotides26,27. To further map which sequence within 34-VSG-UTR can be recognized by TbRAP1 RRM, we performed NMR titration using 16-VSG-UTR, an oligo that contains only the 16-mer consensus sequence in VSG 3’UTRs (Supplementary Table 2). 16-VSG-UTR and 34-VSG-UTR induced the same pattern of chemical shifts in both TbRAP1639-761 and TbRAP1639-733 (Fig. 2a, b, d, e; Supplementary Fig. 2c–f). Therefore, the 16-mer consensus sequence in VSG 3’UTRs is sufficient to be recognized by TbRAP1 RRM. In addition, 16-VSG-UTR also induced stronger chemical shifts in TbRAP1639-761 than TbRAP1639-733, further validating the supporting role of the DB domain (Supplementary Fig. 2c–f). Furthermore, the magnitude of chemical shifts induced by 16-VSG-UTR for the aromatic residues F655 and F694 in RRM was ~50% lower than those induced by 34-VSG-UTR (Supplementary Fig. 2d, f vs. Fig. 2b, e, respectively). These subtle differences suggest that TbRAP1 RRM may recognize additional sequence motifs in the longer 34-VSG-UTR substrate, which leads to stronger binding and more prominent chemical shifts. Since RRM domains are known to have promiscuous binding activities, it is likely that TbRAP1 RRM can recognize more than one sequence within the VSG RNA.
We also used the fluorescence polarization assay as a biophysical technique to assess the RNA binding activity of TbRAP1. Fluorophore-labeled 16-VSG-UTR was titrated to TbRAP1639-761, TbRAP1639-733, and TbRAP1639-7615A and the estimated binding affinity Kd were ~258, 929, and 969 μM, respectively (Supplementary Fig. 2g–i). These data corroborate our NMR studies to confirm that TbRAP1 RRM recognizes the 16-mer consensus sequence of VSG 3’UTRs. This RNA binding activity requires the two conserved aromatic residues, F655 and F694 in RNP2 and RNP1, respectively, and is enhanced by the DB domain.
We subsequently performed EMSA to validate the TbRAP1 RRM-mediated RNA binding activity. Initially, TrxA-His6 (TH6) or GST-tagged TbRAP1 fragments were used (Supplementary Fig. 3a). TH6-tagged TbRAP1639-761, TbRAP1639-733, and TbRAP1639-7615A (Supplementary Table 3) all bound 170-VSG-UTR, a 170 nt RNA containing the VSG2 3’UTR sequence (Supplementary Table 2), while TH6 alone or TH6–TbRAP1639-7612FA&5A (F655AF694A, 737RKRRR741 mutated to 737AAAAA741, Supplementary Table 3) did not (Supplementary Fig. 3b–d). In addition, GST-TbRAP1414-855 bound this RNA, while GST alone and the GST-tagged duplex telomere DNA-binding TbTRF30 did not (Supplementary Fig. 3b, c, e; Supplementary Table 3).
To examine TbRAP1-specific RNA binding activity without any possible interference by the fusion tag, we cleaved the TH6 tag by 3C and purified tagless TbRAP1 fragments (Supplementary Fig. 3f). Both TbRAP1639-761 and TbRAP1639-733 bound 170-VSG-UTR (Fig. 3a, b) but TbRAP1639-7332FQ (F655QF694Q), TbRAP1639-7332FL (F655LF694L), or TbRAP1639-7332FA (F655AF694A) did not (Fig. 3c, d; Supplementary Fig. 3g; Supplementary Table 3). In addition, more than one TbRAP1639-733 molecule can bind the same 170-VSG-UTR substrate when the protein:RNA ratio is increased (Fig. 3b).
Untagged recombinant TbRAP1639-761 (a, e, j), TbRAP1639-733 (b–d, f–i), TbRAP1639-7332FQ (c, g), and TbRAP1639-7332FL (d, h) were incubated with 170-VSG-UTR (a–d), 170-no-VSG (e–h), 35-VSG-UTR (i), 35-random (i), or 16-VSG-UTR (j) (Supplementary Table 2). The concentration of protein (µM) used in each reaction is indicated on top of each lane. Samples were electrophoresed in 0.8% agarose gels (a–i) or a 1.2% agarose gel (j) in 0.5 x TBE buffer. Source data are provided as a Source Data file.
Unexpectedly, tagless TbRAP1639-761 and TbRAP1639-733 bound 170-no-VSG (Fig. 3e, f; Supplementary Table 2) but none of TbRAP1639-7332FQ, TbRAP1639-7332FL, or TbRAP1639-7332FA did (Fig. 3g, h; Supplementary Fig. 3g). Similarly, TH6-tagged TbRAP1639-761, TbRAP1639-733, and TbRAP1639-7615A also bound 170-no-VSG (Supplementary Fig. 3h, i). TbRAP1639-733 did exhibit higher affinity to 170-VSG-UTR than to 170-no-VSG (Supplementary Fig. 3j), indicating that TbRAP1 RRM prefers the VSG 3’UTR sequence. Nevertheless, the observation that TbRAP1639-733 bound 170-no-VSG (Fig. 3f) seems inconsistent with the fact that TbRAP1639-733 does not bind 35-random in NMR titration (Fig. 2f; Supplementary Fig. 2b). We, therefore, examined whether TbRAP1639-733 binds 35-random in EMSA. 35-VSG-UTR was used as a positive control, which contains both the 9-mer and the 16-mer consensus motifs in VSG 3’UTR (Supplementary Table 2)20. TbRAP1639-733 bound 35-VSG-UTR but not 35-random in EMSA (Fig. 3i), confirming the NMR titration result. RRM domains usually recognize a short RNA sequence of 2-8 nucleotides26,27. It is possible that 170-no-VSG may contain additional sequences that can be recognized by TbRAP1 RRM other than the consensus sequences in VSG 3’UTRs.
We further performed EMSA using the shorter 16-VSG-UTR substrate (Supplementary Table 2), to better explore the sequence specificity of TbRAP1’s RNA binding activity. Interestingly, TbRAP1639-761 clearly bound 16-VSG-UTR (Fig. 3j) but TbRAP1639-733’s binding affinity appears to be too weak to be detected by EMSA. This observation supports our NMR titration results and further validates the importance of the DB domain in the RRM-mediated RNA binding. Additionally, Kd values estimated by EMSA show stronger affinity of TbRAP1639-761 for 35-VSG-UTR than 16-VSG-UTR, which is consistent with our NMR titration results (Supplementary Fig. 2j).
The in vivo TbRAP1-VSG RNA interaction depends on the conserved aromatic residues in RRM
We generated TbRAP1F/mut strains by replacing the WT TbRAP1 allele with various RRM mutants in TbRAP1F/+ cells (Supplementary Fig. 4a, b; Table 1)17. To specifically examine the in vivo RNA binding activities of TbRAP1 mutants, we did RNA CLIP after deleting the loxP-flanked TbRAP1 (the F allele) by Cre, as RRM mutants can interact with WT TbRAP1 through its BRCT domain17. Removal of the TbRAP1 F allele was confirmed by PCR (Supplementary Fig. 4c–g). TbRAP1∆RRM and TbRAP1∆MybL (MybLike deletion)18 were expressed at a subtly lower level than WT TbRAP1 (Supplementary Fig. 4h), while TbRAP1-2FQ, TbRAP1-2FL, TbRAP1-2FA, and TbRAP1-2FA&5A were expressed the same as TbRAP1 (Supplementary Fig. 4i–l). As expected, TbRAP1∆MybL and TbRAP1∆RRM mutants that lack the whole RRM domain lost the TbRAP1-VSG2 RNA interaction (Fig. 4a). Similarly, TbRAP1-2FQ, TbRAP1-2FA, and TbRAP1-2FA&5A did not bind VSG2 RNA, either (Fig. 4a). Interestingly, although TbRAP1-2FL bound significantly lower amount of VSG2 RNA than WT TbRAP1, this mutant appeared to have a smaller RNA binding defect than other mutants (Fig. 4a).
a RNA CLIP experiments were performed in various TbRAP1F/mut strains (expressing VSG2) after a 30-h induction of Cre. The presence of the active VSG2 RNA in the RNA CLIP product was determined by qRT-PCR. The enrichment of VSG2 RNA (CLIP/Input) was calculated for the CLIP experiment using the HA antibody 12CA5 and that using IgG. Relative enrichment was calculated using the enrichment of IgG CLIP as a reference. Average and standard deviation were calculated from two to seventeen independent experiments (the exact number of experiments was indicated in parentheses following each strain name). P values of two-sided unpaired t-tests between the TbRAP1F2H+/- and TbRAP1F/mut are shown on top of corresponding columns. Data for WT TbRAP1 is the same as that in Fig. 1a. b–d ChIP experiments using the HA antibody 12CA5, a TbTRF rabbit antibody30, and IgG were done in TbRAP1F2H+/- cells and Cre-induced (for 30 h) TbRAP1F/2FA&5A (b) TbRAP1F/2FQ (c) and TbRAP1F/2FL (d) cells. Average and standard deviation were calculated from two to five independent experiments (exact number of samples are indicated beneath bottom labels). P values of two-sided unpaired t-tests are shown (ChIP using 12CA5, TbRAP1F/mut vs TbRAP1F2H+/-). Source data are provided as a Source Data file. e IF analyses were done in TbRAP1F2H+/- (top), TbRAP1F/2FQ (middle), and TbRAP1F/2FL (bottom) cells. 12CA5 and a TbTRF chicken antibody15 were used. TbRAP1 genotypes are listed on the left. DNA was stained by DAPI. All images are of the same scale and a size bar is shown in one of the images.
Because TbRAP1 DB enhances the RNA binding activity in vitro, we further examined the effect of DB domain mutations on VSG RNA binding in vivo. We previously reported that TbRAP1∆DB and TbRAP1-5A were expressed at the same level as TbRAP118. Surprisingly, both TbRAP1∆DB and TbRAP1-5A only pulled down background level of VSG2 RNA (Fig. 4a). Since neither mutant is associated with the telomere chromatin18, this observation suggests that being localized at the telomere is a prerequisite for TbRAP1 to bind the active VSG RNA, which has a high concentration only at the active VSG locus.
We also performed Chromatin IP (ChIP) to test whether the RRM domain is necessary for TbRAP1’s localization to the telomere. TbRAP1-2FA&5A did not associate with the telomere chromatin (Fig. 4b), presumably because the 5A mutation already abolished TbRAP1’s DNA binding activities18. In contrast, TbRAP1-2FQ, TbRAP1-2FL, and TbRAP1-2FA still associated with the telomere chromatin (Fig. 4c, d; Supplementary Fig. 4m). Immunofluorescence (IF) analysis further showed that both TbRAP1-2FQ and 2FL were partially colocalized with TbTRF that binds the duplex telomere DNA30 the same way as WT TbRAP1 (Fig. 4e). Hence, in vivo binding of TbRAP1 to the active VSG RNA depends on RRM’s two conserved residues F655 and F694 and the R/K patch within the DB domain. Additionally, the RRM-mediated RNA binding activity is not required for TbRAP1’s association to the telomere chromatin.
TbRAP1’s RNA binding activity is important for VSG MAE and telomere integrity
We examined phenotypes of TbRAP1F/∆RRM, TbRAP1F/2FQ, TbRAP1F/2FL, TbRAP1F/2FA, and TbRAP1F/2FA&5A after a 30–48 h Cre induction (Supplementary Fig. 5a–e). In TbRAP1F/∆RRM, TbRAP1F/2FQ, TbRAP1F/2FA, and TbRAP1F/2FA&5A cells, Cre induction led to an acute growth arrest (Supplementary Fig. 5f–i). However, TbRAP1F/2FL cells showed a slower but not arrested growth phenotype upon Cre induction (Supplementary Fig. 5j), which is consistent with the observation that the 2FL mutant affects the RNA binding less than 2FQ, 2FA, and ∆RRM (Fig. 4a). In addition, substituting an aromatic ring in the phenylalanine residue with a long hydrophobic chain in the leucine residue likely has a weaker effect than substituting it with an alanine.
VSG MAE has two essential aspects: silencing all but one VSGs and a full-level expression of the active VSG. In TbRAP1∆RRM, TbRAP12FQ, and TbRAP12FL mutants, qRT-PCR analysis after the 30-48 h Cre induction detected a significant decrease (~40-60%) in the active VSG2 RNA level, while RNA levels of silent VSGs increased several hundred-fold (Supplementary Fig. 5k–m), indicating that TbRAP1 RRM is essential for both aspects of VSG MAE. The decrease in VSG2 level is particularly striking because the active VSG RNA is ~10,000 fold more abundant than any silent VSG RNA (Fig. 5a)16. Thus, ~50% reduction of the active VSG2 RNA represents a more overwhelming change than the several hundred-fold increase in RNA levels of silent VSGs. This decrease is also in distinct contrast to the phenotype of TbRAP1F/5A and TbRAP1F/∆DB cells that mutated the R/K patch, where silent VSGs were similarly derepressed but the active VSG RNA remained at ~90% of the WT level (Supplementary Fig. 5n, o)18. Interestingly, in TbRAP12FA&5A mutant, the active VSG RNA level was also only decreased to ~87% of the WT level (Supplementary Fig. 5p). Hence, TbRAP1’s RNA binding activity is particularly essential for keeping the active VSG fully transcribed, while mutating the R/K patch leads to a global VSG derepression and renders the TbRAP1-VSG RNA interaction unimportant.
a A diagram illustrating the ~10,000-fold difference between the active VSG RNA amount and any silent VSG RNA amount. Spheres are not drawn to scale. b–e, g qRT-PCR of RNA levels of the active VSG2 (indicated in red), several silent ES-linked VSGs, and chromosome internal TbTERT and tubulin in TbRAP1F/2FQ (b), TbRAP1F/2FL (c), TbRAP1F/2FA (d), TbRAP1F/5A (e), and TbRAP1F/2FA&5A (g) cells. The fold changes in RNA level are shown in the log scale. Average and standard deviation were calculated from two to nine independent experiments (exact number of samples are indicated beneath each column). The change in VSG2 RNA level in these mutants is plotted again in the linear scale in (f). At the 12 h point, derepression of VSG3, 6, and 9 in TbRAP1F/2FQ, TbRAP1F/2FL, TbRAP1F/2FA, and TbRAP1F/2FA&5A cells was compared to that in TbRAP1F/5A by two-sided unpaired student t-tests, and p values of significant differences are indicated on top of corresponding columns in (b–d, g). The changes in the VSG2 RNA level at all time points were compared to that in TbRAP1F/5A cells in the same way. P values of significant differences are indicated on top of corresponding columns in (f). Error bars in (b–g) represent standard deviation. Source data are provided as a Source Data file. h IF analysis of TbRAP1F/2FL cells before and after the Cre induction. Antibodies specifically recognizing VSG6 (green) and VSG3 (red), which were silent in WT cells, were used. DAPI was used to stain DNA. All panels are of the same scale, and a size bar is shown in one of the panels.
We further examined the RNA levels of the active VSG2 at early time points of 12–36 h after Cre induction in TbRAP1F/mut cells, aiming to assess direct effects of TbRAP1 RRM mutations on VSG expression. Western analysis confirmed the decrease of the total TbRAP1 level in these cells (Supplementary Fig. 6a–e). Strikingly, the active VSG2 RNA level showed significant drop by 12 h and continued to decrease over time, dropping to 58%, 68%, and 50% of the WT level by 24 h in TbRAP1F/2FQ, TbRAP1F/2FL, and TbRAP1F/2FA cells, respectively (Fig. 5b–d, f). In contrast, the VSG2 RNA level remained close to the WT level in TbRAP1F/5A (~90% by 30 h after Cre induction) and TbRAP1F/2FA&5A cells (~87% by 36 h after Cre induction) (Fig. 5e–g). Our temporal profiling of the VSG2 RNA level further confirms that TbRAP1 RRM is critical for sustaining full-level expression of the active VSG.
We also examined the RNA levels of several silent VSGs at the time points of 12–36 h after Cre induction in TbRAP1F/mut cells. Notably, derepression of silent VSG 3, 6, and 9 at 12 h after the Cre induction in RRM point mutants was only ~10 fold, significantly milder than the ~100 fold observed in TbRAP1F/5A cells (Fig. 5b–e). The magnitude of depression became similar at later time points of 18, 24, and 36 h (Fig. 5b–e). Nevertheless, both VSG3 and VSG6, two silent VSGs in uninduced TbRAP1F/2FL cells, were expressed simultaneously in individual cells upon Cre induction (Fig. 5h). Overall, these results confirm that disrupting TbRAP1’s RNA binding indeed led to VSG derepression, albeit with a slower kinetic profile compared to mutations in the DB domain.
Subsequently, we examined the transcriptome profiles in TbRAP1F/2FQ and TbRAP1F/2FA&5A cells by RNAseq. ~5,000 genes were up-regulated and 200-1500 genes were down-regulated in the TbRAP1F/2FQ and TbRAP1F/2FA&5A cells (Supplementary Fig. 6f, g). A large number of VSG genes and pseudogenes were up-regulated in both mutants, including all silent VSGs in bloodstream form VSG ESs (Supplementary Figs. 7a, b and 8). GO term analysis indicated that significantly derepressed genes are predominantly involved in host immune response evasion (Supplementary Fig. 7c). We also estimated the VSG2 RNA half-life in TbRAP1F/∆MybL, TbRAP1F/∆RRM, and TbRAP1F/2FQ cells. The VSG2 RNA levels were examined by qRT-PCR after various lengths of time of Actinomycin D treatment, but the half-life of VSG2 RNA did not change in RRM mutants (Supplementary Fig. 9).
We previously showed that TbRAP1 suppresses VSG switching by maintaining genome integrity at the telomere and subtelomere23. Since TbRAP1-2FL is viable, we estimated the VSG switching rate in TbRAP1-/2FL cells (Table 1), which is twice as high as that in WT cells (Fig. 6a), suggesting that the TbRAP1’s RNA binding activity also helps suppress VSG switching. In addition, the level of γH2A, an indicator of DNA damage31, was increased mildly (Fig. 6b), and significantly more γH2A was associated with the telomere chromatin in Cre-induced TbRAP1F/2FL cells (Fig. 6c). TbRAP1∆RRM, TbRAP1-2FQ, TbRAP1-2FA, and TbRAP1-2FA&5A mutants exhibited a strong growth arrest phenotype (Supplementary Fig. 5f–i), which prevented us from determining the VSG switching rate in these mutants. Therefore, we examined the γH2A levels. An increased level of γH2A was observed in Cre-induced TbRAP1F/2FQ, TbRAP1F/2FA&5A, TbRAP1F/2FA, and TbRAP1F/∆RRM cells (Fig. 6d, e; Supplementary Fig. 6h–i), indicating that these mutants had more DNA damage. Particularly, we observed an increased amount of γH2A associated with the telomere chromatin (Fig. 6f) and the active ES (Fig. 6g) in TbRAP1F/2FQ cells after the Cre induction, indicating that the TbRAP1’s RNA binding activity is also critical for telomeric and subtelomeric integrity. Telomeric DNA breaks, particularly those at the active VSG vicinity, can lead to cell death in >80% of parasites32, which can explain why RRM mutants have growth defects.
a TbRAP1-/2FL exhibits an increased VSG switching rate. Average and standard deviation were calculated from three (WT) and four (TbRAP1-/2FL) independent experiments. P values of two-sided unpaired t-tests are shown (TbRAP1-/2FL vs TbRAP1+/+). b, d, e Western analyses to examine the γH2A protein level in WT cells before and after phleomycin treatment (as a positive control) and in TbRAP1F/2FL (b), TbRAP1F/2FQ (d), and TbRAP1F/2FA&5A (e) cells before and after a 30–48 h Cre induction. A γH2A rabbit antibody23 and the tubulin antibody TAT-149 were used. Molecular marker was run on the left lane in each gel and their sizes are indicated on the left. c, f ChIP using the γH2A rabbit antibody and IgG in TbRAP1F/2FL (c) and TbRAP1F/2FQ (f) cells after a 30 h Cre induction followed by Southern blotting using a telomere and a tubulin probe. Blots were exposed to a phosphorimager. Images were quantified using ImageQuant and average and standard deviation were calculated from two (γH2A antibody, (TTAGGG)n probe in TbRAP1F/2FL cells) or three (all other samples) independent experiments in (c) and three independent experiments in (f). P values of two-sided unpaired t-tests (mutant vs. control cells) are shown. g ChIP using a γH2A rabbit antibody and IgG in TbRAP1F/2FQ cells followed by quantitative PCR using primers specific to the indicated active and silent ES loci. SNAP50 is a chromosome internal gene. Average enrichment (ChIP/Input) was calculated from three independent experiments. P values of two-sided unpaired t-tests (γH2A ChIP products, +Cre vs. -Cre) are shown. Error bars in (a, c, f, g) represent standard deviation. Source data are provided as a Source Data file.
TbRAP1 binds DNA and RNA in a mutually exclusive manner
Our NMR structure of TbRAP1639-761 shows that the DB domain forms a long and flexible loop that does not contact the RRM module (Supplementary Fig. 1b). Thus, it is theoretically possible for TbRAP1 to bind DNA and RNA simultaneously. To test this possibility, we conducted EMSA assays using both dsDNA and RNA substrates. We first confirmed that TbRAP1639-761 bound a duplex telomeric DNA probe, 100-ds(TTAGGG) (Fig. 7a; Supplementary Table 2)18. However, when non-radiolabeled 170-VSG-UTR and radiolabeled 100-ds(TTAGGG) were both incubated with TbRAP1639-761, no ternary complex of TbRAP1-RNA-DNA was observed (Fig. 7b). Instead, the amount of TbRAP1-DNA complex gradually decreased in the presence of an increasing amount of 170-VSG-UTR (Fig. 7b). Similarly, when non-radiolabeled 100-ds(TTAGGG) and radiolabeled 170-VSG-UTR were both incubated with TbRAP1639-761, no ternary complex was observed while the amount of TbRAP1-RNA gradually decreased with increasing amount of 100-ds(TTAGGG) (Fig. 7c). EMSA estimated that the Kd values for binding either 100-ds(TTAGGG) or 170-VSG-UTR by TbRAP1639-761 are comparable in the range of ~100–300 nM (Fig. 7d), thus allowing two-way competition. To investigate whether such competition applies to shorter DNA or RNA substrates, we further compared TbRAP1639-761 binding on 80-dsDNA and 81-VSG-UTR (Supplementary Table 2), as the shortest ssDNA and dsDNA that TbRAP1 can bind is ~60 nt and 60 bp, respectively18. TbRAP1639-761 bound 80-dsDNA (Fig. 7e) as expected18. When an increasing amount of 81-VSG-UTR was added to the reaction using radiolabeled 80-dsDNA as the substrate, no ternary complex of TbRAP1-RNA-DNA was observed, but 81-VSG-UTR competed away TbRAP1639-761’s binding on 80-dsDNA (Fig. 7f). Therefore, DNA and RNA bind to TbRAP1 in mutually exclusive and competitive manner due to their overlapping binding site and comparable binding affinities.
EMSA experiments were performed using TbRAP1639–761. Radiolabeled 100-ds(TTAGGG) (a, b), 170-VSG-UTR (c), and 80-dsDNA (e, f) were used as the binding substrates. Non-radiolabeled 170-VSG-UTR (b), 100-ds(TTAGGG) (c), and 81-VSG-UTR (f) were used as competitors. The concentration of proteins (µM) used in each experiment is indicated on top of each lane in (a) and (e). 4.7 µM (b), 2.35 µM (c), and 0.5 µM (f) of TbRAP1639–761 was used in each competition reaction. The molar excess of the competitor is indicated on top of each lane in (b, c, f). Samples were electrophoresed in 0.8% agarose gels in 0.5x TBE buffer. d TbRAP1639-761’s affinities to 100-ds(TTAGGG) and 170-VSG-UTR (Kd values) were estimated by EMSA. Average and standard deviation were calculated from four (for 100-ds(TTAGGG)) or eight (for 170-VSG-UTR) independent experiments. Source data are provided as a Source Data file.
Read more here: Source link