Kinetic cross-linking and analysis of cDNA is a method that allows investigation of the dynamics of protein-RNA interactions in living cells at high temporal resolution. Here the protocol is described in detail, including the growth of yeast cells, UV cross-linking, harvesting, protein purification, and next generation sequencing library preparation steps.
The interaction between RNA-binding proteins (RBPs) and their RNA substrates exhibits fluidity and complexity. Within its lifespan, a single RNA can be bound by many different RBPs that will regulate its production, stability, activity, and degradation. As such, much has been done to understand the dynamics that exist between these two types of molecules. A particularly important breakthrough came with the emergence of ‘cross-linking and immunoprecipitation’ (CLIP). This technique allowed stringent investigation into which RNAs are bound by a particular RBP. In short, the protein of interest is UV cross-linked to its RNA substrates in vivo, purified under highly stringent conditions, and then the RNAs covalently cross-linked to the protein are converted into cDNA libraries and sequenced. Since its conception, many derivative techniques have been developed in order to make CLIP amenable to particular fields of study. However, cross-linking using ultraviolet light is notoriously inefficient. This results in extended exposure times that make the temporal study of RBP-RNA interactions impossible. To overcome this issue, we recently designed and built much-improved UV irradiation and cell harvesting devices. Using these new tools, we developed a protocol for time-resolved analyses of RBP-RNA interactions in living cells at high temporal resolution: Kinetic CRoss-linking and Analysis of cDNAs (χCRAC). We recently used this technique to study the role of yeast RBPs in nutrient stress adaptation. This manuscript provides a detailed overview of the χCRAC method and presents recent results obtained with the Nrd1 RBP.
RNAs often rely on RBPs to exert their function, which has led to great interest in understanding the dynamics between these molecules. Many RBPs have been identified in a wide variety of organisms. However, it has always been notoriously difficult to study RBP-RNA interactions in vivo. A major breakthrough in studying such interactions came with the emergence of CLIP1. This method utilizes ultraviolet (UV, 254 nm) irradiation to induce covalent bonds between RBPs and their directly bound RNAs (i.e., zero-distance cross-linking). Subsequently, the RBP of interest is immunopurified under stringent conditions to ensure that only RNAs covalently cross-linked to the proteins are identified. Bound RNAs are then partially digested with RNases and subsequently converted into cDNA libraries for sequencing. The high purification stringency is important as it greatly increases the specificity of protein and RNA recovery, which is also further enhanced through SDS-PAGE purification of the cross-linked ribonucleoprotein (RNP) complex. CLIP and related methods also provide nucleotide resolution insight into the protein binding site, because during the preparation of the sequencing library, amino acids that cross-linked to the RNA frequently terminate the reverse transcriptase or cause the enzyme to introduce mutations at this site1,2,3.
Since its introduction, the original CLIP protocol has produced a remarkable variety of derivative methodologies. A particularly important breakthrough came with the development of HITS-CLIP (or CLIP-seq), which merges high-throughput sequencing with the CLIP approach3. This has since been adopted by all CLIP-based methodologies. iCLIP introduced improvements in the RNase-mediated trimming and adaptor ligation techniques that facilitate more accurate mapping of the RBP binding sites4. PAR-CLIP combined 4thio-uridine/uracil labeling with cross-linking at 365 nm, making it possible to map cross-linking sites by analyzing T-C substitutions5. CRAC, urea-iCLIP, dCLIP, and uvCLAP introduced denaturing conditions and dual affinity purification steps that further reduce background binding to the affinity resin and further increase the specificity of protein capture2,6,7,8,9. Additionally, CRAC, uvCLAP, and dCLIP introduced tagging the RBP of interest with an affinity tag, thus overcoming the need to generate specific antibodies.
Several optimizations have also been made to expedite the CLIP methodology. The original CLIP protocol utilized radiolabeling of the captured RNAs in order to visualize the RBP-RNA complexes after SDS-PAGE. However, the use of radioactivity can be problematic for laboratories not set up for such work. irCLIP incorporates a fluorophore-coupled adaptor that facilitates visualization through infrared imaging10 and sCLIP utilizes biotinylation of captured RNAs in order to visualize them through streptavidin-conjugated HRP11. Furthermore, eCLIP completely forgoes RNA labeling; instead the protein is excised based solely on its known size12. Streptavidin-based purification has also been used to speed up the process of library preparation in FAST-iCLIP, where a biotinylated 3’ adaptor is ligated onto the RNAs and used to enable purification after reverse transcription and circularization13. Additional enhancements to the iCLIP protocol also greatly increased the complexity of the libraries4.
Finally, CLIP has been modified to enable capture of RBPs from different cellular subcompartments14,15, to visualize newly transcribed RNAs using pulsed induction of photoactivable ribonucleosides5,16,17, to capture methylated RNAs18,19,20, to examine RNA-RNA interactions21,22, and to map 3’ ends23,24.
Despite the great contributions of CLIP-based techniques in aiding our understanding of the interactions between RBPs and RNAs, it has been limited by the inefficiency of UV cross-linking. Although culture cells grown in a monolayer are generally relatively easy to cross-link, this is significantly more challenging in tissues or cells in solution. Tissues can require multiple rounds of UV exposure in order to penetrate to the required cell layers, while microbial cells are often grown in rich mediums that contain aromatic, UV-absorbing compounds25. Indeed, UV irradiation times of up to 30 min have been used to generate sufficient cross-linking between RBPs and their bound RNAs for such samples26,27,28. This extended UV exposure induces stress responses within the cell, such UV-induced DNA damage, which can contaminate the final data in some applications.
The majority of CLIP studies have focused on generating single "snapshots" of specific protein-RNA interactions in a cell. However, protein-RNA interactions are inherently dynamic, particularly when cells are subject to changes in their environment. This can include a sudden reduction in the availability of essential nutrients or rapid changes in temperature. As such, to truly understand the role of an RBP during stress, it is best to perform time-resolved analyses because they can capture the full spectrum of RBP targets during stress and determine at what stage of the stress response the chosen RBP is active. In particular, studies in yeast showed that the first few minutes of adaptation are absolutely crucial for survival and RNA half-lives in bacteria can vary from minutes to seconds29,30,31,32,33. Therefore, such time-resolved analyses should ideally be performed at high temporal resolution. However, the long cross-linking times make the study of early stage adaptive responses particularly challenging.
In order to overcome these issues, we recently developed an improved method that is capable of cross-linking and harvesting cells on minute-long timescales. Our χCRAC method allows quantitative measurement of dynamic changes in RBP-RNA interactions at previously unwitnessed resolution. Crucial to this method was the development of a novel UV irradiation device32 that reduces the required cross-linking time in yeast and bacteria in solution around 10-fold, effectively freezing RBP-RNA interactions instantaneously. In addition, in order to rapidly harvest the cells after UV irradiation, we developed a vacuum filtration device that can harvest exponentially growing yeast in a 0.5 L culture in around 30 s32. These technological innovations allow the study of RBP-RNA dynamics at minute-scale resolution. Additionally, we also introduced several optimizations to the original CRAC protocol2 in order to increase its practicality.
Using χCRAC, we recently studied the targetome of a yeast nuclear RBP, Nab3, in response to glucose deprivation. In Saccharomyces cerevisiae, Nab3 can form a complex with Nrd1, an RBP, and the RNA helicase Sen1 to form the NNS complex. NNS binding to the RNA polymerase and the nascent transcript can trigger transcriptional termination34. This complex is mostly involved in removing cryptic noncoding RNA transcripts but has also been shown to regulate expression of protein-coding genes. The study showed differential targeting of Nab3 to noncoding and coding transcripts after only a minute of stress32. We demonstrated that co-transcriptional termination by Nab3 results in a very transient, pulse-like expression of retrotransposon genes, which would have been difficult to detect using traditional CLIP-based approaches. Additionally, the short UV irradiation times in our UV cross-linker also significantly increased the recovery of short-lived noncoding RNAs32. χCRAC will likely be a crucial tool in elucidating not only how RBPs shape the response to stress on immediate timescales but also their changing roles during the whole lifecycle of a response. This manuscript provides a detailed overview of all the steps in the χCRAC protocol. For illustrative purposes, the method was used to study the yeast Nrd1 protein, which is involved in transcriptional termination and RNA decay35,36, and its RNA targetome in response to glucose deprivation across a multitude of timepoints. Finally, we also demonstrate that our UV irradiation unit can rapidly cross-link RBPs to RNA in HeLa cells, making it possible to also perform high-resolution time-resolved analyses in adherent cells.
TN150 |
50 mM Tris pH 7.8 |
150 mM NaCl |
0.1% NP-40 |
1X protease inhibitor |
TN1000 |
50 mM Tris pH 7.8 |
1M NaCl |
0.1% NP-40 |
NP-PNK |
50 mM Tris-HCl pH 7.8 |
10 mM MgCl2 |
0.1% NP-40 |
5 mM beta-mercaptoethanol |
5 x PNK |
250 mM Tris-HCl pH 7.8 |
50 mM MgCl2 |
50 mM beta-mercaptoethanol |
WB I |
50 mM Tris-HCl pH 7.8 |
300 mM NaCl |
10 mM imidazole |
6M guanidine-HCl |
0.1% NP-40 |
5 mM beta-mercaptoethanol |
WB II |
50 mM Tris-HCl pH 7.8 |
50 mM NaCl |
10 mM imidazole |
0.1% NP-40 |
5 mM beta-mercaptoethanol |
Elution buffer |
50 mM Tris pH 7.8 |
50 mM NaCl |
250 mM imidazole |
0.1% NP-40 |
5 mM beta-mercaptoethanol |
Protease K buffer |
50 mM Tris |
0.1% NP-40 |
5 mM β-mercaptoethanol |
1% SDS |
5 mM EDTA |
50 mM NaCl2 |
Mammalian lysis buffer |
50 mM Tris-HCl pH 8 |
100 mM NaCl |
0.5% v/v Triton X-100 |
0.25% w/v Na-deoxycholate |
0.1% w/v SDS |
5 mM EDTA |
1 mM DTT (added fresh) |
1X protease inhibitor |
Table 1: The buffers required for χCRAC and their compositions.
1. UV cross-linking and lysate production
2. RBP capture
3. Washing the beads and TEV cleavage of the tags
4. On-bead alkaline phosphatase treatment
Component | 1x | 7.5x |
5 x PNK buffer | 12 | 90 |
Alkaline phosphatase | 4 | 30 |
RNase inhibitor | 2 | 15 |
H2O | 42 | 315 |
Final volume | 60 µL | 450 µL |
Table 2: Alkaline phosphatase reaction mixture.
5. On-bead ligation of the App-PE linker to the 3' end of the RNA
Component | 1x | 7.5x |
5 x PNK buffer | 12 | 90 |
App-PE adapter (100 µM) | 0.6 | 4.5 |
T4 RNA ligase 2 truncated K227Q | 3 | 22.5 |
RNase inhibitor | 1.5 | 11.25 |
50% PEG 8000 | 12 | 90 |
H2O | 30.9 | 231.75 |
Final volume | 60 µL | 450 µL |
Table 3: App-PE linker ligation reaction mixture.
Oligonucleotide name | Sequence (5'-3') | |||
L5Aa | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrUrArArGrCrN-OH | |||
L5Ab | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrArUrUrArGrCrN-OH | |||
L5Ac | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrGrCrGrCrArGrCrN-OH | |||
L5Ad | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrCrGrCrUrUrArGrCrN-OH | |||
L5Ba | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrArGrArGrCrN-OH | |||
L5Bb | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrGrUrGrArGrCrN-OH | |||
L5Bc | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrCrArCrUrArGrCrN-OH | |||
L5Bd | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrUrCrUrCrUrArGrCrN-OH | |||
L5Ca | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrCrUrArGrCrN-OH | |||
L5Cb | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrUrGrGrArGrCrN-OH | |||
L5Cc | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrArCrUrCrArGrCrN-OH | |||
L5Cd | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrGrArCrUrUrArGrCrN-OH | |||
L5Da | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrCrGrUrGrArUrN-OH | |||
L5Db | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrGrCrArCrUrArN-OH | |||
L5Dc | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrUrArGrUrGrCrN-OH | |||
L5Dd | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrArUrCrArCrGrN-OH | |||
L5Ea | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrCrArCrUrGrUrN-OH | |||
L5Eb | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrGrUrGrArCrArN-OH | |||
L5Ec | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrUrGrUrCrArCrN-OH | |||
L5Ed | invddT-ACACrGrArCrGrCrUrCrUrUrCrCrGrArUrCrUrNrNrNrArCrArGrUrGrN-OH | |||
App_PE | App-NAGATCGGAAGAGCACACGTCTG-ddC |
Table 4: The sequences of the DNA and RNA adaptors required for ligation onto the 5’ and 3’ ends of captured RNAs. These were purified through RNase-free HPLC.
6. On-bead phosphorylation of the 5’ ends of the RNA
Component | 1x | 7.5x |
5 x PNK buffer | 16 | 120 |
32P-ɣATP (10 µCi/µL) | 3 | 22.5 |
T4 PNK | 3 | 22.5 |
H2O | 58 | 435 |
Final volume | 80 µL | 600 µL |
Table 5: Phosphorylation reaction mixture.
7. On-bead ligation of the 5’ linker
NOTE: The 5’ linkers contain an RNA barcode that is used for identification of each sample after sequencing. Thus, it is absolutely crucial to note which linker is used for which sample.
Component | 1x | 7.5x |
5 x PNK buffer | 16 | 120 |
ATP (10 mM) | 8 | 60 |
RNase inhibitor | 2 | 15 |
T4 RNA ligase | 4 | 30 |
H2O | 48 | 360 |
Final volume | 78 µL | 585 µL |
Table 6: 5’ linker ligation reaction mixture.
8. Elution, SDS-PAGE, and RNA extraction
9. Reverse transcription
Oligonucleotide name | Sequence (5'-3') | |||
P5 forward | AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT | |||
BC1 | CAAGCAGAAGACGGCATACGAGATCGTGATGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT | |||
BC3 | CAAGCAGAAGACGGCATACGAGATGCCTAAGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT | |||
BC4 | CAAGCAGAAGACGGCATACGAGATTGGTCAGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT | |||
BC5 | CAAGCAGAAGACGGCATACGAGATCACTGTGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT | |||
BC7 | CAAGCAGAAGACGGCATACGAGATCAGATCGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT | |||
BC8 | CAAGCAGAAGACGGCATACGAGATTAGCTTGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT | |||
BC9 | CAAGCAGAAGACGGCATACGAGATGATCAGGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT | |||
BC10 | CAAGCAGAAGACGGCATACGAGATATCACGGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT | |||
PE_reverse | CAGACGTGTGCTCTTCCGATCT |
Table 7: The PCR primers (including the barcode sequences) and the reverse transcription primer.
10. qPCR reaction
Component | 1x |
2x qPCR reaction mastermix | 5 |
0.1 µM P5 primer (forward) | 0.8 |
0.1 µM BC primer (reverse) | 0.8 |
cDNA (or water as a negative control) | 1 |
H2O | 2.4 |
Final volume | 10 µL |
Table 8: qPCR reaction mixture.
11. PCR reaction and gel extraction
Component | 1x |
10 x proof-reading polymerase buffer | 5 |
10 µM P5 primer (forward) | 1 |
10 µM BC primer (reverse) | 1 |
5 mM dNTPs | 2.5 |
Proof reading polymerase enzyme | 1 |
cDNA | 5 |
H2O | 34.5 |
Final volume | 50 µL |
Table 9: PCR reaction mixture.
To demonstrate the efficacy of the χCRAC method, a time-course experiment with yeast strains expressing an HTP-tagged Nrd1 protein was performed. A detailed schematic representation describing how the method works is provided in Figure 1. Like Nab3, Nrd1 is involved in nuclear RNA decay of a variety of RNA transcripts37. Previous work from the Corden lab suggested that Nrd1 binding to its RNA targets changes significantly when cells are subjected to glucose starvation28,38. As such, cells growing exponentially in medium containing glucose (SD-TRP) were shifted to the same medium without glucose (S-TRP) over a time-course to monitor dynamic changes in Nrd1-RNA interactions. Samples were taken and cross-linked in the Vari-X-linker chamber (Figure 3A) before the shift and then after 1, 2, 4, 8, 14, and 20 min. The medium used for cell growth was deliberately deficient in tryptophan to reduce UV absorption by this aromatic amino acid. Note that it is best to use synthetic medium that is filter sterilized as because autoclaving the medium can lead to caramelization of the sugars. This then reduces the cross-linking efficiency.
Figure 4A shows a representative autoradiograph from a χCRAC experiment. Note that in this example, the samples were not pooled together. Instead, each was run individually on the gel. This is recommended for initial experimental tests to show that the protein cross-links effectively to RNA at all of the tested timepoints. A particularly intense signal was observed at the expected molecular weight of the RBP, representing the protein bound to very short, radiolabeled RNAs not amenable for sequencing. Therefore, the smeary signal above this band, which is the protein crosslinked to longer RNA fragments, was isolated. The fragment was cut from just above the protein band plus around 30 kDa. Figure 4B shows an autoradiogram following excision, with the protein cross-linked to short RNAs left in the gel and the previously smeary signal now excised.
After reverse transcription, the cDNA library must be amplified using PCR. However, overamplification of the library must be avoided as this can introduce bias towards sequences preferentially amplified by the polymerase and generate PCR artefacts. Overamplified libraries also contain a large number of duplicate sequences that waste reads on the sequencer. In order to calculate the ideal number of PCR cycles for amplification of the final library, an aliquot of the cDNA was amplified through qPCR using the P5 and BC oligonucleotides. The first cycle at which the library reached peak fluorescence was chosen as the PCR cycle count. Figure 4C gives an example of a qPCR from a typical cDNA library, which yielded a peak cycle count of 16. This value was then used for the final χCRAC PCR. In order to process the sequenced data, we used software previously developed in our lab (pyCRAC) and the corresponding pipeline for analysis of the kinetic CRAC data (Nues et al., 2017; https://git.ecdf.ed.ac.uk/sgrannem/pycrac, https://bitbucket.org/sgrann/kinetic_crac_pipeline/src/default/). These open source software tools enable demultiplexing and trimming of the data, removal of PCR duplicates, identification of statistically significant peaks, cluster reads into contiguous sequences, and identify binding motifs39. Further details of how these tools operate are found on their respective webpages.
We also started to develop a χCRAC protocol for mammalian cells. The majority of mammalian cell lines are grown as a monolayer and the tray in our crosslinker with the UV-permeable bag is not suitable for experiments with adherent cells. To overcome this problem, we developed a stage where users can UV irradiate 1–2 Petri dishes (150 mm diameter and 25 mm in depth) with adherent cells (Figure 3B). As a first test, the efficiency of the cross-linker for mammalian cells was measured through cross-linking and capture of stably tagged GFP-RBM7 using anti-GFP antibodies and a traditional CLIP-based purification. As shown in Figure 5A, the cross-linker was able to recover protein-RNA complexes from mammalian cells grown as a monolayer using 254 nm UV irradiation at efficiencies comparable to a widely used UV irradiation device. However, standard cell culture plasticware normally used for UV cross-linking experiments is impenetrable to 254 nm UV. Therefore, in our cross-linker the cells would only receive irradiation from the upper bank of UV lamps. To overcome this, we developed a UV-permeable quartz Petri dish for cell growth and cross-linking. Use of the quartz cultureware displayed robust recovery of protein-RNA complexes with as few as 2 s of UV irradiation (Figure 5B). When combined with RBP capture methods for mammalian cells such as CLIP technologies, these short cross-linking times are amenable with timecourses to recover spatiotemporal RNA-binding profiles of RBPs in response to genotoxic stresses or rapid depletions of protein factors, or in parallel with transcriptional or cell cycle synchronization.
Figure 6 shows several examples of the Nrd1 data processed by the χCRAC pipeline. This figure was prepared using the bedgraph files generated by the pipeline and the python GenomeBrowser package (https://pypi.org/project/GenomeBrowser/1.6.3/), which we designed to simplify making publication-quality genome browser images of the data. The grey rectangles represent genomic regions that expressed noncoding RNAs, such as the cryptic unstable transcript (CUTs), stable uncharacterized transcripts (SUTs)40, and Xrn1-sensitive unstable transcripts (XUTs)41. The data in Figure 6 show that Nrd1 binds to many of these noncoding RNA transcripts, consistent with the idea that this protein is involved in degradation of this class of transcripts42. Figure 6A shows a ~15 kb region on chromosome IV. Here there was a significant increase in binding of Nrd1 to transcripts encoding the high-affinity glucose transporters HXT6 and HXT7, both of which are upregulated during glucose starvation. It is likely that transcription termination by the NNS complex may influence the induction kinetics of these genes during glucose starvation. Figure 6B shows an example of Nrd1 cross-linking to the Imd3 transcript, which is known to be regulated by Nab343. In this case the data demonstrated a significant reduction in binding upon glucose starvation. Previous work showed decreased binding of Nab3 to the Tye7 transcript during glucose starvation44. Consistent with this observation,the χCRAC data suggest that binding of Nrd1 decreased during glucose starvation and Nrd1 cross-linking to Tye7 was at its lowest after 8 min of stress (Figure 4C). However, it appears that this effect was only transient, because after 14 min of glucose starvation, Nrd1 binding went back to starting levels.
Figure 1: Schematic representation of the χCRAC protocol. Tagged strains were grown until the desired density. RBP indicates RNA-binding protein. Afterwards, a reference sample was taken and cross-linked with 254 nm UV light. The remaining cells were harvested by filtration and then rapidly shifted to the stress-inducing medium. For the χCRAC experiment described here, samples were taken and cross-linked 1, 2, 4, 8, 14, and 20 min after the shift (1). The RBP of interest was then purified using a highly stringent two-step affinity purification (2). Next, the captured cross-linked RNAs were partially digested with RNases, radiolabeled at the 5’ end and adapters were ligated onto them (3). The 5’ adapters contained unique “in-read” barcode sequences so that the individual samples could be separated bioinformatically after sequencing. The RBP-RNA complexes were then eluted, pooled, and precipitated together (4), resolved by SDS-PAGE and visualized through autoradiography (5). Subsequently, a single gel slice containing the radioactive signal just above the main band, illustrated with dashed red box in the autoradiography image, was cut from the gel (5). The gel slices were treated with protease K and the RNA was subsequently extracted (6), converted into cDNAs, and amplified through PCR (7). The PCR step introduced additional barcodes (yellow block introduced by P7 oligo) so that many libraries could be multiplexed into a single lane. Please click here to view a larger version of this figure.
Figure 2: Cross-linking and vacuum filtration. (A) The cross-linker. The cell suspension is poured into a funnel located in the top right of the machine (also see Figure 3A for a close-up) and held in a UV-transparent bag located in the middle tray. This bag is flanked by two shutters that remain closed until the user instructs the machine to start the irradiation step. The cells are irradiated with UV light from the trays both above and below. The machine comes supplied with 254 and 365 nm UV lamps, with the latter being applicable for PAR-CLIP experiments. The machine is operated through a touchscreen panel located on the top right which allows one to control UV dosage or exposure time. (B) Following cross-linking, the cells are drained from the left-hand side of the machine. Cell suspensions are recovered through vacuum and drained into a glass flask where they can be subsequently poured into a vacuum filtration device for harvesting. (C) Vacuum filtration devices. These are opened and closed via a clip and a filter is inserted between it. Four filtration devices were used in parallel for very short time series to not lose any time as a result of changing filters. (D) Following filtration, the media supernatant was drained into flasks for subsequent disposal. Valves were installed below the vacuum filtration devices to maintain the vacuum in the system when the filter is removed. Please click here to view a larger version of this figure.
Figure 3: Cross-linking suspended vs. adherent cells. (A) The cross-linker with the Vari-X-linker chamber for suspension cells. The cell culture is poured into the sample inlet (funnel) located in the top right of the tray. (B) Tray that can hold plastic or quartz Petri dishes for cross-linking adherent cells or small volumes of suspension cells. Please click here to view a larger version of this figure.
Figure 4: Library preparation. (A) Example of an autoradiogram from a Nrd1-HTP χCRAC experiment. The strong, concentrated signal represents the protein crosslinked to very short RNAs, while the smear above represents the protein cross-linked to RNAs of sufficient length for sequencing. (B) The smear was excised as shown in an autoradiogram taken after gel excision. (C) A representative qPCR from a χCRAC cDNA library. In this example, maximum amplification of the cDNA was reached at 16 cycles. Thus, 16 cycles were used for the final amplification. The error bar represents the standard deviation of three technical qPCR replicates. (D) Example of a phosphorimage from a cDNA library on a 6% TBE gel. (E) cDNA length and quality analysis from a chip-based capillary electrophoresis. Please click here to view a larger version of this figure.
Figure 5: High RNase test iCLIP experiment to test crosslinking in mammalian cells. Shown are autoradiograms from GFP-RBM7 iCLIP experiments which tested the efficiency of RNP recovery across various cross-linking energies. Immunoprecipitations were performed using anti-GFP antibodies coupled to magnetic beads on cross-linked cells which stably expressed GFP-RBM7. Immunoprecipitates were incubated with high concentrations of RNase I in order to trim associated RNAs to short, uniform lengths. RNPs were visualized by 32P labelling and SDS-PAGE and migrate as a defined band, close to the migration of the non-cross-linked protein. Quantification indicates the results of densitometric analyses of radiolabeled RBM7-RNA signal normalized to the anti-GFP western blot signal. (A) Cross-linking time-course of the commonly used UVP cross-linker versus our cross-linker (Vari-X-linker; VxL). (B) Cross-linking time-course of our cross-linker on quartz (left) and plastic (right) cultureware. Please click here to view a larger version of this figure.
Figure 6: Example genome browser plots showing the power of χCRAC to show differential, temporal binding of Nrd1 to its targets. Each box shows plots for individual genomic regions. The arrows indicate on which strand the genes are encoded (left pointing arrow = minus strand; right pointing arrow = plus strand). The timepoints (min) are indicated by t0, t1, t2, etc on the y-axes of each subplot. Roman numerals indicating the chromosomes and the coordinates are shown. (A) Upon glucose deprivation, Nrd1 binds two high-affinity glucose transporters, HXT6 and HXT7, which are both upregulated in this condition. (B) Nrd1 is observed to bind to Imd3, an already validated target of Nab344, with reducing intensity following glucose starvation. (C) Nrd1 binding of Tye7 exhibits a dynamic and transient nature, decreasing after glucose starvation to a minimum after 8 min of stress. However, binding subsequently returns to basal levels after 14 min. Reads were normalized to “reads per million” (RPM; y-axis). Grey boxes indicate regions encoding noncoding RNAs. Please click here to view a larger version of this figure.
The χCRAC method, combined with the new cross-linking and cell harvesting devices, has great potential because it is applicable to a wide range of model organisms and therefore should be of general interest to the RNA field. There are many areas in which χCRAC can be utilized. For example, the method could be used to measure the hierarchical assembly of proteins into large macromolecular complexes, such as the spliceosome and the ribosome, which often involves dynamic interactions between proteins and RNA molecules. We also now routinely use it to monitor interactions between RNA decay factors and their substrates when cells are subjected to diverse kind of stresses. This enables us to determine at what stage of the adaptive response these factors are most active, what substrates they bind to, and how dynamic these interactions are. Such data should enable researchers to determine the relative contribution of each factor in adaptation to environmental changes.
χCRAC uses dual affinity purification tags (HTF or HTP) to purify the protein under highly stringent and denaturing conditions. This ensures that the copurified RNA is highly enriched for RNAs that were covalently cross-linked to the protein of interest. However, relying on affinity tags has disadvantages. For example, the tag could interfere with protein function, which could give a distorted readout of its RNA-binding interactome. Additionally, for some model organisms it may not always be possible to utilize tags because the genetic tools to integrate DNA fragments into the genome or to transform expression plasmids are not yet available. However, it is straightforward to alter some parts of the χCRAC protocol to make it compatible with CLIP-based protocols that rely on antibodies for purification of the RBP. Indeed, this study showed that it is possible to combine iCLIP-based purifications with our crosslinker. We are now in the process of developing CLIP protocols to study the temporal association of human RNA-binding proteins with nascent RNA transcripts.
When performing χCRAC on a new protein, the UV exposure must be optimized in order to induce maximal cross-linking. This is important because high UV exposures can reduce the recovery of RNA during the purification step. Cells expressing the recombinant RBP were exposed to various UV doses, 100 mJ/cm2, 250 mJ/cm2, 500 mJ/cm2, and 1 J/cm2. The RNPs were then captured and the RNAs were fragmented and radiolabeled. Afterwards, the RNPs were resolved by SDS-PAGE and an autoradiogram was taken in order to deduce which exposure gave the most intense signal (i.e. the maximal cross-linking).
Once the experimental conditions are optimized, several control experiments are recommended when performing χCRAC. First, a UV irradiated, untagged sample can be used to monitor background binding to the purification beads. Second, when applying χCRAC during a shift experiment, a second time series where the cells are shifted back into the original medium enables investigation into whether the filtration of the cells itself induces changes in RNA levels or protein-RNA interactions.
As mentioned in the Introduction, numerous recently published papers suggest a number of optimizations to the CLIP protocol. This includes the use of fluorescently labeled adapters for detecting the protein-RNA complex through infrared scanning10 as well as optimizations to various nucleic acid purification and size selection steps shown to increase the complexity of the resulting libraries12,45. We are currently implementing some of these improvements to further refine the χCRAC protocol. The protocol presented here already contains a number of improvements to the original CRAC and χCRAC protocols that increase the complexity of the data. For example, previously, after resolving the cross-linked, radioactive protein-RNA complexes on SDS-PAGE gels, they were transferred to a nitrocellulose membrane and the cross-linked RNA was isolated from the blot. However, the transfer of the RNP and subsequent RNA extraction can be very inefficient, particularly when dealing with large RBPs such as RNA polymerase subunits. This can result in a significant reduction in the recovery of the cross-linked RNA. In the current protocol, the cross-linked RNA is extracted directly from SDS-PAGE gel slices as illustrated in Figure 1. This increased the recovery of cross-linked RNAs. Additionally, after PCR amplification of the cDNAs the product was originally resolved on 3%, low melting temperature agarose gels, and then 175–300 bp PCR products were extracted from the gel. However, these gels can be easily overloaded, resulting in very poor separation of the DNA. Replacing agarose gels with precast TBE gels resulted in more consistent size separation and better recovery of PCR products.
The authors have nothing to disclose.
This work was supported by grants from the Wellcome Trust (091549 to S.G and 109093/Z/15/A to S.M.), the Wellcome Trust Centre for Cell Biology core grant (092076) and Medical Research Council Non-Clinical Senior Research Fellowship (MR/R008205/1 to S.G.), the European Molecular Biology Organization under a long-term postdoctoral fellowship (ALTF 1070-2017 to R.A.C), and the Independent Research Fund Denmark (T.H.J).
1,4-dithioreitol | Merck | 10708984001 | Buffer component in mammalian cell lysis |
1.5 mL tubes | Eppendorf | 0030 120.086 | General reaction tube |
2 mL tubes | Eppendorf | 0030 123.344 | For holding columns and collection of waste |
32P-yATP | Perkin Elmer | NEG502Z-250 | For radiolabelling the 5' end of the RNA |
4-12% Bis-Tris gel | Invitrogen | NP0321BOX | SDS-PAGE gel |
4X loading buffer | Novex | NP0008 | Protein loading dye concentrate |
50 bp ladder | New England Biolabs | N3236 | Reference ladder for excising region of interest from the amplified cDNA library |
50% PEG | NEB | B100045 | For the L5 linker ligation |
6% TBE gel | Invitrogen | EC6265BOX | For separation and purification of the cDNA library |
Acetone | ACROS Organics | 423245000 | Washing of TCA-precipitated proteins |
anti-FLAG beads | Sigma Aldrich | M8823-1ML | For purifcation of FLAG-tagged RBPs |
ATP (100 mM) | Thermo Fisher Scientific | R0441 | For ligation of the L5 linker onto the 5' end of captured RNAs |
Beta-mercaptoethanol | Sigma Aldrich | M3148-100ML | Buffer component |
Biomax MS intensifying screen | Sigma Aldrich | Z363162-1EA | For intensifying the autoradiogram signal |
Chloroform | Thermo Fisher Scientific | 1010219 | For phenol-chloroform extraction following RNA purification |
cOmplete EDTA-free protease inhibitor cocktail | Roche | 11873580001 | For inhibition of cellular proteases after lysis |
Complete supplement mixture -TRP | Formedium | DCS0149 | For preparation of synthetic defined medium |
Costar Spin-X 0.22 µm filters | Sigma Aldrich | CLS8160 | For isolating the excised cDNAs following gel extraction |
DNase RQ1 | Promega | M6101 | For DNA digest following cell lysis |
dNTPs (10 mM) | Sigma Aldrich | 4638956001 | For reverse transcription and PCR |
Ethanol | Thermo Fisher Scientific | 10041814 | For phenol-chloroform extraction following RNA purification and DNA precipitation |
Ethylenediaminetetraacetic acid | Invitrogen | AM9261 | For protease K buffer |
Exonuclease I | New England Biolabs | M0293 | For degradation of primers following PCR |
Glass microfiber filters | Whatman | 1823-010 | For isolating the excised cDNAs following gel extraction |
Glucose | Formedium | GLU03 | For preparation of glucose-containing, synthetic defined medium |
Glycogen (20 mg/mL) | Roche | 10901393001 | Precipitation of proteins, RNA and DNA |
GST-TEV | Homemade | Construct and purification protocol is available upon request | |
Guanidium hydrochloroide | Thermo Fisher Scientific | 10071503 | Required for pulldown denaturing conditions and washing buffer |
IgG beads | GE Healthcare | 17-0969-01 | For purification of protein A-tagged RBPs |
Imidazole | Sigma Aldrich | I2399-100G | For elution of captured proteins from Nickel beads |
Isoamyl alcohol | Thermo Fisher Scientific | A393-500 | For phenol-chloroform extraction following RNA purification |
Luna Universal One-Step RT-qPCR | NEB | E3005S | For qPCR of the cDNA in order to calculate required number of PCR cycles |
Magnesium chloride | Fluka Analytical | 63020-1L | For PNK buffer |
Membrane filters | Millipore | AAWP09000 for yeast or HAWP09000 for bacteria | For vacuum filtration of cells |
Micro bio-spin columns | Biorad | 732-6204 | For collecting eluate after gel extraction |
Ni-NTA beads | Qiagen | 30210 | For secondary protein capture |
NP-40 | Sigma Aldrich | I8896-100ML | Buffer component |
Pfu polymerase | Promega | M7741 | For amplification of the cDNA library |
Phenol | Sigma Aldrich | P4682-400ML | For phenol-chloroform extraction following RNA purification |
Pierce spin columns | Thermo Fisher Scientific | 69725 | For on-column enzymatic reactions |
Protease K | Roche | 3115887001 | For degradation of the RBP following gel extraction |
Quartz Petri dish | UVO3 | N/A | For cross-linking of adherent cells. Available from https://www.vari-x-link.com for 400 GBP |
Radiography films | Amersham | 28906843 | For autoradiography visualisation |
RNAClean XP beads | Beckmann | A63987 | SPRI beads for clean up of RNAs and cDNAs |
RNase H | New England Biolabs | M0297 | For degradation of RNAs following reverse transcription |
RNase-It | Agilent | 400720 | For RNA digestion |
rRNasin | Promega | N2511 | For inhibition of any contaminating RNases during enzymatic reaction |
Sodium acetate | Sigma Aldrich | S2889-1KG | For phenol-chloroform extraction following RNA purification and DNA precipitation |
Sodium chloride | Thermo Fisher Scientific | 7647-14-5 | Buffer component |
Sodium deoxycholate | Sigma Aldrich | D6750-100G | Buffer component in mammalian cell lysis |
Sodium dodecylsulfate | Sigma Aldrich | L3771-1KG | For protease K buffer |
SUPERase-In | Invitrogen | AM2694 | For inhibition of cellular RNases after lysis |
SuperScript IV | Thermo Fisher Scientific | 18090010 | For reverse transcription |
T4 PNK | New England Biolabs | M0201 | For radiolabelling the 5' end of the RNA |
T4 RNA ligase 1 | New England Biolabs | M0204 | For ligation of the L5 adaptor onto the RNA 5' end |
T4 RNase ligase 2, truncated K222Q | NEB | M0351S | For ligation of the App_PE linker onto the 3' end of captured RNAs |
TBE buffer (10X) | Invitrogen | 15581-028 | For running TBE gels |
TEV protease | Homemade | For eluting captured proteins following FLAG capture | |
Thermosensitive alkaline phosphatase | Promega | M9910 | For 5' and 3' dephosphorylation of RNAs |
Trichloroacetic acid (100%) | Sigma Aldrich | T0699-100ML | For precipitation of RBP-RNA complexes |
Tris hydrochloride | Invitrogen | 15504-020 | Buffer component |
Triton X-100 | Sigma Aldrich | T8787-100ML | Buffer component in mammalian cell lysis |
Vari Filter | UVO3 | N/A | Device for vacuum harvesting cells. Available from https://www.vari-x-link.com for 100 GBP |
Vari-X-Linker | UVO3 | N/A | Cross-linker for cross-linking cells. Available from https://www.vari-x-link.com for 16,000 GBP |
Yeast nitrogen base | Formedium | CYN0410 | For preparation of synthetic defined medium |
Zirconia beads | Thistle | 11079105Z for yeast or 11079101Z for bacteria | For cell lysis via bead beating |