Genome-Wide Gene Expression Profiling of Nucleus Accumbens Neurons Projecting to Ventral Pallidum Using both Microarray and Transcriptome Sequencing

• Published in: Frontiers in neuroscience 2011-01, Vol.5, p.98-98
• Main Authors: Chen, Hao, Liu, Zhimin, Gong, Suzhen, Wu, Xingjun, Taylor, William L, Williams, Robert W, Matta, Shannon G, Sharp, Burt M
• Format: Article
• Language: English
• Subjects: Addictions; Bioinformatics; Brain; DNA microarrays; Exons; GABA; GABAergic; Gene expression; Genomes; Genomics; Laboratory animals; Laser capture microdissectio; Lasers; Microarray; Neurons; Neuroscience; Neurosciences; Nucleus Accumbens; Pallidum (ventral); Reinforcement; Ribonucleic acid; RNA; RNA sequencing
• ISSN: 1662-453X; 1662-4548; 1662-453X
• DOI: 10.3389/fnins.2011.00098

The cellular heterogeneity of brain poses a particularly thorny issue in genome-wide gene expression studies. Because laser capture microdissection (LCM) enables the precise extraction of a small area of tissue, we combined LCM with neuronal track tracing to collect nucleus accumbens shell neurons that project to ventral pallidum, which are of particular interest in the study of reward and addiction. Four independent biological samples of accumbens projection neurons were obtained. Approximately 500 pg of total RNA from each sample was then amplified linearly and subjected to Affymetrix microarray and Applied Biosystems sequencing by oligonucleotide ligation and detection (SOLiD) transcriptome sequencing (RNA-seq). A total of 375 million 50-bp reads were obtained from RNA-seq. Approximately 57% of these reads were mapped to the rat reference genome (Baylor 3.4/rn4). Approximately 11,000 unique RefSeq genes and 100,000 unique exons were identified from each sample. Of the unmapped reads, the quality scores were 4.74 ± 0.42 lower than the mapped reads. When RNA-seq and microarray data from the same samples were compared, Pearson correlations were between 0.764 and 0.798. The variances in data obtained for the four samples by microarray and RNA-seq were similar for medium to high abundance genes, but less among low abundance genes detected by microarray. Analysis of 34 genes by real-time polymerase chain reaction showed higher correlation with RNA-seq (0.66) than with microarray (0.46). Further analysis showed 20-30 million 50-bp reads are sufficient to provide estimates of gene expression levels comparable to those produced by microarray. In summary, this study showed that picogram quantities of total RNA obtained by LCM of ∼700 individual neurons is sufficient to take advantage of the benefits provided by the transcriptome sequencing technology, such as low background noise, high dynamic range, and high precision.