$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
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.