Post-transcriptional modifications of RNA represent an understudied layer of translation regulation that has recently been linked to central nervous system plasticity. Here, sample preparation and liquid chromatography-tandem mass spectrometry approach is described for simultaneous characterization of numerous RNA modifications in single neurons.
Post-transcriptional modifications (PTMs) of RNA represent an understudied mechanism involved in the regulation of translation in the central nervous system (CNS). Recent evidence has linked specific neuronal RNA modifications to learning and memory paradigms. Unfortunately, conventional methods for the detection of these epitranscriptomic features are only capable of characterizing highly abundant RNA modifications in bulk tissues, precluding the assessment of unique PTM profiles that may exist for individual neurons within the activated behavioral circuits. In this protocol, an approach is described—single-neuron RNA modification analysis by mass spectrometry (SNRMA-MS)—to simultaneously detect and quantify numerous modified ribonucleosides in single neurons. The approach is validated using individual neurons of the marine mollusk, Aplysia californica, beginning with surgical isolation and enzymatic treatment of major CNS ganglia to expose neuron cell bodies, followed by manual single-neuron isolation using sharp needles and a micropipette. Next, mechanical and thermal treatment of the sample in a small volume of buffer is done to liberate RNA from an individual cell for subsequent RNA digestion. Modified nucleosides are then identified and quantified using an optimized liquid chromatography-mass spectrometry method. SNRMA-MS is employed to establish RNA modification patterns for single, identified neurons from A. californica that have known morphologies and functions. Examples of qualitative and quantitative SNRMA-MS are presented that highlight the heterogeneous distribution of RNA modifications across individual neurons in neuronal networks.
Modifications in the canonical nucleosides of RNA have been increasingly recognized for their myriad roles in the regulation of protein translation. Over 150 unique RNA modifications have been reported to date that range in complexity from methylation, heteroatom addition, to conjugation with cellular metabolites1,2. This expanded RNA alphabet, also known as the epitranscriptome, is generated by enzymatic writers and is responsible for altering the stability3, folding4, and translation efficiency5,6 of cellular RNAs. Select RNA modifications may also be reversed via enzymatic erasers7,8, whereas others are appended to RNAs substoichiometrically9,10, rendering a complex landscape of modified and unmodified RNA sequences in biological systems.
The importance of RNA modifications in biological function is indicated by the unequal distribution of modifications across different organs and tissues, including subregions of the central nervous system (CNS)11. This chemical heterogeneity has been linked to CNS development12, stress response13, and activity-dependent plasticity14. The CNS subregions further comprise heterogeneous populations of cells in which individual cells exhibit distinct chemical profiles15,16,17,18. Even single cells of the same type may display unique transcriptomes, in part due to tissue microenvironments19 and stochastic gene expression20. However, while characterization of the single-cell transcriptome is somewhat routine, there are no analogous methods for establishing single-cell epitranscriptomes for multiple RNA modifications. New approaches that are capable of profiling the distribution of RNA modifications in individual cells are needed for comprehensive analysis of the cellular heterogeneity and regulatory influence of post-transcriptional modifications (PTMs) in the CNS and other biological systems.
Simultaneous measurement of numerous RNA modifications in bulk cells/tissues is readily accomplished using liquid chromatography-tandem mass spectrometry (LC-MS/MS) techniques. For LC-MS/MS analysis of modified ribonucleosides, RNA is extracted from cells (typically 103-106 cells), purified by precipitation and resuspension, and subsequently digested into nucleosides. The sample mixture consisting of canonical and modified nucleosides is then injected into the LC-MS system for analyte separation and detection, leading to determination of the full complement of RNA modifications in an organism21,22,23. LC-MS/MS was recently used to determine 26 RNA modifications in the CNS of the neurobiological model, Aplysia californica (A. californica). The abundances of some of these epitranscriptomic marks exhibited time- and region-dependent dynamics that correlated with behavioral changes in the animal24. However, it was only possible to detect RNA modifications in bulk samples containing >103 cells due to the limited sensitivity of the method. These larger samples likely concealed unique and functionally important RNA modification profiles of individual cells with population averages. Although careful control of sample preparation conditions has improved the detection limits for RNA modifications in small-volume samples25,26,27,28, there remains a need for analytical methods that can detect and quantify multiple modified ribonucleosides in single cells.
This protocol introduces single-neuron RNA modification analysis by mass spectrometry (SNRMA-MS), which permits the detection of over a dozen RNA modifications in single neurons from the CNS of A. californica29. The approach consists of surgical isolation of single, identified cells from major CNS ganglia followed by an optimized sample preparation workflow involving mechanical cell lysis, RNA denaturation, and enzymatic hydrolysis in an MS-compatible buffer. Identification and quantification of post-transcriptionally modified nucleosides is then accomplished using LC-MS/MS. SNRMA-MS fulfills an unmet need in the field of RNA modification analysis by facilitating the acquisition of post-transcriptional RNA modification profiles for single neurons and has potential for future application to other cell types.
1. Preparation of materials and solutions
2. Isolation of single neuron
3. Cell lysis and RNA digestion for SNRMA-MS
4. Liquid chromatography-tandem mass spectrometry
NOTE: Analyze the single-neuron digests and authentic modified nucleoside standards using an LC-MS/MS system equipped with an electrospray ionization source and six-port diverter valve.
5. Data analysis
SNRMA-MS involves the manual isolation of identified neurons into small sample volumes for lysis, digestion, and LC-MS/MS analysis (Figure 1A). This workflow routinely detected over a dozen RNA modifications in single neurons from the CNS of A. californica (Figure 1B), representing a coverage of nearly half of the known epitranscriptome of this animal24 in a single cell. For example, subjecting the LPl1 neuron (~500 µm diameter) to SNRMA-MS resulted in the detection of 15 ± 1 RNA modifications (n = 3). Modified nucleosides were positively identified on the basis of three attributes: LC retention properties, exact mass, and MS2 fragmentation profiles compared to values provided in the database2. Table 1 shows a list of all the RNA modifications detected in the LPl1 neuron. The high-resolution quadrupole time-of-flight mass spectrometer used for these experiments enabled a mass accuracy of 4 ppm as well as detection of the characteristic MS2 fragment ion at m/z 164 for the modified nucleoside N6,6-dimethyladenosine (m66A, Figure 1B). Combined with the LC separation data, which agrees with findings deposited in the database (m66A elutes after N6-methyladenosine (m6A)), the SNRMA-MS approach demonstrated correct assignment of modified nucleoside identities.
The SNRMA-MS platform can be leveraged for establishing RNA modification profiles of single neurons and investigating their relationship to bulk tissues. R2 neurons are large (~500 µm in diameter) cholinergic cells that reside in the abdominal ganglion. SNRMA-MS was used to analyze RNA modifications in R2 neurons as well as the surrounding bulk abdominal ganglion (Figure 2A). Normalized peak areas for RNA modifications in each neuron/ganglion sample (n = 7) were used as inputs for PCA, revealing that R2 neurons exhibit distinct RNA modification profiles compared to the ganglia in which they reside (Figure 2B). This is evidenced by data points for R2 neurons and abdominal ganglia occupying different regions of the PCA score plot. Further support for unique modified nucleoside patterns was obtained from a separate cohort of animals (n = 7) in which pairwise comparisons were performed for 13 RNA modifications that were commonly detected in both the single neurons and bulk tissue (Figure 2C). Two modified nucleosides, pseudouridine (Ψ) and 2'-O-methylguanosine (Gm), were at significantly higher abundance in the abdominal ganglia compared to R2 neurons. When viewing a subset of the RNA modifications with R2 neuron-ganglion pairs indicated, all of the abdominal ganglia exhibited higher levels of Ψ and Gm, and all but one of the R2 neurons had higher abundances of 2'-O-methyladenosine (Am) than their corresponding ganglion (Figure 2D). Overall, the SNRMA-MS results reveal for the first time that RNA modification profiles of single cells can diverge from bulk cells in the same tissue.
Using SNRMA-MS to investigate the model animal A. californica provides a unique opportunity to characterize RNA modification profiles in identified, functionally distinct neurons. RNA modifications were evaluated by SNRMA-MS in four identified cells: R2 and LPl1 (homologous, cholinergic cells involved in defensive mucous release)32, MCCs (serotonergic modulatory neurons involved in feeding)33, and B2 cells (peptidergic neurons involved in gut motility)34. PCA of six RNA modifications in these identified neurons, either isolated immediately after enzymatic treatment or cultured in a ganglion preparation for 48 h, demonstrated the stability and dynamics of single cell epitranscriptomes. RNA modifications in functionally different cells formed unique clusters in the score plot while homologous R2/LPl1 neurons co-clustered (Figure 3A). The loading plot shows that differences were primarily driven by the abundance of positional isomers of methyladenosine, including 2'-O-methyladenosine (Am) and N1-methyladenosine (m1A) (Figure 3B). In the same analysis, a comparison of freshly isolated cells and cells cultured in situ (i.e., in their respective ganglia) for 48 h was performed. As shown in the PCA score plot, functionally different cells remained distinguishable by their RNA modification profiles.
Quantitative SNRMA-MS can be used for determining absolute amounts of select RNA modifications for which authentic standards are available. External calibration curves were generated for m1A, Ψ, 2'-O-methylcytidine (Cm), Am, and m6A, and the amount of each modified nucleoside in the MCC and R2/LPl1 cell pairs was determined by interpolation (Figure 4A–E). The intracellular quantities of m1A and Ψ in two pairs of symmetrical MCCs appeared to be similar, while larger differences in the amounts of these modifications were observed in three pairs of R2/LPl1 cells. In order to account for differences due to the physical size of the cells studied, RNA modification quantities were normalized by cell volumes calculated from optical measurement of cell diameters to yield intracellular concentrations of modified nucleosides. Significant differences in intracellular concentrations of Cm and Am were observed between MCCs and R2/LPl1 neurons. Overall, SNRMA-MS enables both qualitative and quantitative profiling of RNA modifications in single neurons.
Figure 1: SNRMA-MS workflow and detection of multiple RNA modifications in single neurons by LC-MS/MS. (A) Photographs of desheathed buccal ganglion and single neuron isolation into a sample tube. Scale bar = 200 µm, arrows indicate identified B1 cell. A diagram of the sample preparation procedure for LC-MS/MS analysis is also shown. (B) Overlaid EICs for RNA modifications in a single LPl1 neuron, with inset showing MS1 and MS2 spectra for N6, N6-dimethyladenosine (m66A). See Table 1 for m/z values used to generate EICs for modified nucleosides. This figure has been modified from29. Please click here to view a larger version of this figure.
Figure 2: SNRMA-MS distinguishes RNA modification profiles of single neurons and bulk tissue. (A) Schematic of A. californica CNS and workflow for analyzing R2 neurons and the surrounding abdominal ganglion. The photograph shows the abdominal ganglion and R2 neuron (injected with Fast Green dye for visibility). (B) Relative peak areas from 13 RNA modifications were used to generate the PCA score (top) and loading (bottom) plots. (C) Pairwise comparison of RNA modifications in the abdominal ganglion and R2 neuron. Error bars represent ±1 standard deviation (SD), *p < 0.05, ***p < 5 x 10−4, paired t-test with Bonferroni−Holm correction. (D) Comparison of select RNA modifications from panel C in which R2-abdominal ganglion pairs for each animal are shown with droplines. This figure has been modified from29. Please click here to view a larger version of this figure.
Figure 3: Profiling RNA modifications in identified, functionally different neurons from the A. californica CNS. (A) PCA score plot for MCCs, and B2, R2, and LPl1 cells that were either freshly isolated or isolated following in situ culture for 48 h (denoted by cMCC, cB2, cR2, cLPl1), and (B) loading plots for six RNA modifications commonly detected in these cells. This figure has been modified from29. Please click here to view a larger version of this figure.
Figure 4: Quantitative SNRMA-MS in A. californica neurons. Quantitative SNRMA-MS provides absolute amounts and intracellular concentrations for several modified nucleosides in single, identified A. californica neurons. Photographs of (A) left and right MCCs (LMCC and RMCC, respectively) in the cerebral ganglion and (B) R2 in the abdominal ganglion and LPl1 in the pleural ganglion. Cells are circled to enhance visibility. Linear calibration plots for (C) m1A and (D) Ψ (triangles) used for interpolation of modified nucleoside quantities in single cells (colored dots). Cell pairs from each animal are labeled 1-3. (E) Intracellular concentrations of five modified nucleosides in MCC and R2/LPl1 cell pairs. Thick lines connect cell pairs. *p < 0.05, **p < 0.005, paired t-test with Bonferroni−Holm correction, n = 2 animals (four cells total) for MCCs and n = 3 animals (six cells total) for R2/LPl1. This figure has been modified from29. Please click here to view a larger version of this figure.
RNA modification | Abbreviation | Elution order (C18) | m/z for EIC | m/z for MS2 |
dihydrouridine | D | 1 | 247.09 | 115 |
pseudouridine | Y | 2 | 245.08 | 209/179/155 |
3-methylcytidine | m3C | 3 | 258.11 | 126 |
N1-methyladenosine | m1A | 4 | 282.12 | 150 |
5-methylcytidine | m5C | 5 | 258.11 | 126 |
N7-methylguanosine | m7G | 6 | 298.12 | 166 |
2'-O-methylcytidine | Cm | 7 | 258.11 | 112 |
inosine | I6A | 8 | 269.09 | 137 |
2'-O-methylguanosine | Gm | 9 | 298.12 | 152 |
N2-methylguanosine | m2G | 10 | 298.12 | 166 |
N2,N2,N7-trimethylguanosine | m227G | 11 | 326.15 | 194 |
N2,N2,-dimethylguanosine | m22G | 12 | 312.13 | 180 |
2'-O-methyladenosine | Am | 13 | 282.12 | 136 |
N6-methyladenosine | m6A | 14 | 282.12 | 150 |
N6,N6-dimethyladenosine | m66A | 15 | 296.14 | 164 |
N6-isopentenyladenosine | i6A | 16 | 336.17 | 204/136/148 |
Table 1: RNA modifications detected in single neurons from A. californica. Attributes for characterization of modified nucleosides are provided including LC retention order, m/z for generating EICs, and corresponding CID fragments.
SNRMA-MS leverages an optimized sample preparation approach, resulting in a small, MS-compatible sample volume that can be delivered to the LC-MS platform. The initial enzymatic pre-treatment of CNS ganglia dictates both the ease with which they can be desheathed and the durability of single neurons during isolation. The cerebral ganglion often requires extended enzymatic treatment due to its relatively thick sheath compared to the buccal, pleural, and abdominal ganglia. Individual researchers performing the single neuron isolations may have differing preferences for how durable the sheath should be when using microscissors and fine forceps. However, it is important that the ganglia are not over-digested, as this may lead to a loss in their structural integrity, loss of positional information that is critical for identifying target cells, and/or cell lysis. Following isolation from the ganglion, it is important to ensure that a minimal volume of isolation medium is aspirated when transferring the neuron to the sample tube. The ASW medium contains a high concentration of salts that may interfere with digestive enzymes during RNA hydrolysis and will also dilute the sample.
During the mechanical lysis step, it is common for large neurons (>250 µm diameter) to rupture upon repeatedly passing through a micropipette. Smaller neurons may require additional attention to ensure cell lysis, which typically involves pressing a glass capillary on the cell. In this instance, it is possible that a partial volume of the sample buffer will be drawn into the capillary because of capillary forces. This volume can be delivered back into the sample tube by applying pressure to the end of the glass capillary to ensure that no sample is lost.
The best results are obtained by including a heating step prior to adding enzymes for RNA digestion. This is likely because heating at 95 °C denatures RNA secondary structures that may impede the activity of RNases35 and diminish the quantity of nucleosides released from RNA biopolymers. Control experiments using samples spiked with stable isotope-labeled methionine were previously performed to investigate whether heat-induced RNA modification artifacts were the cause of the increased peak areas observed for RNA modifications relative to a no-heat SNRMA-MS protocol29. No such labeling was observed, indicating that the improved signals for RNA modifications using the optimized SNRMA-MS method were due to superior digestion of RNA.
Conventional methods for isolating total RNA from cells involve liquid-liquid extraction (LLE) with phenol-chloroform and subsequent RNA precipitation, washing, and resuspension. These approaches have proven useful for reverse transcription polymerase chain reaction experiments where the expression of select genes can be readily monitored in identified A. californica neurons36,37. However, LLE methods cannot recover sufficient RNA for the detection of modified ribonucleosides by LC-MS, whereas the SNRMA-MS method described herein enables detection of numerous RNA modifications29. In order to assess modification profiles of specific RNA types (e.g., rRNA, tRNA, mRNA) in bulk tissue/cell samples, anion-exchange solid phase extraction24, hybridization probe-based enrichment38, and chromatographic fractionation39 approaches have been applied, but similar methods are not yet available for single-cell RNA purification. The development of RNA fractionation approaches that are capable of isolating specific RNA types from single cells would provide additional insight toward the function of RNA modifications.
SNRMA-MS revealed previously uncharacterized heterogeneity in the RNA modification landscape of single neurons in A. californica and it is conceivable that similarly distinct PTM profiles exist across mammalian cells. Because mammalian cells are relatively small compared to the large A. californica neurons analyzed in this protocol, improvements in the handling of small sample volumes are needed to facilitate lower detection limits. Though currently SNRMA-MS is limited to volumes of ~5 µL, it is expected that substantial improvements can be attained by incorporating microfluidic liquid handling devices into the workflow. Furthermore, the implementation of automated or semi-automated cell isolations would increase sample throughput and permit single-cell RNA modification analysis of larger cell populations. By interfacing with nanoflow LC separations27, the characterization of epitranscriptomic marks in even smaller cells will be achievable.
The authors have nothing to disclose.
This work was funded by the National Institute on Drug Abuse under Award no. P30DA018310 and the National Human Genome Research Institute under Award no. RM1HG010023. K.D.C. acknowledges support from a Beckman Institute Postdoctoral Fellowship. The content is solely the responsibility of the authors and does not necessarily represent the official views of the funding agencies.
Animals | |||
Aplysia californica | National Resource for Aplysia (Miami, FL) | 150–250 g (adult) | |
Benchtop equipment | |||
Glass capillary puller | Sutter | P-97 | Equipped with online degasser, autosampler, and thermostatted column compartment |
Milli-Q water purification system | Millipore | Equipped with ESI source | |
Minicentrifuge for PCR tubes | LW Scientific | ZS-1 | |
Optical Microscope | Zeiss | Stemi 2000C | |
Thermocycler | Techne | EW-93945-01 | |
HPLC column and consumables | |||
Acclaim RSLC 120 C18 column | Thermo Scientific | 71399 | |
Autosampler vials | Thermo Scientific | C4011-13 | |
LC and MS Instrumentation and Software | |||
DataAnalysis 4.4 software | Bruker | ||
Dionex Ultimate 3000 nanoLC | Thermo Scientific | ||
Impact HD UHR QqTOF mass spectrometer | Bruker | ||
RStudio | RStudio | ||
Microdissection tools | |||
Microscissors extra fine vannas 3.5” | Roboz | RS-5640 | |
Tungsten needles | Roboz | RS-6065 | |
Reagents/Materials | |||
4-(2-hydroxyethyl)-1-piperazineethane-sulfonic acid (HEPES) | Fisher Scientific | H3375 | |
Alkaline phosphatase | Worthington Biochemical Corp. | LS004081 | |
Ammonium acetate | Honeywell | 17836 | |
benzonase (endonuclease from S. marcescens) | EMD Millipore | 70746-4 | |
bovine serum albumin | Sigma-Aldrich | A2153 | |
calcium chloride | Sigma-Aldrich | C4901 | |
gentamycin sulfate | Fisher Scientific | G1264 | |
Magnesium chloride | Sigma-Aldrich | M9272 | |
magnesium sulfate | Sigma-Aldrich | 208094 | |
nucleosides test mix | Sigma-Aldrich | 47310-U | |
penicillin G | Sigma-Aldrich | P7794 | |
pentostatin | Sigma-Aldrich | SML0508 | |
phosphodiesterase I | Worthington Biochemical Corp. | LS003926 | |
protease type XIV from Streptomyces griseus | Sigma-Aldrich | P5147 | |
sodium chloride | Sigma-Aldrich | S9888 | |
Standard glass capillaries | A-M Systems | 626000 | 1 mm o.d., 0.5 mm i.d., 4 in |
streptomycin sulfate | Sigma-Aldrich | S9137 |