RNA localization is an important regulatory layer of gene expression and cell functioning. The protocol guides through the Proximity RNA-seq method, in which RNA molecules are sequenced in their spatial, cellular context to derive RNA co-localization and transcriptome organization. Transcripts in individual subcellular particles from chemically crosslinked cells are tagged with the same, unique DNA barcode in water-in-oil emulsion droplets. First, single DNA barcodes are PCR amplified and immobilized on single, small magnetic beads in droplets. Subsequently, 3' ends of bead-bound barcode copies are tailed with random pentadecamers. Then beads are encapsulated again into droplets together with crosslinked subcellular particles containing RNA. Reverse transcription using random pentadecamers as primers is performed in droplets, which optimally contain one bead and one particle, in order to tag RNAs co-localized to the same particle. Sequencing such cDNA molecules identifies the RNA molecule and the barcode. Subsequent analysis of transcripts that share the same barcode, i.e., co-barcoding, reveals RNA co-localization and interactions. The technique is not restricted to pairs of RNAs but can as well detect groups of transcripts and estimates local RNA density or connectivity for individual transcripts. We provide here a detailed protocol to perform and analyze Proximity RNA-seq on cell nuclei to study spatial, nuclear RNA organization.