Isolation of Specific Neurons from C. elegans Larvae for Gene Expression Profiling

Abstract

Background

The simple and well-described structure of the C. elegans nervous system offers an unprecedented opportunity to identify the genetic programs that define the connectivity and function of individual neurons and their circuits. A correspondingly precise gene expression map of C. elegans neurons would facilitate the application of genetic methods toward this goal. Here we describe a powerful new approach, SeqCeL (RNA-Seq of C. elegans cells) for producing gene expression profiles of specific larval C. elegans neurons.

Methods and Results

We have exploited available GFP reporter lines for FACS isolation of specific larval C. elegans neurons for RNA-Seq analysis. Our analysis showed that diverse classes of neurons are accessible to this approach. To demonstrate the applicability of this strategy to rare neuron types, we generated RNA-Seq profiles of the NSM serotonergic neurons that occur as a single bilateral pair of cells in the C. elegans pharynx. These data detected >1,000 NSM enriched transcripts, including the majority of previously known NSM-expressed genes.

Significance

This work offers a simple and robust protocol for expression profiling studies of post-embryonic C. elegans neurons and thus provides an important new method for identifying candidate genes for key roles in neuron-specific development and function.

About the Authors:

W. Clay Spencer

Current address: Department of Neurosciences, Case Western Reserve University, Cleveland, Ohio, United States of America

Affiliation: Department of Cell and Developmental Biology, Vanderbilt University, Nashville, Tennessee, United States of America

Rebecca McWhirter

Affiliation: Department of Cell and Developmental Biology, Vanderbilt University, Nashville, Tennessee, United States of America

Tyne Miller

Affiliation: Program in Neuroscience, Vanderbilt University, Nashville, Tennessee, United States of America

Pnina Strasbourger

Affiliation: Department of Genome Sciences, University of Washington School of Medicine, Seattle, Washington, United States of America

Owen Thompson

Affiliation: Department of Genome Sciences, University of Washington School of Medicine, Seattle, Washington, United States of America

LaDeana W. Hillier

Affiliation: Department of Genome Sciences, University of Washington School of Medicine, Seattle, Washington, United States of America

Robert H. Waterston

Affiliation: Department of Genome Sciences, University of Washington School of Medicine, Seattle, Washington, United States of America

David M. Miller III

* E-mail: [email protected]

Affiliations Department of Cell and Developmental Biology, Vanderbilt University, Nashville, Tennessee, United States of America, Program in Neuroscience, Vanderbilt University, Nashville, Tennessee, United States of America

Introduction

With its well-defined, compact nervous system and facile genetics, C. elegans is widely exploited for studies of neural development and function. The morphology and connectivity of the C. elegans nervous system is catalogued in comprehensive wiring diagrams [1], [2], [3] that facilitate functional analysis [4], [5]. The unrivaled precision of this nervous system model is complemented by the complete sequence of the C. elegans genome and its extensive annotation derived from direct RNA-Seq analysis [6], [7]. Cell-specific profiling experiments have identified subsets of genes that are highly expressed in particular neurons or that may be regulated by transcription factors with key roles in neuron-specific differentiation [8], [9], [10], [11], [12], [13], [14]. Expression profiling substantially narrows the list of candidate genes for tests of function in vivo and therefore offers an efficient strategy for identifying critical determinants of neuron differentiation and activity. Fluorescent reporter transgenic lines have been generated for thousands of individual C. elegans genes and their expression in specific neurons has been documented [15]. For embryos, cells can be dissociated and neurons that are marked with GFP reporters can be readily isolated by Fluorescence-Activated-Cell-Sorting (FACS) [9], [16], [17]. For example, profiling data generated for BAG sensory neurons by this approach led to the identification of a guanylate cyclase receptor that detects CO₂ and a conserved ETS transcription factor that regulates BAG neuron fate [18], [19]. For post-embryonic animals, the mRNA-tagging method has been extensively utilized to profile larval and adult neurons [10], [11], [13], [20], [21], [22], [23], [24]. In this strategy, an epitope-tagged mRNA binding protein is selectively expressed in target neurons for immuno-precipitation of cell-specific transcripts [20]. Although useful, this approach requires custom-built transgenic lines and the biochemical preparation may include significant background RNA that limits specificity [21].

The recent development of a simple protocol by Jeff Kuhn’s laboratory for generating dissociated populations of viable cells from C. elegans larvae offers the potential alternative to the mRNA tagging method of profiling postembryonic cells isolated by FACS. However, neurons were reportedly under-represented in these preparations [25]. This apparent limitation would restrict ready access to specific types of neurons, the majority of which are rare since they are defined by either a single cell or by bilateral pairs of similar neurons in each animal [1]. But our investigations demonstrate that larval neurons are readily released by the Kuhn cell dissociation protocol. The discrepancy is explained by a requirement in the Kuhn method for rapid adherence of cells to the culture dish. In the first instance, we used FACS to show that GFP-labeled cells comprise at least 30% of viable cells obtained from a transgenic line in which all neurons are marked with a GFP reporter. Specific classes of sensory and motor neurons were also isolated by FACS at a fraction predicted by their relative abundance in vivo. To test the applicability of this approach to profiling a specific neuron type, we used FACS to isolate NSM neurons from L1 stage larvae. The NSM neurons consist of two morphologically similar serotonergic neurons located in the pharynx [26], [27]. RNA-Seq profiles confirmed expression of known NSM-specific genes and also detected >1,000 additional transcripts that are enriched in NSM relative to other cell types. We therefore conclude that C. elegans larval neurons are readily accessible to isolation by FACS for gene expression profiling and predict that this simple approach will be highly useful for studies of neural development and function in this model organism.

Results

Viable C. elegans larval neurons can be readily isolated by FACS (Fluorescence Activated Cell Sorting)

In previous expression profiling studies of embryonic cell-types, we determined that the approximate fraction of a specific cell-type in a FACS (Fluorescence-Activated-Cell-Sorting) profile of all viable cells is correlated with the relative abundance of these target cells in the intact embryo [9], [11], [28]. To test this prediction for larval neurons, we generated primary cultures from L1 larvae labeled with the pan-neural marker, F25B3.3::GFP (Fig. 1A, B) [29]. Morphologically distinct, GFP-labeled neurons were well represented in these cultures within 36 hours after plating (Fig. 1C). A FACS profile of the initial dissociated preparation of L1 larval cells detected ∼32% of cells as GFP-labeled (Fig. 1D, E), which compares favorably to the fraction of neurons (222/558 = 40%) in newly hatched L1 larva. Our results stand in contrast to an earlier report that neurons comprise less than 10% of cultured cells derived from L1 larvae. In this previous experiment, the culture medium was replaced within 24 hours of plating and before GFP-labeled cells were counted [25]. We have observed that neurons are only loosely attached at this stage (data not shown) and therefore largely discarded when media are changed after 1 day in culture. In contrast, body muscle cells adhere more quickly, which likely accounts for the previous observation that 81% of cultured L1 cells were derived from muscle, although this fraction comprises only 15% (81/558) of newly-hatched L1 cells [25]. We found that viable neurons from F25B3.3::GFP could be isolated by FACS to ∼90% purity by directly counting GFP-labeled cells immediately after sorting (Fig. 1F, G).

[Figure omitted. See PDF.]

Figure 1. Isolation of C. elegans neurons from larval animals by Fluorescence Activated Cell Sorting (FACS).

(A) Expression of the pan-neural GFP marker, evIs111, in an L2 stage larval animal. Anterior is to left. (B) Procedure for generating dissociated suspensions of larval cells for primary culture and for isolation of GFP-marked cells (green) by FACS. (C) Primary culture of L1-stage cells 36 hr after plating. Arrows point to elongated processes extending from GFP-labeled neurons. Scatterplot of FACS profile for cells dissociated from the wild-type (N2) reference strain (D) and from the pan-neural evIs111 line (E). Propidium Iodide (PI) marks dead cells (blue). GFP-labeled cells (green) were isolated by FACS (outlined with box). The majority (∼90%) of FACS-derived cells (F) are marked with GFP (G) in primary cultures examined within 2 hr of FACS isolation. Note that neurons are loosely attached which likely accounts for the displacement of DIC vs GFP images of individual cells in these micrographs. Scale bars = 10 micron.

https://doi.org/10.1371/journal.pone.0112102.g001

We generated dissociated preparations of larval cells labeled with fluorescent markers for sensory and motor neurons to determine if specific neuron types were also accessible to isolation by FACS. Touch neurons (PLM and ALM) and serotonergic sensory neurons (ADF) from L1 larvae (Fig. 2A, B) and GABA motor neurons from L4 stage animals (Fig. 2C) were readily detected in primary cultures. To confirm that specific neuron classes could be purified by FACS, we isolated A-class (Fig. 2D) and VB type motor neurons from L2 stage larvae (Fig. 2E, F). Detection of ADF neurons in culture, which comprise only 0.3% of all L1 cells (2/558), suggested that rare types of neurons could also be isolated by this approach for gene expression profiling.

[Figure omitted. See PDF.]

Figure 2. Specific sensory and motor neurons are accessible to isolation by FACS from multiple larval stages.

Primary cultures of L1 stage cells 24 hr after dissociating from transgenic lines expressing (A) mec-4::mCherry to mark ALM and PLM neurons (red) and (B) srh-142p::dsRed to label ADF sensory neurons (red). (C) GABA motor neuron marked with unc-47::mCherry (red) and cultured for 24 hour after dissociating from L4 larval animals. (D) Primary culture of A-class motor neurons marked with unc-4::GFP (green) and isolated by FACS from L2 stage larvae. (E) del-1::GFP labels VB motor neurons (green) in the ventral nerve cord of an L2 stage larva. Anterior to left. (F) FACS profile of cells dissociated from del-1::GFP L2 larvae. Propidium Iodide (PI) marks dead cells (blue). Viable GFP-labeled cells (green) are outlined with the box. Scale bars are 10 microns.

https://doi.org/10.1371/journal.pone.0112102.g002

To test the applicability of this strategy to rare neuron classes, we chose the NSM serotonergic neurons for FACS isolation and cell-specific profiling experiments. Two NSM neurons are generated as a left-right pair (NSML, NSMR) in the embryo [30] and function as serotonerigic neurosecretory neurons in the pharynx (Fig. 3A) [31], [32]. To mark these cells, we used a tph-1::GFP transgene that exclusively labels NSML and NSMR with a strong GFP signal in L1 larvae [33] (Fig. 3B) (See Experimental Procedures). The lateral branches that are characteristic of adult NSM neurons (Fig. 3A) emerge later in development at the L4 stage [1]. The FACS profile shows that bright tph-1::GFP-labeled neurons were readily dissociated from L1 larvae and represented ∼0.3% of all cells, a fraction in good agreement with the in vivo proportion (2/558 = 0.36%) of NSM neurons. To confirm viability, FACS-isolated tph-1::GFP labeled cells were cultured in vitro where they attached to the substrate and initiated process outgrowth (Fig. 3D–F).

[Figure omitted. See PDF.]

Figure 3. Isolation of NSM serotonergic neurons by FACS.

(A) Head region of an adult hermaphrodite depicting the NSML neuron (orange) in the pharynx (green) (WormAtlas). Anterior is to left. (B) Confocal image of NSML in L1 larval stage. Anterior is to left. (C) FACS profile of cells dissociated from L1 stage larvae expressing tph-1::GFP in NSM neurons. Propidium Iodide (PI) marks dead cells (blue). Viable GFP-labeled cells (green) are outlined with the box. Images of FACS-isolated NSM neurons (D) 24 hr, (E), 48 hr and (F) 72 hr after plating. Scale bars are 10 microns.

https://doi.org/10.1371/journal.pone.0112102.g003

Generating RNA-Seq profiles from small amounts of total RNA

Preparations of dissociated cells were generated from approximately 3 million L1 larvae. We obtained 5–10 ng of total RNA from 30,000–50,000 FACS-isolated NSM neurons (∼0.25 pg total RNA/cell). Samples were amplified using the NuGEN Ovation V2 protocol to produce cDNA for library construction. Independent replicates were generated from two separate samples of FACS-isolated NSM neurons. Total RNA was also obtained from whole L1 larvae to produce duplicate reference data sets of transcripts expressed in all larval cells at this developmental stage. RNA-Seq data were collected for each sample in a single lane of an Illumina Hi-Seq 2000 to produce >150 million of either PE-75 (paired-end 75) or PE-100 (paired-end 100) reads for each independent library. Each data set yielded at least 8 million reads that map to non-ribosomal RNA genes and detected >10,000 coding genes at a threshold value of ≥1 FPKM (Table S1). This initial list was filtered to exclude transcripts that likely derived from contaminating cells in the FACS-isolated NSM preparation to yield an conservative estimate of 6,200 coding sequence transcripts that are expressed in the NSM neurons (See Experimental Procedures) (File S1). Fig. 4B, C shows average 5′ to 3′ sequence coverage and reads mapping to selected protein coding genes. Fig. 4D features an example of reads mapping to a noncoding RNA gene as expected since the amplification method uses both poly dT and random primers for cDNA synthesis [34].

[Figure omitted. See PDF.]

Figure 4. 5′-to-3′ coverage.

(A) The average relative coverage for each library is shown for normalized transcript length. The Y-axis denotes the average coverage of each position in the transcript normalized to the total number of mapped reads. Coverage for coding sequence genes (B) tph-1 and (C) bas-1 and for (D) a small nuclear RNA (snoRNA) gene, ZK643.9. Lighter vs darker shades of color depict results for independent NSM (blue) and reference (green) samples.

https://doi.org/10.1371/journal.pone.0112102.g004

RNA-Seq profiles of larval NSM neurons detect highly expressed serotonergic transcripts

Duplicate RNA-Seq data sets for the NSM neurons are well correlated as are the reference samples (Fig. 5A, B) but they show significant differences when compared to each other (e.g., NSM vs reference) (Fig. 5C). We detected a total of 1,073 transcripts with significantly elevated expression (≥2.4 fold, ≤0.05 q-value) in the NSM data set vs all L1 larval cells (see File S1). To validate this result, we compared the data set of NSM-enriched transcripts to a list of genes annotated in WormBase as previously assigned to NSM by direct observation. The overlap is highly significant with 21 of 36 known NSM genes represented in the NSM-enriched data set (p = 1.6×e⁻¹⁹) and an additional 13 genes that are detected at ≥1 FPKM (Table S2). Of particular note in the NSM-enriched list are genes with established roles in NSM differentiation or function. For example, canonical serotonergic genes are highly expressed in NSM. These include enzymes tph-1 (tryptophan hydroxylase/TPH) and bas-1 (aromatic amino acid decarboxylase/AAADC) required for serotonin (5-HT) synthesis, a transporter for selective uptake of monamine neurotranmitters into synaptic vesicles, cat-1 (vesicular monamine transporter/VMAT) and the synaptically localized transporter, mod-5 (serotonin reuptake transporter/SERT). mRNAs for the POU and LIM homeodomain transcription factors, unc-86 and ttx-3, respectively, which are known to activate expression of these serotonergic genes [31], [35] are also enriched in the NSM profile (Fig. 5C–E). The specificity of this data set is also underscored by the depletion of transcripts that are known to be highly expressed in other tissues or in different developmental periods e.g., tnt-3 and unc-89 (body muscle), dpf-6 (pharynx, intestine) and cht-1 (embryo) (Fig. 5C) [11], [36], [37]. As an additional test of the reliability of the NSM-enriched data set, we scored in vivo expression of promoter-GFP reporters for additional genes on this list that have not been previously assigned to NSM. The transcripts selected for this experiment display a broad range of statistical ranks (42–591) and half (2/4) of the corresponding GFP reporters were positively identified as expressed in NSM (Fig. 6) (Table 1). Together, these results strongly support the conclusion that these data sets provide an accurate representation of transcripts that are highly expressed in NSM neurons in vivo.

[Figure omitted. See PDF.]

Figure 5. Differential RNA-Seq analysis detects transcripts that are highly expressed in NSM neurons.

Pairwise comparisons of (A) NSM and (B) Reference data sets with expression values represented as log₂(FPKM) and Spearman rank-order correlation coefficients, r_s. (C) Volcano scatter plot of log₁₀(p-value) vs log₂(fold change) of transcript expression in NSM neurons relative to the reference sample derived from all L1 larval cells. Significantly enriched or depleted transcripts (≥2.4 fold, p<0.012) are indicated in red. (D) Schematic of NSM serotonergic presynaptic terminus depicting genes that are highly enriched (E) in the NSM RNA-Seq profile. See text for additional information about specific genes.

https://doi.org/10.1371/journal.pone.0112102.g005

[Figure omitted. See PDF.]

Figure 6. GFP reporters validate NSM-enriched RNA-Seq data set.

Transgenic animals expressing promoter::GFP reporter genes for NSM-enriched transcripts. Confocal images of adult hermaphrodites, anterior to left, ventral down. A–B, Bright tph-1::GFP (vsIs45) expression is limited to NSML and NSMR; note dim tph-1::GFP signal in ADF neurons (arrows). C, Cartoon depicting NSM neurons marked with tph-1::GFP in the anterior pharyngeal bulb. Y39B6A.19::GFP (D, E) and F41E7.3::GFP (F, G) are specifically expressed in neurons, including NSM, head neurons and ventral cord motor neurons (asterisks). Y39B6A.19::GFP is also expressed in sensory neuron, PDE (E). Scale bar is 25 microns.

https://doi.org/10.1371/journal.pone.0112102.g006

[Figure omitted. See PDF.]

Table 1. Expression of promoter-GFP reporters for transcripts enriched in the NSM RNA-Seq data set.

https://doi.org/10.1371/journal.pone.0112102.t001

Additional RNA-Seq data sets were generated from these RNA samples after treatment with the DSN (Duplex-Specific Nuclease) method to deplete rRNA sequences [38]. As expected, most of the DSN-treated samples include a higher fraction of reads derived from non-rRNA sequences (Table S1) and show the strongest correlation with the corresponding RNA-Seq data sets obtained from Total RNA (e.g., NSM-DSN vs NSM total RNA, Fig. S1). Finally, the total RNA samples were also assayed on tiling arrays [11] (File S3). NSM-enriched transcripts were calculated from the DSN-treated RNA-Seq data sets and from the tiling array results to compile the union of NSM-enriched transcripts (879 genes) showing elevated expression in at least one of the three profiles (i.e., NSM Total RNA, NSM-DSN or NSM tiling array) (See Experimental Procedures) (File S4).

Discussion

Here we describe a new protocol, SeqCeL (RNA-Seq of C. elegans cells), for gene expression profiling of specific larval C. elegans cells. Our approach exploits a recently developed method for dissociating larval cells for primary culture [25], [39]. Although this earlier work reported that neurons are poorly preserved in these preparations, we have now demonstrated that diverse neuron types are readily released at fractions predicted by their representation in vivo. We showed that GFP-positive neurons constitute approximately 30% of all cells obtained from transgenic L1 stage larval animals expressing a pan-neural reporter; this fraction is comparable to the relative number of neurons in the newly hatched L1 larva (222 neurons/558 total cells; 40%). We have now exploited this finding to demonstrate that viable larval neurons can be readily isolated by FACS for microarray and RNA-Seq analysis. To illustrate the feasibility of this cell isolation procedure for profiling rare neuron types, we produced RNA-Seq profiles for the NSM serotonergic neurons, which occur as a single bilateral pair of cells in the C. elegans pharynx [1], [26]. These RNA-Seq data sets identified >6,200 coding sequence transcripts (“Expressed” tab, File S1) in the FACS-isolated cells (See Experimental Procedures). Of these genes, 1,070 unique transcripts are also differentially expressed in the NSM neurons in comparison to all L1 stage larval cells. This NSM-enriched data set was validated by robust elevation of known NSM mRNAs such as genes encoding serotonin-related functions (e.g., tph-1, bas-1, mod-5) and the transcription factors that regulate their expression (e.g., unc-86, ttx-3) [27], [31]. Thus, these data sets provide a detailed list of NSM-expressed transcripts that are now available for experimental tests to reveal specific roles in NSM development and function.

This new SeqCeL method of using FACS to isolate neurons from C. elegans larvae for gene expression profiling complements the previously described application of this approach to embryonic cells (i.e., MAPCeL) [8], [9], [14]. Thus, neurons at multiple developmental stages, from embryos to adults (RLM and DMM, unpublished data), are now accessible for RNA-Seq analysis. What are the advantages and disadvantages of this technique in comparison to other approaches? We have previously used the mRNA tagging method, for example, to produce expression profiles of larval neurons. Although this technique is effective, it depends upon the production, in each case, of a custom-built transgenic line in which mRNA can be specifically immuno-precipitated from each target neuron [11], [13], [20], [21]. Gene expression profiles have also been generated from larval cells after microdissection [40] and by methods designed to isolate cell-specific nuclei [41], [42]. To date, however, these approaches have not been applied to neurons. In contrast, the FACS method of cell isolation that we have described here takes advantage of the large number of fluorescent live-cell markers that are already available for most neurons types in C. elegans [15] and extensively annotated in WormBase. Furthermore, for neurons that may not be specifically marked with a single reporter gene, pairs of different-colored labels (e.g., GFP + mCherry) that are uniquely co-expressed in the target neuron can be readily utilized for FACS isolation [11]; this combinatorial strategy greatly enhances the immediate applicability of our approach [43]. Although we have successfully generated robust neuron-specific profiles, these RNA-seq data sets are dominated by reads derived from ribosomal RNA (rRNA). We have confirmed however, that rRNA sequences can be relatively depleted by the DSN method. Recently, we have also successfully used alternative protocols for excluding rRNA [Ribo-Zero (Invitrogen) and SMARTer (Clontech)] before library construction to produce RNA-Seq data sets from small amounts (<10 ng) of total RNA (RLM, DMM and RHW, unpublished data). These approaches offer the important advantage of reducing the net sequencing bandwidth for each sample. Thus, it should be possible in the future to produce representative neuron-specific profiles of the non-rRNA transcriptome at significantly lower cost. We note the possibility that the cell isolation procedure could potentially alter gene expression. However, our analysis of the NSM data set (Fig. 5D) and of other neuron-specific profiles that we have produced by this method (RLM and DMM, data not shown) has confirmed that known cell-specific transcripts are faithfully represented in these results.

Conclusion

We have reported a new method, SeqCeL (RNA-Seq of C. elegans cells), that uses FACS to purify fluorescently-labeled neurons from C. elegans larvae for gene expression profiling by RNA-Seq analysis. We validated this approach by demonstrating that multiple different neuron types are available for FACS isolation and by confirming that an RNA-Seq data set obtained from the NSM serotonergic neurons includes most known NSM-expressed genes. Given the large number and ready availability of neuron-specific markers (e.g., GFP, mCherry, etc.) for C. elegans it should now be possible to extend this approach a wide array of different neuron classes.

Experimental Procedures

Nematode strains

C. elegans cultures were maintained as described [42]. N2 was used as the wild-type strain. Additional strains used in this study:

1. NW1229 [dpy-20(e1362) IV; evIs111 (F25B3.3::GFP + dpy-20(+))] [29]

2. oyIs51(srh-142p::dsRed) [44], LX837 [vsIs45 (tph-1p::GFP)] [33], NC197 [dpy-20(e1362) IV; wdIs4 (unc-4::GFP + dpy-20(+)] [45], NC2537 [unc-119(ed3); wdEx848 (mec-4::mCherry + unc-119(+))] [12]; NC138[dpy-20 (e1282) wdIs3(X)] [46]; XE1374 [wpIs39 (unc-47::mCherry)]

3. BC14934 [dpy-5(e907); sEx14934 (rCesZC317.3::GFP + pCeh361)]

4. BC12468 [dpy-5(e907); sEx12468 (rCesF56F3.6::GFP + pCeh361)]

5. BC13337 [dpy-5(e907); sEx13337 (rCesY39B6A.19::GFP + pCeh361)]

6. BC12693 [dpy-5(e907); sEx12693 (rCesF41E7.3::GFP + pCeh361)]

Preparation of dissociated larval cells and primary cultures

Worm strains were initially grown on 8P nutrient agar 150 mm plates (see recipe below) seeded with E. coli strain NA22 to produce a thick bacterial lawn [21]. For synchronized cultures of L1 larvae, embryos produced by hypochlorite treatment of adult hermaphrodites were allowed to hatch overnight in M9 buffer at 20°C. Approximately 3 million synchronized L1 larvae were typically used for each preparation (Determined by visually counting a sample aliquot on a glass microscope slide). L2 and L4 larval animals were obtained after growth of synchronized L1 larvae on NA22-seeded plates. For isolation of NSM neurons, synchronized L1 larvae were generated from vsIs45 (tph-1p::GFP) worms grown on twenty 8P nutrient agar plates (150 mm) seeded with NA22. Bright GFP expression is limited to NSM neurons in vsIs45 L1 larvae (Fig. 3B, Fig. 6A, B) [33] ADF neurons show dim GFP expression in vsIs45 but are readily excluded by FACS from the NSM preparation (Fig. 3C). Preparations of dissociated larval cells were generated as previously described [23], [34] with minor changes. L2 and L4 larvae were allowed to settle in cold M9 buffer on ice for 2×30 minute periods to remove bacteria before generating the larval cell prep. Cells were passed through a 5 micron filter at the final step prior to sorting [9]. Primary cell cultures were maintained as previously described [25], [28]. Recipe for 1 liter of 8P nutrient agar: 25 g agar, 20 g bactopeptone, 3 g NaCl [47].

FACS analysis

Sorting experiments were performed as previously described [7], [9] on a BD FACSAria equipped with a 70 micron diameter nozzle. Propidium iodide or DAPI was included to mark damaged cells. Profiles of GFP and mCherry marker strains were compared to an N2 standard to exclude auto-fluorescent cells. Sorted cells were collected in L-15-10 cell culture medium for primary cultures or in Trizol LS for RNA extraction. Yields of target neurons were calculated as the fraction of potentially available neurons that were isolated by FACS (Table S3). For example, we obtained 30,000–50,000 FACS-isolated NSM neurons from approximately 3 million L1 larvae. With two NSM neurons per animal, 3 million L1 larvae should contain 6 million NSM neurons for an overall yield of <0.85% of purified NSM neurons. FACS yields of <1% were also obtained for other classes of neurons (Table S3). Although yields of sorted neurons are low, this approach is not limiting because large numbers of L1 larvae are easily generated by standard culture conditions (see above).

RNA extraction and sample preparation for tiling arrays and RNA-Seq analysis

Total RNA was isolated and amplified for application to Affymetrix tiling arrays [9]. RNA integrity and final concentrations were determined in an Agilent Bioanalyzer. For RNA-Seq analysis, purified total RNA was amplified with Ovation RNA-Seq System V2 (NuGEN) and libraries were sequenced using the HiSeq 2000 system (Illumina). The DSN protocol was used to deplete rRNA sequences [38].

Microarray analysis

A custom chip definition file was generated using probes mapped to gene models for expression quantification as previously described [9]. Tiling array data were quantile normalized and median polished using RMA from the Affy package (v1.43.2) in Bioconductor (v2.14) [48]. A linear model and moderated t-statistic were used to determine differentially expressed genes as implemented by the limma package (v3.21.4) [49].

RNA sequencing analysis

Sequencing reads were mapped to the Caenorhabditis elegans genome and transcriptome (UCSC ce10, Wormbase WS220) using annotation supplied by Illumina (ftp://igenome:[email protected]/Caenorhabditis_elegans/UCSC/ce10/Caenorhabditis_elegans_UCSC_ce10.tar.gz). Reads were mapped using the splice-aware aligner tophat2 v2.0.11 [50]. Gene expression quantification and differential expression was analyzed using cufflinks v2.2.1 [51]. FPKM expression values were normalized by upper-quartile normalization [52]; FPKMs were scaled by the ratio of 75^th percentile fragment counts to the mean 75^th percentile value across all NSM and reference libraries; reads derived from rRNA genes were effectively excluded from this calculation. Gene body coverage analysis was performed using RSeQC v2.3.9. Genomic alignments were visualized using the UCSC genome browser with raw read counts normalized to Reads Per Million mapped reads (RPM). All RNA-Seq data sets used for this work are available at the NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra). An independent analysis of these RNA-Seq data sets is described in a separate publication [53].

Hypergeometric probability was calculated with publicly available software at http://www.geneprof.org/GeneProf/tools/hypergeometric.jsp.

Transcripts showing significant expression in NSM neurons were filtered by the test NSM FPKM < reference FPKM to remove false positives that are likely due to contamination from non-GFP cells in the preparation of FACS-isolated NSM neurons [11], [21]. For the union of all three NSM enriched data sets (i.e., NSM Total RNA, NSM-DSN or NSM tiling array), a Bonferroni-style correction was applied to the FDR threshold to limit the accumulation of false positives [54].

Microscopy

DIC and epifluorescence images of dissociated cells in culture were collected on a Zeiss Axiovert microscope with 40X and 63X objectives. Confocal images were obtained on Leica TCS SP5 and Nikon Eclipse TI confocal microscopes.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Heat-map depicting Spearman correlation coefficients for pair-wise comparisons of RNA-Seq data sets.

https://doi.org/10.1371/journal.pone.0112102.s001

(PDF)

Table S1.

Summary results for RNA-Seq data sets.

https://doi.org/10.1371/journal.pone.0112102.s002

(XLSX)

Table S2.

Known NSM-expressed genes are highly represented in the NSM RNA-Seq data set.

https://doi.org/10.1371/journal.pone.0112102.s003

(DOCX)

Table S3.

Yields for larval neurons isolated by FACS.

https://doi.org/10.1371/journal.pone.0112102.s004

(DOCX)

File S1.

Transcripts detected by RNA-Seq of total RNA from L1 Reference samples and NSM neurons. File contains worksheets for differential expression results. “All” contains all cuffdiff output. “Significant” contains all genes called differentially expressed ≥2.4X and q-value≤0.05. “Enriched” contains only genes significantly enriched in the NSM neurons versus the reference. “Expressed” lists genes that have a FPKM value >1 in either the NSM samples or reference samples.

https://doi.org/10.1371/journal.pone.0112102.s005

(XLSX)

File S2.

Transcripts detected by RNA-Seq of DSN-treated total RNA from L1 Reference sample and NSM neurons. This file contains worksheets for differential expression results. “All” contains all cuffdiff output. “Significant” contains all genes called differentially expressed ≥3.2X and q-value≤0.05. “Enriched” lists genes significantly enriched in the DSN-treated NSM neuron samples versus the DSN-treated reference. “Expressed” lists genes with FPKM >1 in either the DSN-treated NSM or DSN-treated reference samples.

https://doi.org/10.1371/journal.pone.0112102.s006

(XLSX)

File S3.

Transcripts detected by tiling array analysis of total RNA from L1 Reference sample and NSM neurons. This file contains worksheets for differential expression results of tiling array analysis. “All” contains all limma output. “Enriched” contains only genes significantly enriched in the NSM neurons versus the reference ≥1.5X and ≤5% FDR.

https://doi.org/10.1371/journal.pone.0112102.s007

(XLSX)

File S4.

Union of NSM-enriched transcripts from Files S1, S2 and S3.

https://doi.org/10.1371/journal.pone.0112102.s008

(XLSX)

Acknowledgments

We thank P. Sengupta (Brandeis University), J. Culotti (Mount Sinai Hospital, Toronto) D. Colon-Ramos (Yale University), M. Hammarlund (Yale University) and the C. elegans Genetics Center (University of Minnesota) for strains, D. Colon-Ramos and S. Jang (Yale University) for the confocal image of NSM in Fig. 3, Vanderbilt Flow Cytometry Shared Resource for FACS analysis and Vanderbilt Technologies for Advanced Genomics (VANTAGE) for generating microarray and RNA-Seq data sets.

Author Contributions

Conceived and designed the experiments: WCS DMM RHW. Performed the experiments: WCS RM TM PS. Analyzed the data: WCS OT LWH DMM. Contributed reagents/materials/analysis tools: RM PS OT WCS. Contributed to the writing of the manuscript: RHW LWH WCS DMM.

Citation: Spencer WC, McWhirter R, Miller T, Strasbourger P, Thompson O, Hillier LW, et al. (2014) Isolation of Specific Neurons from C. elegans Larvae for Gene Expression Profiling. PLoS ONE 9(11): e112102. https://doi.org/10.1371/journal.pone.0112102

References

1. White JG, Southgate E, Thomson JN, Brenner S (1986) The structure of the nervous system of the nematode Caenorhabditis elegans. Philosophical transactions of the Royal Society of London Series B, Biological sciences 314: 1–340.

2. Jarrell TA, Wang Y, Bloniarz AE, Brittin CA, Xu M, et al. (2012) The connectome of a decision-making neural network. Science 337: 437–444.

3. Varshney LR, Chen BL, Paniagua E, Hall DH, Chklovskii DB (2011) Structural properties of the Caenorhabditis elegans neuronal network. PLoS computational biology 7: e1001066.

4. Shipley FB, Clark CM, Alkema MJ, Leifer AM (2014) Simultaneous optogenetic manipulation and calcium imaging in freely moving C. elegans. Frontiers in neural circuits 8: 28.

5. Husson SJ, Costa WS, Wabnig S, Stirman JN, Watson JD, et al. (2012) Optogenetic analysis of a nociceptor neuron and network reveals ion channels acting downstream of primary sensors. Current biology: CB 22: 743–752.

6. Gerstein MB, Lu ZJ, Van Nostrand EL, Cheng C, Arshinoff BI, et al. (2010) Integrative analysis of the Caenorhabditis elegans genome by the modENCODE project. Science 330: 1775–1787.

7. Hillier LW, Reinke V, Green P, Hirst M, Marra MA, et al.. (2009) Massively parallel sequencing of the poly-adenylated transcriptome of C. elegans. Genome Res.

8. Zhang Y, Ma C, Delohery T, Nasipak B, Foat BC, et al. (2002) Identification of genes expressed in C. elegans touch receptor neurons. Nature 418: 331–335.

9. Fox RM, Von Stetina SE, Barlow SJ, Shaffer C, Olszewski KL, et al. (2005) A gene expression fingerprint of C. elegans embryonic motor neurons. BMC Genomics 6: 42.

10. Petersen SC, Watson JD, Richmond JE, Sarov M, Walthall WW, et al. (2011) A Transcriptional Program Promotes Remodeling of GABAergic Synapses in Caenorhabditis elegans. The Journal of neuroscience: the official journal of the Society for Neuroscience 31: 15362–15375.

11. Spencer WC, Zeller G, Watson JD, Henz SR, Watkins KL, et al. (2011) A spatial and temporal map of C. elegans gene expression. Genome Res 21: 325–341.

12. Smith CJ, O’Brien T, Chatzigeorgiou M, Spencer WC, Feingold-Link E, et al. (2013) Sensory Neuron Fates Are Distinguished by a Transcriptional Switch that Regulates Dendrite Branch Stabilization. Neuron 79: 266–280.

13. Smith CJ, Watson JD, Spencer WC, O’Brien T, Cha B, et al. (2010) Time-lapse imaging and cell-specific expression profiling reveal dynamic branching and molecular determinants of a multi-dendritic nociceptor in C. elegans. Dev Biol 345: 18–33.

14. McKay SJ, Johnsen R, Khattra J, Asano J, Baillie DL, et al. (2003) Gene expression profiling of cells, tissues, and developmental stages of the nematode C. elegans. Cold Spring Harb Symp Quant Biol 68: 159–169.

15. Hunt-Newbury R, Viveiros R, Johnsen R, Mah A, Anastas D, et al. (2007) High-throughput in vivo analysis of gene expression in Caenorhabditis elegans. PLoS Biol 5: e237.

16. Cinar H, Keles S, Jin Y (2005) Expression profiling of GABAergic motor neurons in Caenorhabditis elegans. Curr Biol 15: 340–346.

17. Colosimo ME, Brown A, Mukhopadhyay S, Gabel C, Lanjuin AE, et al. (2004) Identification of thermosensory and olfactory neuron-specific genes via expression profiling of single neuron types. Curr Biol 14: 2245–2251.

18. Guillermin ML, Castelletto ML, Hallem EA (2011) Differentiation of carbon dioxide-sensing neurons in Caenorhabditis elegans requires the ETS-5 transcription factor. Genetics 189: 1327–1339.

19. Hallem EA, Spencer WC, McWhirter RD, Zeller G, Henz SR, et al. (2011) Receptor-type guanylate cyclase is required for carbon dioxide sensation by Caenorhabditis elegans. Proc Natl Acad Sci U S A 108: 254–259.

20. Roy PJ, Stuart JM, Lund J, Kim SK (2002) Chromosomal clustering of muscle-expressed genes in Caenorhabditis elegans. Nature 418: 975–979.

21. Von Stetina S, Watson J, Fox R, Olszewski K, Spencer W, et al. (2007) Cell-specific microarray profiling experiments reveal a comprehensive picture of gene expression in the C. elegans nervous system. Genome Biol 8: R135.

22. Watson JD, Wang S, Von Stetina SE, Spencer WC, Levy S, et al. (2008) Complementary RNA amplification methods enhance microarray identification of transcripts expressed in the C. elegans nervous system. BMC Genomics 9: 84.

23. Von Stetina SE, Fox RM, Watkins KL, Starich TA, Shaw JE, et al. (2007) UNC-4 represses CEH-12/HB9 to specify synaptic inputs to VA motor neurons in C. elegans. Genes Dev 21: 332–346.

24. Chatzigeorgiou M, Yoo S, Watson JD, Lee WH, Spencer WC, et al. (2010) Specific roles for DEG/ENaC and TRP channels in touch and thermosensation in C. elegans nociceptors. Nature neuroscience 13: 861–868.

25. Zhang S, Banerjee D, Kuhn JR (2011) Isolation and culture of larval cells from C. elegans. PloS one 6: e19505.

26. Albertson DG, Thomson JN (1976) The pharynx of Caenorhabditis elegans. Philosophical transactions of the Royal Society of London Series B, Biological sciences 275: 299–325.

27. Zhang F, Bhattacharya A, Nelson JC, Abe N, Gordon P, et al. (2014) The LIM and POU homeobox genes ttx-3 and unc-86 act as terminal selectors in distinct cholinergic and serotonergic neuron types. Development 141: 422–435.

28. Christensen M, Estevez A, Yin X, Fox R, Morrison R, et al. (2002) A primary culture system for functional analysis of C. elegans neurons and muscle cells. Neuron 33: 503–514.

29. Altun-Gultekin Z, Andachi Y, Tsalik EL, Pilgrim D, Kohara Y, et al. (2001) A regulatory cascade of three homeobox genes, ceh-10, ttx-3 and ceh-23, controls cell fate specification of a defined interneuron class in C. elegans. Development 128: 1951–1969.

30. Sulston J, Schierenberg E, White J, Thomson J (1983) The embryonic cell lineage of the nematode Caenorhabditis elegans. Dev Biol 100: 64–119.

31. Sze JY, Zhang S, Li J, Ruvkun G (2002) The C. elegans POU-domain transcription factor UNC-86 regulates the tph-1 tryptophan hydroxylase gene and neurite outgrowth in specific serotonergic neurons. Development 129: 3901–3911.

32. Horvitz HR, Chalfie M, Trent C, Sulston JE, Evans PD (1982) Serotonin and octopamine in the nematode Caenorhabditis elegans. Science 216: 1012–1014.

33. Nelson JC, Colon-Ramos DA (2013) Serotonergic neurosecretory synapse targeting is controlled by netrin-releasing guidepost neurons in Caenorhabditis elegans. The Journal of neuroscience: the official journal of the Society for Neuroscience 33: 1366–1376.

34. Adiconis X, Borges-Rivera D, Satija R, DeLuca DS, Busby MA, et al. (2013) Comparative analysis of RNA sequencing methods for degraded or low-input samples. Nature methods 10: 623–629.

35. Zhang S, Sokolchik I, Blanco G, Sze JY (2004) Caenorhabditis elegans TRPV ion channel regulates 5HT biosynthesis in chemosensory neurons. Development 131: 1629–1638.

36. Meissner B, Warner A, Wong K, Dube N, Lorch A, et al. (2009) An integrated strategy to study muscle development and myofilament structure in Caenorhabditis elegans. PLoS genetics 5: e1000537.

37. Shen C, Nettleton D, Jiang M, Kim SK, Powell-Coffman JA (2005) Roles of the HIF-1 hypoxia-inducible factor during hypoxia response in Caenorhabditis elegans. The Journal of biological chemistry 280: 20580–20588.

38. Shagina I, Bogdanova E, Mamedov IZ, Lebedev Y, Lukyanov S, et al. (2010) Normalization of genomic DNA using duplex-specific nuclease. BioTechniques 48: 455–459.

39. Zhang S, Kuhn JR (2013) Cell isolation and culture. WormBook: the online review of C elegans biology: 1–39.

40. Schwarz EM, Kato M, Sternberg PW (2012) Functional transcriptomics of a migrating cell in Caenorhabditis elegans. Proceedings of the National Academy of Sciences of the United States of America 109: 16246–16251.

41. Haenni S, Ji Z, Hoque M, Rust N, Sharpe H, et al. (2012) Analysis of C. elegans intestinal gene expression and polyadenylation by fluorescence-activated nuclei sorting and 3′-end-seq. Nucleic acids research 40: 6304–6318.

42. Steiner FA, Talbert PB, Kasinathan S, Deal RB, Henikoff S (2012) Cell-type-specific nuclei purification from whole animals for genome-wide expression and chromatin profiling. Genome research 22: 766–777.

43. Zhang S, Ma C, Chalfie M (2004) Combinatorial marking of cells and organelles with reconstituted fluorescent proteins. Cell 119: 137–144.

44. Kim K, Kim R, Sengupta P (2010) The HMX/NKX homeodomain protein MLS-2 specifies the identity of the AWC sensory neuron type via regulation of the ceh-36 Otx gene in C. elegans. Development 137: 963–974.

45. Pflugrad A, Meir JY, Barnes TM, Miller DM 3rd (1997) The Groucho-like transcription factor UNC-37 functions with the neural specificity gene unc-4 to govern motor neuron identity in C. elegans. Development 124: 1699–1709.

46. Winnier AR, Meir JY, Ross JM, Tavernarakis N, Driscoll M, et al. (1999) UNC-4/UNC-37-dependent repression of motor neuron-specific genes controls synaptic choice in Caenorhabditis elegans. Genes Dev 13: 2774–2786.

47. Schachat F, Garcea RL, Epstein HF (1978) Myosins exist as homodimers of heavy chains: demonstration with specific antibody purified by nematode mutant myosin affinity chromatography. Cell 15: 405–411.

48. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, et al. (2004) Bioconductor: open software development for computational biology and bioinformatics. Genome biology 5: R80.

49. Smyth GK, Michaud J, Scott HS (2005) Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics 21: 2067–2075.

50. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, et al. (2013) TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome biology 14: R36.

51. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, et al. (2012) Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nature protocols 7: 562–578.

52. Bullard JH, Purdom E, Hansen KD, Dudoit S (2010) Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC bioinformatics 11: 94.

53. Gerstein MB RJ, Yan K, Wang D, Cheng C, Brown JB, et al.. (2014) Comparative analysis of the transcriptome across distant species. Nature: in press.

54. Bonferroni E (1935) Il calcolo delle assicurazioni su gruppi di teste. Rome, Italy.