Studying how macromolecules fold and interact with each other is the key to understanding gene regulation in the cell. While much effort has been focused in the past decade on understanding how DNA and proteins contribute to gene regulation, relatively less is known about post-transcriptional regulation of gene expression. RNA carries information in both its linear sequence and in its secondary and tertiary structure1. Its ability to base pair with itself and with others is important for its function in vivo. Recent advances in high throughput RNA secondary structure probing has provided valuable insights into the locations of double and single stranded regions in the transcriptome2,3,4,5,6,7,8, however information on the pairing interaction partners is still largely missing. To determine which RNA sequence is interacting with another RNA region in the transcriptome, we need global pair-wise information.
Mapping pair-wise RNA interactions in a global, unbiased manner has traditionally been a major challenge. While previous approaches, such as CLASH9, hiCLIP10 and RAP11, are used to identify RNA interactions in a large scale manner, these techniques typically map RNA base pairing for a subset of RNAs that either interact with a particular protein or RNA species. Recent developments in studying global RNA interactions include the method RPL12, which does not stabilize RNA interactions in vivo and hence may only capture a subset of in vivo interactions. To overcome these challenges, we and others developed genome-wide, unbiased strategies to map RNA interactomes in vivo, using modified versions of the crosslinker psoralen13,14,15. In this protocol, we describe the details for performing Sequencing of Psoralen crosslinked, Ligated, and Selected Hybrids (SPLASH), which utilizes biotinylated psoralen to crosslink base pairing RNAs in vivo, followed by proximity ligation and high throughput sequencing to identify RNA base-pairing partners genome-wide (Figure 1)15.
In this manuscript, we describe the steps to perform SPLASH using cultured adherent cells, in this case HeLa cells. The same protocol can be easily adapted to suspension mammalian cells and to yeast and bacteria cells. Briefly, HeLa cells are treated with biotinylated psoralen and irradiated at 365 nm to crosslink interacting RNA base pairs in vivo. The RNAs are then extracted from the cells, fragmented and enriched for crosslinking regions using streptavidin beads. Interacting RNA fragments are then ligated together using proximity ligation and made into a cDNA library for deep sequencing. Upon sequencing, the chimeric RNAs are mapped onto the transcriptome/genome to identify the RNA interacting regions that are paired to each other. We have successfully utilized SPLASH to identify thousands of RNA interactions in vivo in yeast and different human cells, including intramolecular and intermolecular RNA base pairing in diverse classes of RNAs, such as snoRNAs, lncRNAs and mRNAs, to glimpse into the structural organization and interaction patterns of RNAs in the cell.