Xiaojun Ren1,2, Ruijie Deng1, Kaixiang Zhang1, Yupeng Sun1, Xucong Teng1, Jinghong Li1. 1. Department of Chemistry, Key Laboratory of Bioorganic Phosphorus Chemistry & Chemical Biology, Tsinghua University, Beijing 100084, China. 2. School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China.
Abstract
Immune cell heterogeneity due to the differential expression of RNA splicing variants still remains unexplored. This is mainly because single-cell imaging technologies of splicing variants with precise sequence or base resolution are now not readily available. Herein, we design a splice-junction anchored padlock-probe-mediated rolling circle amplification assay (SpliceRCA) for single-cell imaging of splice isoforms of essential regulatory immune gene (CD45) upon T-cell activation. Two recognition regions in the padlock probe can target the splice-junction sequence, resulting in a close proximity for triggering in situ one-target-one-amplicon amplification. With the read length of ∼30 nucleotides, this method allows discrimination of isoforms with single-base precision and quantification of isoforms with single-molecule resolution. We applied SpliceRCA to single-cell image splice variants of essential regulatory immune gene (CD45) upon T-cell activation. It is found that CD45RO isoform presents a distal nuclear spatial distribution and is coregulated with CD45RB upon activation. Our strategy provides a single-cell analysis platform to investigate the mechanism of complex immune responses and may further guide immunotherapy.
Immune cell heterogeneity due to the differential expression of RNA splicing variants still remains unexplored. This is mainly because single-cell imaging technologies of splicing variants with precise sequence or base resolution are now not readily available. Herein, we design a splice-junction anchored padlock-probe-mediated rolling circle amplification assay (SpliceRCA) for single-cell imaging of splice isoforms of essential regulatory immune gene (CD45) upon T-cell activation. Two recognition regions in the padlock probe can target the splice-junction sequence, resulting in a close proximity for triggering in situ one-target-one-amplicon amplification. With the read length of ∼30 nucleotides, this method allows discrimination of isoforms with single-base precision and quantification of isoforms with single-molecule resolution. We applied SpliceRCA to single-cell image splice variants of essential regulatory immune gene (CD45) upon T-cell activation. It is found that CD45RO isoform presents a distal nuclear spatial distribution and is coregulated with CD45RB upon activation. Our strategy provides a single-cell analysis platform to investigate the mechanism of complex immune responses and may further guide immunotherapy.
Alternative splicing is a fundamental
regulatory process of gene
expression that allows generation of multiple mRNA isoforms from single
genes,[1,2] thereby increasing transcriptome complexity.[3,4] Splicing variability among individual cells accounts for a large
part of gene expression heterogeneity, which plays a significant role
in the immune system for the efficient battling rich variability of
pathogens.[5,6] Isoforms with diversified functions were
evolved in the generation of the immune response.[7] CD45, the prototypic receptor-like protein tyrosine phosphatase
gene, acts as an essential regulator of signal transduction pathways
in immune cells.[8] The transition from naïve
to activated T cells is marked by CD45 pre-mRNA alternative splicing.[9] Abnormal CD45 splice variant expressions are
associated with autoimmune and infectious diseases.[10,11] However, the roles of differential CD45 splicing variants, such
as their relative expression level, intracellular localization, expression
heterogeneity, and regulation in eliciting immune responses, are underexplored.
Thus, visualizing[12] intracellular mRNA
variants in single cells is essential for understanding differential
splicing-mediated immune cell heterogeneity and resolving the complexity
and heterogeneity of RNA variant-related immune diseases.Currently,
only a limited number of methods have been explored
to detect RNA splicing variants at the single-cell level. Single-molecule
fluorescence in situ hybridization (smFISH) is a
facile and powerful method to in situ image RNA splice
variants in single cells, and brings a significant advance of RNA
splicing in single-cell study.[13,14] However, the FISH method
can only detect RNA sequences with a minimum length of 600 nucleotides
(nt) (usually need >30 different hybridization probes with ∼20
nt to ensure that the labeling signals are distinguished from the
background).[15,16] However, the human exon is with
the average length of only 320 nt; thus, large amounts of exons can
hardly be detected by smFISH.[17,18] Recently, a newly proposed
plasmonic in situ hybridization (plasmonic ISH) based
on gold nanoparticle (AuNP) labeling is masterly adapted to analyze
RNA splicing and is able to differentiate splice variants with ∼300
nt sequence.[19] However, limitations on
AuNPs remain currently unsolvable, such as their large size for impeding
the delivery efficiency, alteration in a real cell state, and easy
aggregation, etc.[20,21] The visualization of short exon
mRNA variants still remains a challenge. Meanwhile, a large amount
of short exon mRNA variants are involved in crucial biological processes,[22] like CD45 isoforms, of which the average length
of the alternative exon is less than 200 nt.[23] There is still an urgent need for developing a versatile RNA imaging
method with high sequence resolution capable of in situ analyzing splice isoforms with different lengths.To address
these issues, we develop a high base-resolution strategy
termed as splice-junction anchored padlockprobe-mediated rolling circle
amplification[24] (SpliceRCA), enabling single-cell
imaging of CD45 splicing variants. Rolling circle amplification (RCA)
can achieve localized isothermal amplification, converting the target
sequence into a long single-stranded DNA or RNA product with thousands
of tandem repeats.[25] Attributed to its
one-target-one-amplicon amplification process, RCA can achieve target
RNA/DNA detection or imaging at the single-molecule level.[26] Furthermore, the padlockprobe-based RCA method
has the capability to target short RNAs and discriminate highly similar
sequences to genotype RNAs with single-nucleotide variations,[27] which prompted us to explore the potential of
RCA in RNA splicing variant detection. In this study, two recognition
regions in the padlock probe are specifically hybridized to a newly
formed splice-junction sequence, resulting in a close proximity for
triggering in situ one-target-one-amplicon amplification,[28] achieving shortening the read length of the
imaging method to ∼30 nt. This method allows discrimination
of isoforms with single-base precision and quantification of isoforms
with single-molecule[29] resolution. With
this method, we measure the isoform expression and spatial distribution
of three critical variants of CD45, and analyze how the isoforms covariation
between individual cells changes in T-cell activation. The method
provides a single-cell-level RNA splicing analysis platform to explore
the differential expression of splicevariant-mediated immune cell
heterogeneity in a quantitative manner.
Results and Discussion
Overview
of SpliceRCA
Scheme illustrates the procedures of the direct
visualization of RNA splice variants in single cells by SpliceRCA.
Alternative splicing of exon 4 (198 nt), 5 (141 nt), and 6 (144 nt)
in CD45 pre-mRNA is strictly regulated in T-cell activation: naïve
T cells express various larger isoforms whereas memory T cells tend
to produce small isoforms.[30] We chose three
critical alternative splicing variants of the CD45 gene involved in
the T-cell activation process which include the largest isoform CD45RA
(all exons retained), middle isoform CD45RB (exon skipping, exon 4
deleted), and smallest isoform CD45RO (exon skipping, exon 4, 5, and
6 deleted) (Scheme A). SpliceRCA is performed with a splice-junction anchored padlock
probe composed of four modules: the recognition of exon junction sites
(Rx, Ry), universal priming region (P), and tag motif (T) modules
(Scheme B). The splice
junction in the target splice isoform brings close proximity between Rx and Ry in the padlock probe for circularizing,
following triggering in situ RCA upon the hybridization
of primer with the P. By tuning the sequence of T corresponding to different fluorophores, three splice isoforms
can be simultaneously differentiated, and visualized with single-molecule
resolution attributed to the in situ one-target-one-amplicon
amplification method. Benefiting from the recognition of the target
splice-junction sequence, the read length of the imaging method can
be shortened to ∼30 nt. Moreover, the ligation reaction confers
high-specificity for discrimination between perfectly matched and
mismatched oligonucleotides, which can realize the imaging of RNA
isoforms with even single-base resolution.[31]
Scheme 1
Schematic Diagram of Multiplex Detection of mRNA Variants in Single
Cells by SpliceRCA
(A) Alternative
splicing patterns
of CD45 during T-cell activation. Isoforms (CD45RA, CD45RB, CD45R0)
with decreasing exon inclusion were expressed upon T-cell activation.
(B) The procedures of SpliceRCA for detecting splice variants in single
cells. The splice-junction anchored padlock probe is composed of four
modules: the recognition of exon junction sites (Rx, Ry), universal
priming region (P), and tag motif (T) modules. The newly formed splice
junction in the target splice isoform brings close proximity between
Rx and Ry in the padlock probe for circularizing, following primer
hybridized with the P, triggering in situ RCA. Upon
tuning of the sequence of T corresponding to different fluorophores,
the three RNA splicing variants can thus be simultaneously differentiated,
and visualized with single-molecule resolution attributed to the in situ one-target-one-amplicon amplification method.
Schematic Diagram of Multiplex Detection of mRNA Variants in Single
Cells by SpliceRCA
(A) Alternative
splicing patterns
of CD45 during T-cell activation. Isoforms (CD45RA, CD45RB, CD45R0)
with decreasing exon inclusion were expressed upon T-cell activation.
(B) The procedures of SpliceRCA for detecting splice variants in single
cells. The splice-junction anchored padlock probe is composed of four
modules: the recognition of exon junction sites (Rx, Ry), universal
priming region (P), and tag motif (T) modules. The newly formed splice
junction in the target splice isoform brings close proximity between
Rx and Ry in the padlock probe for circularizing, following primer
hybridized with the P, triggering in situ RCA. Upon
tuning of the sequence of T corresponding to different fluorophores,
the three RNA splicing variants can thus be simultaneously differentiated,
and visualized with single-molecule resolution attributed to the in situ one-target-one-amplicon amplification method.
Multiplex in Situ Imaging
of CD45 mRNA Splicing
Variants in Single Cells by SpliceRCA
The specific scheme
of in situ RCA is shown in Figure A. The splice-junction sequence was specifically
recognized by the padlock probes and as the template for ligating
the padlock probe. Then, RCA is initiated with the help of an additional
DNA primer, resulting in a long DNA amplicon with hundreds of copies
of the padlock probe.[32,33] The RCA amplicon forms as a nanoclew
that can be visible as a diffraction-limited fluorescent spot upon
hybridization with different detection probes. To demonstrate the
feasibility of this principle for detection of RNA splice isoforms,
we first performed the assay in vitro to amplify
the target synthetic sequence of CD45RO, CD45RA, and CD45RB (Table S1, Supporting Information, SI). This method
can effectively discriminate isoforms and exhibit a high-sensitivity
performance for in vitro detection (Figures S1 and S2).
Figure 1
Simultaneous imaging of CD45 splicing variants
in single Jurkat
T cells. (A) Schematic diagram of SpliceRCA for detection of CD45
RNA isoforms using in situ RCA. (B) Fluorescent image
of CD45RO, CD45RA, and CD45RB isoforms visualized by SpliceRCA in
Jurkat T cells. Inset: the fluorescent image of the boxed subregion
with colored labels indicating potential mRNA isoforms. The outline
of the Jurkat T cell is marked with a gray dot line. Scale bars: 10
μm (overview images) and 500 nm (insets). (C) Dynamic light
scattering (DLS) analysis for RCA amplicons. (D) Intensity plots of
the lines in panel B show that mRNA spot signals are distinguished
from background. (E) Frequency histogram of amplicons of three splicing
isoforms in single cells (cell number >100). (F) Comparison of
SpliceRCA
results in simultaneous detection with separate detection (n > 100; NS, not significant). (G) Quantification of
the
average numbers of CD45RB amplicons detected in Jurkat cells by using
Mis-0, Mis-1, and Mis-2 splice-junction anchored padlock probes (n > 100; **P < 0.001).
Simultaneous imaging of CD45 splicing variants
in single Jurkat
T cells. (A) Schematic diagram of SpliceRCA for detection of CD45
RNA isoforms using in situ RCA. (B) Fluorescent image
of CD45RO, CD45RA, and CD45RB isoforms visualized by SpliceRCA in
Jurkat T cells. Inset: the fluorescent image of the boxed subregion
with colored labels indicating potential mRNA isoforms. The outline
of the Jurkat T cell is marked with a gray dot line. Scale bars: 10
μm (overview images) and 500 nm (insets). (C) Dynamic light
scattering (DLS) analysis for RCA amplicons. (D) Intensity plots of
the lines in panel B show that mRNA spot signals are distinguished
from background. (E) Frequency histogram of amplicons of three splicing
isoforms in single cells (cell number >100). (F) Comparison of
SpliceRCA
results in simultaneous detection with separate detection (n > 100; NS, not significant). (G) Quantification of
the
average numbers of CD45RB amplicons detected in Jurkat cells by using
Mis-0, Mis-1, and Mis-2 splice-junction anchored padlock probes (n > 100; **P < 0.001).Next, we explored the potential of the SpliceRCA
method for imaging
CD45 RNA splicing variants in Jurkat T cells. As illustrated in Figure B, the generated
superbright dots amplified from the target splice isoform could be
clearly distinguished from the background inside cell. These bright
dot signals were generated from one-target-one-amplicon in
situ amplification, thus presumably corresponding to separated
mRNAs. Characterization of RCA amplicons by dynamic light scattering
(DLS) and transmitted electron microscopy (TEM) demonstrated that
the RCA amplicons are monodisperse particles with a size of ∼300
nm (Figure C and Figure S3). The size of the RCA amplicon is larger
than the resolution of optical imaging, thus conferring the single-molecular
RCA amplicons resolvable by confocal microscopy.[16,34] The signal intensities of the RCA amplicon amplified from individual
splice isoforms were much larger than the background (Figure D). Next, for verification
that signals were amplified from the target, several experiments were
performed as control. As shown in Figure S4, no bright spot was observed when no padlock probe or trigger primer
was used. There was only a rare fluorescence signal when random probes
were added. Additionally, to confirm that the signal resulted from
the target sequence, we blocked the binding sites with an unlabeled
complementary probe before carrying out SpliceRCA. Only less than
0.2 amplicons per cell could be seen, thus suggesting that the bright
spots came from the amplification of the target RNA sequence. Furthermore,
we also performed a siRNA knockdown experiment to suppress the expression
of CD45; the signals of CD45RO, CD45RA, and CD45RB isoforms decreased
24.59%, 31.27%, and 27.01%, respectively, with the knockdown of CD45
(Table S3). The knockdown efficiency was
verified by RT-qPCR assay (Figure S5 and Table S4). All of these results suggest that the bright dot signals
come from target mRNA splicing isoforms.As shown in Figure E, CD45 isoforms
were presented at different levels. The average
numbers of amplicons for CD45RO, CD45RA, and CD45RB were 3.70, 5.98,
and 11.07 per cell. The in situ detection efficiency
of mRNA splice isoform-initiated RCA was estimated to be 10–20%
on the basis of a comparison to RT-qPCR data (Table S5 and discussed in the Supporting Information). The
obvious variability in the copy numbers of CD45 isoforms indicates
that significant cell-to-cell variation in isoform expression would
be exhibited even in the same batch cells. In addition, SpliceRCA
has been successfully applied to image splicing isoforms of BRCA1,
breast cancer susceptibility gene 1.[19,35] The alternative
splicing of BRCA1 was closely related with the transformation of malignant
breast cancer (Figure S6).[36,37] As a validation, we performed the reverse transcription quantitative
PCR (RT-qPCR) assay for expression comparison. The results of SpliceRCA
are in good accordance with the RT-qPCR results in general (Figures S7 and S8). The single-molecule resolution
of the SpliceRCA method allows precise quantification of isoform abundances.Moreover, the assurance of the discrimination among the various
splicing variants, especially which have sharing sequences in splice
junction, highly depends on the stringent specificity of the RNA profiling
method. To assess the accuracy of SpliceRCA for profiling RNA splicing
variants inside cells, we compared the isoform expression level measured
by SpliceRCA in simultaneous and separate detection. The quantification
of the three isoforms by simultaneous and separate detection shows
negligible difference (Figure F and Figure S9, P < 0.001), indicating minimal crosstalk among the targets simultaneously
detected. Thus, SpliceRCA confers low crosstalk for multiplex in situ mRNA isoforms detection. Further, we test the base
precision of the SpliceRCA method, in which the splice-junction padlock
probe used for imaging CD45RB was altered by one or two bases. The
copy number of amplicons per cell decreased sharply after the introduction
of a one- or two-base mismatch (Figure S10). The average numbers of amplicons when using one-base mismatched
(Mis-1) and two-base mismatched (Mis-2) padlock probes were 0.88 per
cell and 0.20 per cell, respectively, much less than the 11.02 per
cell measured with the matched probe (Mis-0) (Figure G). Thus, the SpliceRCA presents high specificity
conferring single-nucleotide resolution, ensuring the precise recognition
of mRNA isoforms in situ.
Quantification of Isoform
Variability in Jurkat Cells upon T-Cell
Activation
Determining the alternative splicing changes induced
by T-cell activation is crucial for understanding the cellular outcomes
of the antigen challenge. To characterize the extent of expression
variability on CD45 mRNA isoforms and decipher its functional implications,
we applied SpliceRCA to profile the isoform expression in the Jurkat
T-cell response to phorbol-12-myristate13-acetate (PMA) stimulation.
Immunofluorescence imaging of HnRNPLL, one marker of T-cell maturation,[38,39] shows a skewed increase in expression level after PMA stimulation
(Figure S11), indicating that Jurkat T
cells were efficiently stimulated. Fluorescent images reveal a different
isoform pattern after stimulation, and the copy numbers of three mRNA
isoforms all present a shift upon T-cell activation (Figure A,B). Moreover, the expressions
of CD45RO and CD45RB isoforms show 70% and 32% increase (Figure C). Then, RT-qPCR
analysis further confirms the isoform expression patterns and the
splicing shift in PMA stimulation (Figure S12). These results indicate that T-cell activation leads to the skipping
of three variable exons in CD45 pre-mRNA (Figure D). The skipping process would reduce the
phosphatase activity of the CD45 molecule and preserve the normal
status to prevent autoimmune disease.[30,40] Interestingly,
cell-to-cell variability in isoform ratios differed between the two
cell states before and after PMA stimulation. After stimulus, the
three isoforms showed subtle lower cell-to-cell variability [coefficients
of variation (CV, SD/mean) of CD45RO, CD45RA, and CD45RB were 0.328,
0.361, and 0.283, respectively] than resting T cells (isoform CV >
0.4) (Table S6). The variability in single-cell
isoform expression may reflect functional differences in the stimulated
Jurkat T cell population.
Figure 2
Cell-to-cell variation in isoform expression
upon PMA stimulation.
(A) Fluorescent images of three CD45 splicing isoforms visualized
by SpliceRCA before and after PMA stimulation. (B) Histograms of the
copy numbers of three mRNA isoforms in Jurkat T cells measured by
SpliceRCA before and after stimulus. The arrow shows the expression
shift; arrow to the right (left) means increased (decreased) expression
upon PMA stimulation. (C) Quantification of average expression for
three CD45 isoforms before and after stimulus. (D) Schematic diagram
of the alternative splicing upon T-cell activation.
Cell-to-cell variation in isoform expression
upon PMA stimulation.
(A) Fluorescent images of three CD45 splicing isoforms visualized
by SpliceRCA before and after PMA stimulation. (B) Histograms of the
copy numbers of three mRNA isoforms in Jurkat T cells measured by
SpliceRCA before and after stimulus. The arrow shows the expression
shift; arrow to the right (left) means increased (decreased) expression
upon PMA stimulation. (C) Quantification of average expression for
three CD45 isoforms before and after stimulus. (D) Schematic diagram
of the alternative splicing upon T-cell activation.
Mapping the Spatial Distribution of CD45
mRNA Variants in Single
Cells
The location distribution of splicing variants from
the same gene has not been well-characterized in mammalian cells,
let alone suspension cells such as T cells. We then exploit the in situ RNA visualization ability of SpliceRCA to investigate
the spatial distributions of multiple splice variants simultaneously.
From the visual inspection of each mRNA variant in resting Jurkat
T cells, isoform CD45RO appeared enriched near the cell periphery,
whereas CD45RA and CD45RB appeared in a random distribution throughout
the cell (Figure A).
Quantitative analysis of the distances between each mRNA splicing
variant and the cell nucleus or the cell periphery further confirmed
the visual impression (Figure B). The average distances to the nucleus for the mRNA variants
CD45RO, CD45RA, and CD45RB were 3.25, 3.32, and 4.15 μm, respectively.
The average distances to cell edge for mRNA variants CD45RA, CD45RB,
and CD45RO were 2.30, 3.24, and 2.47 μm, respectively. Furthermore,
the spatial distribution of splice isoforms after T-cell activation
was investigated. As shown in Figure C, the average distances to the nucleus for CD45RO,
CD45RA, and CD45RB increased by 0.07, 0.08, and 0.67 μm, respectively,
after T-cell activation. The average distances to the cell edge for
CD45RO, CD45RA, and CD45RB decreased by 0.02, 0.08, and 0.54 μm,
respectively, after T-cell activation. The results indicate that CD45RO
and CD45RA showed a modest change in the distance to the nucleus and
cell edge upon T-cell activation, while the CD45RB presented an obvious
change to be more likely expressed in the cell periphery after T-cell
activation. The different spatial patterns may result from the expression
of some RNA binding proteins as splice mRNA localization is closely
related with regulation of many interaction proteins.[41,42] The ability for spatial mapping of splice isoforms in single cells
may offer us new ways to predict and study the rich diversity in splicing
functions.
Figure 3
Mapping the spatial distribution of CD45 isoforms in the Jurkat
T cells. (A) Example images of the spatial distributions observed
for CD45 isoforms before and after PMA stimulation. The RCA amplicons
of CD45RO, CD45RA, and CD45RB are marked with green, blue, and red
circles, respectively, and the outlines of cell are marked by a gray
dotted line. Scale bars: 10 μm. (B) Single-cell profiles for
distance to cell nucleus (yellow dots) and the edge (black dots). N is the number of cell dots per plot. (C) Quantification
of average distances for three CD45 isoforms to the cell nucleus or
the edge before and after PMA stimulation.
Mapping the spatial distribution of CD45 isoforms in the Jurkat
T cells. (A) Example images of the spatial distributions observed
for CD45 isoforms before and after PMA stimulation. The RCA amplicons
of CD45RO, CD45RA, and CD45RB are marked with green, blue, and red
circles, respectively, and the outlines of cell are marked by a gray
dotted line. Scale bars: 10 μm. (B) Single-cell profiles for
distance to cell nucleus (yellow dots) and the edge (black dots). N is the number of cell dots per plot. (C) Quantification
of average distances for three CD45 isoforms to the cell nucleus or
the edge before and after PMA stimulation.
Analysis of Expression Covariation among Different Isoforms
We next examined the expression covariation among different isoforms
that might arise from differential activity of immune activation.
Pairwise correlation coefficients for the RNA species were calculated
as shown in Figure A. Upon T-cell activation, the correlation coefficients of CD45RO
and CD45RB presented a pronounced increase whereas pair correlations
of CD45RO and CD45RA, and CD45RA and CD45RB, presented a modest change.
The single-cell correlation coefficients of CD45RO and CD45RB vary
from 0.281 to 0.740 upon stimulus (Figure B). Covariation across single cells between
splicing isoforms of a gene upon T-cell activation would represent
a potential regulatory interaction between alternative splicing and
T-cell activation. Isoforms CD45RO and CD45RB encode the transcript
mRNA which both skip exon 4. T-cell activation would lead to the up-regulation
of some regulators, such as HnRNP LL, SR proteins, or other RNA binding
protein which can alter alternative splice sites by binding to different
exons.[43] We suggest that the expression
of some RNA binding proteins involved in exon 4 skipping was up-regulated
with T-cell activation, which thus may lead to the simultaneous up-regulation
of CD45RO and CD45RB, resulting in the increase of their correlation
coefficients.[44] The normal functioning
of the immune system may be associated with the complex interplay
of regulators that mediate the appropriate splicing of CD45 exon 4.
Analysis of covariations in the expression levels of different isoforms
on the single-cell level could reveal coregulated isoforms and help
to elucidate splicing regulatory circuits in the T-cell maturation.
Figure 4
Cell-to-cell
pairwise correlations for the RNA species upon T-cell
activation. (A) Matrix of the pairwise correlation coefficients of
the cell-to-cell variation in expression for CD45 splice isoforms
before and after PMA stimulation. (B) CD45RO–CD45RB expression
correlation analysis in a single cell before (blue) and after (red)
PMA stimulation. (P value <0.001). The Pearson
correlation coefficient changes from 0.281 to 0.740 upon PMA stimulation
(P value <0.001).
Cell-to-cell
pairwise correlations for the RNA species upon T-cell
activation. (A) Matrix of the pairwise correlation coefficients of
the cell-to-cell variation in expression for CD45 splice isoforms
before and after PMA stimulation. (B) CD45RO–CD45RB expression
correlation analysis in a single cell before (blue) and after (red)
PMA stimulation. (P value <0.001). The Pearson
correlation coefficient changes from 0.281 to 0.740 upon PMA stimulation
(P value <0.001).
Conclusions
In summary, we develop a facile and high
base-resolution strategy
for multiplexed profiling of RNA splicing variants in single cells.
This imaging method, termed SpliceRCA, has several advantages: (1)
Benefitting from recognition of the target splice-junction sequence,
the read length of the imaging method can be shortened to ∼30
nt, which resolved the limitations in read length of current imaging
technologies. (2) The ligation process and one-target-one-amplicon
amplification involved in SpliceRCA confer high specificity with single-base
precision and quantification of splice isoforms with single-molecule
resolution. (3) Insight into RNA splicing variants on expression level,
spatial distribution, and cell heterogeneity was simultaneously acquired
by SplicRCA. The method was used to image splice variants of the essential
regulatory immune gene (CD45) upon T-cell activation. We have identified
that CD45RO shows a distal nuclear spatial distribution, which is
different from other isoforms. After T-cell activation, CD45RB presented
an obvious change in spatial distribution, which appeared enriched
near the cell periphery. Furthermore, it is found that CD45RO and
CD45RB coregulate upon PMA stimulation, which has not been observed
by previous cell-population-based RNA quantification methods. Our
study demonstrates the promise of single-cell alternative splicing
analysis in deciphering rich diversity in functional communities and
understanding the essential features of immune responses for assisting
clinical monitoring.
Materials and Methods
Oligonucleotide Sequences
The DNA sequences (Tables S1 and S2 in
the Supporting Information)
were purchased from Shanghai Sangon Biological Engineering Technology
& Services Co., Ltd. (Shanghai, China). The sequences modified
with Alexa488, Alexa555, and Cy5 were purchased from Thermo Fisher
Scientific (Beijing, China), and were purified by HPLC. RNA sequences
were created by in vitro T7 transcription reactions
with ordering single-stranded DNA used as templates.
Cell Culture
and Stimulations
The Jurkat cells were
maintained in a standard RPMI 1640 medium supplemented with 15% fetal
bovine serum, 100 units mL–1 of penicillin and 100
μg mL–1 of streptomycin. The cells were grown
at 37 °C, 5% CO2, and 95% air humidity. For stimulations,
cells were diluted to 3 × 106 cells mL–1 and
incubated in medium with PMA plus ionomycin (100 ng mL–1 and 1 μM, respectively).
Cell Fixation and Permeabilization
Cells were fixed
at a proper density on a 22 mm × 22 mm gelatin-coated cover glass
(VWR, Radnor, PA) enclosed by a PDMS with a chamber (5 mm in diameter)
by being maintained in 4% (w/v) paraformaldehyde in phosphate buffered
saline (PBS) for 15 min at room temperature (20–25 °C),
washed twice with 1× DEPC-treated PBS (DEPC-PBS). Then, the cells
were permeabilized for 5 min with 0.5% v/v Triton-X100 in 1×
PBS at room temperature, and washed twice with DEPC-PBS.
In
Situ Visualization of Splicing Variants
by SpliceRCA
The hybridization of the exon-junction padlock
probe with the target mRNA splicing variants was conducted in a volume
of 20 μL containing 2 μL 20× saline-sodium citrate
buffer (SSC) (Ambion, AM9763), 2 μL of each phosphorylated padlock
probe (10 μM), 1 μL of DTT (100 mM), 2 μL of yeast
tRNA (10 mg mL–1), and 0.5 μL of RiboLock
RNase inhibitor (40 U μL–1) overnight at 37
°C. The sample was then washed twice using PBS-T (DEPC-PBS with
0.05% Tween-20) for 3 min at room temperature. The ligation reaction
was carried out in 10 μL of circularization reaction mixture
[1 μL 10× ligase reaction buffer, 1 μL of 20×
SSC, 1 μL of T4 DNA ligase (5 U μL–1), 0.25 μL of RiboLock RNase inhibitor (40 U μL–1), 6.75 μL of RNase-free water] at 37 °C for 2 h. After
ligation, a 20 μL mixture containing 2 μL of 20×
SSC, 1 μL of each primer (4 μM), 2 μL of formamide,
1 μL of DTT (100 mM), 13.5 μL of RNase-free water, and
0.5 μL of RiboLock RNase inhibitor (40 U μL–1) were added to the sample and incubated for 60 min at 37 °C,
following a wash using PBS-T. RCA was then conducted with a 10 μL
mixture containing 1 μL of 10× phi29 DNA polymerase reaction
buffer, 0.5 μL of phi29 DNA polymerase, 3 μL of dNTPs
(10 mM for each of dATP, dGTP, dCTP, and dTTP), 5.25 μL of RNase-free
water, and 0.25 μL of RiboLock RNase inhibitor (40 U μL–1) for 120 min at 37 °C. The incubation was followed
by a wash in PBS-T. Then, the hybridization of amplicons with detection
probes was conducted in a 20 μL mixture of 100 nM fluorophore-labeled
detection probes, 2× SSC, 15% formamide, and 10 ng μL–1 salmon sperm DNA for 30 min at 37 °C, following
two washes using PBS-T. After mounting with Fluoromount-G [containing
4′,6-diamidino-2-phenylindole (DAPI) to counterstain the cell
nuclei], the slides were ready for imaging.[45,46]
Image Acquisition and Analysis
Fluorescence imaging
was performed using a Zeiss LSM 710 META confocal microscope (Zeiss).
The cellular images were acquired with a 63× (oil) objective.
Argon/2 (488 nm) was used as excitation source for Alexa488-labeled
probe, and a 500–535 nm bandpass filter was used for fluorescence
detection. The Alexa555 dye was excited with HeNe1 (561 nm) laser
and detected with a 580–620 nm bandpass filter. The Cy5 dye
was excited with HeNe2 (633 nm) laser and detected with a 650–750
nm bandpass filter. The DAPI dye was excited with a 405 nm diode laser
and detected with a 430–550 nm bandpass filter. Z-stacks were collected at step size of 0.4 μm for 20 slices
to image the entire cell. The images were processed by ImageJ version
1.46r software. The outline of cells was determined by the bright
field images. The amplicons which had bright fluorescent signals were
distinguished from background by setting the threshold value.[16,47] For determination of the copy number of amplicons per cell, the
number of bright pixels was counted by particle analysis in ImageJ
software.[48]
Total RNA from Jurkat cells was harvested
by using TransZol according to the included protocol. The reverse
transcription was performed using TransScript one-step gDNA removal
and cDNA synthesis kit. In brief, a total volume of 20 μL of
mixture containing 10 μL of 2× TS reaction mixture, 2 μL
of the total RNA (from 50 ng to 5 μg), 1 μL of RT primer
(0.5 μg μL–1), 1 μL of TransScript
RT/RI enzyme mix, 1 μL of gDNA remover, and 5 μL of RNase-free
water was incubated at 42 °C for 15 min, following heat inactivation
of reverse transcriptase at 85 °C for 5 s. The produced cDNA
samples were then stored at −80 °C for future use. Upstream
primers used for the reverse transcription reaction are listed in Table S1.For qPCR analysis, following
the manufacturer’s instructions on a Bio-Rad C1000TM (Bio-Rad)
instrument, the 20 μL reaction solution contained 10 μL
of 2× SYBR Select master mix, 2 μL of cDNA sample, 2 μL
of forward primer (5 μM), 2 μL of reverse primer (5 μM),
and 4 μL of RNase-free water. The qPCR was done by staying at
50 °C for 2 min for the hot start, annealing at 95 °C for
2 min, then followed by 40 cycles of 15 s at 95 °C, 60 s at 60
°C, and 5 min at 60 °C. Ct values were converted into absolute
GAPDH copy numbers using a standard curve from a control RNA (humanGAPDH mRNA in RevertAid First Strand cDNA synthesis kit). The experiments
were performed in triplicate. The copy numbers of target mRNA splicing
variants CD45RO, CD45RA, and CD45RB were evaluated by referring to
the expression of GAPDH mRNA using the 2−ΔΔCt method.
Calculations of mRNA copy numbers were based on the number of counted
cells at harvest. No unexpected or unusually high safety hazards were
encountered in the experiments.
Authors: Zuopeng Wu; Xinying Jia; Laura de la Cruz; Xun-Cheng Su; Bruz Marzolf; Pamela Troisch; Daniel Zak; Adam Hamilton; Belinda Whittle; Di Yu; Daniel Sheahan; Edward Bertram; Alan Aderem; Gottfried Otting; Christopher C Goodnow; Gerard F Hoyne Journal: Immunity Date: 2008-12-19 Impact factor: 31.745